EES Program Listing and Output

Program Listing:

{adjustable parameters}
P=pressure(R600,T=Tp,x=0)
Tp=T1+adj
adj=.1
T1=315
T3=275
T11=375
{Qe1=1000}
m1=1

{Limiting Temperatures}
Tla=temperature(AMMONIA,P=P,x=0)
Tlb=temperature(R600,P=P,x=0)

{State 1:  Strong Butane}

CALL raoultab(T1a,P:xa1,xb1,xma1,xmb1,ya1,yb1,yma1,ymb1,Pa1,Pb1)
h1al=enthalpy(AMMONIA,T=T1,x=0)
h1bl=enthalpy(R600,T=T1,x=0)
u1al=intenergy(AMMONIA,T=T1,x=0)
u1bl=intenergy(R600,T=T1,x=0)
v1al=volume(AMMONIA,T=T1,x=0)
v1bl=volume(R600,T=T1,x=0)
s1al=entropy(AMMONIA,T=T1,x=0)
s1bl=entropy(R600,T=T1,x=0)
s1mixl=-xma1*Ra*ln(xa1)-xmb1*Rb*ln(xb1)
u1=xma1*u1al+xmb1*u1bl
v1=xma1*v1al+xmb1*v1bl
h1=u1+P*v1
s1=xma1*s1al+xmb1*s1bl+s1mixl

{State 2:  Strong Butane}
T2=T4
xma1=xma2
xmb1=xmb2
s1mixl=s2mixl
u2al=intenergy(AMMONIA,T=T2,x=0)
u2bl=intenergy(R600,T=T2,x=0)
v2al=volume(AMMONIA,T=T2,x=0)
v2bl=volume(R600,T=T2,x=0)
s2al=entropy(AMMONIA,T=T2,x=0)
s2bl=entropy(R600,T=T2,x=0)
v2=xma2*v2al+xmb2*v2bl
u2=xma2*u2al+xmb2*u2bl
h2=u2+P*v2
s2=xma2*s2al+xmb2*s2bl+s2mixl

{State 3:  Ammonia-Butane Mixture}
CALL raoultab(T3,P:xa3,xb3,xma3,xmb3,ya3,yb3,yma3,ymb3,Pa3,Pb3)
h3al=enthalpy(AMMONIA,T=T3,x=0)
h3bl=enthalpy(R600,T=T3,x=0)
u3al=intenergy(AMMONIA,T=T3,x=0)
u3bl=intenergy(R600,T=T3,x=0)
v3al=volume(AMMONIA,T=T3,x=0)
v3bl=volume(R600,T=T3,x=0)
s3al=entropy(AMMONIA,T=T3,x=0)
s3bl=entropy(R600,T=T3,x=0)
s3mixl=-xma3*Ra*ln(xa3)-xmb3*Rb*ln(xb3)
u3=xma3*u3al+xmb3*u3bl
v3=xma3*v3al+xmb3*v3bl
h3=u3+P*v3
s3=xma3*s3al+xmb3*s3bl+s3mixl

{State 4: Ammonia-Water}
CALL raoultaw(T4,P:xa4,xw4,xma4,xmw4,ya4,yw4,yma4,ymw4,Pa4,Pw4)
{liquid}
u4al=intenergy(AMMONIA,T=T4,x=0)
u4wl=intenergy(WATER,T=T4,x=0)
v4al=volume(AMMONIA,T=T4,x=0)
v4wl=volume(WATER,T=T4,x=0)
s4al=entropy(AMMONIA,T=T4,x=0)
s4wl=entropy(WATER,T=T4,x=0)
s4mixl=-xma4*Ra*ln(xa4)-xmw4*Rw*ln(xw4)
u4=xma4*u4al+xmw4*u4wl
v4=xma4*v4al+xmw4*v4wl
h4l=u4+P*v4
s4l=xma4*s4al+xmw4*s4wl+s4mixl

{vapor}
h4av=enthalpy(AMMONIA,T=T4,P=Pa4)
h4wv=enthalpy(WATER,T=T4,P=Pw4)
s4av=entropy(AMMONIA,T=T4,P=Pa4)
s4wv=entropy(WATER,T=T4,P=Pw4)
s4mixv=-yma4*Ra*ln(ya4)-ymw4*Rw*ln(yw4)
h4v=yma4*h4av+ymw4*h4wv
s4v=yma4*s4av+ymw4*s4wv

{Saturated V-L Mixture}
yma13=Q4*yma4+(1-Q4)*xma4
h4=Q4*h4v+(1-Q4)*h4l
s4=Q4*s4v+(1-Q4)*s4l
m4a=(Q4*yma4+(1-Q4)*xma4)*m13
m4w=(Q4*ymw4+(1-Q4)*xmw4)*m13
m4v=Q4*m13
m4l=(1-Q4)*m13
m4=m4v

{State 5:  100% Ammonia}
T5=T3
h5=enthalpy(AMMONIA,T=T5,P=P)
s5=entropy(AMMONIA,T=T5,P=P)

{State 6:  Ammonia-Butane mixture}
T6=T3
CALL raoultab(T6,P:xa6,xb6,xma6,xmb6,ya6,yb6,yma6,ymb6,Pa6,Pb6)
h6av=enthalpy(AMMONIA,T=T6,P=Pa6)
h6bv=enthalpy(R600,T=T6,P=Pb6)
u6av=intenergy(AMMONIA,T=T6,P=Pa6)
u6bv=intenergy(R600,T=T6,P=Pb6)
s6av=entropy(AMMONIA,T=T6,P=Pa6)
s6bv=entropy(R600,T=T6,P=Pb6)
s6mixv=-yma6*Ra*ln(ya6)-ymb6*Rb*ln(yb6)
h6=yma6*h6av+ymb6*h6bv
u6=yma6*u6av+ymb6*u6bv
s6=yma6*s6av+ymb6*s6bv

{State 7: Ammonia-Butane mixture}
CALL raoultab(T7,P:xa7,xb7,xma7,xmb7,ya7,yb7,yma7,ymb7,Pa7,Pb7)
h7av=enthalpy(AMMONIA,T=T7,P=Pa7)
h7bv=enthalpy(R600,T=T7,P=Pb7)
u7av=intenergy(AMMONIA,T=T7,P=Pa7)
u7bv=intenergy(R600,T=T7,P=Pb7)
s7av=entropy(AMMONIA,T=T7,P=Pa7)
s7bv=entropy(R600,T=T7,P=Pb7)
s7mixv=-yma7*Ra*ln(ya7)-ymb7*Rb*ln(yb7)
s7mixl=-xma7*Ra*ln(xa7)-xmb7*Rb*ln(xb7)
h7=yma7*h7av+ymb7*h7bv
u7=yma7*u7av+ymb7*u7bv
s7=yma7*s7av+ymb7*s7bv

{State 8: superheated ammonia-butane}
T8=T13
h8ash=enthalpy(AMMONIA,P=Pa7,T=T8)
s8ash=entropy(AMMONIA,P=Pa7,T=T8)
h8bsh=enthalpy(R600,P=Pb7,T=T8)
s8bsh=entropy(R600,P=Pb7,T=T8)
h8=yma7*h8ash+ymb7*h8bsh
s8=yma7*s8ash+ymb7*s8bsh

{The Evaporator}

{mass balances}
{upper}
m1+m6=m7+m3                             {overall}
m1*xmb1+m6*ymb6=m3*xmb3+m7*ymb7         {butane}
mchke1=m6*yma6-m7*yma7-m3*xma3-m1*xma1  {ammonia}

{lower}
m4+m3=m6                 {overall}
m4+m3*xma3=m6*yma6       {ammonia}
mchke2=m3*xmb3-m6*ymb6   {butane}

{heat transfer & entropy generation}
Qe1=m7*h7-m1*h2-m4*h4v
Qe2=m6*h6-m4*h5-m3*h3
Qe1=Qe2

{the pre-cooler/rectifier}
m1*h2+m7*h8+m13*h4=m1*h1+m7*h7+m13*h13

{State 13: strong ammonia-water sat. vapor}
T13=T9
CALL raoultaw(T13,P:xa13,xw13,xma13,xmw13,ya13,yw13,yma13,ymw13,Pa13,Pw13)
h13av=enthalpy(AMMONIA,T=T13,P=Pa13)
h13wv=enthalpy(WATER,T=T13,P=Pw13)
s13av=entropy(AMMONIA,T=T13,P=Pa13)
s13wv=entropy(WATER,T=T13,P=Pw13)
h13=yma13*h13av+ymw13*h13wv
s13=yma13*s13av+ymw13*s13wv

{State 9:  Ammonia-Water Strong liquid solution}
T9=T1
CALL raoultaw(T9,Pa7:xa9,xw9,xma9,xmw9,ya9,yw9,yma9,ymw9,Pa9,Pw9)
u9al=intenergy(AMMONIA,T=T9,x=0)
u9wl=intenergy(WATER,T=T9,x=0)
v9al=volume(AMMONIA,T=T9,x=0)
v9wl=volume(WATER,T=T9,x=0)
s9al=entropy(AMMONIA,T=T9,x=0)
s9wl=entropy(WATER,T=T9,x=0)
s9mixl=-xma9*Ra*ln(xa9)-xmw9*Rw*ln(xw9)
u9=xma9*u9al+xmw9*u9wl
v9=xma9*v9al+xmw9*v9wl
h9=u9+Pa7*v9
s9=xma9*s9al+xmw9*s9wl+s9mixl

{State 9i (internal): Liquid}
CALL raoultaw(T9i,P:xa9i,xw9i,xma9i,xmw9i,ya9i,yw9i,yma9i,ymw9i,Pa9i,Pw9i)
u9ial=intenergy(AMMONIA,T=T9i,x=0)
u9iwl=intenergy(WATER,T=T9i,x=0)
v9ial=volume(AMMONIA,T=T9i,x=0)
v9iwl=volume(WATER,T=T9i,x=0)
s9ial=entropy(AMMONIA,T=T9i,x=0)
s9iwl=entropy(WATER,T=T9i,x=0)
s9imixl=-xma9i*Ra*ln(xa9i)-xmw9i*Rw*ln(xw9i)
u9i=xma9i*u9ial+xmw9i*u9iwl
v9i=xma9i*v9ial+xmw9i*v9iwl
h9il=u9i+P*v9i
s9il=xma9i*s9ial+xmw9i*s9iwl+s9imixl

{State 9:  Ammonia-Water Strong vapor solution}
h9av=enthalpy(AMMONIA,T=T9,P=Pa9)
h9wv=enthalpy(WATER,T=T9,P=Pw9)
s9av=entropy(AMMONIA,T=T9,P=Pa9)
s9wv=entropy(WATER,T=T9,P=Pw9)
h9v=yma9*h9av+ymw9*h9wv
s9v=yma9*s9av+ymw9*s9wv

{State 9i (internal): Vapor}

h9iav=enthalpy(AMMONIA,T=T9i,P=Pa9i)
h9iwv=enthalpy(WATER,T=T9i,P=Pw9i)
s9iav=entropy(AMMONIA,T=T9i,P=Pa9i)
s9iwv=entropy(WATER,T=T9i,P=Pw9i)
h9iv=yma9i*h9iav+ymw9i*h9iwv
s9iv=yma9i*s9iav+ymw9i*s9iwv

{State 11:  Ammonia-Water Weak liquid solution}

CALL raoultaw(T11,P:xa11,xw11,xma11,xmw11,ya11,yw11,yma11,ymw11,Pa11,Pw11)
u11al=intenergy(AMMONIA,T=T11,x=0)
u11wl=intenergy(WATER,T=T11,x=0)
v11al=volume(AMMONIA,T=T11,x=0)
v11wl=volume(WATER,T=T11,x=0)
s11al=entropy(AMMONIA,T=T11,x=0)
s11wl=entropy(WATER,T=T11,x=0)
s11mixl=-xma11*Ra*ln(xa11)-xmw11*Rw*ln(xw11)
u11=xma11*u11al+xmw11*u11wl
v11=xma11*v11al+xmw11*v11wl
h11=u11+P*v11
s11=xma11*s11al+xmw11*s11wl+s11mixl

{State 12:  Ammonia-Water Weak sub-cooled liquid solution}
T12=T9
xma11=xma12
xmw11=xmw12
s11mixl=s12mixl
u12al=intenergy(AMMONIA,T=T12,x=0)
u12wl=intenergy(WATER,T=T12,x=0)
v12al=volume(AMMONIA,T=T12,x=0)
v12wl=volume(WATER,T=T12,x=0)
s12al=entropy(AMMONIA,T=T12,x=0)
s12wl=entropy(WATER,T=T12,x=0)
v12=xma12*v12al+xmw12*v12wl
u12=xma12*u12al+xmw12*u12wl
h12=u12+P*v12
s12=xma12*s12al+xmw12*s12wl+s12mixl

{SHX}
Qshx=m11*(h11-h12)
Qshx2=m9v*h9v+m9il*h9il-m9*h9-m9iv*h9iv
m9v+m9il=m9+m9iv
m9v*yma9+m9il*xma9i=m9*xma9+m9iv*yma9i
m9v=m13
m9v*h9v+m9il*h9il+m11*h12=m9*h9+m9iv*h9iv+m11*h11

{Generator}
{mass conservation}

m9=m13+m11                           {overall}
xma9*m9=m13*yma13+xma11*m11          {ammonia}
mchkg1=xmw9*m9-xmw11*m11-m13*ymw13   {water}

{Energy & Entropy}
Tg=(T11+T9i)/2
Qg=m13*h13+m11*h12-m9*h9
COP=Qe1/Qg

{Condenser}
CALL raoultab(T9,P:xac,xbc,xmac,xmbc,yac,ybc,ymac,ymbc,Pac,Pbc)
xmb1=ymbc
Qc=m9*h9+m1*h1-m11*h12-m7*h8-m4l*h4l
Qchk=Qc+Qe1+Qg
mchkc1=m11+m4l+m7-m9-m1
mchkc2=m4l*xma4+m11*xma11+m7*yma7-m9*xma9-m1*xma1
Tc=T9

{2nd Law}
COPideal=((Tg-Tc)/(Tg))/((Tc-T3)/T3)
Sgen=-(Qg/Tg)-(Qe1/T3)-(Qc/Tc)
COPact=COPideal-((T3*Tc)/(Tc-T3))*Sgen/Qg
 

Output: