Mathematical Modelling of Cacl Aqueous Solutions Thermophysical Properties
Mathematical Modelling of Cacl Aqueous Solutions Thermophysical Properties
417 – 426
Chimie şi Inginerie Chimică, Biotehnologii, Industrie Alimentară
INTRODUCTION
Calcium chloride is a nontoxic inorganic salt which can be easily obtained at accessible
prices. It is characterized by a high hygroscopicity and an important sorption capacity. It
has also a good thermal conductivity, a high fusion latent heat and a low volume
modification in phase transition [1].
Its properties make it suitable for uses in various fields such as chemistry, food
products, cosmetics etc. In chemical industry, pure calcium chloride is used for
producing different compounds such as calcium carbonate, calcium phosphate or
calcium chromate [2, 3]. In food industry, this salt can be incorporated in coating
materials [4 – 7]; it can reduce the non-desirable storage effects (softening,
discoloration etc.) of fruits and vegetable by offering a good protection against enzymes
action [5, 6, 8, 9]; it enhance meat products properties [10, 11] etc. In cosmetic
production processes, calcium chloride, mixed with other ingredients, can be used for
example for bigels manufacturing [12]. In pharmaceutics development, calcium chloride
presents synergic effects on in vitro drugs skin permeation [13] or it can serve to
formulate products aimed to improve oral delivery of low soluble drugs [14]. Due to its
cooling properties, calcium chloride is also used to obtain selective water composite
adsorbents [15, 16], liquid desiccants [17 – 19] or cooling agents [1].
In many of the above mentioned utilizations examples, calcium chloride is often
employed as aqueous solutions whose thermophysical properties (namely density,
dynamic viscosity, thermal conductivity etc.) present high importance in the appropriate
transport pipes dimensions and industrial equipment choosing and / or designing
processes.
As consequence, over the time, several studies were aimed to establish values for the
cited calcium chloride properties and to correlate them with various parameters such as
pressure, temperature, solute concentration different data being nowadays available [20
– 22]. These data are frequently expressed in tabular or graphical form fact that makes
their employment rather complicate. A more elegant, adequate and easy to implement
method consist in transforming the existing data in mathematical equations.
Considered a suitable alternative in different other areas such as juice production [23],
study of flour and dough rheological properties [24], study of dough dynamic viscosity
evolution [25] etc., this data interpretation way was successfully applied also to other
fluids thermodynamic properties studies like aqueous solutions of sodium chloride [26],
glucose [27] or glycerol [28, 29] etc.
Therefore, the present work was focused on the study of the behavior of calcium
chloride aqueous solutions density and dynamic viscosity both known as significant in
characterizing the behavior when passing through pipes as well as in heat and mass
transfer in many processing processes.
Data from existing publications were used for establish mathematical correlations for
these properties when influenced by parameters such as temperature (both in positive
and negative ranges) and solutions concentration.
The obtained equations have a high similarity with the experimental data fact that
confirms the accuracy of the mathematical models.
Table 1. Variation of CaCl2 aqueous solutions density with temperature and concentration [30]
Concentration,
Density, ρ [kg·m-3]
C% [%, w/w]
Temperature, T [K]
268.15 273.15 293.15 303.15 313.15 333.15 353.15 373.15 393.15 413.15
2 - 1017.1 1014.8 1012 1008.4 999.4 988.1 974.8 959.6 942.8
4 - 1034.6 1031.6 1028.6 1024.9 1015.8 1004.6 991.5 976.5 960.1
8 1070.8 1070.3 1065.9 1062.6 1058.6 1049.2 1038.2 1025.7 1011.1 995.4
12 1108.3 1107.2 1101.5 1097.8 1093.7 1084.0 1073.0 1061.0 1046.6 1031.7
16 1147.1 1145.4 1138.6 1134.5 1130.1 1120.2 1109.2 1097.3 1083.5 1069.1
20 1187.4 1185.3 1177.5 1173.0 1168.4 1158.1 1147.1 1135.2 1121.9 1108.0
25 - 1237.6 1228.4 1223.6 1218.6 1207.9 1196.5 1184.6 - -
30 - 1292.2 1281.6 1276.4 1270.9 1259.7 1247.8 1235.9 - -
35 - - 1337.3 1331.6 1325.5 1313.7 1301.3 1289.3 - -
40 - - 1395.7 1389.5 1382.6 1370 1357.1 1345 - -
Table 2. Variation of CaCl2 aqueous solutions density with temperature and concentration [30]
Density, ρ [kg·m-3]
Concentration,
Temperature, T [K]
C% [%, w/w]
288.15 273.15 263.15 253.15 243.15 233.15
15 1132 1137 1140 - - -
16 1142 1147 1150 - - -
17 1151 1157 1160 - - -
18 1161 1167 1170 - - -
19 1171 1177 1180 - - -
20 1181 1187 1190 - - -
21 1191 1197 1201 1205 - -
22 1201 1207 1211 1215 - -
23 1211 1218 1222 1226 - -
24 1222 1228 1233 1237 - -
25 1232 1239 1244 1248 - -
26 1243 1250 1254 1259 1263 -
27 1252 1261 1266 1270 1275 -
28 1264 1272 1277 1282 1287 -
29 1275 1283 1288 1293 1298 1303
30 1286 1294 1298 1304 1310 1315
Microsoft Excel™ 2013, CurveExpert® and TableCurve 3D® v.4 software’s were
employed for mathematical models development and for data 2D and 3D plotting.
Table 3. Variation of CaCl2 aqueous solutions dynamic viscosity with temperature and
concentration [30]
Concentration,
Temperature, T [K]
293.15 283.15 273.15 268.15 263.15 258.15 253.15 248.15 243.15 238.15 233.15 228.15 223.15
9.4 12.36 15.49 21.58 25.48 - - - - - - - - -
14.7 14.91 18.63 25.6 30.4 40.6 - - - - - - - -
18.9 18.04 22.36 29.91 34.32 46.68 61.49 - - - - - - -
20.9 20.01 24.52 32.75 38.25 50.7 65.9 - - - - - - -
23.8 23.54 28.73 38.15 44.13 59.23 75.51 94.73 115.72 - - - - -
25.7 26.28 32.17 42.56 - 66.78 83.65 105.62 129.15 148.08 - - - -
27.5 29.32 36.09 48.05 - 75.22 - 118.66 147.1 171.61 215.74 - - -
28.4 31.38 - 51.19 - 80.22 - 126.51 159.85 188.29 245.17 304.01 - -
29.4 34.03 - 54.92 - 86.3 - 138.27 - 212.8 254.97 323.62 402.07 490.33
29.9 35.11 - 56.88 - 90.43 - 144.16 - 225.55 284.39 353.04 431.49 509.95
Density
(Equation 2) for A2 and A3 coefficients with good regression coefficients (Table 5) have
been established.
Coefficient a1 a 2 C % a 3 C %2 (2)
Coefficient a1 a 2 C % (3)
Equations 1, 2 and 3 were combined in order to obtain the general model expressed by
Equation 4:
( a1 a2 C% a3 C% 2 ) A ( a1 a2 C% ) A T ( a1 a2 C% )A T 2
1 2 3
(4)
A parallel between data calculated with the mathematical model and the experimental
ones was realized by the mean of relative error equation (5).
Dataexperimental Datacalculated
100[%] (5)
Datacalculated
An average of 0.630% (in absolute value) was obtained. The regression coefficient R2 of
the proposed model was 0.9946 for densities values presented in Table 1 and higher
(0.9999) for those given in Table 2. In both cases the absolute value of the relative error
was of 0.524%.
The ANOVA (Two-Factor with Replication) test was used also to compare tabular and
calculated density values. The results revealed sample P-values of 0.900567 and
0.941336, greater than the targeted alpha 0.05, and Fcrit values of 3.89364 superior to
F-test values (0.015655 and 0.005431) for calcium chloride aqueous solutions density
values variation in positive (Table 1) and negative (Table 2) temperature ranges. In this
case, the null hypothesis cannot be rejected and no statistical difference between tabular
and calculated data was registered.
Density data were introduced in TableCurve 3D® v.4 software which served to generate
two different polynomial equations (6 and 7): one classified Rank 49 (Eqn. 313) with a
precision of R2 = 0.9999966601, FitSdErr = 0.2326931809, Fstat. = 2.46181E+06 for
values shown in Table 1 and one of Rank 16 (Eqn. 310) with R2 = 0.9999560422,
FitSdErr = 0.3701979907, Fstat. = 139015.96133 for values exposed in Table 2.
b1 b2 C% b3 ln(T ) b4 C%2 b5 (ln T ) 2 b 6 C% ln(T ) b7 T 3
2
(6)
b8 (ln T ) 3 b9 C % (ln T ) 2 b10 C % ln T
b1 b2 C% b3 T b4 C%2 b5 T 2 b6 C% T b 7 C% 3 b8 T 3
2
(7)
b9 C % T 2 b10 C % T
Coefficients of the equations 6 and 7 are presented in Tables 9 and 10.
1300
Density, [kg.m-3]
1300
Density, [kg.m-3 ]
Density, [kg.m-3]
1300
1250 1250
1200 1250 1250
1200
1150 1150
1100 1200 1200
1100
1050 1050
1000 1150
1000 1150
950 950
900 1100
900 1100
Co 35 0 27 .5 5 24
n ce 3 5 30 Co 2 25 0
2 0 nce
nt r
atio 20 5 35 ntr 22 .5 0 26 0
0 ] atio 2 0 [K ]
n, C 1 T [K n, C 17 .
5 27 0 re ,
T
% 10 40 r e, 28 at u
[% 5 0 atu % 15 0 pe r
,w 0 450 p er [%
,w .5 29 m
/w] T em /w ] 12 0 Te
Figure 1. CaCl2 aqueous solutions density values plotted in TableCurve 3D and fitted
with polynomial type (Equations 6 and 7) and its residuals
Dynamic viscosity
The same calculus algorithm was followed for calcium chloride aqueous solutions
dynamic viscosity data. In this case, 10 quadratic correlations expressed by equation 8
have been established by plotting the logarithmic values of dynamic viscosity μ [mPa·s]
vs temperature T [K] at fixed CaCl2 concentrations.
log( ) A1 A2T A3T (8)
Table 10 depicts A1, A2 and A3 coefficients values.
A1, A2 and A3 coefficients were correlated with CaCl2 concentrations C% [% w/w], and
various mathematical models (1st, 2nd and 3rd degree polynomial equations, “vapour
pressure” model, “heat capacity” model etc.) were generated.
Among them the best fit was given by the quadratic equation 9 whose coefficients values
are shown in Table 11.
Coefficien t a1 a 2 C % a 3 C %2 (9)
The last mentioned two equations (8 and 9) were combined leading to the following final model:
log( ) ( a1 a2 C% a3 C% ) A1 ( a1 a2 C% a3 C% ) A2 T
2 2
2
(10)
( a1 a2 C% a3 C% )A3 T 2
The error of this model was of 0.897% (2.488% in absolute value) and its regression
coefficient was of 0.996.
The ANOVA test revealed a sample P-value of 0.9647 (superior to targeted alpha of
0.05) and an Fcrit value of 3.881505 (larger than the F-test value of 0.001956). The
same conclusion as for density values analysis was drawn, namely that there is no
statistical difference between tabular and calculated data.
Multiple equations were generated also with TableCurve 3D® v.4 software but only one
was selected due to its high precision (Rank 33, Eqn. 254122972). This equation has a
correlation coefficient of 0.9990438298 and its FitSdErr and Fstat. values were
6922.0578888 and 3.2745062671 respectively. The established equation (11) has the
coefficients presented in Table 14 and the graphical form pictured in Figure 2.
ln( ) b1 b2 C%2 ln( C% ) b3 C% b4 C% b5 T b6 T 0.5 ln( T )
2.5 3
(11)
b7 ln( T )2 b8 T 0.5 b9 ln( T )
500
Kinematic viscosity
Using the data obtained for density and dynamic viscosity of calcium chloride aqueous
solutions, with the help of equation 12, one can calculate also the kinematic viscosity (ν).
[m2·s-1] (12)
Compared to results presented in other studies on the same thematic [20, 21], the
mathematical models developed in this work present a higher degree of accuracy. In
some cases, they are also simpler and therefore easier to use.
CONCLUSION
This paper was able to establish mathematical models relating density and dynamic
viscosity of calcium chloride aqueous solutions with their temperature and concentration.
The results released with the help of Excel Microsoft software showed that a second
order polynomial equation can be considered a suitable mathematical expression of
density while the dynamic viscosity can be presented as a linear model. More complex
equations were obtained with CurveExpert and TableCurve 3D softwares. Correlation
coefficient, relative error and ANOVA values revealed no significant differences
between the experimental and calculated data. The high degree of accuracy registered
for all the mathematical models developed in this study recommends their use for data
base enrichment but also for equipment and processes sizing and optimization.
REFERENCES
1. N’Tsoukpoe, K.E., Rammelberg, H.U., Lele, A.F., Korhammer, K., Watts, B.A., Schmidt, T.,
Ruck, W.K.L.: A Review on the Use of Calcium Chloride in Applied Thermal Engineering,
Applied Thermal Engineering, 2015, 75, 531-531;
2. Harja, M., Ciocîntă, R.C., Bărbuţă, M., Rusu, L., Simion, A.I., Bistriceanu, I.L.: CaCO3
Controllable Synthesis by Double Exchange Method Using CaCl2 Residual Solutions,
Environmental Engineering and Management Journal, 2010, 9 (11), 1571-1577;
3. Fakeev, A.A., Vasil’eva, L.V., Sukhanovskaya, A.I.: A Study of Calicum Chloride Purification,
Russian Journal of Applied Chemistry, 2003, 76 (2), 167-170;
4. Chuang, L., Panyoyaim N., Katopo, L., Shanks, R., Kasapis, S.: Calcium Chloride Effects on the
Glass Transition of Condensed Systems of Potato Starch, Food Chemistry, 2016, 199, 791-798;
5. Khaliq, G., Mahomed, M.T.M., Ali, A., Ding, P., Ghazali, H.M.: Effect of Gum Arabic Coating
Combined with Calcium Chloride on Physico-Chemical and Qualitative Properties o Mango
(Mangifera indica L.) Fruit During Low Temperature Storage, Scientia Horticulturae, 2015, 190,
187-194;
6. Khaliq, G., Mahomed, M.T.M., Ghazali, H.M., Ding, P., Ali, A.: Influence of Gum Arabic Coating
Enriched with Calcium Chloride on Physiological, Biochemical and Quality Responses of Mango
(Mangifera infica L.) Fruit Stored Under Low Temperature Stress, Postharvest Biology and
Technology, 2016, 111, 362-369;
7. Liu, X., Ren, J., Zhu, Y., Han, W., Xuan, H., Ge, L.: The Preservation Effect of Ascorbic Acid and
Calcium Chloride Modified Chitosan Coating on Fresh-Cut Apples at Room Temperature,
Colloids and Surfaces A: Physicochem. Eng. Aspects, 2016, 502, 102-106;
8. Rahman, M.U., Sajid, M., Rab, A., Ali, S., Shahid, M.O., Alam, A., Israr, M., Ahmad, I.: Impact
of Calcium Chloride Concentrations and Storage Duration on Quality Attributes of Peach (Prunus
persica), Russian Agricultural Sciences, 2016, 42 (2), 130-136;
9. Teoh, L.S., Lasekan, O., Adzahan, N.M., Hashim, N.: The Effect of Combinations of UV-C
Exposure with Ascorbate and Calcium Chloride Dips on the Enzymatic Activities and Total
Phenolic Content of Minimally Processed Yam Slices, Postharvest Biology and Technology, 2016,
120, 138-144;
10. Pan, T., Guo, H., Li, Y., Song. J., Ren, F.: The Effects of Calcium Chloride on the Gel Properties
of Porcine Mycosin-k-Carrageenan Mixtures, Food Hydrocolloids, 2017, 63, 467-477;
11. Bunmee, T., Jaturasitha, S., Kreuzer, M., Wicke, M.: Can Calcium Chloride Injection Facilitate the
Ageing-Derived Improvement in the Quality of Meat from Culled Dairy Cows?, Meat Science,
2014, 96 (4), 1440-1445;
12. Lupi, F.R., Shakeel, A., Greco, V., Rosi, C.O., Baldino, M., Gabriele, D.: A Rheological and
Microstructural Characterisation of Bigels for Cosmetic and Pharmaceutical Uses, Materials
Science and Engineering C, 2016, 69, 358-365;
13. Tokudome, Y., Sugibayashi, K.: Mechanism of the Synergic Effects of Calcium Chloride and
Electroporation on the in Vitro Enhanced Skin Permeation of Drugs, Journal of Controlled
Release, 2004, 95 (2), 267-274;
14. Yasmin, R., Rao, S., Bremmell, K., Prestidge, C.: Synergistic Role of Solid Lipid and Porous
Silica in Improving the Oral Delivery of Weakly Basic Poorly Water Soluble Drugs, European
Journal of Pharmaceutical Sciences, 2017, 96, 508-514;
15. Saha, B.B.; Chakraborty, A., Koyama, S., Aristov, Y.I.: A New Generation Cooling Device
Employing CaCl2-in-Silica Gel-Water System, International Journal of Heat and Mass Transfer,
2009, 52 (1-2), 516-524;
16. Miyazaki, T., Akisawa, A., Saha, B.B., El-Sharkawy, I.I., Chakraborty, A.: A New Cycle Time
Allocation for Enhancing the Performance of Two-Bed Adsorption Chillers, International Journal
of Refrigeration, 2009, 32 (5), 846-853;
17. Bansal, P., Jain, S., Moon, C.: Performance Comparison of an Adiabatic and an Internally Cooled
Structured Packed-Bed Dehumidifier, Applied Thermal Engineering, 2011, 31 (1), 14-19;
18. Rafique, M.M., Gandhidasan, P., Bahaidarah, H.M.S.: Liquid Desiccant Materials and
Dehumidifiers – A Review, Renewable and Sustanable Energy Reviews, 2016, 56, 179-195;
19. Ghiasi, M.M., Bahadori, A., Zendehboudi, S.: Estimation of the Water Content of Natural Gas
Dried by Solid Calcium Chloride Dehydrator Units, Fuel, 2014, 117, 33-42;
20. Abdulagatov, I.M., Azizov, N.D.: Viscosity of Aqueous Calcium Chloride Solutions at High
Temperatures and High Pressures, Fluid Phase Equilibria, 2006, 240 (2), 204-219;
21. Kodejs, Z., Slama, I., Novak, J.: Viscosity of the Concentrated Aqueous Solutions of Calcium
Chloride and Nitrate, Chemicke Zvesti, 1976, 30 (4), 439-445;
22. Akhmedova-Azizova, L.A., Abdulagatov, I.M.: Thermal Conductivity of Aqueous CaCl2
Solutions at High Temperatures and High Pressures, Journal of Solution Chemistry, 2014, 43, 421-
444;
23. Simion, A.I., Dobrovici, P.E., Rusu, L., Gavrilă, L.: Modelling of the Thermo-Physical Properties
of Grapes Juice III. Viscosity and Heat Capacity, Scientific Study and Research: Chemistry and
Chemical Engineering, Biotechnology, Food Industry, 2011, 12 (4), 409-420;
24. Temea (Moroi), A.M., Pircu (Vartolomei), N., Simion, A.I., Grigoraş, C.-G., Ungureanu
(Cărbune), R.E., Alexe, P.: Improvement of Flour and Dough Rheological Properties by
Maturation Process, Romanian Biotechnological Letters, 2016, 21 (2), 11381-11392;
25. Simion, A.I., Grigoraş, C.-G., Moroi, A., Vartolomei, N.: Mathematical Modelling of Pasta Dough
Dynamic Viscosity. Thermal Conductivity and Diffusivity, The Annals of the University Dunărea
de Jos of Gala i, 2015, 39 (1), 81-92;
26. Simion, A.I., Grigoraş, C.-G., Roşu, A.M., Gavrilă, L.: Mathematical Modelling of Density and
Viscosity of NaCl Aqueous Solutions, Journal of Agroalimentary Processing and Technologies,
2015, 21 (1), 41-52;
27. Simion, A.I., Grigoraş, C.-G., Rusu, L., Gavrilă, L.: Modelling of the Thermophysical Properties
of Liquids Involved in Food Processing, Environmental Engineering and Management Journal,
2012, 11 (11), 2035-2045;
28. Simion, A.I., Grigoraş, C.-G., Gavrilă, L.: Modeling Thermodynamic Properties of Glycerol
Involved in Cosmetic and Pharmaceutical Industries. Boiling Point, Heat Capacity and Thermal
Conductivity, The 5th IEEE International Conference on E-Health and Bioengineering – EHB
2015, Grigore T. Popa University of Medicine and Pharmacy, Iaşi, Romania, November 19-21,
2015, DOI: 10.1109/EHB.2015.7391366;
29. Simion, A.I., Grigoraş, C.-G., Gavrilă, L.: Mathematical Modeling of Glycerol Aqueous Solutions
Thermophysical Properties. Case Study of Density and Viscosity, Journal of Agroalimentary
Processing and Technologies, 2015, 21 (2), 125-135;
30. Macovei, V.M., Culegere de Caracteristici Termofizice pentru Biotehnologie şi Industrie
Alimentară, Ed. Alma, Galaţi, 2000.