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

105MW=sum(Y.*M);R=Ru/MW;cv=cp-R;Qconv1=hcoeff*(pi*b^2/2+4*V/b)*(y(2)-Tw);Qconv2=hcoeff*(pi*b^2/2+4*V/b)*(Tv-Tw);Const1=1/omega/mass;yprime(3)=p*dVdtheta;yprime(4)=Const1*mass*(Qconv1-Qconv2);yprime(5)=Cblowby*mass/omega*h;Uchangeptheta=-yprime(4)-yprime(3)-yprime(5)-Cblowby*masswaterin*XSteam('h_pT',Pr*patm,Tin)*1000/omega;yprime(2)=Uchangeptheta/(mass*cv);yprime(1)=yprime(2)/Tb-p*dVdtheta/V;% Two-zone Thermodynamic Modelload('ecpdata','dlvlTb','dlvlpb','cpb','Rb','ub','vb');load('Wdata','dlvlTw','dlvlpw');% mass of waterdi=thetaWa*180/pi;mwaterT=eps1*phi*Xv*18.05/di;mwpmf=mwaterT/mfuel;masswaterT=massfuel*mwpmf;RW=Ru/18.05;QconvW=hcoeff*(pi*b^2/2+4*V/b)*(Tw-Tv);% heat lossQconvb=hcoeff*(pi*b^2/2+4*V/b)*(Tw-y(2));% heat lossvw=XSteam('v_pT',Pr*patm,Tin); % Specific volume of watercpw=XSteam('Cp_pT',Pr*patm,Tin)*1000*(1/274.15); % Specific isobaric heatcapacity of wateruw=XSteam('u_pT',Pr*patm,Tin)*1000; % Specific internal energy of waterhw=XSteam('h_pT',Pr*patm,Tin)*1000; % Specific enthapy of water% Equally VolumeRPW=0;%0.1*yprime(1)*omega;TemW1V=Tv*RPW*((1-(dlvlTw+dlvlpw)))/(Pr*patm*101.325*1000)+QconvW/(masswaterT*RW);TemW2V=(1+(cpw/RW-dlvlTw));RTWV=TemW1V/TemW2V;MasB1V=Qconvb;MasB2V=(y(1)+Pr*patm*101.325*1000)*V*((cpb/Rb-dlvlTb)+(dlvlTb+dlvlpb));MasB3V=ub-uw-RW*Tv-(cpb*y(2)-(y(1)+Pr*patm*101.325*1000)*vb*dlvlTb)+(p)*vw*((cpb/RbdlvlTb)+(dlvlTb+dlvlpb));RMBV=(MasB1V+MasB2V)/MasB3V;RPWV=0;RPBV=yprime(1)*omega;%dVdtheta*omegaTemB1V=RMBV*(hw-ub)+Qconvb-V*RPBV*(dlvlTb+dlvlpb);TemB2V=mass1*cpb-mass1*Rb*dlvlTb;RTBV=TemB1V/TemB2V;