Chapman Enskog Viscosity Model#
The ChapmanEnskogLennardJones
pure component gas dynamic viscosity model can be imported from idaes.models.properties.modular_properties.pure
,
while additional rules and utility functions can be imported from idaes.models.property_models.modular_properties.pure.ChapmanEnskog
.
Formulation#
The description of Chapman Enskog theory here is derived from that in The Properties of Gases and Liquids, Section 9-3. Pure component gas viscosity is given by the formula
in which \(\mu\) is the dynamic viscosity, \(M\) is the molar mass, \(\sigma\) is the collision diameter of the molecule (equal to the molecule’s diameter for hard spheres and typically taken to be the Lennard-Jones parameter \(\sigma\) for real molecules), and \(\Omega_v(T)\) is a dimensionless quantity known as the collision integral (equal to 1 for hard spheres but a function of temperature for real molecules). \(C\) is a constant with value \(26.69\;\mu\text{PÅ}^2/\sqrt{(\text{g}/\text{mol})\text{K}}\).
The collision integral is specific for the property of viscosity (e.g., there is a different collision integral for diffusivity calculations) and is given in terms of the dimensionless temperature
in which \(k_B\) is the Boltzmann constant and \(\varepsilon\) is the Lennard-Jones well depth. For good viscosity predictions, it is vital that \(\sigma\) and \(\varepsilon\) come from the same source, preferably from fitting experimental viscosity data (see Reichenberg (1973) for more details). Several different forms of the collision integral, of varying complexity, are available to use. If no callback is specified, the form proposed by Neufeld is used.
Defined Collision Integral Callbacks#
These callbacks provide expressions for the viscosity collision integral.
- idaes.models.properties.modular_properties.pure.ChapmanEnskog.collision_integral_neufeld_callback(T_dim)[source]#
Equation 9.4.3 from Properties of Gases and Liquids.
\[\Omega(T^*) = A(T^*)^{-B} + C\exp(-DT^*) + E\exp(-FT^*)\]in which \(A=1.16145\), \(B=0.14874\), \(C=0.52487\), \(D=0.77320\), \(E=2.16178\), and \(F=2.43787\).
List of Parameters#
Parameter Name |
Description |
Units |
---|---|---|
|
Molecular weight \(M\) |
Mass/Amount |
|
Lennard-Jones ‘particle size’ \(\sigma\) |
Length |
|
Reduced Lennard-Jones well depth \(\varepsilon /k_B\) |
Temperature |
|
Callback to use for viscosity integral |
n/a |
Example#
The code snippet below demonstrates how to specify use of the ChapmanEnskogLennardJones
model for pure component vapor viscosity as part of the modular property
framework. Note that if you specify visc_d_phase_comp
for one phase, you must specify it for all phases, even if only to pass None
as the
method.
from idaes.models.properties.modular_properties.pure import ChapmanEnskogLennardJones
from idaes.models.properties.modular_properties.pure.ChapmanEnskog import collision_integral_neufeld_callback, collision_integral_kim_ross_callback
configuration = {
"components":{
"H2O": {
"type": Component,
"valid_phase_types": [PhaseType.vaporPhase, PhaseType.liquidPhase],
...
"visc_d_phase_comp": {"Vap": ChapmanEnskogLennardJones, "Liq": None},
"viscosity_collision_integral_callback": collision_integral_neufeld_callback,
"parameter_data": {
"mw": (0.01801528, pyunits.kg / pyunits.mol),
"lennard_jones_sigma": (2.641, pyunits.angstrom),
"lennard_jones_epsilon_reduced": (809.1, pyunits.K),
}
...
}
...
}
...
}
References#
Poling, Bruce, E. et al. The Properties of Gases and Liquids. 5th ed. New York: NcGraw-Hill, 2001.
Reichenberg, Daniel. “The indeterminacy of the values of potential parameters as derived from transport and virial coefficients.” AIChE Journal 19.4 (1973): 854-856.