Skip to content

Commit 8f1b219

Browse files
dcordiercnrsdcordiercnrs
committed
PCSAFT starting point
0 parents  commit 8f1b219

27 files changed

+25731
-0
lines changed

‎BINARY_RESULTS_EXAMPLES/binary_diag_N2CH4_L_110K.dat

Lines changed: 3009 additions & 0 deletions
Large diffs are not rendered by default.

‎BINARY_RESULTS_EXAMPLES/binary_diag_N2CH4_L_115K.dat

Lines changed: 3009 additions & 0 deletions
Large diffs are not rendered by default.

‎BINARY_RESULTS_EXAMPLES/binary_diag_N2CH4_L_120K.dat

Lines changed: 3009 additions & 0 deletions
Large diffs are not rendered by default.

‎BINARY_RESULTS_EXAMPLES/binary_diag_N2CH4_V_110K.dat

Lines changed: 3009 additions & 0 deletions
Large diffs are not rendered by default.

‎BINARY_RESULTS_EXAMPLES/binary_diag_N2CH4_V_115K.dat

Lines changed: 3009 additions & 0 deletions
Large diffs are not rendered by default.

‎BINARY_RESULTS_EXAMPLES/binary_diag_N2CH4_V_120K.dat

Lines changed: 3009 additions & 0 deletions
Large diffs are not rendered by default.

‎DATABASE_PCSAFT/COMPOUNDS_DATA_PC-SAFT-PARAM_INTERAC.data

Lines changed: 175 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,175 @@
1+
# ----------------------------------------------------------------------------------------------------------------------------------
2+
# D. Cordier -- Institut UTINAM, Besançon, 13 novembre 2013
3+
# Database of interaction parameters "K_ij" used in PC-SAFT theory.
4+
# ----------------------------------------------------------------------------------------------------------------------------------
5+
#
6+
# ----------------------------------------------------------------------------------------------------------------------------------
7+
# On commence par les couples d'espces n'ayant pas un 'kij' qui dpend de la temprature (coefficient 'b' valant zro), on lit
8+
# la Table 3 de Tan et al. (2013) par colonne. Ces valeurs sont celles apparaissant en roman (ie par en italique) dans la table,
9+
# ceci ventuellement entre parenthses.
10+
# Comp. 1 | Comp. 2 | a | b | Comment
11+
N2 CH4 0.0300 0.0000 Tan et al. (2013)
12+
#N2 CH4 0.0000 0.0000 Test pour alcanofère
13+
# Maintenant ceux que j'ajoute pour calculer la dissolution des solides :
14+
N2 HCN 0.0000 0.0000 Moi
15+
N2 C4H10 0.0000 0.0000 Moi
16+
N2 CH3CN 0.0000 0.0000 Moi
17+
N2 CO2 0.0000 0.0000 Moi
18+
#N2 CO2 0.8000 0.0000 Moi
19+
#N2 CO2 0.6000 0.0000 Moi
20+
#N2 CO2 -0.0200 0.0000 Moi
21+
N2 C6H6 0.0000 0.0000 Moi
22+
#
23+
CO N2 0.0080 0.0000 Tan et al. (2013)
24+
#
25+
Ar CO 0.0000 0.0000 Tan et al. (2013)
26+
Ar H2 -0.1530 0.0000 Tan et al. (2013)
27+
#
28+
C2H6 N2 0.0400 0.0000 Tan et al. (2013)
29+
#C2H6 N2 0.0000 0.0000 Valeur de test pour l'influence de PC-SAFT sur le gradient de l'alcanofère
30+
#C2H6 N2 0.0400 0.0000 Tan et al. (2013)
31+
#C2H6 N2 0.0600 0.0000 Dt. le 13 mars 2017 partir des donnes de Malaska et al. (2017).
32+
#C2H6 N2 0.0700 0.0000 Valeur de travail pur TiLSS, 6 Nov. 2018, JPL.
33+
#
34+
#
35+
#C2H6 N2 0.0500 0.0000 Moi en fittant les donnes de Gabis (1991) fournies par Glein & Shock (2013)
36+
C2H6 CH4 0.0000 0.0000 Tan et al. (2013)
37+
C2H6 CO 0.0100 0.0000 Tan et al. (2013)
38+
# Maintenant ceux que j'ajoute pour calculer la dissolution des solides :
39+
C2H6 HCN 0.0000 0.0000 Moi
40+
C2H6 C4H10 0.0000 0.0000 Moi
41+
C2H6 CH3CN 0.0000 0.0000 Moi
42+
C2H6 CO2 0.1300 0.0000 Moi
43+
C2H6 C6H6 -0.1388 15.0699 Moi: ajuste pour fitter la valeur de Diez-y-Riega et al. (2014)
44+
#
45+
C2H2 N2 0.0800 0.0000 Tan et al. (2013)
46+
#C2H2 CH4 0.1000 0.0000 Tan et al. (2013)
47+
C2H2 CH4 0.1150 0.0000 Cordier, juillet 2014, par comparaison avec les data exp. de Neumann & Mann (1969)
48+
#C2H2 CH4 -0.0735 0.0000
49+
#C2H2 CH4 -0.0805 0.0000 Test pour Sunny (24/08/2015)
50+
C2H2 CO 0.0400 0.0000 Tan et al. (2013)
51+
C2H2 C2H6 0.1050 0.0000 Tan et al. (2013)
52+
#C2H2 C2H6 -0.1500 0.0000 Test pour Sunny (24/08/2015)
53+
#
54+
C3H8 N2 0.0730 0.0000 Tan et al. (2013)
55+
C3H8 CH4 0.0000 0.0000 Tan et al. (2013)
56+
C3H8 C2H6 0.0050 0.0000 Tan et al. (2013)
57+
C3H8 C2H2 0.1000 0.0000 Tan et al. (2013)
58+
#
59+
C2H4 N2 0.0620 0.0000 Tan et al. (2013)
60+
C2H4 CH4 0.0120 0.0000 Tan et al. (2013)
61+
#C2H4 CH4 0.0000 0.0000 Test pour Ellen Czaplinski (Arkansas), 27/01/2020
62+
#C2H4 CH4 -0.0340 0.0000 Test pour calculs pour Sunny (24/08/2015).
63+
C2H4 CO 0.0200 0.0000 Tan et al. (2013)
64+
#C2H4 Ar 0.0000 0.0000 Tan et al. (2013)
65+
#C2H4 Ar 0.1500 0.0000 pour le travail avec Carlos Manzanares
66+
C2H4 Ar 0.1970 0.0000 pour le travail avec Carlos Manzanares
67+
C2H4 C2H6 0.0100 0.0000 Tan et al. (2013)
68+
#C2H4 C2H6 0.0450 0.0000 Test pour Ellen Czaplinski (Arkansas), 27/01/2020
69+
#C2H4 C2H6 0.0120 0.0000 Test pour Sunny (24/08/2015).
70+
C2H4 C2H2 0.0450 0.0000 Tan et al. (2013)
71+
C2H4 C3H8 0.0200 0.0000 Tan et al. (2013)
72+
#
73+
# ----------------------------------------------------------------------------------------------------------------------------------
74+
# Paramtre d'interaction des soluts avec C3H8 (11/12/2014) :
75+
C3H8 HCN 0.0000 0.0000 Moi : supposition
76+
C3H8 C4H10 0.0000 0.0000 Moi : supposition
77+
C3H8 CH3CN 0.0000 0.0000 Moi : supposition
78+
C3H8 CO2 0.0000 0.0000 Moi : supposition
79+
C3H8 C6H6 0.0000 0.0000 Moi : supposition
80+
#
81+
# ----------------------------------------------------------------------------------------------------------------------------------
82+
# On passe aux couples d'espces ayant un 'kij' ayant une dpendance avec T, on lit toujours la Table 3 de Tan et al. (2013) par
83+
# colonne. Les "a" sont les valeurs au-dessus de la diagonnale, en italique et soulignes; les "b" sont en dessous de la diagonnale,
84+
# en italique.
85+
# Comp. 1 | Comp. 2 | a | b | Comment
86+
N2 H2 0.2894 -48.0300 Tan et al. (2013)
87+
N2 Ar 0.0152 -1.8100 Tan et al. (2013)
88+
#
89+
CH4 H2 0.1700 -24.6100 Tan et al. (2013)
90+
CH4 CO 0.0300 -1.2240 Tan et al. (2013)
91+
CH4 Ar 0.0527 -2.5000 Tan et al. (2013)
92+
# Maintenant ceux que j'ajoute pour calculer la dissolution des solides :
93+
CH4 HCN 0.0000 0.0000 Moi
94+
CH4 C4H10 0.0220 0.0000 Gross & Sadowski (2001)
95+
CH4 CH3CN 0.0000 0.0000 Moi
96+
#CH4 CO2 0.0650 0.0000 Gross & Sadowski (2001)
97+
CH4 CO2 0.0850 0.0000 Moi
98+
CH4 C6H6 0.0370 0.0000 Gross & Sadowski (2001)
99+
#CH4 C6H6 0.0000 0.0000 TEST
100+
#
101+
H2 CO 0.1066 -24.4400 Tan et al. (2013)
102+
H2 C2H6 0.1212 -9.8300 Tan et al. (2013)
103+
H2 C2H2 0.2070 -17.0500 Tan et al. (2013)
104+
H2 C3H8 0.0857 -5.5900 Tan et al. (2013)
105+
H2 C2H4 0.1641 -13.4400 Tan et al. (2013)
106+
#
107+
CO C3H8 0.0649 -0.9300 Tan et al. (2013)
108+
#
109+
Ar C2H6 0.0162 2.5500 Tan et al. (2013)
110+
Ar C2H2 -0.0160 -2.5500 Tan et al. (2013)
111+
Ar C3H8 0.0150 2.5200 Tan et al. (2013)
112+
# ----------------------------------------------------------------------------------------------------------------------------------
113+
# Données de Kr et Xe ajoutées le 18/11/2020 pour l'article sur l'alcanofère :
114+
Kr N2 0.0000 0.0000 DC 18/11/2020 - Dans l'ignorance du kij
115+
Kr CH4 0.0000 0.0000 DC 18/11/2020 - Dans l'ignorance du kij
116+
Kr C2H6 0.0000 0.0000 DC 18/11/2020 - Dans l'ignorance du kij
117+
Xe N2 0.0000 0.0000 DC 18/11/2020 - Dans l'ignorance du kij
118+
Xe CH4 0.0000 0.0000 DC 18/11/2020 - Dans l'ignorance du kij
119+
Xe C2H6 0.0000 0.0000 DC 18/11/2020 - Dans l'ignorance du kij
120+
#
121+
# ----------------------------------------------------------------------------------------------------------------------------------
122+
# On passe aux couples d'espces ayant un 'kij' dans le Maity (2003)
123+
# Comp. 1 | Comp. 2 | a | b | Comment
124+
n-octadecane n-heptane 0.0003 0.0000 Maity (2003) p. 38.
125+
n-dotriacontanen-heptane 0.0004 0.0000 Maity (2003) p. 39.
126+
biphenyl benzene 0.0000 0.0000 Maity (2003) p. 40.
127+
e-caprolactone toluene 0.0000 0.0000 Maity (2003) p. 41.
128+
#
129+
# **********************************************************************************************************************************
130+
# Maintenant ceux que j'ajoute pour calculer la dissolution des solides :
131+
HCN C4H10 0.0000 0.0000 Moi
132+
HCN C6H6 0.0000 0.0000 Moi
133+
HCN C2H2 0.0000 0.0000 Moi
134+
HCN CH3CN 0.0000 0.0000 Moi
135+
HCN CO2 0.0000 0.0000 Moi
136+
#
137+
C4H10 C2H2 0.0000 0.0000 Moi
138+
C4H10 CH3CN 0.0000 0.0000 Moi
139+
C4H10 CO2 0.0000 0.0000 Moi
140+
C4H10 C6H6 0.0000 0.0000 Moi
141+
#
142+
C2H2 CH3CN 0.0000 0.0000 Moi
143+
C2H2 CO2 0.0000 0.0000 Moi
144+
C2H2 C6H6 0.0000 0.0000 Moi
145+
#
146+
CH3CN CO2 0.0000 0.0000 Moi
147+
CH3CN C6H6 0.0000 0.0000 Moi
148+
#
149+
CO2 C6H6 0.0000 0.0000 Moi
150+
#
151+
naphthalene CH4 0.0000 0.0000 Moi
152+
naphthalene C2H6 -0.0494 0.0000 Moi
153+
naphthalene N2 0.0000 0.0000 Moi
154+
#
155+
biphenyl CH4 0.0000 0.0000 Moi
156+
biphenyl C2H6 -0.0030 0.0000 Moi
157+
biphenyl N2 0.0000 0.0000 Moi
158+
#
159+
# **********************************************************************************************************************************
160+
# Cas de l'eau :
161+
H2O CH4 0.0000 0.0000 Moi
162+
H2O C2H6 0.0000 0.0000 Moi
163+
H2O N2 0.0000 0.0000 Moi
164+
#
165+
# **********************************************************************************************************************************
166+
# Espces ajoues pour le travail sur l'absorption vers 2.4 microns dans l'atmosphre de Titan (10 mai 2016) :
167+
C6H12 N2 0.0000 0.0000 Moi
168+
n-pentane N2 0.0000 0.0000 Moi
169+
#
170+
# **********************************************************************************************************************************
171+
# Espces ajoues pour le travail avec Carlos Manzanares (14 juin 2016) :
172+
propyne Ar 0.1430 0.0000 Det. pour travail avec C. Manzanares 20/02/2017
173+
2methyl2butene Ar 0.1187 0.0000 Moi
174+
2m13butadiene Ar 0.0960 0.0000 Moi
175+

‎DATABASE_PCSAFT/COMPOUNDS_DATA_PC-SAFT.data

Lines changed: 83 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,83 @@
1+
# ----------------------------------------------------------------------------------------------------------------------------------
2+
# D. Cordier -- Institut UTINAM, Besançon, 1er novembre 2013
3+
# Database of parameters used in PC-SAFT theory.
4+
# ----------------------------------------------------------------------------------------------------------------------------------
5+
#
6+
# Comp | Molar mass | Num. seg. | Seg. diam. | Seg. engergy | Ref.
7+
# | (kg.mol^-1) | "m" | "sig" (Ang) | e/kB (K) |
8+
N2 28.01340E-03 1.2414 3.2992 89.2230 NIST Used by Tan et al. (2013)
9+
CH4 16.04250E-03 1.0000 3.7039 150.0300 NIST Used by Tan et al. (2013), repris sur G&S (2001) qui eux mêmes ...
10+
H2 2.01588E-03 1.0000 2.9860 19.2775 - Used by Tan et al. (2013)
11+
CO 28.01010E-03 1.3680 3.2035 89.4037 NIST Used by Tan et al. (2013)
12+
#CO 28.01010E-03 1.3097 3.2507 92.1500 Gross & Sadowski (2001) qui reprennent ... cf. leur papier.
13+
Ar 39.94800E-03 1.0000 3.3768 117.8538 NIST Used by Tan et al. (2013)
14+
#Ar 39.94800E-03 0.9200 3.4784 122.2300 Gross & Sadowski (2001) qui reprennent ... cf. leur papier.
15+
C2H6 30.06900E-03 1.6114 3.5245 190.9926 NIST Used by Tan et al. (2013)
16+
#C2H6 30.06900E-03 1.6069 3.5206 191.4200 Gross & Sadowski (2001) qui reprennent ... cf. leur papier.
17+
C2H2 26.03730E-03 2.1569 2.9064 168.5506 Din (1962) Used by Tan et al. (2013)
18+
C3H8 44.09560E-03 2.0022 3.6209 207.7721 NIST Used by Tan et al. (2013)
19+
C2H4 28.05320E-03 1.6018 3.4029 176.5698 NIST Used by Tan et al. (2013)
20+
#
21+
#HCN 27.02530E-03 1.3320 3.1890 259.1600 Formules p. 118 thèse Belkadi (2008)
22+
HCN 27.02530E-03 2.4340 3.2929 248.4800 Moi, Cordier et al. (2014 ou 2015)
23+
#
24+
#
25+
#C4H10 58.12220E-03 8.0000 3.5100 190.9000 TEST TEST
26+
#
27+
#
28+
C4H10 58.12220E-03 2.6300 3.5100 190.9000 Tab. III-1 p. 76 de Belkadi (2008), repris sur Tamouza et al. (2004)
29+
#C4H10 58.12220E-03 2.3316 3.7086 222.8800 Gross & Sadowski (2001) qui reprennent ... cf. leur papier.
30+
#CH3CN 41.05190E-03 1.4500 3.7000 268.0000 Belkadi (2008), Tab. V-4 p. 117
31+
CH3CN 41.05190E-03 2.2661 3.3587 313.0400 Spuhl et al. (2004)
32+
C6H6 78.11180E-03 2.4653 3.6478 287.3500 Maity (2003), Tab. 5.1 p. 35, qui reprend Gross & Sadowski (2001) qui reprennent eux mêmes ...
33+
CO2 44.00950E-03 2.0729 2.7852 169.2100 Gross & Sadowski (2001) qui reprennent ... cf. leur papier.
34+
#
35+
C6H12 84.14700E-03 2.5303 3.8499 278.1100 Gross & Sadowski (2001) qui reprennent ... cf. leur papier.
36+
# ----------------------------------------------------------------------------------------------------------------------------------
37+
# Dans la suite : Table 5.1 p. 35 de Maity (2003) en vue des tests de validation sur les données expérimentales utilisées par
38+
# Maity (2003).
39+
#
40+
# | (kg.mol^-1) | "m" | "sig" (Ang) | e/kB (K) |
41+
# ------- Solutes ---------
42+
n-dodecane 170.33480E-03 5.3060 3.8959 249.210 Maity (2003)
43+
n-hexadecane 226.44000E-03 6.6485 3.9592 254.700 Maity (2003)
44+
n-octadecane 254.49430E-03 7.3271 3.9668 256.200 Maity (2003)
45+
n-octacosane 394.76010E-03 10.3622 4.0217 252.000 Maity (2003)
46+
n-dotriacontane 450.86640E-03 11.8350 4.0217 252.000 Maity (2003)
47+
biphenyl 154.20780E-03 3.8877 3.8151 227.420 Maity (2003)
48+
e-caprolactone 114.14240E-03 2.8149 3.9902 255.050 Maity (2003)
49+
# ------- Solvents ---------
50+
n-pentane 72.14880E-03 2.6896 3.7729 231.200 Maity (2003)
51+
n-heptane 100.20190E-03 3.4831 3.8049 238.400 Maity (2003)
52+
n-decane 142.28170E-03 4.6627 3.8384 243.870 Maity (2003)
53+
# Rq. : il y a déjà le benzene au-dessus sous le nom C6H6 ---
54+
benzene 78.11180E-03 2.4653 3.6478 287.35 Maity (2003)
55+
toluene 92.13840E-03 2.8149 3.7169 285.69 Maity (2003)
56+
#
57+
# ----------------------------------------------------------------------------------------------------------------------------------
58+
# Collectées pour écrire des papiers :
59+
#
60+
# Comp | Molar mass | Num. seg. | Seg. diam. | Seg. engergy | Ref.
61+
# | (kg.mol^-1) | "m" | "sig" (Ang) | e/kB (K) |
62+
naphthalene 128.17000E-03 3.0915 3.8333 348.4000 Ting, P.D., PhD, Rice University (2003).
63+
#
64+
# ----------------------------------------------------------------------------------------------------------------------------------
65+
# Cas de l'eau
66+
#
67+
# Comp | Molar mass | Num. seg. | Seg. diam. | Seg. engergy | Asso? | Ass. energy | Vol. asso | Ref.
68+
# | (kg.mol^-1) | "m" | "sig" (Ang) | e/kB (K) | 0/1 | e_ass/kB (K) | kappa (AA^3) |
69+
H2O 18.01528E-03 2.6500 2.0700 151.3000 1 1515.0000 0.4320 Belkadi p. 81
70+
#
71+
# ----------------------------------------------------------------------------------------------------------------------------------
72+
# Pour faire des tests de calcul de Cp :
73+
# Comp | Molar mass | Num. seg. | Seg. diam. | Seg. engergy | Asso? | Ass. energy | Vol. asso | Ref.
74+
# | (kg.mol^-1) | "m" | "sig" (Ang) | e/kB (K) | 0/1 | e_ass/kB (K) | kappa (AA^3) |
75+
octane 114.23100E-03 3.8176 3.8373 242.7800 0 0.0000 0.0000 Gross & Sadowski (2001) Tab. 2 p. 1248
76+
#
77+
# ----------------------------------------------------------------------------------------------------------------------------------
78+
# 14 juin 2016 : estimations faites avec la méthode de "Group Contributions" de Tihic et al. (2008) pour les calculs faits pour Carlos Manzanares
79+
# Comp | Molar mass | Num. seg. | Seg. diam. | Seg. engergy | Asso? | Ass. energy | Vol. asso | Ref.
80+
# | (kg.mol^-1) | "m" | "sig" (Ang) | e/kB (K) | 0/1 | e_ass/kB (K) | kappa (AA^3) |
81+
propyne 40.06390E-03 1.8167 3.4710 229.6611 0 0.0000 0.0000 Mon estimation thèse Tihic (different du papier, plus de choses)
82+
2methyl2butene 70.13300E-03 2.4470 3.7800 243.8960 0 0.0000 0.0000 Mon estimation à partir de la GC de Tihic et al. (2008)
83+
2m13butadiene 68.11700E-03 2.4160 3.7470 242.1880 0 0.0000 0.0000 Mon estimation à partir de la GC de Tihic et al. (2008)

0 commit comments

Comments
 (0)