550

F-PROPATH: Binary Mixture by General Equations

2.2 CSD Equation
Thermodynamic properties are evaluated by Carnahan–Starling–DeSantis equation[1]. No provision is madefor transport properties at this moment.

2.2.1 The Names of Subroutine, Library File and Single Shot Program
Library File for UNIX:
Library File for DOS,Windows95/NT:
Single Shot Program for UNIX:
Single Shot Program for DOS,Windows95/NT:

libjcsdmx.a
JCSDMX.LIB
csdmx--ss
CSDMX–SS.EXE

2.2.2 Reference State
By user’s selecton of KSTAN, which is the ﬁrst ARGUMENT of SUBROUTINE STNKAS, the reference stateis determined as shown below.
Table V-2.2-1 Reference State Speciﬁed by KSTAN
KSTAN
0
1
othersReference Statevalue of enthalpy and entropy
298.15K, 1bar : ideal gas state of pure component
0.0kJ/kmol, 0.0kJ/(kmol·K)
273.15K : saturated liquid of pure component
200kJ/kg, 1.0kJ/(kg·K)
same with KSTAN=0

2.2.3 Missing SUBPROGRAM
Of the SUBPROGRAMs listed in 1.6, the next is missing in F–PROPATH by CSD equation.
CRPM, SUBCRT, SUBPAR, PBT, PDT, QMIX, TBP, TDP

2.2.4 Initialization of F–PROPATH by CSD equation
You can initialize F–PROPATH by writing three lines shown below. See 1.1 and Sample Program 1 of 2.1.6.

CALL KPAMES(KPA, MESS)CALL STNKAS(KSTAN, KAS)CALL START1(J,KOMBI)ARGUMENTs to the ﬁrst two SUBROUTINEs are explained in 1.3. The last, START1, works as follows.
START1SUBROUTINE START1(J,KOMBI)
input
INTEGERKOMBI : user selects one from the table below.

outputINTEGERJ−2: error detection code
: invalid ARGUMENT(s)Table V-2.2-2 Combination of Substances Speciﬁed by KOMBI
KOMBI
1
2
3
4
5

component 1
R22
R13B1
R22
R23
R13

component 2
R123
R152a
R12
R13
R12

KOMBI
6
7
8
9
10

component 1 component 2
R12
R152a
R22
R114
R23
R12
R22
R11
not deﬁned

CSD Equation

551

2.2.5 Example of User’s PROGRAM UNIT
Sample Program 1
EX1.for :
c
c
c
c
c
c

-----------------------------------------------CSD Equation
This program calculates the VLE properties of
R22-R123 mixture at p=1.0[bar],T=150[K]-----------------------------------------------REAL T,P,X,Y,VL,VV,HL,HV,SL,SVINTEGER J

c
c
c
c
c
c

c
c
c
c
c

c
c
c
c
c
c
c
c
c
c

-------------------------------------------------T=[K],p=[Pa]standard state=(bibliography [1] of chapter 10)
amount of substance=[kmol]-------------------------------------------------CALL KPAMES(0,1)
CALL STNKAS(0,0)
CALL START1(J,1)------------------------------------------------Calculation of property on saturated state andoutput results.
------------------------------------------------P=1.0E5
T=290.0
CALL SUBXY(J,T,P,X,Y,VL,VV,HL,HV,SL,SV)
IF (J.NE.0) THEN
WRITE(*,*) ’ ERROR IN SUBXY J=’,J
STOP
ENDIF
HL=HL*1.0E-6
HV=HV*1.0E-6
SL=SL*1.0E-3
SV=SV*1.0E-3
WRITE(*,*) ’ X[kmol/kmol]=’,X, ’
y[kmol/kmol]
WRITE(*,*) ’ vL[m^3/kmol]=’,VL,’ vV[m^3/kmol]
WRITE(*,*) ’ hL[MJ/kmol]=’,HL,’ hV[MJ/kmol]
WRITE(*,*) ’ sL[MJ/(kmol K)] =’,SL,’ sV[MJ/(kmol K)]
STOP
END

=’,Y
=’,VV
=’,HV
=’,SV

------------------------------------------------------------------Results------------------------------------------------------------------X[kmol/kmol]vL[m^3/kmol]hL[MJ/kmol]sL[MJ/(kmol K)]

=
=
=
=

4.420045E-02
1.019272E-01-31.306840-58.291830

y[kmol/kmol]vV[m^3/kmol]hV[MJ/kmol]sV[MJ/(kmol K)]

=
=
=
=

3.520642E-01
23.330370-9.433830E-01
2.713856

2.2.6 Sample Output of Single Shot Program========================================================================|
Single Shot Program for Peng-Robinson Equation,||
an Application Program by F-PROPATH Ver.9.1||
Thermophysical Properties of Binary Mixtures|

552

F-PROPATH: Binary Mixture by General Equations|----------------------------------------------------------------------|| ## First Menu ##|| No.
Current|| 1 --> Go to Second Menu|
| 2 --> Select Mixture(Component1 - Component2)||[ R22 - R152a ]|| 3 --> Read Help|
| 4 --> Set Logfile[OFF]|| 0 --> Quit|========================================================================
Input No.
3
This is help for this program (CSDMX-SS.EXE)
Nomenclature==== First Charactor ====
H : Enthalpy
P : Pressure
S : Entropy
T : Temperature
V : Volume
X : Component1 Composition of Liquid
Y : Component1 Composition of Vapor
Z : Total Composition of Component1==== Second Charactor ====
B : Bubble Point
D : Dew Point
L : Liquid
V : Vapor----- Hit RETURN Key ----========================================================================|
Single Shot Program for Peng-Robinson Equation,||
an Application Program by F-PROPATH Ver.9.1||
Thermophysical Properties of Binary Mixtures||----------------------------------------------------------------------|| ## First Menu ##|| No.
Current|| 1 --> Go to Second Menu|
| 2 --> Select Mixture(Component1 - Component2)||[ R22 - R152a ]|| 3 --> Read Help|
| 4 --> Set Logfile[OFF]|| 0 --> Quit|========================================================================
Input No.
4
1: Logfile ON
2: Logfile OFF
Input No.
1
Input Filename
csdmix.log========================================================================|
Single Shot Program for Peng-Robinson Equation,||
an Application Program by F-PROPATH Ver.9.1||
Thermophysical Properties of Binary Mixtures||----------------------------------------------------------------------|| ## First Menu ##|| No.
Current|| 1 --> Go to Second Menu|
| 2 --> Select Mixture(Component1 - Component2)||[ R22 - R152a ]|| 3 --> Read Help|
| 4 --> Set Logfile[ON] csdmix.log|| 0 --> Quit|========================================================================
Input No.
2

CSD Equation=====================================| No.| Component1 - Component2||----|------------------------------|| 1 | [ R22 - R123 ]|| 2 | [ R13 - R152a ]|| 3 | [ R22 - R12 ]|| 4 | [ R23 - R13 ]|| 5 | [ R13 - R12 ]|| 6 | [ R12 - R152a ]|| 7 | [ R22 - R114 ]|| 8 | [ R23 - R12 ]|| 9 | [ R22 - R11 ]|=====================================
Input No.
1========================================================================|
Single Shot Program for Peng-Robinson Equation,||
an Application Program by F-PROPATH Ver.9.1||
Thermophysical Properties of Binary Mixtures||----------------------------------------------------------------------|| ## First Menu ##|| No.
Current|| 1 --> Go to Second Menu|
| 2 --> Select Mixture(Component1 - Component2)||[ R22 - R123 ]|| 3 --> Read Help|
| 4 --> Set Logfile[ON] csdmix.log|| 0 --> Quit|========================================================================
Input No.
1===============================================================| ## Second Menu ##|| <Single and Two Phase States>|| No. Input
Output|| 1 --(P,T,Z)-->H HL HV S SL SV V VL VV X Y QUALITY || 2 --(P,Z,H)--> T
HL HV S SL SV V VL VV X Y QUALITY || 3 --(P,Z,S)--> T H HL HVSL SV V VL VV X Y QUALITY || 4 --(P,Z,V)--> T H HL HV S SL SV
VL VV X Y QUALITY ||-------------------------------------------------------------|| <Two Phase State>|| 5 --(P,Z)----> X=Z : TB HB SB VB
Y=Z : TD HD SD VD|| 6 --(T,Z)----> X=Z : PB HB SB VB
Y=Z : PD HD SD VD|| 7 --(P,T)----> Coexisting Phases: X Y HL HV SL SV VL VV||-------------------------------------------------------------|| 8 -----------> Fundamental Constants|| 9 -----------> Conversion of Composition and Relative||
Molcular Mass of Mixture||10 -----------> Change System of Unit, and Standard Values of||Enthalpy and Entropy||99 -----------> Return to First Menu|
| 0 -----------> Quit|===============================================================
Input No.
10=======================================================
No.
current=======================================================
1 --> Pressure and Temperature[bar, C]2 --> Amount of Substance[ kg ]
3 --> Standard Values of h and s [2] IIR
0 --> Return to Second Menu=======================================================
Input No.
1===============================
No. <Pressure> <Temperature>===============================
1 -->
Pa
K

553

554

F-PROPATH: Binary Mixture by General Equations

2 -->
bar
C
3 -->
bar
K
4 -->
Pa
C
===============================
Input No.
1=======================================================
No.
current=======================================================
1 --> Pressure and Temperature[ Pa, K]2 --> Amount of Substance[ kg ]
3 --> Standard Values of h and s [2] IIR
0 --> Return to Second Menu=======================================================
Input No.
2===========================No.
Amount of Substance
1 -->
kmol
2 -->
kg
===========================
Input No.
1=======================================================
No.
current=======================================================
1 --> Pressure and Temperature[ Pa, K]2 --> Amount of Substance[kmol]
3 --> Standard Values of h and s [2] IIR
0 --> Return to Second Menu=======================================================
Input No.
3================================================
No.
Standard States and Values of h and s
1 --> Ideal Gas State of Pure Components
h=0.0[kJ/kg], s=0.0[kJ/(kg*K)]
Pure Ideal Gases at 0[C], 1.0[bar]2 --> Convention of International Institute ofRefrigeration(IIR):
h=200[kJ/kg], s=1.0[kJ/(kg*K)]
Saturated Pure Liquids at 0[C]================================================
Input No.
1=======================================================
No.
current=======================================================
1 --> Pressure and Temperature[ Pa, K]2 --> Amount of Substance[kmol]
3 --> Standard Values of h and s [1] Ideal Gas State0 --> Return to Second Menu=======================================================
Input No.
0===============================================================| ## Second Menu ##|contents not shown===============================================================
Input No.
9==================================================================
Conversion of Composition and Relative Molecular Mass of Mixture==================================================================

CSD Equation555

1: [kmol/kmol] --> [kg/kg]
and Relative Molecular Mass of Mixture
2: [kg/kg]--> [kmol/kmol] and Relative Molecular Mass of Mixture
0: Return to Second Menu
Input No.
1
Input Z[kmol/kmol]
0.5
Z =.50000000[kmol/kmol]
Z =.36118780[kg/kg]
Relative Molecular Mass of Mixture =119.69950000[kg/kmol]
Relative Molecular Mass of [ R22 ]86.46800000[kg/kmol]
Relative Molecular Mass of [ R123 ]152.93100000[kg/kmol]

1: [kmol/kmol] --> [kg/kg]
and Relative Molecular Mass of Mixture
2: [kg/kg]--> [kmol/kmol] and Relative Molecular Mass of Mixture
0: Return to Second Menu
Input No.
2
Input Z[kg/kg]
0.5
Z =.50000000[kg/kg]
Z =.63881220[kmol/kmol]
Relative Molecular Mass of Mixture =110.47360000[kg/kmol]
Relative Molecular Mass of [ R22 ]86.46800000[kg/kmol]
Relative Molecular Mass of [ R123 ]152.93100000[kg/kmol]
1: [kmol/kmol] --> [kg/kg]
and Relative Molecular Mass of Mixture
2: [kg/kg]--> [kmol/kmol] and Relative Molecular Mass of Mixture
0: Return to Second Menu
Input No.
0===============================================================| ## Second Menu ##|contents not shown===============================================================
Input No.
8=======================
Fundamental Constants=======================<Constants>
Relative Molecular Mass
Gas Constants
Critical Temperature
Critical Pressure
Critical Volume[ R22 ]
86.46800000
96.15707000
369.30000000
96.15001000
49.71000000.16861500E+00.19500280E-02[ R123 ]152.93100000
54.36772000
456.86000000
183.71000000
36.66000000.27527600E+00.18000010E-02----- Hit RETURN Key ----===============================================================| ## Second Menu ##|contents not shown===============================================================
Input No.
1======================================(P,T,Z) --> H, S, V, QUALITY
X, HL, SL, VL
Y, HV, SV, VV======================================
Input P[Pa](P=0 : Second Menu)
100000.0[kg/kmol][J/(kg*K)]
[K]
[C]
[bar]
[m**3/kmol][m**3/kg]

556

Input T[K]280.0
Input Z[kmol/kmol]
0.5
T =280.00000000[K]P =Z =.50000000[kmol/kmol] H =S =-10.61071000[kJ/(kmol*K)]V =X =.09671523[kmol/kmol] Y =HL = -32138.17000000[kJ/kmol]
HV =
SL =
-61.69294000[kJ/(kmol*K)]SV =VL =.98700080E-01[m**3/kmol] VV =QUALITY =.81843440[kmol/kmol]
Input P[Pa](P=0 : Second Menu)
100000.0
Input T[K]300.0
Input Z[kmol/kmol]
0.5
T =300.00000000[K]P =Z =.50000000[kmol/kmol] H =S =5.99217600[kJ/(kmol*K)]V =Vapor Region

F-PROPATH: Binary Mixture by General Equations

100.00000000[kPa]-7103.10700000[kJ/kmol].18527250E+02[m**3/kmol].58946680[kmol/kmol]-1549.20300000[kJ/kmol]
721.63350000[J/(kmol*K)].22615530E+02[m**3/kmol]

100.00000000[kPa]-57.29559000[kJ/kmol].24316680E+02[m**3/kmol]

Input P[Pa](P=0 : Second Menu)
0===============================================================| ## Second Menu ##|contents not shown===============================================================
Input No.
2======================================(P,Z,H) --> T, S, V, QUALITY
X, HL, SL, VL
Y, HV, SV, VV======================================
Input P[Pa](P=0 : Second Menu)
100000.0
Input Z[kmol/kmol]
0.5
Input H[J/kmol]-57295.59
T =300.00000000[K]P =Z =.50000000[kmol/kmol] H =S =5.99217600[kJ/(kmol*K)]V =Vapor Region

100.00000000[kPa]-57.29559000[kJ/kmol].24316680E+02[m**3/kmol]

Input P[Pa](P=0 : Second Menu)
0===============================================================| ## Second Menu ##|contents not shown===============================================================
Input No.
3======================================(P,Z,S) --> T, H, V, QUALITY
X, HL, SL, VL
Y, HV, SV, VV======================================
Input P[Pa](P=0 : Second Menu)
100000.0
Input Z[kmol/kmol]
0.5
Input S[J/(kmol*K)]
5992.176
T =300.00000000[K]P =Z =.50000000[kmol/kmol] H =100.00000000[kPa]-57.29567000[kJ/kmol]

CSD Equation

S =5.99217600[kJ/(kmol*K)]V
Vapor Region

557=.24316680E+02[m**3/kmol]

Input P[Pa](P=0 : Second Menu)
0===============================================================| ## Second Menu ##|contents not shown===============================================================
Input No.
4========================================(P,Z,V) --> T, H, S, QUALITY
X, HL, SL, VL
Y, HV, SV, VV========================================
Input P[Pa](P=0 : Second Menu)
100000.0
Input Z[kmol/kmol]
0.5
Input V[m**3/kmol]
24.31668
T =300.00010000[K]P =Z =.50000000[kmol/kmol] H =S =5.99218900[kJ/(kmol*K)]V =Vapor Region

100.00000000[kPa]-57.29181000[kJ/kmol].24316680E+02[m**3/kmol]

Input P[Pa](P=0 : Second Menu)
0===============================================================| ## Second Menu ##|contents not shown===============================================================
Input No.
5===============================================(P,Z) --> X=Z : TB, HB, SB, VB
Y=Z : TD, HD, SD, VD===============================================
Input P[Pa](P=0 : Second Menu)
100000.0
Input Z[kmol/kmol]
0.5
P =100.00000000[kPa]Z =.50000000[kmol/kmol]
TB =
246.39250000[K]TD =284.17860000[K]HB = -31560.63000000[kJ/kmol]HD =-1315.70200000[kJ/kmol]
SB =-67.17770000[kJ/(kmol*K)]SD =1.72309500[kJ/(kmol*K)]
VB =.80985700E-01[m**3/kmol] VD =.22921670E+02[m**3/kmol]
Input P[Pa](P=0 : Second Menu)
0===============================================================| ## Second Menu ##|contents not shown===============================================================
Input No.
6===============================================(T,Z) --> X=Z : PB, HB, SB, VB
Y=Z : PD, HD, SD, VD===============================================
Input Z[kmol/kmol](Z < 0 : Second Menu)
0.5
Input T[K]273.15T =273.15000000[K]Z =.50000000[kmol/kmol]
PB =263.49430000[kPa]PD =62.22648000[kPa]
HB = -28486.09000000[kJ/kmol]
HD =-2080.65200000[kJ/kmol]

558

F-PROPATH: Binary Mixture by General Equations

SB =
VB =

-53.59378000[kJ/(kmol*K)]SD =.85064980E-01[m**3/kmol] VD =-1.19924600[kJ/(kmol*K)].35736690E+02[m**3/kmol]

Input Z[kmol/kmol](Z < 0 : Second Menu)-1===============================================================| ## Second Menu ##|contents not shown===============================================================
Input No.
7==============================================(P,T) --> Coexisting Phases
X, HL, SL, VL
Y, HV, SV, VV==============================================
Input P[Pa](P=0 : Second Menu)
100000.0
Input T[K]273.15T =273.15000000[K]P =100.00000000[kPa]X =.14236480[kmol/kmol] Y =.70970110[kmol/kmol]
HL = -32552.08000000[kJ/kmol]
HV =-1889.64500000[kJ/kmol]
SL =-63.78323000[kJ/(kmol*K)]SV =-1.18756400[kJ/(kmol*K)]
VL =.96228840E-01[m**3/kmol] VV =.22094880E+02[m**3/kmol]
Input P[Pa](P=0 : Second Menu)
0===============================================================| ## Second Menu ##|contents not shown===============================================================
Input No.
0***** See You Again ! *****
Stop - Program terminated.

References[1] G. Morrison, M. O. McLinden, Application of a hard sphere equation of state to refrigerants and refrigerant
mixtures, NBS Technical Note 1226, 1986