{!Dane}
n_dot_g=0,7[kmol/s]
beta=24
T[1]=285[K]
p[1]=101[kPa]
T[4]=1560[K]
phi_0=0,88
eta_iS=0,9
eta_iT=0,94
p_od=101,325[kPa]
T_od=292[K]

{!Stechiometria}

z1_CH4=1 "Paliwo"
n1_C=z1_CH4 "ilość substancji palnej w paliwie"
n1_H2=2*z1_CH4

z1_O2as=0,20946
z1_N2as=1-z1_O2as
z1_H2O=phi_0*p_sat0/(p_od-phi_0*p_sat0)
p_sat0=p_sat(Steam_IAPWS;T=T[1])

z1_a=1+z1_H2O
z1_O2=z1_O2as/z1_a
z1_H2O=z1_H2Oas/z1_a
z1_N2=1-z1_O2-z1_H2O

{Minimalna ilość pow. do spalania}

n1_O2min=n1_C+0,5*n1_H2
n1_amin=n1_O2min/z1_O2
n1_a=lambda*n1_amin

{Skład spalin}

n2_CO2=n1_C
n2_H2O=n1_H2+n1_a*z1_H2O
n2_N2=n1_a*z1_N2
n2_O2=(1-lambda)*n1_O2min

{Udziały molowe składników w spalinach}
n2_s=n2_CO2+n2_O2+n2_H2O+n2_N2
z2_CO2=n2_CO2/n2_s
z2_O2=n2_O2/n2_s
z2_H2O=n2_H2O/n2_s
z2_N2=n2_N2/n2_s

n_dot_a=n1_a*n_dot_g
n_dot_sp=n2_s*n_dot_g

{!Entalpie odniesienia}

i_0g=enthalpy(CH4;T=T_od) "Paliwo"
s_0g=entropy(CH4;T=T_od;P=P_od)

i_0a=z1_O2*enthalpy(O2;T=T_od)+z1_N2*enthalpy(N2;T=T_od)+z1_H2O*enthalpy(H2O;T=T_od) "Powietrze"
s_0a=z1_O2*entropy(O2;T=T_od;P=P_od)+z1_N2*entropy(N2;T=T_od;P=P_od)+z1_H2O*entropy(H2O;T=T_od;P=P_od)

i_0sp=z2_O2*enthalpy(O2;T=T_od)+z2_N2*enthalpy(N2;T=T_od)+z2_H2O*enthalpy(H2O;T=T_od)+z2_CO2*enthalpy(CO2;T=T_od) "Spaliny"
s_0sp=z2_O2*entropy(O2;T=T_od;P=P_od)+z2_N2*entropy(N2;T=T_od;P=P_od)+z2_H2O*entropy(H2O;T=T_od;P=P_od)+z2_CO2*entropy(CO2;T=T_od;P=P_od)


{!Obieg gazowy}
Duplicate i=1;2

i[i]=z1_O2*enthalpy(O2;T=T[i])+z1_N2*enthalpy(N2;T=T[i])+z1_H2O*enthalpy(H2O;T=T[i])-i_0a
s[i]=z1_O2*entropy(O2;T=T[i];P=P[i])+z1_N2*entropy(N2;T=T[i];P=P[i])+z1_H2O*entropy(H2O;T=T[i];P=P[i])-s_0a
End
Duplicate j=4;5

i[j]=z2_O2*enthalpy(O2;T=T[j])+z2_N2*enthalpy(N2;T=T[j]+z2_H2O*enthalpy(H2O;T=T[j]))-i_0sp

s[j]=z2_O2*entropy(O2;T=T[j];P=P[j])+z2_N2*entropy(N2;T=T[j];P=P[j])+z2_H2O*entropy(H2O;T=T[j];P=P[j])-s_0sp
End

{!Bilans sprezarki}
p[2]=p[1]*beta
eta_iS=(i_2s-i[1])/(i[2]-i[1])
i_2s=(z2_O2*enthalpy(O2;s=entropy(O2;T=t[1];P=P[1]);P=P[2]) +z2_H2O*enthalpy(H2O;s=entropy(H2O;T=T[1];P=P[1]);P=P[2]) + z2_N2*enthalpy(N2;s=entropy(N2;T=T[1];P=P[1]);P=P[2]))-i_0a


{!Bilans komory spalania}

p[2]=p[4]
T[3]=T[1]
p[3]=p[2]
I_dot_a=n_dot_a*i[2]
I_dot_sp=n_dot_sp*i[2]
I_dot_g=n_dot_g*(enthalpy(CH4;T=T[3])-i_0g)
I_dot_=I_dot_g+I_dot_a