Calculation of thermal conductivity in FPROPS is in development. This development is driven by an application requiring transport properties of carbon dioxide, so the first correlations being implemented will be for that. Although textbooks such as Incropera and DeWitt, Holman and Cengel use for thermal conductivity, most publications on thermophysical properties use the convention of , and that convention will be used on this page.
Correlations typically specify conductivity in terms of ideal (zero-density limit) part , residual part and critical region enhancement function , as follows:
Here, is the isobaric heat capacity divided by the gas constant (either both specific or both molar), is length scaling parameter (in units of nm, which is messy), a constant particular to the fluid in question, and is relative molecular mass of the fluid.
The reduced temperature is calculated as
where the constants would be tabulated for each fluid separately.
Essentially, the ideal part of the thermal conductivity, then, is a ratio of two power series in terms of : one derived from experimentally-determined , and the other a reduced effective cross-section derived experimentally from low-density conductivity data, . To accommodate the many different ways that authors may publish their results, we can allow arbitrary power series on both numerator and denominator, and we can allow the already-defined function to be referenced/used as the denominator if desired.
The calculation of the effective cross section can be generalised a little better by defined it as
with , and is a scaling temperature provided as part of the correlation data. The specific correlation above would be writen with hence .
(Note that the reduced effective cross section is dimensionless, and is reduced from the 'real' effective cross-section area using .)
Note that one publication provides the function in the form
and this form does not comply with the above ratio of power series, since zero-density viscosity is, in turn, defined as a ratio of power series (see FPROPS/Viscosity). This form could be accommodated by 'flagging' one of the terms in the numerator power series in a special way, and providing a pointer to the appropriate data/function, perhaps.
The residual part is correlated by Lemmon and Jacobsen as:
with , , , all given as correlation data for the particular fluid, and defined as zero when is zero, or one otherwise.
Lemmon and Jacobsen define and ; but for flexibility, we can define arbitrarily-scaled temperature and density, with and provided as part of the correlation data,
Vesovic et al find that they can represent the residual part for carbon dioxide using power series in terms only of , as follows:
This form can be accommodated by the form of Lemmon and Jacobsen by setting , , and and to zero, and and to one.
The approach attributed to Olchowy and Sengers is used in both the carbon dioxide data and nitrogen data, although the equations presented in each of these references appear somewhat different (Specifically, in Lemmon and Jacobsen, has replaced from Vesovic et al, and the function has been divided by ). Using the formulae of Lemmon and Jacobsen, the following parameters are required as correlation data:
|symbol||value and units (for CO2)||name|
There are also several 'universal' theoretical constants in the Olchowy and Sengers equations:
- = 1.01
- = 0.630
- = 1.2415
We define the 'reduced symmetrized compressibility' as
then the critical enhancement to thermal conductivity is calculate from
- V Vesovic, W A Wakeham, G A Olchowy, J V Sengers, J T R Watson and J Millat, 1990. The Transport Properties of Carbon Dioxide, J Phys Chem Ref Data 19, 763. doi:10.1063/1.555875.
- E W Lemmon and R T Jacobsen, 2004. Viscosity and Thermal Conductivity Equations for Nitrogen, Oxygen, Argon, and Air Int J Thermophys 25. doi:10.1023/B:IJOT.0000022327.04529.f3.
- V Vesovic, 1994. "On Correlating the Transport Properties of Supercritical Fluids", in Supercritical Fluids: Fundamentals for Application, Springer, pp 273-283. doi:10.1007/978-94-015-8295-7_10. (Note the error in equation (6) in this reference.)