CoolProp  4.2.5
An open-source fluid property and humid air property database
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Brine.cpp
Go to the documentation of this file.
1 
2 /*
3 Brine data, by Ian Bell, 2012
4 
5 Coefficients taken from
6 
7 Ake Melinder, "Properties of Secondary Working Fluids for Indirect Systems", 2010. Published by IIR.
8 
9 Erratum:
10  --> In the tables of coefficients, the range for y should always be between tF(x) < y < ymax
11  --> There is a typo in the specific heat for Ethylene Glycol, coefficients give a result in kJ/kg
12  --> Watch out for some values with commas that should be decimals
13 
14 */
15 #if defined(_MSC_VER)
16 #define _CRTDBG_MAP_ALLOC
17 #define _CRT_SECURE_NO_WARNINGS
18 #include <stdlib.h>
19 #include <crtdbg.h>
20 #else
21 #include <stdlib.h>
22 #endif
23 
24 #include "math.h"
25 #include "stdio.h"
26 #include <string.h>
27 #include "Brine.h"
28 #include "CoolProp.h"
29 
30 // From Melinder, 2010
31 static double PG[18][5]={
32  {-13.25,1018,3882,0.4513,0.6837},
33  {-0.0000382,-0.5406,2.699,0.0007955,-0.03045},
34  {0.0000007865,-0.002666,-0.001659,0.00000003482,0.0002525},
35  {-0.000000001733,0.00001347,-0.00001032,-0.000000005966,-0.000001399},
36  {-0.6631,0.7604,-13.04,-0.004795,0.03328},
37  {0.000006774,-0.00945,0.0507,-0.00001678,-0.0003984},
38  {-0.00000006242,0.00005541,-0.00004752,0.00000008941,0.000004332},
39  {-0.0000000007819,-0.0000001343,0.000001522,0.0000000001493,-0.0000000186},
40  {-0.01094,-0.002498,-0.1598,0.00002076,0.00005453},
41  {0.00000005332,0.000027,0.00009534,0.0000001563,-0.000000086},
42  {-0.000000004169,-0.0000004018,0.00001167,-0.000000004615,-0.00000001593},
43  {0.00000000003288,0.000000003376,-0.0000000487,0.000000000009897,-0.00000000004465},
44  {-0.0002283,-0.000155,0.0003539,-0.00000009083,-0.0000039},
45  {-0.00000001131,0.000002829,0.00003102,-0.000000002518,0.0000001054},
46  {0.0000000001918,-0.000000007175,-0.000000295,0.00000000006543,-0.000000001589},
47  {-0.000003409,-0.000001131,0.00005,-0.0000000005952,-0.00000001587},
48  {0.00000000008035,-0.00000002221,-0.0000007135,-0.00000000003605,0.0000000004475},
49  {0.00000001465,0.00000002342,-0.0000004959,0.00000000002104,0.000000003564}
50 };
51 
52 // From Melinder 2010
53 static double EG[18][5]={
54  {-15.25,1034,3.737,0.472,0.4705},
55  {-0.000001566,-0.4781,0.00293,0.0008903,-0.0255},
56  {-0.0000002278,-0.002692,-0.000004675,-0.000001058,0.0001782},
57  {0.000000002169,0.000004725,-0.00000001389,-0.000000002789,-0.0000007669},
58  {-0.808,1.311,-0.01799,-0.004286,0.02471},
59  {-0.000001339,-0.006876,0.0001046,-0.00001473,-0.0001171},
60  {0.00000002047,0.00004805,-0.0000004147,0.0000001059,0.000001052},
61  {-0.00000000002717,0.0000000169,0.0000000001847,-0.0000000001142,-0.00000001634},
62  {-0.01334,0.0000749,-0.00009933,0.00001747,0.000003328},
63  {0.00000006322,0.00007855,0.0000003516,0.00000006814,0.000001086},
64  {0.0000000002373,-0.0000003995,0.000000005109,-0.000000003612,0.00000001051},
65  {-0.000000000002183,0.000000004982,-0.00000000007138,0.000000000002365,-0.0000000006475},
66  {-0.00007293,-0.0001062,0.00000261,0.00000003017,0.000001659},
67  {0.000000001764,0.000001229,-0.000000001189,-0.000000002412,0.000000003157},
68  {-0.00000000002442,-0.00000001153,-0.0000000001643,0.00000000004004,0.0000000004063},
69  {0.000001006,-0.0000009623,0.00000001537,-0.000000001322,0.00000003089},
70  {-0.00000000007662,-0.00000007211,-0.0000000004272,0.00000000002555,0.0000000001831},
71  {0.00000000114,0.00000004891,-0.000000001618,0.00000000002678,-0.000000001865}
72 };
73 
74 // From Melinder 2010
75 static double EA[18][5]={
76  {-19.41,961.9,4204,0.4067,1.474},
77  {-0.0003668,-0.5222,2.319,0.0006775,-0.04745},
78  {-0.00004005,-0.003281,-0.03042,0.0000003105,0.0004314},
79  {0.000001524,0.00001569,0.000686,-0.00000002,-0.000003023},
80  {-0.954,-1.433,-21.02,-0.005008,0.01565},
81  {-0.00001209,-0.01989,0.4927,-0.00002377,-0.00004106},
82  {0.000002877,0.000187,-0.003072,-0.00000003216,-0.000005135},
83  {-0.00000004394,-0.0000009154,-0.0000569,0.00000000008362,0.00000007004},
84  {-0.002648,-0.0226,-0.3714,0.00002801,-0.0008435},
85  {-0.0000003173,0.0002281,-0.002335,0.0000002669,0.0000164},
86  {0.000000008652,-0.00000008581,-0.0000196,-0.000000003606,-0.0000001091},
87  {-0.0000000003717,0.000000004056,0.0000007461,0.00000000001552,-0.000000001967},
88  {0.0003851,-0.000169,0.01743,-0.00000002009,0.000007552},
89  {0.0000000134,0.000008594,-0.0002969,-0.000000006813,-0.0000001118},
90  {-0.000000002091,-0.00000009607,0.000001901,0.0000000001429,0.000000001899},
91  {-0.0000002858,0.00001291,-0.00006292,-0.000000001506,0.0000001529},
92  {0.0000000009312,-0.000000159,0.000005353,0.0000000001167,-0.0000000009481},
93  {-0.000000167,-0.00000008318,-0.00000829,-0.00000000001653,-0.00000000413}
94 };
95 
96 // From Melinder 2010 (Methanol = Methyl Alcohol = MA)
97 static double MA[18][5]={
98  {-26.29,958.1,3887,0.4175,1.153},
99  {-0.000002575,-0.4151,7.201,0.0007271,-0.03866},
100  {-0.000006732,-0.002261,-0.08979,0.0000002823,0.0002779},
101  {0.000000163,0.0000002998,-0.000439,0.000000009718,-0.000001543},
102  {-1.187,-1.391,-18.5,-0.004421,0.005448},
103  {-0.00001609,-0.0151,0.2984,-0.00002952,0.0001008},
104  {0.000000342,0.0001113,-0.001865,0.00000007336,-0.000002809},
105  {0.0000000005687,-0.0000003264,-0.00001718,0.0000000004328,0.000000009811},
106  {-0.01218,-0.01105,-0.03769,0.00002044,-0.0005552},
107  {0.0000003865,0.0001828,-0.01196,0.0000003413,0.000008384},
108  {0.000000008768,-0.000001641,0.00009801,-0.000000003665,-0.00000003997},
109  {-0.0000000002095,0.0000000151,0.000000666,-0.00000000002791,-0.0000000003466},
110  {-0.00006823,-0.0001208,-0.003776,0.0000002943,0.000003038},
111  {0.00000002137,0.000002992,-0.00005611,-0.0000000009646,-0.00000007435},
112  {-0.0000000004271,0.000000001455,-0.0000007811,0.00000000003174,0.0000000007442},
113  {0.0000001297,0.000004927,-0.0001504,-0.0000000008666,0.00000006669},
114  {-0.0000000005407,-0.0000001325,0.000007373,-0.0000000000004573,-0.0000000009105},
115  {0.00000002363,-0.00000007727,0.000006433,-0.0000000002033,-0.0000000008472}
116 
117 };
118 
119 // From Melinder 2010
120 static double Glycerol[18][5]={
121  {-13,1093,3486,0.4532,1.52},
122  {-0.0008638,-0.3624,3.766,0.0009782,-0.03729},
123  {0.000006895,-0.002451,-0.0001222,-0.0000001196,0.0003572},
124  {0.0000005229,0.00001547,0.00003219,-0.000000008778,-0.000003648},
125  {-0.5742,2.74,-22.5,-0.003223,0.0445},
126  {-0.000007991,-0.008373,0.1811,-0.00001951,-0.0002688},
127  {-0.0000007515,0.00009596,0.0003766,0.0000001178,0.0000008876},
128  {0.0000000171,-0.0000006999,-0.0000173,0.0000000001048,-0.00000002209},
129  {-0.009119,0.004081,-0.03258,0.000005539,0.0003633},
130  {0.000002973,0.00001808,0.0007249,-0.00000006878,-0.000004088},
131  {0.000000002379,-0.0000008516,-0.00002133,-0.000000002587,0.00000006219},
132  {-0.000000001237,0.00000001441,0.0000004907,0.00000000002262,0.0000000006331},
133  {-0.0001641,-0.00004744,0.002922,-0.0000002073,0.000001069},
134  {-0.000000005313,0.000001833,-0.00005346,-0.0000000002235,-0.0000001248},
135  {0.0000000004546,-0.00000003077,0.00000023,-0.00000000001421,0.000000003019},
136  {-0.000002408,-0.000001415,0.00002238,0.000000003,0.00000005729},
137  {-0.000000001682,0.00000001863,-0.0000005383,0.0000000001604,-0.00000000178},
138  {-0.000000007734,-0.000000006097,-0.0000005944,0.0000000002221,0.000000002116}
139 };
140 
141 // From Melinder 2010
142 static double NH3[18][5]={
143  {-25.76,944.5,4233,0.4551,0.9255},
144  {-0.0001817,-0.2743,-1.618,0.001673,-0.03439},
145  {0.00001204,-0.003113,0.0161,-0.000002214,0.0003217},
146  {0.0000005567,0.000003349,0.00001662,0.0000001228,-0.000004544},
147  {-2.385,-2.914,1.145,-0.005216,0.01327},
148  {0.00002315,-0.02307,0.02715,0.000003544,0.0001856},
149  {0.0000001415,0.0001341,0.001072,0.000001057,-0.00001646},
150  {-0.00000004244,-0.0000005151,-0.00005266,-0.00000003474,0.0000003004},
151  {-0.07636,0.02262,-0.001965,0.00008909,-0.0005979},
152  {-0.00000229,0.00006645,0.003472,0.000003526,-0.00002184},
153  {-0.000000262,-0.0000087,-0.00009051,-0.0000001782,0.000001297},
154  {0.0000000001786,0.00000008999,0.000002106,0.000000001858,-0.00000001141},
155  {-0.00261,-0.0006685,0.002131,0.000006807,-0.0001097},
156  {-0.000000376,-0.000001002,0.00004117,-0.0000003394,0.00000167},
157  {0.0000000136,0.0000003309,0.0000004446,0.000000008315,-0.00000003377},
158  {-0.000073,0.000001635,0.0002136,-0.0000001898,0.000003065},
159  {0.00000003524,0.0000004465,-0.00001354,0.000000006304,-0.00000006166},
160  {-0.000001075,0.0000006298,-0.000008551,-0.00000001361,0.0000003244}
161 };
162 
163 // From Melinder 2010
164 static double K2CO3[18][5]={
165  {-10.3,1216,3217,0.5622,0.8063},
166  {-0.0001575,-0.4114,1.492,0.001656,-0.02362},
167  {0.000003598,-0.003383,-0.001539,0.000002108,0.0001851},
168  {0.000000004324,0.0000299,-0.00003546,-0.00000004482,-0.000002372},
169  {-0.7786,10.64,-37.33,-0.000892,0.03624},
170  {0.00001112,-0.007614,-0.01775,0.000002031,-0.00001262},
171  {-0.0000003479,0.00005214,0.0005416,-0.0000002616,-0.0000003022},
172  {-0.000000001244,-0.0000008087,0.00000659,0.0000000007965,-0.000000007761},
173  {-0.02766,0.04413,0.2023,-0.0000005233,0.0006659},
174  {0.000001616,0.00007806,0.004553,-0.0000002962,-0.00001611},
175  {-0.00000001681,-0.000001173,0.00003587,-0.000000009174,0.000000153},
176  {0.00000000003847,0.00000005658,-0.0000003707,0.0000000001027,0.000000001061},
177  {-0.0008226,-0.0001333,0.01971,-0.0000009283,0.00001077},
178  {-0.00000004913,-0.000002381,0.0001367,-0.00000001814,-0.00000009796},
179  {0.000000001395,0.0000001696,-0.000003667,0.0000000008767,0.00000000307},
180  {-0.000002372,-0.00001246,0.0003405,-0.00000001011,-0.0000001097},
181  {-0.000000002886,0.0000002708,-0.00001676,0.0000000008471,0.000000007825},
182  {0.0000003251,0.0000002035,-0.00003488,0.000000001311,-0.000000008453}
183 };
184 
185 // From Melinder 2010
186 static double CaCl2[18][5]={
187  {-16.21,1171,3133,0.558,0.8939},
188  {-0.0001344,-0.1463,2.81,0.00146,-0.02647},
189  {0.000005073,-0.001302,-0.01563,0.0000003861,0.0001718},
190  {-0.0000000482,-0.0001871,-0.00001233,0.00000001307,-0.0000007918},
191  {-1.555,9.847,-44.8,-0.00115,0.04389},
192  {0.00002146,-0.02488,0.03271,-0.00001008,0.0002102},
193  {-0.0000015,-0.000553,-0.001205,-0.0000000654,-0.0000008688},
194  {0.00000002219,0.00001665,0.000009346,0.0000000004728,-0.00000004353},
195  {-0.05496,0.03389,0.9511,-0.00001784,0.0009121},
196  {0.0000009415,-0.002302,-0.005191,0.0000003496,0.000003993},
197  {0.00000006185,0.00003526,0.0002282,-0.00000000484,0.0000003785},
198  {-0.000000001723,0.0000002788,-0.000000929,-0.0000000002011,-0.000000009979},
199  {-0.002624,0.001062,0.01472,-0.0000004415,0.000001963},
200  {-0.0000001082,0.00006291,0.0001615,-0.000000003462,-0.0000004892},
201  {0.000000003036,0.000001806,0.000005073,0.0000000003922,0.00000000001526},
202  {-0.0001726,0.00002785,-0.001346,0.00000004699,0.0000002997},
203  {-0.000000004396,0.000006423,-0.000009754,0.0000000006179,-0.00000003631},
204  {-0.000004494,-0.000001811,-0.00006674,0.000000002929,0.00000003435}
205 };
206 
207 // From Melinder 2010
208 static double MgCl2[18][5]={
209  {-15.12,1124,3365,0.5461,0.9573},
210  {-0.0004843,-0.3072,2.229,0.001784,-0.03065},
211  {0.00001113,-0.003295,-0.004627,-0.0000008171,0.0001115},
212  {0.0000001858,0.00001015,0.00009186,-0.00000006594,-0.000002923},
213  {-1.885,9.071,-52.22,-0.00273,0.04968},
214  {-0.00005461,-0.006513,0.1162,-0.00001483,0.0001559},
215  {0.000003579,0.00004664,0.001249,0.000000385,-0.00001796},
216  {-0.00000003999,0.000002287,0.000002421,-0.000000005556,0.0000003051},
217  {-0.05455,0.02449,0.6202,0.000008675,-0.002722},
218  {0.00001887,0.00003574,0.002337,-0.000001489,-0.00001753},
219  {-0.0000003171,0.000004337,0.0000724,-0.00000003882,0.000002021},
220  {-0.000000006246,0.0000006044,-0.000003613,0.0000000009282,0.00000002614},
221  {-0.0007257,0.003402,0.01052,0.0000008651,0.00009351},
222  {0.0000007588,-0.0001409,-0.000459,0.0000001992,-0.000008353},
223  {-0.00000004102,0.000001857,-0.00002477,-0.000000001196,0.0000001901},
224  {-0.0003208,0.0003344,-0.001067,-0.0000004779,0.00007364},
225  {-0.0000000492,-0.00000683,-0.0001048,0.00000001797,-0.0000004014},
226  {-0.00001794,0.000007239,-0.0000696,-0.00000002503,0.000003717}
227 };
228 
229 // From Melinder 2010
230 static double NaCl[18][5]={
231  {-9.383,1099,3593,0.5736,0.4369},
232  {-0.00002581,-0.3758,1.669,0.001595,-0.02666},
233  {0.000001423,-0.002023,-0.02019,-0.0000004003,0.0002035},
234  {0,0,0,0,0},
235  {-0.9039,7.723,-32.48,-0.0009383,0.02346},
236  {0.000002578,-0.01426,-0.03331,-0.00001248,-0.00005368},
237  {-0.00000003318,0.0001535,-0.001164,0.0000003353,0.000002871},
238  {0,0,0,0,0},
239  {-0.02204,0.02567,0.6453,-0.00001057,0.0004276},
240  {0.0000001192,0.0003994,-0.009314,0.0000004158,-0.000004526},
241  {-0.000000008993,-0.000007281,0.0002236,-0.00000002032,0.0000001838},
242  {0,0,0,0,0},
243  {-0.0004827,0.0001108,-0.01629,-0.0000004853,0.000007386},
244  {-0.00000001467,0.0000003522,0.0007927,-0.00000001587,0.0000005437},
245  {0,0,0,0,0},
246  {0.000002247,-0.00001686,0.0002331,-0.000000004654,0.0000004688},
247  {0,0,0,0,0},
248  {0,0,0,0,0}
249 };
250 
251 // From Melinder 2010
252 static double KAC[18][5]={
253  {-17.04,1138,3327,0.4958,1.042},
254  {-0.0001082,-0.3565,1.806,0.00134,-0.03071},
255  {0.000006892,-0.00202,-0.001766,0.00000006837,0.0002819},
256  {-0.0000001397,0.000004205,0.00004357,0.000000002637,-0.00000219},
257  {-1.228,5.796,-28.95,-0.002931,0.03405},
258  {0.0000002302,-0.0079,0.04131,-0.00001477,-0.0001511},
259  {-0.0000008711,0.00002753,0.0004429,0.00000004659,0.000001172},
260  {0.00000002016,0.0000000514,0.00001125,0.0000000002886,-0.00000002379},
261  {-0.03862,0.01306,0.04663,0.00001032,0.0005017},
262  {0.000001565,0.00006845,0.0007775,-0.0000002396,-0.000007779},
263  {0.000000007565,-0.00000113,0.00003463,-0.000000004352,0.00000009125},
264  {-0.0000000003063,0.00000002433,-0.0000007261,-0.0000000000223,-0.000000001888},
265  {-0.0004571,-0.001427,-0.001249,-0.0000002024,0.000005637},
266  {-0.00000003734,0.0000008304,0.00005115,-0.000000005541,0.00000002534},
267  {0.000000001268,0.00000001303,-0.000002987,0.00000000008984,0.000000001596},
268  {0.00002969,0.000009353,0.0001659,-0.000000002371,-0.0000002922},
269  {-0.000000002817,0.00000002322,-0.000005193,0.0000000005573,0.000000004601},
270  {0.000000869,0.000002285,-0.0000004612,0.0000000005515,-0.00000000796}
271 };
272 
273 // From Melinder 2010
274 static double KFO[18][5]={
275  {-20.19,1189,3144,0.5253,0.8088},
276  {-0.0001703,-0.3515,1.698,0.001241,-0.02556},
277  {-0.000007478,-0.001918,-0.001303,-0.00000003799,0.0002195},
278  {0.0000003761,0.00003132,0.00005177,-0.000000002951,-0.000001667},
279  {-1.106,7.044,-29.94,-0.001972,0.01758},
280  {-0.000004203,-0.00656,0.0229,-0.000006322,0.00008603},
281  {-0.0000005737,0.00007018,0.000003898,0.00000002654,0.000002498},
282  {0.00000001474,-0.000001459,0.000007391,-0.0000000007324,-0.00000003569},
283  {-0.021,0.01889,0.2944,-0.00002596,-0.00008372},
284  {-0.0000003802,0.0001034,-0.002482,-0.00000004693,-0.000001601},
285  {0.00000006504,-0.000003889,0.000033,-0.000000004362,0.000000004815},
286  {-0.000000001877,0.000000009225,-0.0000004871,0.00000000002899,-0.000000001861},
287  {-0.0002638,-0.0002262,0.001161,-0.0000005062,-0.000001184},
288  {0.00000004027,0.0000003692,0.00006758,-0.000000005229,-0.0000001331},
289  {-0.0000000004789,0.00000006609,-0.000001277,0.0000000001035,-0.000000005489},
290  {0.0000008327,0.00002368,-0.0001429,0.00000001501,0.000001088},
291  {0.0000000009507,0.00000008766,0.0000009949,0.0000000005376,-0.000000007003},
292  {0.00000006345,0.000001061,-0.000003221,0.0000000005562,0.00000003098}
293 };
294 
295 // From Melinder 2010
296 static double LiCl[18][5]={
297  {-23.29,1088,3383,0.5362,1.013},
298  {0.0006555,-0.1772,3.958,0.001454,-0.03062},
299  {-0.0001208,-0.002619,-0.0003035,-0.0000000326,0.000294},
300  {0.000002616,0.000006209,-0.000003477,-0.0000000142,-0.000002719},
301  {-3.051,6.056,-50.36,-0.001855,0.0392},
302  {-0.0003972,-0.008588,0.4415,-0.00001405,0.00006246},
303  {0.00003674,0.0001567,-0.0002609,-0.000000005424,-0.000001752},
304  {-0.0000005569,-0.000001847,0.000003651,0.0000000009821,0.00000008346},
305  {-0.179,0.02556,0.6298,0.00001017,0.000332},
306  {0.00001391,0.00007194,-0.004384,0.0000006821,0.000000784},
307  {-0.000001997,-0.00001053,0.0001039,-0.000000008674,-0.00000031},
308  {0.00000002931,0.00000009716,-0.000001076,0.0000000001095,0.00000001381},
309  {-0.002917,0.0009407,0.04544,0.0000007084,0.000002206},
310  {0.00000149,-0.000007253,-0.0008787,-0.00000007434,-0.0000006011},
311  {-0.00000006904,0.0000003144,-0.000008457,0.0000000006988,0.000000004023},
312  {0.0005715,-0.00008105,0.002527,-0.0000001273,0.000001745},
313  {0.0000001186,0.000001072,-0.0000358,-0.000000003058,-0.00000007094},
314  {0.00002757,-0.000003974,0.00004058,-0.000000009124,0.00000006699}
315 };
316 
317 static int a[18][2]={
318  {0,0},
319  {0,1},
320  {0,2},
321  {0,3},
322  {1,0},
323  {1,1},
324  {1,2},
325  {1,3},
326  {2,0},
327  {2,1},
328  {2,2},
329  {2,3},
330  {3,0},
331  {3,1},
332  {3,2},
333  {4,0},
334  {4,1},
335  {5,0}
336 };
337 double SecFluidsSI(std::string Output, double T, double p, std::string Ref)
338 {
339  double p_kSI = convert_from_SI_to_unit_system(iP,p,UNIT_SYSTEM_KSI);
340 
341  double out = SecFluids(Output, T, p_kSI, Ref);
342 
343  return convert_from_unit_system_to_SI(get_param_index(Output),out,get_standard_unit_system());
344 }
345 double SecFluids(std::string Output, double T, double p, std::string Ref)
346 {
347  double Tfreeze,Tmax,rho,cp,k,mu,u,s,TC,C_gly;
348  // Temperature and Pressure are the inputs
349 
350  long iOutput = get_param_index(Output);
351  // Check for trivial output values
352  if (!Output.compare("T"))
353  return T;
354  else if (!Output.compare("P"))
355  return p;
356 
357  if (!Ref.compare("HC-10"))
358  {
359  // Curve fits generated from Microsoft Excel fit of manufacturer data
360  // Input temperature is in deg C
361  TC=T-273.15;
362  switch(iOutput)
363  {
364  case iD:
365  return -4.52609118E-01*TC + 1.19919457E+03;
366  case iC:
367  return 2.4797494E-03*TC + 3.2708330E+00;
368  case iL:
369  return 1.00000E-06*TC + 5.04400E-04;
370  case iV:
371  return 2.937072591E-12*TC*TC*TC*TC - 1.713938938E-09*TC*TC*TC + 3.826311605E-07*TC*TC - 4.253611683E-05*TC + 2.509839772E-03;
372  case iTfreeze:
373  return 263.15;
374  case iTmax:
375  return 390;
376  case iH:
377  return (2.4797494E-03*TC*TC/2.0 + 3.2708330E+00*TC)/1000+p/(-4.52609118E-01*TC + 1.19919457E+03);
378  case iS:
379  return 2.4797494E-03*(T-298.15) + 3.2708330E+00*log(T/298.15);
380  default:
381  return _HUGE;
382  }
383  }
384  // Split into fluid, concentration pair
385  std::vector<std::string> fluid_concentration = strsplit(Ref,'-');
386 
387  // Check it worked
388  if (fluid_concentration.size() != 2){
389  throw ValueError(format("Format of brine string [%s] is invalid, should be like \"EG-20%\" ", Ref.c_str()) );
390  }
391 
392  // Convert the concentration into a string, discarding anything to the right of the numbers
393  C_gly = strtod(fluid_concentration[1].c_str(),NULL);
394 
395  // Set the temperature to within the band to avoid zero
396  if (iOutput == iTmax || iOutput == iTfreeze){ T=300;}
397 
398  // Melinder correlations are in degree C, so need to convert the temp
399  Brine((char *)fluid_concentration[0].c_str(), T - 273.15, C_gly, &Tfreeze, &Tmax, &rho, &cp, &k, &mu,&u,&s);
400 
401  switch(iOutput)
402  {
403  case iD:
404  return rho;
405  case iC:
406  return cp/1000; // J/kg -> kJ/kg
407  case iL:
408  return k/1000; // W/m/K -> kW/m/K
409  case iV:
410  return mu / 1000.0; // mPa-s --> Pa-s
411  case iTfreeze:
412  return Tfreeze+273.15; // C --> K
413  case iTmax:
414  return Tmax+273.15; // C --> K
415  case iH:
416  return u/1000+p/rho; // J/kg -> kJ/kg
417  case iS:
418  return s/1000; // J/kg/K -> kJ/kg/K
419  default:
420  return _HUGE;
421  }
422 }
423 
424 // inputs in C, percent mass,
425 // outputs in kg/m^3, J/kg-K, mW/m-K, Pa-s, J/kg, J/kg/K
426 int Brine(const char * Mix, double T, double C, /*in --- out */double *Tfreeze, double *Tmax, double *rho, double *cp, double *k, double *visc, double *u, double *s)
427 {
428 
429  double (*A)[18][5];
430  double xm,ym,x,y;
431  double f_Tfreeze=0.0, f_rho=0.0,f_cp=0.0,f_k=0.0,f_visc=0.0,f_u=0.0,f_s=0.0;
432  double Cmin,Cmax;
433  int i;
434 
435  if (!strcmp(Mix,"PG"))
436  { A=&PG; xm=30.7031; ym=32.7083; Cmin=0.0; Cmax=60.0; *Tmax=100;}
437  else if (!strcmp(Mix,"EG"))
438  { A=&EG; xm=30.8462; ym=31.728; Cmin=0.0; Cmax=60.0; *Tmax=100;}
439  else if (!strcmp(Mix,"EA"))
440  { A=&EA; xm=29.2361; ym=8.1578; Cmin=0.0; Cmax=60.0; *Tmax=40;}
441  else if (!strcmp(Mix,"Methanol") || !strcmp(Mix,"MethylAlcohol") || !strcmp(Mix,"MA"))
442  { A=&MA; xm=30.5128; ym=3.5359; Cmin=0; Cmax=60; *Tmax=40;}
443  else if (!strcmp(Mix,"Glycerol"))
444  { A=&Glycerol; xm=36.1905; ym=8.911; Cmin=0; Cmax=60; *Tmax=40;}
445  else if (!strcmp(Mix,"NH3") || !strcmp(Mix,"NH3/H2O"))
446  { A=&NH3; xm=16.0784; ym=-4.649; Cmin=0; Cmax=30; *Tmax=30;}
447  else if (!strcmp(Mix,"K2CO3"))
448  { A=&K2CO3; xm=22.0833; ym=11.2422; Cmin=0.0; Cmax=40.0; *Tmax=40;}
449  else if (!strcmp(Mix,"CaCl2"))
450  { A=&CaCl2; xm=18.7414; ym=7.5257; Cmin=0.0; Cmax=30.0; *Tmax=40;}
451  else if (!strcmp(Mix,"MgCl2"))
452  { A=&MgCl2; xm=14.1327; ym=9.3163; Cmin=0.0; Cmax=30.0; *Tmax=40;}
453  else if (!strcmp(Mix,"NaCl"))
454  { A=&NaCl; xm=13.3897; ym=12.6179; Cmin=0.0; Cmax=23.0; *Tmax=40;}
455  else if (!strcmp(Mix,"KAC"))
456  { A=&KAC; xm=25.6757; ym=6.7757; Cmin=0.0; Cmax=45.0; *Tmax=40;}
457  else if (!strcmp(Mix,"KFO"))
458  { A=&KFO; xm=29.1447; ym=5.8908; Cmin=0.0; Cmax=48.0; *Tmax=40;}
459  else if (!strcmp(Mix,"LiCl"))
460  { A=&LiCl; xm=14.8000; ym=1.4895; Cmin=0.0; Cmax=24.0; *Tmax=40;}
461  else{
462  throw ValueError(format("The name of your Brine diluent [%s] is invalid",Mix));
463  }
464 
465  if (C<Cmin || C>Cmax)
466  {
467  throw ValueError(format("Brine conc. [%g %%] out of range [%g,%g]",C,Cmin,Cmax));
468  *Tfreeze=-1;
469  *Tmax=-1;
470  *rho=-1;
471  *cp=-1;
472  *k=-1;
473  *visc=-1;
474  return -1;
475  }
476 
477  for(i=0;i<18;i++)
478  {
479  x=C; y=T;
480  f_Tfreeze += (*A)[i][0] * pow((x-xm),a[i][0]) * pow((y-ym),a[i][1]);
481  f_rho += (*A)[i][1] * pow((x-xm),a[i][0]) * pow((y-ym),a[i][1]);
482  f_cp += (*A)[i][2] * pow((x-xm),a[i][0]) * pow((y-ym),a[i][1]);
483  f_k += (*A)[i][3] * pow((x-xm),a[i][0]) * pow((y-ym),a[i][1]);
484  f_visc += (*A)[i][4] * pow((x-xm),a[i][0]) * pow((y-ym),a[i][1]);
485  f_u += (*A)[i][2] * pow((x-xm),a[i][0]) * pow((y-ym),a[i][1]+1)/(a[i][1]+1);
486  int n = a[i][1];
487  if (n == 0){
488  f_s += (*A)[i][2] * pow((x-xm),a[i][0]) * log((y+273.15)/298.0);
489  }
490  else if (n == 1){
491  f_s += (*A)[i][2] * pow((x-xm),a[i][0]) * ((y+273.15)-(ym+273.15)*log(y+273.15));
492  }
493  else if (n == 2){
494  f_s += (*A)[i][2] * pow((x-xm),a[i][0]) * (pow(ym+273.15,2)*log(y+273.15)+0.5*(y+273.15)*((y+273.15)-4*(ym+273.15)));
495  }
496  else if (n == 3){
497  double _y = y+273.15, b = ym+273.15;
498  f_s += (*A)[i][2] * pow((x-xm),a[i][0]) * (-pow(b,3)*log(_y)-1.0/6.0*(b-_y)*(11*pow(b,2)-7*b*_y+2*pow(_y,2)));
499  }
500  }
501 
502  if (T < f_Tfreeze || T > *Tmax)
503  {
504  throw ValueError(format("Brine Temp.(%g C) out of range [%g, %g] for fluid %s\n",T,f_Tfreeze,*Tmax,Mix));
505  *Tfreeze=f_Tfreeze;
506  *rho=-2;
507  *cp=-2;
508  *k=-2;
509  *visc=-2;
510  *u=-2;
511  *s=-2;
512  return -1;
513  }
514  else
515  {
516  // There is a typo in the specific heat for Ethylene Glycol, coefficients give a result in kJ/kg
517  if (!strcmp(Mix,"EG")){
518  f_cp *= 1000;
519  f_u *= 1000;
520  f_s *= 1000;
521  }
522  *Tfreeze=f_Tfreeze;
523  *rho=f_rho;
524  *cp=f_cp;
525  *k=f_k;
526  *visc=exp(f_visc);
527  *u=f_u;
528  *s=f_s;
529  return 0;
530  }
531 }
CoolProp.h
iC
Definition: GlobalConstants.h:13
Brine.h
strsplit
std::vector< std::string > strsplit(std::string s, char del)
Definition: CoolPropTools.cpp:41
get_standard_unit_system
EXPORT_CODE int CONVENTION get_standard_unit_system(void)
Definition: CoolPropDLL.cpp:248
format
std::string format(const char *fmt,...)
Definition: CoolPropTools.cpp:22
convert_from_SI_to_unit_system
double convert_from_SI_to_unit_system(long iInput, double value, int new_system)
Definition: Units.cpp:123
iH
Definition: GlobalConstants.h:13
iS
Definition: GlobalConstants.h:13
x
std::vector< double > x(ncmax, 0)
iD
Definition: GlobalConstants.h:13
SecFluids
double SecFluids(std::string Output, double T, double p, std::string Ref)
Definition: Brine.cpp:345
iL
Definition: GlobalConstants.h:13
Brine
int Brine(const char *Mix, double T, double C, double *Tfreeze, double *Tmax, double *rho, double *cp, double *k, double *visc, double *u, double *s)
Definition: Brine.cpp:426
iV
Definition: GlobalConstants.h:13
SecFluidsSI
double SecFluidsSI(std::string Output, double T, double p, std::string Ref)
Definition: Brine.cpp:337
get_param_index
long get_param_index(std::string param)
Definition: CoolProp.cpp:288
iTmax
Definition: GlobalConstants.h:15
convert_from_unit_system_to_SI
double convert_from_unit_system_to_SI(long iInput, double value, int old_system)
Definition: Units.cpp:5
iP
Definition: GlobalConstants.h:13
iTfreeze
Definition: GlobalConstants.h:13
UNIT_SYSTEM_KSI
Definition: GlobalConstants.h:39
Generated on Sun Aug 10 2014 16:02:39 for CoolProp by   doxygen 1.8.7