Analytical Determination of the Nucleation-Prone, Low-Density Fraction of Subcooled Water

Hellmuth, Olaf; Feistel, Rainer

doi:10.3390/e22090933

Open AccessArticle

Analytical Determination of the Nucleation-Prone, Low-Density Fraction of Subcooled Water

by

Olaf Hellmuth

^1,*

and

Rainer Feistel

²

¹

Leibniz Institute for Tropospheric Research (TROPOS), Permoserstraße 15, D-04318 Leipzig, Germany

²

Leibniz Institute for Baltic Research (IOW), Seestraße 15, D-18119 Rostock-Warnemünde, Germany

^*

Author to whom correspondence should be addressed.

Entropy 2020, 22(9), 933; https://doi.org/10.3390/e22090933

Submission received: 19 June 2020 / Revised: 12 August 2020 / Accepted: 17 August 2020 / Published: 25 August 2020

(This article belongs to the Special Issue Crystallization Thermodynamics)

Download

Browse Figure

Review Reports Versions Notes

Abstract

:

Subcooled water is the primordial matrix for ice embryo formation by homogeneous and heterogeneous nucleation. The knowledge of the specific Gibbs free energy and other thermodynamic quantities of subcooled water is one of the basic prerequisites of the theoretical analysis of ice crystallization in terms of classical nucleation theory. The most advanced equation of state of subcooled water is the IAPWS G12-15 formulation. The determination of the thermodynamic quantities of subcooled water on the basis of this equation of state requires the iterative determination of the fraction of low-density water in the two-state mixture of low-density and high-density subcooled water from a transcendental equation. For applications such as microscopic nucleation simulation models requiring highly frequent calls of the IAPWS G12-15 calculus, a new two-step predictor-corrector method for the approximative determination of the low-density water fraction has been developed. The new solution method allows a sufficiently accurate determination of the specific Gibbs energy and of all other thermodynamic quantities of subcooled water at given pressure and temperature, such as specific volume and mass density, specific entropy, isothermal compressibility, thermal expansion coefficient, specific isobaric and isochoric heat capacities, and speed of sound. The misfit of this new approximate analytical solution against the exact numerical solution was demonstrated to be smaller than or equal to the misprediction of the original IAPWS G12-15 formulation with respect to experimental values.

Keywords:

IAPWS G12-15; equation of state of subcooled water; low-density and high-density water fraction

1. Introduction

The knowledge of the properties of supercooled or subcooled water (the notion used here) is of outstanding importance for a variety of natural and technical processes related to phase transitions of water such as the microphysical evolution of atmospheric clouds (see, e.g., Meyers et al. [1], Khvorostyanov and Sassen [2], Lohmann and Kärcher [3], Lohmann et al. [4], Pruppacher and Klett [5], Heymsfield et al. [6], Jensen and Ackerman [7], Barahona and Nenes [8], Jensen et al. [9], Zasetsky et al. [10], Khvorostyanov and Curry [11], Khvorostyanov and Curry [12], Hellmuth et al. [13], Khvorostyanov and Curry [14], and Lohmann et al. [15]), the cryopreservation of organelles, cells, tissues, extracellular matrices, organs, and foods (see, e.g., Pegg [16] and Espinosa et al. [17,18]); and water vitrification (see, e.g., Debenedetti and Stanley [19], Bhat et al. [20], and Zobrist et al. [21]). Subcooled water is the primordial matrix for ice embryo formation by homogeneous and heterogeneous nucleation. Recently, Hellmuth et al. [22] studied the variation of the homogeneous freezing rate of subcooled water in dependence on pressure and temperature. Within the framework of classical nucleation theory both the thermodynamic driving force of freezing and the surface tension of the water–ice system were scaled in terms of fundamental thermodynamic properties of the macrophases of subcooled water and hexagonal ice, such as the Gibbs energy, entropy, enthalpy, mass density, isobaric heat capacity, isothermal compressibility, and isobaric thermal expansivity. These quantities appeared the generating ones for the determination of the nucleation metrics.

For applications to the aforementioned questions of interest, several equations of state (EoS) for subcooled water became available (for a brief review see Holten et al. [23]). In their study on ice nucleation, Hellmuth et al. [22] employed the Thermodynamic Equation Of Seawater 2010 (TEOS-10), which is an international standard for the thermodynamic properties of seawater, ice, and humid air, adopted in 2009 by the Intergovernmental Oceanographic Commission (IOC) of the UNESCO and in 2011 by the International Union of Geodesy and Geophysics (IUGG) (Feistel et al. [24,25]; Wright et al. [26]; IOC, SCOR, and IAPSO [27]; Feistel [28]; Feistel [29]; TEOS-10 website http://www.teos-10.org/). This seawater standard is based on thermodynamic potentials for fluid water, hexagonal ice, seasalt dissolved in water, and dry air. The fluid water part is described according to Wagner and Pruß [30] and IAPWS R6-95 [31], known as “IAPWS-95”, and took into account quality-assured experimental data which were available by the middle of the year 1994. For details of TEOS-10 and IAPWS-95 the reader is referred to the Supplementary Material (SM), Section SM-1.

In order to better represent the thermodynamic properties of subcooled liquid water, Holten et al. [23,32,33,34] and Holten and Anisimov [35] proposed a new EoS of subcooled water, which considers the experimental data of liquid water in the metastable region of the phase diagram above the homogeneous nucleation temperature. The thermodynamic formulation of this EoS is a mean-field version of an EoS developed by Holten and Anisimov [35]. It is based on the assumption that liquid water is a mixture of a high-density structure A and a low-density structure B of water (“two-state model”, cf. Section 3.1). Such an assumption is justified by empirical findings (see references in Holten and Anisimov [35]). This new EoS was implemented in the IAPWS guideline IAPWS G12-15 [36], given in form of the dimensionless Gibbs energy of subcooled water as a function of temperature, pressure, and the equilibrium mole fraction of the low-density structure of water in the two-state mixture. Having at ones disposal the Gibbs energy of subcooled water, one can determine the mass density, the specific entropy, the isothermal compressibility, the isothermal expansivity, the specific isobaric and isochoric heat capacities, and the speed of sound of subcooled water. While temperature and pressure serve as the independent variables, the equilibrium mole fraction of low-density water is a dependent variable, which enters the determination of the Gibbs energy in form of a function of pressure and temperature. The determination of the low-density water mole fraction requires the numerical solution of a transcendental equation, which serves as a physical constraint for the water composition at prescribed pressure and temperature.

For applications which require frequent calls of the EoS such as for computational fluid dynamics (CFD) modeling, the iterative determination of the low-density water mole fraction might be too time-consuming and should be circumvented. Alternatively, one can either precalculate look-up tables of the low-density water mole fraction with sufficiently high resolution in pressure and temperature, from which the sought after value can be determined by interpolation in two dimensions, or one can try to find a sufficiently accurate approximative solution, which avoids the numerical solution of a transcendental equation. In the present study, we will follow the second attempt.

The aim of the present study is twofold: (i) derivation of an approximative solution of the transcendental equation for the low-density water mole fraction, which allows a highly accurate analytical determination of this quantity as a function of pressure and temperature, and (ii) determination of the deviations of the resulting approximative values of the aforementioned thermodynamic variables from their numerical (“exact”) counterparts throughout the validity range of the underlying EoS. The paper is organized as follows. In Section 2, the material and the method of the present analysis is briefly described. Section 3 addresses the theoretical background of the analysis. In Section 4, the new calculus is applied to the IAPWS G12-15 [36] formulation for subcooled water, and its deviation from the traditional numerical method is analyzed. Section 5 contains the discussion of the results. Finally, the conclusions are drawn in Section 6. The paper is completed by Appendix A containing the algorithm of the calculus in a compact form and by supplementary material, which includes additional explanations, details of the calculus, and a comprehensive table part displaying the verification results. The enumeration of sections, equations, and tables which are presented in the supplementary material, is specified with the prefix SM.

2. Materials and Methods

The starting point of the present analysis is the IAPWS G12-15 [36] formulation for subcooled water, the calculus of which is presented in Section SM-3. The subject of interest of the present study is the governing equation for the mole fraction of low-density water therein, given by Equation (SM-3.8), which is a transcendental equation and ordinarily requires a numerical iteration to find the roots. However, by means of a Taylor expansion up to the quadratic term (quadratization), this equation is brought into the form of a cubic equation, which in turn can be solved analytically by means of Cardano’s method presented in the supplementary material, Section SM-4. To define a benchmark for the subsequent verification of the new analytical solution, at first a numerical solver was implemented and evaluated (cf. Section SM-5). On this basis the differences between the traditional numerical and the new analytical solution of the mole fraction of subcooled water together with the thermodynamic deliveries of the IAPWS G12-15 [36] formulation are determined as functions of pressure and temperature, such as the specific Gibbs energy, the specific volume and mass density, the specific entropy, the isothermal compressibility, the thermal expansion coefficient, the specific isobaric and isochoric heat capacities, and speed of sound (cf. Sections SM-6 and SM-7). In order to assess the goodness of the analytical solution of the EoS, the misprediction of the analytical form with respect to the observations is compared with those of the numerical form. For the comparison, the experimental data of Holten et al. [23] are used.

3. Theoretical Background

3.1. On the Two-State Character of Subcooled Water

The observed anomalous behavior of the compressibility, the thermal expansivity, and viscosity of water and the existence of mass density maximum already prompted Röntgen [37] to view subcooled water as a two-state molecular mixture. According to Röntgen’s perception, molecules of the first kind are ice molecules, which go over into molecules of the second kind upon heat supply. If in turn heat is extracted from the liquid, then a corresponding part of ice molecules will reappear. Subcooled water is considered a saturated solution of ice molecules; the saturated fraction of which depends on temperature. The saturated concentration of ice molecules increases with decreasing temperature. Röntgen [37] further assumed that the transition of molecules of the first kind into the second kind is accompanied by a decrease of the volume of the mixture, as can be observed during melting of ice. In its modern form, this perception is the basis for the IAPWS G12-15 [36] EoS, the rationale of which will be briefly outlined hereafter. According to Holten et al. [34], subcooled water is assumed to be a mixture of two interconvertible states or structures of water: a high-density state A and a low-density state B. The fraction of molecules in state B is denoted by x, and is controlled by the following equilibrium pseudo-reaction (Holten et al. [34], Equation (1) therein):

A ⇌ B .

(1)

As argued by those authors, the states A and B could correspond to different arrangements of the hydrogen-bonded network or to two kinds of local coordination of the water molecules. The molar Gibbs energy of a two-state mixture reads (Holten et al. [34], Equation (2) therein):

G = G^{A} + x G^{BA} + R T [x ln x + (1 - x) ln (1 - x)] + G^{E} .

(2)

Here, R denotes the universal gas constant, T is the temperature,

G^{A}

is the molar Gibbs energy of the pure structure A, defines the background of the properties, and can be approximated by a well-defined polynomial as a function of T and p, i.e.,

G^{A} = G^{A} (T, p)

. The quantity

G^{BA} \equiv G^{B} - G^{A}

is the difference in the molar Gibbs energies between the pure configurations A and B, and

G^{E}

is the excess Gibbs energy of mixing. The difference

G^{BA}

is related to the equilibrium constant

K (T, p)

of the pseudo-reaction given by Equation (1), which reads (Holten et al. [34], Equation (3) therein):

ln K (T, p) = \frac{G^{BA}}{R T} .

(3)

From Equation (3), it follows that

G^{BA} = G^{BA} (T, p)

. The excess Gibbs energy

G^{E}

originates from the non-ideality of the mixture and is given by the sum of contributions of the enthalpy of mixing

H^{E}

and of the excess entropy

S^{E}

(Holten et al. [34], Equation (6) therein):

G^{E} = H^{E} - T S^{E} .

(4)

Subcooled water is approximated by an athermal solution, i.e., the molecules of the two liquids are allowed to differ appreciably in size, and strong orienting forces caused by dipoles or hydrogen bonding are allowed to occur. In such case,

H^{E} = 0

and

S^{E} = - R ω x (1 - x)

, where

ω (p)

denotes the pressure-dependent molecular interaction parameter, which accounts for the non-ideality of the solution, resulting in

G^{E} = G^{E} (x, T, p)

. Knowing x, the EoS of subcooled water delivers the molar Gibbs energies of the pure configurations A and B,

G^{A} (T, p)

and

G^{B} (T, p)

, as well as that of the mixture,

G (x, T, p)

.

The two-state character of subcooled water has implications for ice crystallization which arise from an empirical judgment of Ostwald [38] (pp. 306–309 therein), which has later become known as Ostwald’s rule of stages. The implications of this rule are discussed in Section SM-2.

3.2. Analytical Determination of the Fraction of Low-Density Water in the Two-State Water Mixture

For the sake of completeness, the full calculus of the IAPWS G12-15 [36] formulation of the properties of subcooled water is presented in Section SM-3. From this section we excerpt here only those parts of the calculus which are directly related to the determination of the low-density water mole fraction.

Introducing the reduced temperature

τ

and the reduced pressure

π

(cf. Section SM-3.1, Equation (SM-3.2)):

τ = \frac{T}{T_{LL}} - 1, π = \frac{p}{{\hat{ϱ}}_{0} R_{W} T_{LL}},

(5)

with

T_{LL}

,

{\hat{ϱ}}_{0}

, and

R_{W}

, respectively, denoting the temperature of the liquid-liquid critical point, the reference mass density of subcooled water, and the specific gas constant of water vapor, the equilibrium mole fraction of low-density water,

x_{e}

, was demonstrated to obey the following transcendental equation (cf. Section SM-3.1, Equation (SM-3.8)):

F (x_{e}) = L (τ, π) + ln (\frac{x_{e}}{1 - x_{e}}) + ω (π) (1 - 2 x_{e}) = 0 .

(6)

The dimensionless quantities

L (τ, π)

and

ω (π)

denote the ordering field and the interaction parameter, respectively, and are well-defined functions of

τ

and

π

(cf. Section SM-3.1, Equations (SM-3.5) and (SM-3.6)). The mixed state of subcooled water is characterized by the condition

0 < x_{e} < 1

. Knowing

x_{e} = x_{e} (τ, π)

from the numerical solution of Equation (6), the dimensionless Gibbs free energy

ψ (τ, π)

is fully determined and given in analytical form, from which all other thermodynamic parameters can be analytically determined (cf. Section SM 3.2).

To derive an approximate analytical solution of Equation (6) we introduce the auxiliary parameters

α = 2 ω, β = - (L + ω),

(7)

by means of which Equation (6) can be rearranged as follows (the subscript “e” for the equilibrium value of x will be omitted hereafter):

ln (\frac{x}{1 - x}) = α x + β ⇝ \frac{x}{1 - x} = e^{β} e^{α x} .

(8)

In order to obtain a first approximation of the mole fraction of low-density water,

x \approx \bar{x}

, the following quadratization by expansion of the exponential term in a Taylor series will be employed:

e^{α \bar{x}} \approx 1 + α \bar{x} + \frac{1}{2} α^{2} {\bar{x}}^{2}, \bar{x} ≪ 1 .

(9)

Therewith, Equation (8) can be rearranged into a cubic equation for the determination of

\bar{x}

:

\begin{matrix} \bar{A} {\bar{x}}^{3} + \bar{B} {\bar{x}}^{2} + \bar{C} \bar{x} + \bar{D} = 0, \\ \bar{A} = - \frac{α^{2}}{2} e^{β} = - \frac{α^{2}}{2} \bar{D}, \\ \bar{B} = α (\frac{α}{2} - 1) e^{β} = α (\frac{α}{2} - 1) \bar{D}, \\ \bar{C} = e^{β} (α - 1) - 1 = \bar{D} (α - 1) - 1, \\ \bar{D} = e^{β} . \end{matrix}

(10)

The parameters

\bar{A}

,

\bar{B}

,

\bar{C}

, and

\bar{D}

are well-defined functions of

τ

and

π

. Equation (10) can be solved using Cardano’s method and has three real and/or complex solutions,

{\bar{x}}_{i} = {\bar{x}}_{i} (\bar{A}, \bar{B}, \bar{C}, \bar{D})

,

i = 1, 2, 3

, representing the first guess of the mole fraction of low-density water. We are only interested in the real solution of Equation (10), which reads (cf. Section SM-4, Equation (SM-4.10))

\begin{matrix} \bar{x} & = & \sqrt[3]{- \frac{\bar{q}}{2} + \sqrt{\bar{Δ}}} + \sqrt[3]{- \frac{\bar{q}}{2} - \sqrt{\bar{Δ}}}, \\ \bar{Δ} & = & {(\frac{\bar{q}}{2})}^{2} + {(\frac{\bar{p}}{3})}^{3} \geq 0, \\ \bar{p} & = & \bar{s} - \frac{{\bar{r}}^{2}}{3}, \bar{q} = \frac{2 {\bar{r}}^{3}}{27} - \frac{\bar{s} \bar{r}}{3} + \bar{t}, \\ \bar{r} & = & \frac{\bar{B}}{\bar{A}}, \bar{s} = \frac{\bar{C}}{\bar{A}}, \bar{t} = \frac{\bar{D}}{\bar{A}} . \end{matrix}

(11)

A prerequisite for the existence of a real solution in the form given by Equation (11) is the non-negativity of the radicand

\bar{Δ}

, the fulfillment of which is ensured throughout the

p - T

range analyzed here.

To obtain the second guess of x, the mole fraction of low-density water is decomposed into the first guess

\bar{x}

and an increment

x^{'}

:

x ≅ \bar{x} + x^{'} .

(12)

Inserting Equation (12) into Equation (8) one obtains

\frac{\bar{x} + x^{'}}{1 - \bar{x} - x^{'}} = e^{α \bar{x} + β} e^{α x^{'}} .

(13)

Employing again the quadratization

e^{α x^{'}} \approx 1 + α x^{'} + \frac{1}{2} α^{2} x^{' 2}, x^{'} ≪ 1,

(14)

Equation (13) can be rearranged into a cubic equation for the determination of

x^{'}

:

\begin{matrix} A^{'} x^{' 3} + B^{'} x^{' 2} + C^{'} x^{'} + D^{'} = 0, \\ A^{'} = - \frac{α^{2}}{2} e^{α \bar{x} + β} = \bar{A} e^{α \bar{x}}, \\ B^{'} = α [\frac{α}{2} (1 - \bar{x}) - 1] e^{α \bar{x} + β} = A^{'} (\frac{\bar{B}}{\bar{A}} + \frac{1}{α} ln \frac{A^{'}}{\bar{A}}), \\ C^{'} = [α (1 - \bar{x}) - 1] e^{α \bar{x} + β} - 1 = \frac{A^{'}}{\bar{A}} (1 + \bar{C} - \bar{D} ln \frac{A^{'}}{\bar{A}}) - 1, \\ D^{'} = (1 - \bar{x}) e^{α \bar{x} + β} - \bar{x} = - \frac{1}{α} (1 + \frac{\bar{D} A^{'}}{\bar{A}}) ln \frac{A^{'}}{\bar{A}} + \frac{\bar{D} A^{'}}{\bar{A}} . \end{matrix}

(15)

The parameters

A^{'}

,

B^{'}

,

C^{'}

, and

D^{'}

are well-defined functions of

τ

and

π

. Equation (15) is of the same structure as Equation (10) and has three real and/or complex solutions,

x_{i}^{'} = x_{i}^{'} (A^{'}, B^{'}, C^{'}, D^{'})

,

i = 1, 2, 3

, representing the residual value of the low-density water mole fraction. Again, we are only interested in the real solution of Equation (15), which reads (cf. Section SM-4, Equation (SM-4.10)):

\begin{matrix} x^{'} & = & \sqrt[3]{- \frac{q^{'}}{2} + \sqrt{Δ^{'}}} + \sqrt[3]{- \frac{q^{'}}{2} - \sqrt{Δ^{'}}}, \\ Δ^{'} & = & {(\frac{q^{'}}{2})}^{2} + {(\frac{p^{'}}{3})}^{3} \geq 0, \\ p^{'} & = & s^{'} - \frac{r^{' 2}}{3}, q^{'} = \frac{2 r^{' 3}}{27} - \frac{s^{'} r^{'}}{3} + t^{'}, \\ r^{'} & = & \frac{B^{'}}{A^{'}}, s^{'} = \frac{C^{'}}{A^{'}}, t^{'} = \frac{D^{'}}{A^{'}} . \end{matrix}

(16)

A prerequisite for the existence of a real solution in the form given by Equation (16) is the non-negativity of the radicand

Δ^{'}

, the fulfillment of which is also ensured throughout the

p - T

range analyzed here.

Knowing the first guess and the incremental value of the low-density water mole fraction,

\bar{x}

and

x^{'}

, the second guess of x is given by Equation (12). The analytical solution given by this equation can be considered as a two-step predictor-corrector approximation.

3.3. Analyses of the Pure States of Water (Limiting Cases)

The limiting cases of unary water, i.e., of water consisting either of pure low-density water with mole fraction x = 1 or of pure high-density water with mole fraction x = 0, set constraints on the ordering field

L (τ, π)

and on the interaction parameter

ω (π)

. In order to determine whether these constraints can be achieved in the definition range of the EoS, we start our consideration with Equations (7) and (8):

f (x) = ln (\frac{x}{1 - x}) = α x + β, α = 2 ω, β = - (L + ω) .

(17)

The transition to the high-density state of water requires the fulfillment of the following constraint:

$lim_{x \to 0} f (x) = - (L + ω) = - \infty .$

(18)

Because of the finiteness of the interaction parameter, which is restricted to the interval $2 \leq ω (π) \leq ω_{\max} (π_{\max}) ≪ \infty$ , the fulfillment of the constraint given by Equation (18) requires the ordering field approaching $L (τ, π) \to \infty$ .
Analogously, the transition to the low-density state of water requires the fulfillment of the following constraint:

$lim_{x \to 1} f (x) = - (L - ω) = \infty .$

(19)

Because of the restriction of the interaction parameter to the interval $2 \leq ω (π) \leq ω_{\max} (π_{\max}) ≪ \infty$ , the fulfillment of the constraint given by Equation (19) requires the ordering field approaching $L (τ, π) \to - \infty$ .

Figure 1 displays the ordering field

L (τ, π)

according to the supplementary material, Section SM-3, Equation (SM-3.5) in the definition range

0 \leq τ \leq 10

and

0 \leq π \leq 10

imbedding the validity subrange of the EoS.

It can be seen that

L (τ, π)

is a positive definite and finite quantity (no singularities) in the EoS validity range. As a consequence, the pure states of either low-density or high-density water corresponding to

x = 0

or

x = 1

are thermodynamically not accessible in the validity range of the EoS.

4. Application to the IAPWS G12-15 Formulation

4.1. Computer Program Verification

The IAPWS G12-15 [36] formulation, as described in Section SM-3, serves as a reference calculus. To ensure the correct implementation of the calculus, we have compared our results with the table values provided by IAPWS G12-15 [36] (Table 5 therein) for computer program verification, presented here in Section SM-3.5, Table SM-3.5.1. Equation (6) was numerically solved with consideration of the constraints on the subintervals for the numerical determination of the equilibrium mole fraction of low-density water,

x_{e}

, specified in Section SM-3.4.

For the numerical root finding of the transcendental Equation (6) we used the standard FORTRAN routine SUBROUTINE zbrak(fx,x1,x2,n,xb1,xb2,nb) from Press et al. [39] (Section 9.1 therein) for root bracketing employing the intermediate value theorem. The quantity fx(x) is given by the function

F (x_{e})

defined by Equation (6) on the interval from x1 to x2 which is subdivided in n equally spaced segments. The root-finding algorithm searches for the zero crossings of the function fx(x) = 0. The integer nb is an input variable and specifies the maximum number of roots sought, and is reset to the number of bracketing pairs xb1(1:nb), xb2(1:nb) that are found. Depending on the parameters

ω (π)

and

L (τ, π)

, the search for the root was executed separately in the three different intervals [x1,x2] defined in Section SM-3.4, Table SM-3.4.1. Each of these three search intervals was subdivided into n subintervals. In the employed standard routine, the variable n is declared as of type INTEGER with an allowed maximum value of

2.147.483.647

. In the present case, we set n =

10^{9}

for the reference solution, allowing

x_{e}

to be accurate to

10^{- 9}

.

Section SM-5, Table SM-5.1 shows the deviation of the numerically determined mass density

\hat{ϱ}

, the thermal expansion coefficient

α_{p}

, the isothermal compressibility

κ_{T}

, the isobaric heat capacity

c_{p}

, the speed of sound w, the equilibrium mole fraction of low-density water

x_{e}

, and the ordering field L using the numerical root finder SUBROUTINE zbrak(fx,x1,x2,n,xb1,xb2,nb) with n =

10^{9}

from the corresponding IAPWS G12-15 [36] reference values (subscript 🟉) provided in Section SM-3.5, Table SM-3.5.1. Except for the signed quantity

α_{p}

all differences are given as relative deviations in units of parts per billion (ppb). The absolute values of the relative deviations for the equilibrium mole fraction of low-density water,

x_{e}

, and for the mass density,

\hat{ϱ}

, are rather small (

< 4 ppb

). The quantity L (last column) is a function of

τ

and

π

only but not of

x_{e}

, i.e., the recognized relative deviation of L from its references value originates from the specific way of implementation. Note, to reproduce the values in Section SM-3.5, Table SM-3.5.1 to the number of digits given, the solution for the fraction

x_{e}

should be accurate to

10^{- 10}

. However, the recognized deviations from the reference values are extremely small and originate from error propagation as a consequence of the somewhat lower accuracy of the calculated

x_{e}

values (employing n =

10^{9}

). Therefore, the implementation of the IAPWS G12-15 [36] formulation is correct. For comparison, Table SM-5.1 also displays the results for n =

10^{k}

,

k

=

8, \dots, 1

. In general, upon decreasing n the deviations remarkably increase. For n

= 10

the relative deviation of

x_{e}

increases to ≈21%, and of

\hat{ϱ}

to ≈1.2% at atmospheric pressure. On the base of the results presented in Section SM-5, Table SM-5.1, one can state the correct implementation of the numerical solver of the IAPWS G12-15 [36] formulation (for benchmark).

4.2. Deviations of the Analytical from the IAPWS G12-15 Reference Formulation

Section SM-6, Table SM-6.1 displays the same type of information as Section SM-5, Table SM-5.1 but now referring to the deviation of the analytically calculated thermodynamic values from the corresponding IAPWS G12-15 [36] reference values shown in Section SM-3, Table SM-3.5.1. The presented results are for n =

10^{9}

. The analysis revealed the fulfillment of the constraints

\bar{Δ}

≥ 0 and

Δ^{'}

≥ 0 in the

T - p

range, i.e., the analytical solution is given by Equation (12) with Cardano’s solutions (11) and (16). The deviations were found to be rather small except for deeply subcooled water at T =

235.15 K

and p =

0.101325 MPa

, where the largest deviation from the IAPWS G12-15 [36] reference values occurs. At these conditions the absolute values of the relative deviations amount

\approx 1.6 ‰

for the low-density water fraction,

\approx 2.6 ‰

for the isothermal compressibility,

\approx 2.3 ‰

for the isobaric heat capacity, and

\approx 0.06 ‰

for the speed of sound, respectively. The absolute deviation for thermal expansion coefficient amounts

0.15 \cdot 10^{- 4} K^{- 1}

at

T = 225.15 K

, resulting in a relative deviation of

\approx 5.1 ‰

by virtue of a thermal expansion coefficient of

\approx - 29.634 \cdot 10^{- 4} K^{- 1}

(cf. Section SM-3.5, Table SM-3.5.1).

4.3. Performance of the Analytical form of the EoS

In the following, the IAPWS G12-15 [36] formulation with use of the root finder SUBROUTINE zbrak(fx,x1,x2,n,xb1,xb2,nb) for the determination of the low-density water fraction will be called “numerical form of the EoS” (thermodynamic properties subscripted with “num”), and the IAPWS G12-15 [36] formulation using the analytical solution for

x_{e}

presented in Section 3.2 will be called “analytical form of the EoS” (values without subscript). The results are presented in form of table values for the deviations of the analytical form from the numerical form in different ranges of p and T. These ranges were chosen according to the availability of experimental data of the different thermodynamic quantities predicted by the IAPWS G12-15 [36] formulation. A comprehensive comparison of the IAPWS G12-15 [36] predictions with experimental data has been performed previously by Holten et al. [23]. The aim of the present study is to assess the accuracy of the analytical form by intercomparison of the differences between the analytical and numerical forms of the EoS with the differences between the numerical form and the experimental data analyzed by Holten et al. [23]. If the absolute value of the former difference is lower than or equal to the absolute value of the latter one, then the approximative form can be considered acceptable.

The results of the analysis on the performance of the approximative form are presented in Section SM-7 and comprise the following tables.

Section SM-7.1: Mass density
-
Table SM-7.1.1: Relative deviation, $(\hat{ϱ} - {\hat{ϱ}}_{num}) / {\hat{ϱ}}_{num}$ in parts per billion $(ppb)$ , of the mass density $\hat{ϱ}$ using the analytically determined low-density water fraction from the mass density ${\hat{ϱ}}_{num}$ using the numerically determined low-density water fraction at p = $0.101325 MPa$ .
-
Table SM-7.1.2: As Table SM-7.1.1, but for $0.1$ ≤ $p / MPa \leq 400$ and $253.15 \leq T / K \leq 303.15$ .
-
Table SM-7.1.3: As Table SM-7.1.1, but for $0.1 \leq p / MPa \leq 1000$ and $235.15 \leq T / K \leq 303.15$ .
-
Table SM-7.1.4: As Table SM-7.1.1, but for $200 \leq p / hPa \leq 1000$ and $235.15 \leq T / K \leq 300.15$ .
Section SM-7.2: Thermal expansivity
-
Table SM-7.2.1: Deviation, $(α_{p} - α_{p, num}) / (10^{- 4} K^{- 1})$ , of the thermal expansivity $α_{p}$ using the analytically determined low-density water fraction from the thermal expansivity $α_{p, num}$ using the numerically determined low-density water fraction for $0.1 \leq p / MPa \leq 600$ and $245.5 \leq T / K \leq 288.0$ .
Section SM-7.3: Isothermal compressibility
-
Table SM-7.3.1: Relative deviation, $(κ_{T} - κ_{T, num}) / κ_{T, num}$ in parts per billion $(ppb)$ , of the isothermal compressibility $κ_{T}$ using the analytically determined low-density water fraction from the thermal expansivity $κ_{T, num}$ using the numerically determined low-density water fraction for $0.1 \leq p / MPa \leq 190$ and $235.15 \leq T / K \leq 300.15$ .
Section SM-7.4: Isobaric heat capacity
-
Table SM-7.4.1: Relative deviation, $({\hat{c}}_{p} - {\hat{c}}_{p, num}) / {\hat{c}}_{p, num}$ in parts per billion $(ppb)$ , of the isobaric heat capacity ${\hat{c}}_{p}$ using the analytically determined low-density water fraction from the isobaric heat capacity ${\hat{c}}_{p, num}$ using the numerically determined low-density water fraction for $0.1 \leq p / MPa \leq 190$ and $235.15 \leq T / K \leq 300.15$ .
-
Table SM-7.4.2: As Table SM-7.4.1, but for $200 \leq p / hPa \leq 1000$ and $235.15 \leq T / K \leq 300.15$ .
Section SM-7.5: Sound speed
-
Table SM-7.5.1: Relative deviation, $(w - w_{num}) / w_{num}$ in parts per billion $(ppb)$ , of the sound speed w using the analytically determined low-density water fraction from the sound speed $w_{num}$ using the numerically determined low-density water fraction for $0.1 \leq p / MPa \leq 400$ and $253.15 \leq T / K \leq 303.15$ .
-
Table SM-7.5.2: As in Table SM-7.5.1, but for $0.1 \leq p / MPa \leq 1000$ and $273.15 \leq T / K \leq 303.15$ .
Section SM-7.6: Gibbs energy and entropy
-
Table SM-7.6.1: Relative deviation, $(\hat{g} - {\hat{g}}_{num}) / {\hat{g}}_{num}$ in units of parts per billion $(ppb)$ , of the specific Gibbs energy $\hat{g}$ using the analytically determined low-density water fraction (Section 3.2) from the specific Gibbs energy ${\hat{g}}_{num}$ using the numerically determined low-density water fraction for $200 \leq p / hPa \leq 1000$ and $235.15 \leq T / K \leq 300.15$ .
-
Table SM-7.6.2: Relative deviation, $(\hat{s} - {\hat{s}}_{num}) / {\hat{s}}_{num}$ in parts per billion $(ppb)$ , of the specific entropy $\hat{s}$ using the analytically determined low-density water fraction from the specific entropy ${\hat{s}}_{num}$ using the numerically determined low-density water fraction for $200 \leq p / hPa \leq 1000$ and $235.15 \leq T / K \leq 300.15$ .

5. Discussion

Table 1 displays the summary of the evaluation of the table values in Section SM-7. Column C1 contains the running number of the table entry, columns C2–C3 show the specification and symbol of the quantity of interest; column C4 displays the corresponding table number in the supplementary material; columns C5–C7 and C8–C10 present the sample interval and resolution of pressure and temperature, respectively; columns C11–C12 display the maxima of the relative and absolute deviations, respectively, of the analytically determined quantity X from its numerical reference value

X_{num}

; and columns C13–C14 contain the pressure and temperature at which the maximum deviations occur. As a general tendency, the deviation was found to increase upon decreasing pressure and temperature. The largest relative deviation over all quantities and throughout the analyzed pressure and temperature ranges occurs for the isothermal compressibility

κ_{T}

, amounting

- 2.58 ‰

at p≈

0.1 MPa

and

T = 235.15 K

. Starting from the values at the maximum deviation, the deviation rapidly decreases to the parts-per-million range upon increasing pressure and temperature.

In order to assess the applicability of the two-step predictor-corrector scheme proposed here, the deviations reported in Table 1 (misfit of the analytical form of the EoS) are compared with the deviations of experimental data from the values calculated using the IAPWS G12-15 [36] formulation (misprediction of the numerical form of the EoS) as reported in Holten et al. [23]. The analysis delivered the following results.

Mass density:
-
Table 1, entry (1) (atmospheric pressure): The maximum relative deviation of the analytically determined mass density from the numerical one amounts to $0.172 ‰$ at p= $0.101325 MPa$ and T= $235.15 K$ . Holten et al. [23] (Figure 6 therein) found the experimental density data deviating from the IAPWS G12-15 [36] formulation by $\pm 0.2 ‰$ at p= $0.101325 MPa$ and T≤ $255 K$ . Hence, for the mass density the misfit of the analytical form of the EoS remains still smaller than the misprediction of the numerical form of the EoS.
-
Table 1, entry (2) (high pressure): The maximum of the relative deviation of the analytical EoS form from the numerical one amounts to $221 ppb$ at p= $0.1 MPa$ and T= $253.15 K$ . Holten et al. [23] (Figure 7 therein) found the experimental density data deviating from the IAPWS G12-15 [36] formulation by more than $- 0.1 %$ at T= $253.15 K$ and p≥ 200 $MPa$ . The deviation at p= $0 MPa$ cannot be resolved at the chosen scale (percent). However, one can expect that in the specified pressure and temperature range the misfit of the analytical form of the EoS does not relevantly exceed the misprediction of the numerical form of the EoS.
-
Table 1, entry (3) (very high pressure): In the specified pressure and temperature range, the maximum of the relative deviation of the analytically determined mass density from the numerical one amounts $0.79 ppb$ at p= $0.1 MPa$ and $T = 273.15 K$ . Throughout most of the $(T, p)$ space depicted in Holten et al. [23] (Figure 21 therein) this value is much smaller than the deviations of the experimental density from the IAPWS G12-15 [36] formulation. At atmospheric pressure, the latter, however, cannot be resolved in Holten et al. [23] (Figure 21 therein) at the chosen scale (percent). Thus, also for the very high pressure range one can safely conclude that the misfit of the analytical form of the EoS does not exceed the misprediction of the numerical form of the EoS.
Thermal expansivity:
-
Table 1, entry (5): The maximum of the deviation of the analytically determined expansivity from the numerical one amounts $0.11 \cdot 10^{- 6} K^{- 1}$ at p= $0.1 MPa$ and T= $245.5 K$ . At these conditions the deviation of the experimental values from the IAPWS G12-15 [36] formulation reported in Holten et al. [23] (Figure 9 therein) amounts to $- 0.2 \cdot 10^{- 4} K^{- 1}$ , i.e., the misfit of the analytical form of the EoS is much smaller than the misprediction of the numerical form of the EoS (see also Holten et al. [23], Figure 10 therein).
Isothermal compressibility:
-
Table 1, entry (6): The maximum of the deviation of the analytically determined compressibility from the numerical one amounts to $- 2.58 ‰$ at p= $0.101325 MPa$ and $T = 235.15 K$ . At these conditions the deviation of the experimental values from the IAPWS G12-15 [36] formulation reported in Holten et al. [23] (Figure 15 therein) amounts to $\approx 6 %$ , i.e., the misfit of the analytical form of the EoS is much smaller than the misprediction of the numerical form of the EoS (see also Holten et al. [23], Figure 16 therein).
Isobaric heat capacity:
-
Table 1, entry (7): The maximum of the deviation of the analytically determined heat capacity from the numerical one amounts to $- 2.35 ‰$ at p= $0.101325 MPa$ and $T = 235.15 K$ . At these conditions, the deviation of the experimental values from the IAPWS G12-15 [36] formulation reported in Holten et al. [23] (Figure 19 therein) amounts up to $- 8 %$ , i.e., the misfit of the analytical form of the EoS is much smaller than the misprediction of the numerical form of the EoS.
Sound speed:
-
Table 1, entry (9) (high pressure): The maximum of the deviation of the analytically determined sound speed from the numerical one over the specified pressure and temperature ranges amounts to $821 ppb$ at p= $0.1 MPa$ and T= $253.15 K$ . At these conditions the deviation of the experimental values from the IAPWS G12-15 [36] formulation reported in Holten et al. [23] (Figure 18 therein) amounts up to $- 0.2 %$ , i.e., the misfit of the analytical form of the EoS is much smaller than the misprediction of the numerical form of the EoS (see also Holten et al. [23], Figure 17 therein).
-
Table 1, entry (10) (very high pressure): The maximum of the deviation of the analytically determined sound speed from the numerical one over the specified pressure and temperature ranges amounts to $4 ppb$ at p= $0.1 MPa$ and T= $273.15 K$ . At these conditions, the deviation of the experimental values from the IAPWS G12-15 [36] formulation reported in Holten et al. [23] (Figure 22 therein) amounts to $\leq 0.1 %$ , i.e., the misfit of the analytical form of the EoS remains much smaller than the misprediction of the numerical form of the EoS.

For the Gibbs free energy and the entropy (entries (11) and (12)) no experimental data for direct comparison are available.

6. Conclusions

Possible implications of the two-state nature of subcooled water for ice crystallization have been discussed. A new two-step predictor-corrector method for the approximative determination of the low-density water fraction of subcooled water has been developed which allows a highly accurate determination of the specific Gibbs energy of subcooled water at given pressure and temperature. The misfit of this analytical solution against the exact numerical solution was demonstrated to be smaller than the misprediction of the original IAPWS G12-15 [36] formulation. The solution proposed here is intended to be used for highly frequent calls of the IAPWS G12-15 [36] formulation in CFD applications. Appendix A contains the algorithm for the calculation of the equilibrium mole fraction of low-density water in subcooled water and the thermodynamic reference values for the check of the correct computer implementation.

The calculus employed here is based on a polynomialization of the underlying transdendental equation (Equations (6) and (8)) by expansion of its exponential term into a Taylor series up to the quadratic term (Equation (9)). This resulted in a cubic equation (Equation (10)), for which an analytical solution formula exists. According to the Galois theory, the existence of a general solution formula for root finding (employing the four basic arithmetic operations and root extraction only) is restricted to polynomials with a degree of up to four, i.e., quintic and higher-degree equations must be solved numerically. Therefore, if the underlying transcendental equation can be taylorized into a polynomial with a degree of lower than or equal to four one can derive analytical solutions. The approach presented here is restricted to cubic equations, i.e., it can be applied to any other EoS of water (or other fluids) provided this equation can be taylorized into a cubic polynomial. For a quaternary polynomial a more generalized solution formula is available.

Supplementary Materials

The document “Supplementary Material to ‘Analytical Determination of the Nucleation-Prone, Low-Density Fraction of Subcooled Water’ ” (SM) is available online at https://www.mdpi.com/1099-4300/22/9/933/s1.

Author Contributions

The authors contributed to the present work as follows: conceptualization by O.H. and R.F; software by O.H. and R.F.; formal analysis and investigation by O.H. and R.F.; original draft preparation by O.H.; idea outlined in Section SM-2 by R.F.; review and editing, O.H. and R.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The publication of this article was funded by the Open Access Fund of the Leibniz Association.

Acknowledgments

The contribution of O. Hellmuth was provided within the framework of the research theme 1 “Aerosols: Process studies at small temporal and spatial scales” of TROPOS Leibniz Institute for Tropospheric Research, Leipzig. Special thanks go to the three anonymous referees for their fair evaluation and helpful comments to improve the manuscript. This paper contributes to the tasks of the IAPWS/SCOR/IAPSO Joint Committee on Seawater (JCS).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Algorithm for the Calculation of the Equilibrium Mole Fraction of Low-Density Water in Subcooled Water and Thermodynamic Reference Values for the Check of the Correct Computer Implementation

Determination of the reduced temperature $τ$ and the reduced pressure $π$ for given temperature and pressure p according to Equation (5):

$τ (T) = \frac{T}{T_{LL}} - 1, π (p) = \frac{p}{{\hat{ϱ}}_{0} R_{W} T_{LL}} .$

(A1)

Here, $T_{LL}$ denotes the temperature of the liquid–liquid critical point, $R_{W}$ denotes the specific gas constant of water, and ${\hat{ϱ}}_{0}$ denotes the reference mass density of subcooled water. The values of $T_{LL}$ , ${\hat{ϱ}}_{0}$ , and $R_{W}$ are given in Section SM-3.1, Table SM-3.1.1.
Determination of the order field $L (τ, π)$ and of the dimensionless interaction parameter $ω (π)$ according to Section SM-3.1, Equations (SM-3.5) and (SM-3.6):

$\begin{matrix} L (τ, π) & = & L_{0} (\frac{K_{2}}{2 k_{1} k_{2}}) [1 + k_{0} k_{2} + k_{1} (π + k_{2} τ) - K_{1} (τ, π))], \\ K_{1} (τ, π) & = & \sqrt{{[1 + k_{0} k_{2} + k_{1} (π - k_{2} τ)]}^{2} - 4 k_{0} k_{1} k_{2} (π - k_{2} τ)}, \\ K_{2} & = & \sqrt{1 + k_{2}^{2}}, \end{matrix}$

(A2)

$ω (π) = 2 + ω_{0} π .$

(A3)

The values of $L_{0}$ , $k_{0}$ , $k_{1}$ , $k_{2}$ , and $ω_{0}$ are given in Section SM-3.1, Table SM-3.1.1.
Determination of the first guess (predictor value) of the low-density mole fraction of water in subcooled water, ${\bar{x}}_{e}$ :

$\begin{matrix} α (π) & = & 2 ω (π), β (τ, π) = - (L (τ, π) + ω (π)), \\ \bar{D} & = & e^{β}, \bar{A} = - \frac{α^{2}}{2} \bar{D}, \bar{B} = α (\frac{α}{2} - 1) \bar{D}, \bar{C} = \bar{D} (α - 1) - 1, \\ \bar{r} & = & \frac{\bar{B}}{\bar{A}}, \bar{s} = \frac{\bar{C}}{\bar{A}}, \bar{t} = \frac{\bar{D}}{\bar{A}}, \\ \bar{p} & = & \bar{s} - \frac{{\bar{r}}^{2}}{3}, \bar{q} = \frac{2 {\bar{r}}^{3}}{27} - \frac{\bar{s} \bar{r}}{3} + \bar{t}, \\ \bar{Δ} & = & {(\frac{\bar{q}}{2})}^{2} + {(\frac{\bar{p}}{3})}^{3} \geq 0, \\ {\bar{x}}_{e} & = & \sqrt[3]{- \frac{\bar{q}}{2} + \sqrt{\bar{Δ}}} + \sqrt[3]{- \frac{\bar{q}}{2} - \sqrt{\bar{Δ}}} . \end{matrix}$

(A4)
Determination of the corrector value of the low-density mole fraction of water in subcooled water, $x_{e}^{'}$ :

$\begin{matrix} A^{'} & = & \bar{A} e^{α {\bar{x}}_{e}}, \\ B^{'} & = & A^{'} (\frac{\bar{B}}{\bar{A}} + \frac{1}{α} ln \frac{A^{'}}{\bar{A}}), \\ C^{'} & = & \frac{A^{'}}{\bar{A}} (1 + \bar{C} - \bar{D} ln \frac{A^{'}}{\bar{A}}) - 1, \\ D^{'} & = & - \frac{1}{α} (1 + \frac{\bar{D} A^{'}}{\bar{A}}) ln \frac{A^{'}}{\bar{A}} + \frac{\bar{D} A^{'}}{\bar{A}}, \\ r^{'} & = & \frac{B^{'}}{A^{'}}, s^{'} = \frac{C^{'}}{A^{'}}, t^{'} = \frac{D^{'}}{A^{'}}, \\ p^{'} & = & s^{'} - \frac{r^{' 2}}{3}, q^{'} = \frac{2 r^{' 3}}{27} - \frac{s^{'} r^{'}}{3} + t^{'}, \\ Δ^{'} & = & {(\frac{q^{'}}{2})}^{2} + {(\frac{p^{'}}{3})}^{3} \geq 0, \\ x_{e}^{'} & = & \sqrt[3]{- \frac{q^{'}}{2} + \sqrt{Δ^{'}}} + \sqrt[3]{- \frac{q^{'}}{2} - \sqrt{Δ^{'}}} . \end{matrix}$

(A5)
Determination of the second guess of the low-density mole fraction of water in subcooled water, $x_{e}$ :

$x_{e} ≅ {\bar{x}}_{e} + x_{e}^{'} .$

(A6)
Determination of the thermodynamic properties according to Section SM-3.2.

Table A1 contains the thermodynamic reference values for the check of the correct computer implementation of the analytical solution.

Table A1. Thermodynamic reference values for the check of the correct computer implementation of the calculus presented in Equations (A1)–(A6).

Quantity	Unit	$(\frac{T}{K}; \frac{p}{MPa})$
Quantity	Unit	$(273.15; 0.101325)$	$(235.15; 0.101325)$	$(250.0; 200.0)$	$(200.0; 400.0)$	$(250.0; 400.0)$
$α$	1	$4.0009272104686762$	$4.0009272104686762$	$5.8301711693577012$	$7.6603423387154033$	$7.6603423387154033$
$β$	1	$- 2.6216683426000236$	$- 2.0922272810837832$	$- 3.6388563972913954$	$- 4.9855677030124106$	$- 5.2646856527278016$
${\bar{x}}_{e}$	1	$9.5712687098620797 \cdot 10^{- 2}$	$0.21690329030435798$	$3.0400304327358713 \cdot 10^{- 2}$	$7.1698813902603409 \cdot 10^{- 3}$	$5.3587744916567814 \cdot 10^{- 3}$
$x_{e}^{'}$	1	$9.4202720002584117 \cdot 10^{- 4}$	$3.7792425779540201 \cdot 10^{- 2}$	$2.8962408084282076 \cdot 10^{- 5}$	$1.9953324550958662 \cdot 10^{- 7}$	$6.2126093464121013 \cdot 10^{- 8}$
$x_{e}$	1	$9.6654714298646638 \cdot 10^{- 2}$	$0.25469571608389818$	$3.0429266735442995 \cdot 10^{- 2}$	$7.1700809235058505 \cdot 10^{- 3}$	$5.3588366177502456 \cdot 10^{- 3}$
$\hat{ϱ}$	$kg m^{- 3}$	$999.84229220373140$	$968.26695400651727$	$1090.4567679510335$	$1185.0279968723719$	$1151.7151726673426$
$α_{p}$	$K^{- 1}$	$- 6.8304230971614192 \cdot 10^{- 5}$	$- 2.9486489045737296 \cdot 10^{- 3}$	$3.2677681131880254 \cdot 10^{- 4}$	$6.7160085822465037 \cdot 10^{- 4}$	$4.9299273113608325 \cdot 10^{- 4}$
$κ_{T}$	${Pa}^{- 1}$	$5.0884994182031361 \cdot 10^{- 10}$	$1.1550865210524601 \cdot 10^{- 9}$	$3.3613108478638421 \cdot 10^{- 10}$	$2.5672366521284537 \cdot 10^{- 10}$	$2.2770292496267704 \cdot 10^{- 10}$
${\hat{c}}_{p}$	$J {kg}^{- 1} K^{- 1}$	$4218.3001572845242$	$5983.4804702453448$	$3708.3901593987844$	$3338.5249744098869$	$3757.2144360129914$
w	$m s^{- 1}$	$1402.3885702265868$	$1134.6519596636447$	$1668.2019636972875$	$1899.3293996688064$	$2015.8781588494110$
$\hat{g}$	$J {kg}^{- 1}$	$101.34282192560295$	$- 12017.427396385123$	$187027.03827476871$	$323271.03039053129$	$365199.08446200017$
$\hat{s}$	$J {kg}^{- 1} K^{- 1}$	$- 0.14754456408819105$	$- 685.64897594084539$	$- 370.50319201333610$	$- 1244.6317196037739$	$- 446.48930388307758$

References

Meyers, M.P.; DeMott, P.J.; Cotton, W.R. New primary ice-nucleation parameterizations in an explicit cloud model. J. Appl. Meteorol. 1992, 31, 708–721. [Google Scholar] [CrossRef] [Green Version]
Khvorostyanov, V.I.; Sassen, K. Toward the theory of homogeneous ice nucleation and its parameterization for cloud models. Geophys. Res. Lett. 1998, 25, 3155–3158. [Google Scholar] [CrossRef]
Lohmann, U.; Kärcher, B. First interactive simulations of cirrus clouds formed by homogeneous freezing in the ECHAM general circulation model. J. Geophys. Res. 2002, 107, AAC 8-1–AAC 8-13. [Google Scholar] [CrossRef]
Lohmann, U.; Kärcher, B.; Timmreck, C. Impact of the Mount Pinatubo eruption on cirrus clouds formed by homogeneous freezing in the ECHAM4 GCM. J. Geophys. Res. 2003, 108. [Google Scholar] [CrossRef] [Green Version]
Pruppacher, H.R.; Klett, J.D. Microphysics of Clouds and Precipitation; Kluwer Academic Publishers: Dordrecht, The Netherlands; Boston, MA, USA; London, UK, 2004; p. 954. [Google Scholar]
Heymsfield, A.J.; Miloshevich, L.M.; Schmitt, C.; Bansemer, A.; Twohy, C.; Poellot, M.R.; Fridlind, A.; Gerber, H. Homogeneous ice nucleation in suptropical and tropical convection and its influence on cirrus anvil microphysics. J. Atmos. Sci. 2005, 62, 41–64. [Google Scholar] [CrossRef] [Green Version]
Jensen, E.J.; Ackerman, A.S. Homogeneous aerosol freezing in the tops of high-altitude tropical cumulonimbus clouds. Geophys. Res. Lett. 2006, 33, L08802. [Google Scholar] [CrossRef]
Barahona, D.; Nenes, A. Parameterization of cirrus cloud formation in large-scale models: Homogeneous nucleation. J. Geophys. Res. 2008, 113, D11211. [Google Scholar] [CrossRef]
Jensen, E.J.; Pfister, L.; Bui, T.V.; Lawson, P.; Baker, B.; Mo, Q.; Baumgardner, D.; Weinstock, E.M.; Smith, J.B.; Moyer, E.J.; et al. Formation of large (≃100 μm) ice crystals near the tropical tropopause. Atmos. Chem. Phys. 2008, 8, 1621–1633. [Google Scholar] [CrossRef] [Green Version]
Zasetsky, A.Y.; Petelina, S.V.; Svishchev, I.M. Thermodynamics of homogeneous nucleation of ice particles in the polar summer mesosphere. Atmos. Chem. Phys. 2009, 9, 965–971. [Google Scholar] [CrossRef] [Green Version]
Khvorostyanov, V.I.; Curry, J.A. Critical humidities of homogeneous and heterogeneous ice nucleation: Inferences from extended classical nucleation theory. J. Geophys. Res. 2009, 114, D04307. [Google Scholar] [CrossRef]
Khvorostyanov, V.I.; Curry, J.A. Parameterization of homogeneous ice nucleation for cloud and climate models based on classical nucleation theory. Atmos. Chem. Phys. 2012, 12, 9275–9302. [Google Scholar] [CrossRef] [Green Version]
Hellmuth, O.; Khvorostyanov, V.I.; Curry, J.A.; Shchekin, A.K.; Schmelzer, J.W.P.; Feistel, R.; Djikaev, Y.S.; Baidakov, V.G. Selected aspects of atmospheric ice and salt crystallisation. In Nucleation Theory and Applications; Schmelzer, J.W.P., Hellmuth, O., Eds.; Special Issues; JINR Joint Institute for Nuclear Research, Bogoliubov Laboratory of Theoretical Physics: Dubna, Russia, 2013; Volume 1, p. 513. ISBN 978-5-9530-0349-0. Available online: http://theor.jinr.ru/meetings/2013/nta/ (accessed on 19 June 2020).
Khvorostyanov, V.I.; Curry, J.A. Thermodynamics, Kinetics, and Microphysics of Clouds, 1st ed.; Cambridge University Press: Cambridge, UK, 2014; p. 782. [Google Scholar]
Lohmann, U.; Lüönd, F.; Mahrt, F. An Introduction to Clouds. From Microscale to Climate; Cambridge University Press: Cambridge, UK, 2016. [Google Scholar]
Pegg, D.E. Principles of Cryopreservation. In Cryopreservation and Freeze-Drying Protocols. Methods in Molecular Biology; Day, J.G., Stacey, G.N., Eds.; Humana Press: Totowa, NJ, USA, 2007; Volume 368. [Google Scholar]
Espinosa, J.R.; Sanz, E.; Valeriani, C.; Vega, C. Homogeneous ice nucleation evaluated for several water models. J. Chem. Phys. 2014, 141, 180529. [Google Scholar] [CrossRef] [Green Version]
Espinosa, J.R.; Zaragoza, A.; Rosales-Pelaez, P.; Navarro, C.; Valeriani, C.; Vega, C.; Sanz, E. Interfacial free energy as the key to the pressure-induced deceleration of ice nucleation. Phys. Rev. Lett. 2016, 117, 135702. [Google Scholar] [CrossRef] [PubMed]
Debenedetti, P.G.; Stanley, H.E. Supercooled and glassy water. Phys. Today 2003, 56, 40–46. [Google Scholar] [CrossRef]
Bhat, S.N.; Sharma, A.; Bhat, S.V. Vitrification and Glass Transition of Water: Insights from Spin Probe ESR. Phys. Rev. Lett. 2005, 95, 235702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Zobrist, B.; Marcolli, C.; Pedernera, D.A.; Koop, T. Do atmospheric aerosols form glasses? Atmos. Chem. Phys. 2008, 8, 5221–5244. [Google Scholar] [CrossRef] [Green Version]
Hellmuth, O.; Schmelzer, J.; Feistel, R. Ice-crystal nucleation in water: Thermodynamic driving force and surface tension. Part I: Theoretical foundation. Entropy 2020, 22, 50. [Google Scholar] [CrossRef] [Green Version]
Holten, V.; Sengers, J.V.; Anisimov, M.A. Equation of state for supercooled water at pressures up to 400 MPa. J. Phys. Chem. Ref. Data 2014, 43, 043101. [Google Scholar] [CrossRef] [Green Version]
Feistel, R.; Wright, D.G.; Kretzschmar, H.-J.; Hagen, E.; Herrmann, S.; Span, R. Thermodynamic properties of sea air. Ocean Sci. 2010, 6, 91–141. [Google Scholar] [CrossRef] [Green Version]
Feistel, R.; Wright, D.G.; Jackett, D.R.; Miyagawa, K.; Reissmann, J.H.; Wagner, W.; Overhoff, U.; Guder, C.; Feistel, A.; Marion, G.M. Numerical implementation and oceanographic application of the thermodynamic potentials of liquid water, water vapour, ice, seawater and humid air—Part 1: Background and equations. Ocean Sci. 2010, 6, 633–677. [Google Scholar] [CrossRef] [Green Version]
Wright, D.G.; Feistel, R.; Reissmann, J.H.; Miyagawa, K.; Jackett, D.R.; Wagner, W.; Overhoff, U.; Guder, C.; Feistel, A.; Marion, G.M. Numerical implementation and oceanographic application of the thermodynamic potentials of liquid water, water vapour, ice, seawater and humid air—Part 2: The library routines. Ocean Sci. 2010, 6, 695–718. [Google Scholar] [CrossRef] [Green Version]
IOC, SCOR, and IAPSO. The International Thermodynamic Equation of Seawater—2010: Calculation and Use of Thermodynamic Properties. Technical Report, Intergovernmental Oceanographic Commission, Manuals and Guides No. 56; UNESCO (English), 2010; 196p. Available online: http://www.teos-10.org (accessed on 19 June 2020).
Feistel, R. TEOS-10: A new international oceanographic standard for seawater, ice, fluid water and humid air. Int. J. Thermophys. 2012, 33, 1335–1351. [Google Scholar] [CrossRef]
Feistel, R. Thermodynamic properties of seawater, ice and humid air: TEOS-10, before and beyond. Ocean Sci. 2018, 14, 471–502. [Google Scholar] [CrossRef] [Green Version]
Wagner, W.; Pruß, A. The IAPWS formulation 1995 for the thermodynamic properties of ordinary water substance for general and scientific use. J. Phys. Chem. Ref. Data 2002, 31, 387–535. [Google Scholar] [CrossRef] [Green Version]
IAPWS R6-95. Revised Release on the IAPWS Formulation 1995 for the Thermodynamic Properties of Ordinary Water Substance for General and Scientific Use; Technical Report; The International Association for the Properties of Water and Steam: Dresden, Germany, September 2016; Available online: http://www.iapws.org (accessed on 19 June 2020).
Holten, V.; Bertrand, C.E.; Anisimov, M.A.; Sengers, J.V. Thermodynamic modeling of supercooled water. In Technical Report for the International Association for the Properties of Water and Steam (IAPWS) (September 2011); Technical Report; Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland: College Park, MD, USA, 2011; p. 43. Available online: http://www.iapws.org (accessed on 19 June 2020).
Holten, V.; Bertrand, C.E.; Anisimov, M.A.; Sengers, J.V. Thermodynamics of supercooled water. J. Chem. Phys. 2012, 136, 094507. [Google Scholar] [CrossRef] [Green Version]
Holten, V.; Limmer, D.T.; Molinero, V.; Anisimov, A.M. Nature of the anomalies in the supercooled liquid state of the mW model of water. J. Chem. Phys. 2013, 138, 174501. [Google Scholar] [CrossRef]
Holten, V.; Anisimov, M. Entropy-driven liquid–liquid separation in supercooled water. Sci. Rep. 2012, 2, 713. [Google Scholar] [CrossRef] [Green Version]
IAPWS G12-15. Guideline on Thermodynamic Properties of Supercooled Water; Technical Report; The International Association for the Properties of Water and Steam: Stockholm, Sweden, July 2015; Available online: http://www.iapws.org (accessed on 19 June 2020).
Röntgen, W.C. Ueber die Constitution des flüssigen Wassers. Ann. Der Phys. 1892, 281, 91–97. [Google Scholar] [CrossRef] [Green Version]
Ostwald, W. Studien über die Bildung und Umwandlung fester Körper. 1. Abhandlung: Übersättigung und Überkaltung. Z. für Phys. Chem. 1897, 22, 289–330. [Google Scholar] [CrossRef]
Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P. Numerical Recipes in Fortran 77. The Art of Scientific Computing. Second Edition. Volume 1 of Fortran Numerical Recipes; Cambridge University Press: New York, NY, USA, 1996. [Google Scholar]

Figure 1. Ordering field

L (x

=

τ,

y =

π)

according to Section SM-3, Equation (SM-3.5) in the interval

0 \leq τ \leq 10

and

0 \leq π \leq 10

. Graphic plot using Wolfram Alpha LLC, 2020, https://www.wolframalpha.com/examples/mathematics/plotting-and-graphics/ (accessed on 22 July 2020).

Figure 1. Ordering field

L (x

=

τ,

y =

π)

according to Section SM-3, Equation (SM-3.5) in the interval

0 \leq τ \leq 10

and

0 \leq π \leq 10

. Graphic plot using Wolfram Alpha LLC, 2020, https://www.wolframalpha.com/examples/mathematics/plotting-and-graphics/ (accessed on 22 July 2020).

Table 1. Summary of the evaluation of the tables presented in Section SM-7. Column C1 contains the running number of the table entry, columns C2–C3 show the quantity of interest and its corresponding symbol; column C4 displays the corresponding table number in the supplementary document; columns C5–C7 and C8–C10 present the sample interval and the resolution of the pressure and temperature, respectively; columns C11–C12 display the maxima of the relative and absolute deviations, respectively, of the analytically determined quantity X from its numerical reference value

X_{num}

(misfit); and columns C13–C14 contain the pressure and temperature at which the maximum deviations occur.

Table 1. Summary of the evaluation of the tables presented in Section SM-7. Column C1 contains the running number of the table entry, columns C2–C3 show the quantity of interest and its corresponding symbol; column C4 displays the corresponding table number in the supplementary document; columns C5–C7 and C8–C10 present the sample interval and the resolution of the pressure and temperature, respectively; columns C11–C12 display the maxima of the relative and absolute deviations, respectively, of the analytically determined quantity X from its numerical reference value

X_{num}

(misfit); and columns C13–C14 contain the pressure and temperature at which the maximum deviations occur.

C1	C2	C3	C4	C5	C6	C7	C8	C9	C10	C11	C12	C13	C14
Entry	Quantity	Symbol $X$	Table	Pressure Range			Temperature Range			${(\frac{X - X_{num}}{X_{num}})}_{\max}$	${(\frac{X - X_{num}}{10^{- 4} K^{- 1}})}_{\max}$	$(p, T)$ at max. dev.
Entry	Quantity	Symbol $X$		$\frac{p_{\min}}{MPa}$	$\frac{p_{\max}}{MPa}$	$\frac{Δ p}{MPa}$	$\frac{T_{\min}}{K}$	$\frac{T_{\max}}{K}$	$\frac{Δ T}{K}$	ppb		$\frac{p}{MPa}$	$\frac{T}{K}$
(1)	Mass density	$\hat{ϱ}$	SM-7.1.1	$0.101325$			$235.15$	$300.15$	5	$0.172 \cdot 10^{6}$		$0.101325$	$235.15$
(2)			SM-7.1.2	$0.1$	400	5	$253.15$	$303.15$	$\approx 2 - 5$	221		$0.1$	$253.15$
(3)			SM-7.1.3	$0.1$	1000	50	$273.15$	$303.15$	$\approx 5 - 20$	$0.79$		$0.1$	$273.15$
(4)			SM-7.1.4	$0.02$	$0.1$	$0.01$	$235.15$	$300.15$	5	$0.174 \cdot 10^{6}$		$0.02$	$235.15$
(5)	Expansivity	$α_{p}$	SM-7.2.1	$0.1$	600	10	$245.5$	$288.0$	$5 - 10$		$0.11 \cdot 10^{- 2}$	$0.1$	$245.5$
(6)	Compressibility	$κ_{T}$	SM-7.3.1	$0.1$	190	$\approx 10 - 50$	$235.15$	$300.15$	5	$- 2.58 \cdot 10^{6}$		$0.101325$	$235.15$
(7)	Heat capacity	${\hat{c}}_{p}$	SM-7.4.1	$0.1$	190	$\approx 10 - 50$	$235.15$	$300.15$	5	$- 2.35 \cdot 10^{6}$		$0.101325$	$235.15$
(8)			SM-7.4.2	$0.02$	$0.1$	$0.01$	$235.15$	$300.15$	5	$- 2.37 \cdot 10^{6}$		$0.02$	$235.15$
(9)	Sound speed	w	SM-7.5.1	$0.1$	400	$\approx 10 - 50$	$253.15$	$303.15$	$\approx 2 - 5$	821		$0.1$	$253.15$
(10)			SM-7.5.2	$0.1$	1000	$\approx 50$	$273.15$	$303.15$	10	4		$0.1$	$273.15$
(11)	Gibbs energy	$\hat{g}$	SM-7.6.1	$0.02$	$0.1$	$0.01$	$235.15$	$300.15$	5	$- 955$		$0.02$	$235.15$
(12)	Entropy	$\hat{s}$	SM-7.6.2	$0.02$	$0.1$	$0.01$	$235.15$	$300.15$	5	$- 0.86 \cdot 10^{6}$		$0.02$	$235.15$

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Hellmuth, O.; Feistel, R. Analytical Determination of the Nucleation-Prone, Low-Density Fraction of Subcooled Water. Entropy 2020, 22, 933. https://doi.org/10.3390/e22090933

AMA Style

Hellmuth O, Feistel R. Analytical Determination of the Nucleation-Prone, Low-Density Fraction of Subcooled Water. Entropy. 2020; 22(9):933. https://doi.org/10.3390/e22090933

Chicago/Turabian Style

Hellmuth, Olaf, and Rainer Feistel. 2020. "Analytical Determination of the Nucleation-Prone, Low-Density Fraction of Subcooled Water" Entropy 22, no. 9: 933. https://doi.org/10.3390/e22090933

APA Style

Hellmuth, O., & Feistel, R. (2020). Analytical Determination of the Nucleation-Prone, Low-Density Fraction of Subcooled Water. Entropy, 22(9), 933. https://doi.org/10.3390/e22090933

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Analytical Determination of the Nucleation-Prone, Low-Density Fraction of Subcooled Water

Abstract

1. Introduction

2. Materials and Methods

3. Theoretical Background

3.1. On the Two-State Character of Subcooled Water

3.2. Analytical Determination of the Fraction of Low-Density Water in the Two-State Water Mixture

3.3. Analyses of the Pure States of Water (Limiting Cases)

4. Application to the IAPWS G12-15 Formulation

4.1. Computer Program Verification

4.2. Deviations of the Analytical from the IAPWS G12-15 Reference Formulation

4.3. Performance of the Analytical form of the EoS

5. Discussion

6. Conclusions

Supplementary Materials

Author Contributions

Funding

Acknowledgments

Conflicts of Interest

Appendix A. Algorithm for the Calculation of the Equilibrium Mole Fraction of Low-Density Water in Subcooled Water and Thermodynamic Reference Values for the Check of the Correct Computer Implementation

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI