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

100Iter=Iter+1;[hb,u,v,s,Y,cp,dlvlT,dlvlp]=ecpw(X,p,Tb,phi,fueltype,airscheme);DeltaT=(hu-hb)/cp;Tb=Tb+DeltaT;endif Iter>=MaxIterwarning('convergence failure in adiabatic flame temperature loop');endD-8. M-file for RatesComp.mfunction yprime=RatesComp(theta,y,flag);global b stroke eps r Cblowby f fueltype airscheme phi thetas thetab omega ...heattransferlaw hcu Tw theta1 Vtdc mass1 ...p=y(1);Tu=y(2);yprime=zeros(5,1);% mass in cylinder accounting for blowby:mass=mass1*exp(-Cblowby*(theta-theta1)/omega);% volume of cylinder:V=Vtdc*(1+(r-1)/2*(1-cos(theta)+ 1/eps*(1-(1-eps^2*sin(theta).^2).^0.5)));% derivate of volume:dVdtheta=Vtdc*(r-1)/2*(sin(theta)+eps/2*sin(2*theta)./sqrt(1-eps^2*sin(theta).^2));switch heattransferlawcase 'constant'hcoeff=hcu;case 'Woschni'upmean=omega*stroke/pi; % mean piston velocityC1=2.28;hcoeff=hcu*130*b^(-0.2)*Tu^(-0.53)*(p/100e3)^(0.8)*C1*upmean;endA=1/mass*(dVdtheta+V*Cblowby/omega);Qconv=hcoeff*(pi*b^2/2+4*V/b)*(Tu-Tw);Const1=1/omega/mass;X=0;[h,u,v,s,Y,cp,dlvlT,dlvlp]=fargw(X,p,Tu,phi,f,fueltype,airscheme);B=Const1*v/cp*dlvlT*Qconv/Tu; % note typo on p174, eq. 4.76C=0;D=0;E=v^2/cp/Tu*dlvlT^2+v/p*dlvlp;yprime(1)=(A+B+C)/(D+E);yprime(2)=-Const1/cp*Qconv+v/cp*dlvlT*yprime(1);yprime(3)=p*dVdtheta;yprime(4)=Qconv/omega;yprime(5)=Cblowby*mass/omega*h;D-9. M-file for RatesComb.mfunction yprime=RatesComb(theta,y,flag);