There are many equations that describe the thermodynamic state and properties of steam. While some of these equations are accurate in generating property tables, they are not suitable for fast CFD computations. Therefore, FLUENT uses a simpler form of the thermodynamic state equations [ 410] for efficient CFD calculations that are accurate over a wide range of temperatures and pressures. These equations are described below.
Equation of State
The steam equation of state used in the solver, which relates the pressure to the vapor density and the temperature, is given by [ 410]:
where , and are the second and the third virial coefficients given by the following empirical functions:
where is given in /kg, = with given in Kelvin, = 10000.0, = 0.0015, = -0.000942, and = -0.0004882.
where is given in / , = with given in Kelvin, = 0.8978, =11.16, = 1.772, and = 1.5E-06.
The two empirical functions that define the virial coefficients and cover the temperature range from 273 K to 1073 K.
The vapor isobaric specific heat capacity is given by:
The vapor specific enthalpy, is given by:
The vapor specific entropy, is given by:
The isobaric specific heat at zero pressure is defined by the following empirical equation:
where is in KJ/kg K, = 46.0, = 1.47276, = 8.38930E-04, = -2.19989E-07, = 2.46619E-10, and = -9.70466E-14.
= , = , = , and = .
Both and are functions of temperature and they are defined by:
where and are arbitrary constants.
The vapor dynamic viscosity and thermal conductivity are also functions of temperature and were obtained from [ 409].
Saturated Vapor Line
The saturation pressure equation as a function of temperature was obtained from [ 302]. The example provided in Section 23.13.5 contains a function called wetst_satP() that represents the formulation for the saturation pressure.
Saturated Liquid Line
At the saturated liquid-line, the liquid density, surface tension, specific heat , dynamic viscosity, and thermal conductivity must be defined. The equation for liquid density, , was obtained from [ 302]. The liquid surface tension equation was obtained from [ 409]. While the values of , and were curve fit using published data from [ 93] and then written in polynomial forms. The example provided in Section 23.13.5 contains functions called wetst_cpl(), wetst_mul(), and wetst_ktl() that represent formulations for , and .
The mixture properties are related to vapor and liquid properties via the wetness factor using the following mixing law:
where represents any of the following thermodynamic properties: , , , , or .