{mass balance}

mLdot=rho1*xA1*Q1-rho2*xA2*Q2

{Wagner Equation}

xA2=(exp(((-7.45514*J5)+(1.202*(J5^1.5))+(-2.43926*(J5^3))+(-3.3559*(J5^6)))/(T2/508.1))*47)

J5=1-(T2/508.1)

{Thermodynamic Balance}

mLNdot*hLN+mLNdot*EnthN1=mLdot*hA+mLdot*EnthL+(rho1*Q1)*EnthAN*(T3-T2)

{Correlations}

rho1=P*MWT1/(R*T3)

rho2=P*MWT2/(R*T2)

Q2=(m2dot/rho2)

MWT1=xA1*MWA+xN1*MWN

MWT2=xA2*MWA+xN2*MWN

xA1+xN1=1

xA2+xN2=1

m2dot=m1dot-mLdot

m1dot=rho1*Q1

{Specific Heats}

CpN1=(31.2+-0.0136 * (TaveN1)+0.0000268 * (TaveN1)^2+ -0.0000000117 * (TaveN1)^3) *(1/MWN)*(1000) {g-mol/kg-mol}

CpN2=(31.2+-0.0136*(TaveN2)+0.0000268*(TaveN2)^2+-0.0000000117*(TaveN2)^3) *(1/MWN)*(1000) {Reid 226}

TaveN1=(T2+T1)/2

TaveN2=(T2+T3)/2

CpL=(2.3388-0.0038942*TaveA+1.1221e-05*(TaveA)^2)*1000 {J/kg K} {correlation with experimental data}

CpA=(6.301+0.261*TaveA-0.000125*(TaveA)^2+0.0000000204*(TaveA)^3)*(1/MWA)*(1000)

TaveA=(T3+T2)/2

CpAN=((xA1+xA2)/2)*CpA+((xN1+xN2)/2)*CpN2

{Enthalpies -- Integrated form of Specific Heats}

EnthN1=(31.2*(T2-T1)+-0.0136/2 * (T2^2-T1^2)^2+0.0000268 /3* (T2^3-T1^3)+ -0.0000000117/4* (T2^4-T1^4)) *(1/MWN)*(1000) {g-mol/kg-mol}

EnthN2=(31.2*(T3-T2)+-0.0136/2*(T3^2-T2^2)+0.0000268/3*(T3^3-T2^3)+-0.0000000117/4* (T3^4-T2^4)) *(1/MWN)*(1000) {Reid 226}

EnthL=(2.3388*(T3-T2)-0.0038942/2*(T3^2-T2^2)+1.1221e-05/3*(T3^3-T2^3))*1000

{J/kg K} {correlation with experimental data}

EnthA=(6.301*(T3-T2)+0.261/2*(T3^2-T2^2)-0.000125/3*(T3^3-T2^3)+0.0000000204/4* (T3^4-T2^4))*(1/MWA)*(1000)

EnthAN=((xA1+xA2)/2)*EnthA+((xN1+xN2)/2)*EnthN2

{Parameters}

P=101325 {N/m^2}

R=8314 {N m/K kg-mol}

T3=294 {K}

T1=77 {K}

MWA=58.1 {kg/kg-mol}

MWN=28.02 {kg/kg-mol}

Q1=(.4{Lpm}/(1000*60)) {m³/sec}

hA=524000*((1-TaveA/508.1)/(1-329.2/508.1))^0.27 {J/Kg}{Reid pg.228}

hLN=199106 {J/Kg}

xA1=0.25