Scribd
Upload
311 views

Matlab Programs For Ignition Delay

The document contains 4 Matlab programs that model diesel fuel properties and combustion. 1) Calculates ignition delay of diesel-ethanol-biodiesel blends using the Hardenberg and Hase method. 2) Models variation of ignition delay with changing cetane number. 3) Models variation of cetane number with changing ethanol volume fraction in diesel-ethanol-biodiesel blends. 4) Models variation of cetane number with changing biodiesel volume fraction in diesel-ethanol-biodiesel blends.

Uploaded by

PravinJadhav
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
311 views

Matlab Programs For Ignition Delay

The document contains 4 Matlab programs that model diesel fuel properties and combustion. 1) Calculates ignition delay of diesel-ethanol-biodiesel blends using the Hardenberg and Hase method. 2) Models variation of ignition delay with changing cetane number. 3) Models variation of cetane number with changing ethanol volume fraction in diesel-ethanol-biodiesel blends. 4) Models variation of cetane number with changing biodiesel volume fraction in diesel-ethanol-biodiesel blends.

Uploaded by

PravinJadhav
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 6

Matlab Programs –

1) Calculation of ignition delay by Hardenberg and Hase


method:-
function t = ignitiondelaycn(CNd,ne,nbd,L,N,T,P)

% CNd = 52; Input 'CNd is cetane number of diesel'


% ne = 0.05; Input 'ne is ethanol ratio in the blend'
% nbd = 0.1; Input 'nbd is biodiesel ratio in the
blend'
% L = 0.1; Input 'L is stroke length in metres'
% N = 4000; Input 'N is speed in rpm'
% P = 45; Input 'P is charge pressure in atm'
% T = 500; Input 'T is charge temperature in Kelvin'

R = 8.314;
%------Calculation of mean piston speed--------
Sp = 2*L*(N/60);

%------Calculation of cetane number of diesel ethanol


biodiesel blend------
CN = CNd - 0.59*ne + 0.55*nbd;

%------Calculation of activation energy--------


Ea = 618840 / (CN + 25);

%------Calculation of ignition delay using hardenberg


and Hase method------
tca = (0.36+0.22*Sp)*exp(Ea*(1/(R*T)-1/17190)*(21.2/(P-
12.4))^0.63);
t = tca / (0.006 * N);
2) Variation of ignition delay with cetane number –
function t = ignitiondelaycn1(CNd,L,N,T,P)

% CNd = 51; Input 'CNd is cetane number of diesel'


% ne = 0.05; Input 'ne is ethanol ratio in the blend'
% nbd = 0.1; Input 'nbd is biodiesel ratio in the
blend'
% L = 0.1; Input 'L is stroke length in metres'
% N = 4000; Input 'N is speed in rpm'
% P = 45; Input 'P is charge pressure in atm'
% T = 500; Input 'T is charge temperature in
Kelvin'
R = 8.314;
%------Calculation of mean piston speed--------
Sp = 2*L*(N/60);

%------Calculation of cetane number of diesel ethanol


biodiesel blend------
ne = 0:2:20;
nbd = 0:2:20;
CN = CNd - 0.59*ne + 0.55*nbd;

%------Calculation of activation energy--------


Ea = 618840 ./ (CN + 25);

%------Calculation of ignition delay using hardenberg


and Hase method------
tca = (0.36+0.22*Sp)*exp(Ea*(1/(R*T)-1/17190)*(21.2/(P-
12.4))^0.63);
t = tca / (0.006 * N);

figure(1);
plot(CN,t);
xlabel('cetane number');
ylabel('ignition delay');
3) Variation of cetane number with volume fraction of
ethanol –

function CNbd = linearbd(CNd)

% CNd = 52; Input ‘cetane number of diesel’

%-----Volume Fraction of Biodiesel = 5%-------


nbd = 0.05;

ne = 0:0.05:0.4;

CNbd = CNd - 0.59 * ne + 0.55 * nbd;

plot(ne,CNbd);
xlabel('Volume fraction of ethanol');
ylabel('Cetane Number');
title('Variation of cetane number with blend % of
ethanol');
hold on; grid on;

%-----Volume Fraction of Biodiesel = 10%--------


nbd1 = 0.1;

ne = 0:0.05:0.4;

CNbd = CNd - 0.59 * ne + 0.55 * nbd1;

plot(ne,CNbd,'m');
xlabel('Volume fraction of ethanol');
ylabel('Cetane Number');
title('Variation of cetane number with blend % of
ethanol');
%-----Volume Fraction of Biodiesel = 20%-----
nbd2 = 0.2;

ne = 0:0.05:0.4;

CNbd = CNd - 0.59 * ne + 0.55 * nbd2;

plot(ne,CNbd,'k');
xlabel('Volume fraction of ethanol');
ylabel('Cetane Number');
legend('Biodiesel = 5%','Biodiesel = 10%','Biodiesel =
20%');
title('Variation of cetane number with blend % of
ethanol');
4) Variation of cetane number with volume fraction of
biodiesel –
function CNet = linearet(CNd)

% CNd = 52; ‘Input cetane number of diesel’

%-----Volume Fraction of Ethanol = 5%-----


ne = 0.05;

nbd = 0:0.05:0.4;

CNet = CNd - 0.59 * ne + 0.55 * nbd;

plot(nbd,CNet);
xlabel('Volume fraction of biodiesel');
ylabel('Cetane Number');
title('Variation of cetane number with blend % of
biodiesel');
hold on;

%-----Volume Fraction of Ethanol = 10%-----


ne1 = 0.1;

nbd = 0:0.05:0.4;

CNet = CNd - 0.59 * ne1 + 0.55 * nbd;

plot(nbd,CNet,'k');
xlabel('Volume fraction of biodiesel');
ylabel('Cetane Number');
title('Variation of cetane number with blend % of
biodiesel');
hold on; grid on;

%-----Volume Fraction of Ethanol = 20%------


ne2 = 0.2;
nbd = 0:0.05:0.4;

CNet = CNd - 0.59 * ne2 + 0.55 * nbd;

plot(nbd,CNet,'m');
xlabel('Volume fraction of biodiesel');
ylabel('Cetane Number');
legend('Ethanol = 5%','Ethanol = 10%','Ethanol = 20%');
title('Variation of cetane number with blend % of
biodiesel');

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy