analysis of water injection into high-temperature mixture of ...

88Appendix D MATLAB program scripts%% Script file to determine the performance of a fuel inducted engine% based on a (user-specified) arbitrary heat release profile as a% function of crank angle.% Method closely follows that of:% Ferguson, C.R., 1986, "Internal Combustion Engines", Wiley.%********************************************************************% input:% enginedata.m - this is another script file that defines all of the% relevant engine parameters and operating conditions.% output:% ahrind.mat - this file contains all of the variables. For plotting% the results, see the example script file plotresults.m%********************************************************************D-1. M-file for ahrind.mtimestart=cputime;global Xv Tv Pr b stroke eps r Cblowby f fueltype airscheme phi thetas thetabthetaWa thetaperiod RPM omega heattransferlaw hcu hcb Tw theta1 Vtdc Vbdcmass1 p1 T1 V1% load the engine parameters and initial conditionsenginedata[alpha,beta,gamma,delta,Afuel]=fueldata(fueltype);% mole fractions and molecular weight of fuel-aireps1=0.210/(alpha+0.25*beta-0.5*gamma);fuel=eps1*phi/(1+eps1*phi);o2=0.21/(1+eps1*phi);n2=0.79/(1+eps1*phi);Mfa=fuel*(12.01*alpha+1.008*beta+16*gamma+14.01*delta)+32*o2+28.02*n2;Mf=fuel*(12.01*alpha+1.008*beta+16*gamma+14.01*delta);% mass of fuel and air in the cylindermfuel=eps1*phi*(12.01*alpha+1.008*beta+16*gamma+14.01*delta);mair=32*0.21+28.02*0.79;mfpma=mfuel/mair;massair=mass1/(mfpma+1);massfuel=mass1-massair;%%%%%%%%%%%%%%%%%%%% mass of water injectionmwater=eps1*phi*Xv*18.05;mwpmf=mwater/mfuel;masswater=massfuel*mwpmf;%%%%%%%%%%%%%%%%%%%%%% Total mass of water injectionmassT=mass1+masswater;mass1=mass1;massair=massair;massfuel=massfuel;