View source for Helmholtz Energy Equations of State

__FORCETOC__


==Overview==
The accurate prediction of the thermodynamic properties of substances is critical in a broad swath of applications in engineering and scientific research.  An example application would be the simulation of pressure vessel failure due to fire exposure.  The accurate simulation of such a scenario would require not only a sophisticated multiscale model of the metal components, but also an equation of state that can accurately describe both the caloric properties and phase equilibria of the fluid within the vessel.  For a small number of important substances, reference equations have been developed that are capable of representing the best experimental physical property data within their reported uncertainties.  The current state of the art is to use an equation of state explicit in the Helmholtz energy to accomplish this goal.  For a pure component<ref name="name">Span, R. 2000. ''Multiparameter equations of state: an accurate source of thermodynamic property data.'' Berlin: Springer.</ref>:
<br>
<br>
<math>\alpha(\tau,\delta) = \alpha^0(\tau,\delta) + \alpha^r(\tau,\delta)</math>
<br>
<br>
<math>\tau=\frac{T_r}{T}</math>
<br>
<br>
<math>\delta=\frac{\rho}{\rho_r}</math>
<br>
<br>
<math>\alpha=\frac{a}{RT}</math>
<br>
<br>
Where <math>\alpha</math> is the reduced Helmholtz energy, <math>T</math> is the temperature, <math>T_r</math> is often (but not necessarily) the critical temperature, <math>\rho</math> is the density, <math>\rho_r</math> is often (but not necessarily) the critical density, <math>\alpha^0</math> is the ideal gas contribution to the reduced Helmholtz energy, <math>\alpha^r</math> is the residual contribution to the reduced Helmholtz energy, and <math>R</math> is the ideal gas constant.  An equation of state explicit in the Helmholtz energy has the advantage that all other thermodynamic properties may be written as a function of its derivatives.  Because it is easier to obtain analytic derivatives than analytic integrals, this allows for a larger number of terms that may be used in the optimization of the functional form.  This equation of state also allows one to utilize all available thermodynamic data in fitting the equation.  This includes data from properties such as saturated vapor pressure, heat of vaporization, liquid density, vapor density, liquid heat capacity, vapor heat capacity, and speed of sound.  The <math>\alpha^0</math> term is typically obtained from the integration of an equation fitted to ideal gas heat capacity data.  An example optimized functional form for the <math>\alpha^r</math> term is the following<ref name="name" /><ref name="name2">Xiang, H.W., and U.K. Deiters. 2008. "A new generalized corresponding-states equation of state for the extension of the Lee–Kesler equation to fluids consisting of polar and larger nonpolar molecules". ''Chemical Engineering Science.'' 63 (6): 1490-1496.</ref><ref>Kunz O., and Wagner W. 2012. "The GERG-2008 wide-range equation of state for natural gases and other mixtures: An expansion of GERG-2004". ''Journal of Chemical and Engineering Data.'' 57 (11): 3032-3091.</ref>:
<br>
<br>
<math>\alpha^r(\tau,\delta) = \sum_{m=1}^{I_{pol}} c_m \delta^{i_m} \tau^{j_m} + \sum_{m=I_{pol}+1}^{I_{pol}+I_{exp}} c_m \delta^{i_m} \tau^{j_m} e^{\delta^{k_m}}</math>
<br>
<br>
Where <math>I_{pol}</math>, <math>I_{exp}</math>, <math>c_m</math>, <math>i_m</math>, <math>j_m</math>, and <math>k_m</math> are constants.  A typical value for <math>I_{pol}</math> and <math>I_{exp}</math> is six, so an equation of this form would have 42 fitted constants (excluding <math>I_{pol}</math> and <math>I_{exp}</math>).  Fitting a new substance to an equation of state of this form does not necessarily mean fitting several dozen constants however.  If there is an existing reference equation of state for a chemically similar molecule (''e.g.'' a saturated hydrocarbon), then it may only be necessary to fit <math>c_m</math> to obtain an adequate fit to the available experimental data.  The <math>c_m</math> constants may also be obtained by a corresponding states correlation.  Span proposes a correlation for nonpolar substances where <math>c_m=f(m,\omega)</math>, where <math>\omega</math> is the accentric factor<ref name="name" />.  Xiang recommends a correlation that applies to nonpolar and polar substances of the form <math>c_m=f(m,\omega,Z_c)</math>, where <math>Z_c</math> is the critical compressibility<ref name="name2" />.  Use of one of these correlations will result in poorer accuracy than fitting <math>c_m</math>, <math>i_m</math>, <math>j_m</math>, and <math>k_m</math> or just <math>c_m</math>, but the results are superior in accuracy over using a cubic equation of state such as Peng-Robinson<ref name="name" /> or a corresponding states approach such as Lee-Kesler<ref name="name2" />, all of which require a comparable amount of user-input information.
<br>
<br>

==Vapor-Liquid Equilibrium==
In implementing a Helmholtz energy equation of state in a code, the ability to distinguish between single-phase vapor or liquid and a two-phase mixture must be made.  One approach is to calculate the vapor-liquid phase envelope once and save the results in a look-up table.  These values are then used as initial guesses for future subroutine calls in either converging on vapor liquid equilibrium or determining that the system is in a single phase region.  Calculation of the vapor-liquid phase envelope entails the identification of all thermodynamic states that satisfy three conditions for phase equilibria.  The first requirement for phase equilibria in a single component system is chemical potential equilibrium:
<br>
<br>
<math>g_L = g_V</math>
<br>
<br>
Where <math>g_L</math> and <math>g_V</math> are the liquid and vapor Gibbs energy.  The second requirement is mechanical equilibrium:
<br>
<br>
<math>P_L = P_V</math>
<br>
<br>
Where <math>P_L</math> and <math>P_V</math> are the liquid and vapor pressures.  The third requirement is thermal equilibrium:
<br>
<br>
<math>T_L = T_V</math>
<br>
<br>
Where <math>T_L</math> and <math>T_V</math> are the liquid and vapor temperatures.  One approach in calculating the phase envelope is to start by solving for the liquid and vapor densities that satisfy the Gibbs energy and pressure constraints at the triple point temperature.  This solution of the two equations for the two unknowns can be accomplished by Newton-Raphson provided suitable initial guesses are provided.  These initial guesses can come from experimental data or simple corresponding states correlations such as those found in ''Properties of Liquids and Gases''<ref>Reid, Robert C., J. M. Prausnitz, and Thomas K. Sherwood. 1977. ''The Properties of Gases and Liquids.'' New York: McGraw-Hill.</ref>.  Once the saturated liquid and vapor densities at the triple point are converged on, these can be used as initial guesses for the saturation properties at <math>T_{tr}+\Delta T</math>, where <math>\Delta T</math> is some small increment in temperature, perhaps 0.1 Kelvin.  Proceeding in this manner to the critical point, the phase envelope is calculated.
<br>
<br>

==Resources==
Access to property information calculated by reference equations of state for 75 substances of great industrial and scientific importance is provided by NIST online<ref>Thermophysical Properties of Fluid Systems. (2011). Retrieved April 26, 2015, from http://webbook.nist.gov/chemistry/fluid/</ref>.  Citation information for the original journal or book where each equation was published is also provided by NIST.  The book by Span is a must have for anyone interested in learning more about Helmholtz energy equations of state and their implementation in scientific codes<ref name="name" />.
<br>
<br>

==References==
<references/>
[[Category: Multiscale Simulations]]
[[Category: ICME Class]]
[[Category: Macroscale]]