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

90[thetawater,pTbWQlHl2]=ode45('RatesWater',[thetas+thetab+thetaperiod:dtheta:thetas…+thetab+thetaperiod+thetaWa],[pw Tbw Ww Qlw Hlw],options);% specification of initial conditions at start of expansion phase e - start of expansionpe=interp1(thetawater,pTbWQlHl2(:,1),thetas+thetab+thetaperiod+thetaWa);Tbe=interp1(thetawater,pTbWQlHl2(:,2),thetas+thetab+thetaperiod+thetaWa);We=interp1(thetawater,pTbWQlHl2(:,3),thetas+thetab+thetaperiod+thetaWa);Qle=interp1(thetawater,pTbWQlHl2(:,4),thetas+thetab+thetaperiod+thetaWa);Hle=interp1(thetawater,pTbWQlHl2(:,5),thetas+thetab+thetaperiod+thetaWa);% integration during expansion phase%disp(['integrating over the expansion phase']);[thetaexp,pTbWQlHl]=ode45('RatesExp', ...[thetas+thetab+thetaperiod+thetaWa:dtheta:pi],[pe Tbe We Qle Hle],options);% error checksmass4=mass1*exp(-Cblowby*2*pi/omega);p4=interp1(thetaexp,pTbWQlHl(:,1),pi);T4=interp1(thetaexp,pTbWQlHl(:,2),pi);W4=interp1(thetaexp,pTbWQlHl(:,3),pi);Ql4=interp1(thetaexp,pTbWQlHl(:,4),pi);Hl4=interp1(thetaexp,pTbWQlHl(:,5),pi);[h4,u4,v4,s4,Y4,cp4,dlvlT4,dlvlp4]= ...farg(p4,T4,phi,1,fueltype,airscheme);U4=u4*mass4;error1=1-v4*mass4/Vbdc;error2=1+W4/(U4-U1+Ql4+Hl4);% indicated mean effective pressure and thermal efficiencyimep=W4/(pi*b^2/4*stroke);eta=W4/mass1*(1+phi*0.06548*(1-f))/phi/0.06548/(1-f)/47870/1e3;% calcuate the heat flux in W/m^2qcomp=calcq(thetacomp,pTuWQlHl,'comp'); % compressionqcombu=calcq(thetacomb,pTbTuWQlHl1,'combu'); % combustion-unburned zoneqcombb=calcq(thetacomb,pTbTuWQlHl1,'combb'); % combustion-burned zoneqperib=calcq(thetaperi,pTbWQlHl1p,'peri'); % periodjection-burned zoneqwaterb=calcq(thetawater,pTbWQlHl2,'inwater'); % injection-burned zoneqexp=calcq(thetaexp,pTbWQlHl,'exp'); % expansiontimefinish=cputime;timetaken=timefinish-timestart;% save all datasave ahrind.mattime=60/(RPM*360*(pi/180))PressureB=pw;PressureW=Pr*pw;TemperatureB=Tbw;TemperatureW=Tv;load('Volume', 'RTWV','RMBV','RPWV','RPBV','RTBV');TWdetaV=(RTWV*time)MBdetaV=(RMBV*time)