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

96dfdp=zeros(4,1);dYdT=zeros(11,1);dYdp=zeros(11,1);B=zeros(4,1);% check if solid carbon will formeps=0.210/(alpha+0.25*beta-0.5*gamma);if phi>(0.210/eps/(0.5*alpha-0.5*gamma))error('phi too high - c(s) and other species will form');endif nargin==6 % no Yguess so estimate the composition using farg[h,u,v,s,Y,cp,dlvlT,dlvlp]=farg(p,T,phi,1,fueltype,airscheme);Y(7:11)=ones(5,1)*MinMol; % since farg only returns first 6 specieselseY=Yguess;endload('Ydata', 'Ym');Y = Ym;% mole water injectionWY=X*eps*phi;% evaluate constantsTKp=[log(T)-1 (T)/2 (T^2)/6 (T^3)/12 (T^4)/20 -(1/T) 1]';% 7x1K=exp(Kp*TKp);i=find(K