Transport Properties
This section includes definitions of the models available in for thermal conductivity, viscosity, and diffusivity:
Thermal Conductivity
Vapor
Ideal
 
The thermal conductivity correlations is according to:
(6-132)
The pressure correction Δλv,P is calculated from the method of Stiel and Thodos, see Ref. 32, which is applicable for ρr<3, but is less accurate for H2, strongly polar gases, and gases with a high degree of hydrogen bonding, such as H2O and NH3
(6-133)
(6-134)
(6-135)
The mixing rules are as suggested by Yorizane, see Ref. 33:
(6-136)
(6-137)
(6-138)
(6-139)
(6-140)
(6-141)
where the binary constants are
(6-142)
(6-143)
The Stiel and Thodos coefficients are
ρr < 0.5
A1=2.702×108
B1=0.535
C1=-1
0.5 ≤ ρr < 2.0
A2=2.528×108
B2=0.670
C2=-1.069
ρr2.0
A3=0.574×108
B3=1.155
C3= 2.016
However, in order to ensure 0th-order continuity at ρr = 0.5 and ρr = 2.0, the following coefficients for 0.5 ≤ ρr < 2.0 are recalculated from
(6-144)
(6-145)
The vapor thermal conductivity correlation must be available for all species. Also critical volumes, Vc,i, critical temperatures, Tc,i, molecular weights Mi, and acentric factors ωi must be specified for all species.
Kinetic Theory
 
Lindsay and Bromley (see Ref. 34) provided an equation for the interaction parameters of the method of Wassiljewa (see Ref. 35) based on the kinetic theory, to provide mixture thermal conductivity from pure species values
(6-146)
(6-147)
where the pressure correction Δλv,P is calculated from Equation 6-133. Both vapor thermal conductivity correlation λi,v and the vapor viscosity correlation ηi,v must be available for all species. In addition, all normal boiling points Ti,b, molecular weights Mi, critical volumes Vc,i, critical temperatures Tc,i, and acentric factors ωi must be specified.
Water (IAPWS)
 
The International Association of the Properties of Water and Steam recommend an equation Ref. 36 which is valid in the following range:
(6-148)
(6-149)
(6-150)
Liquid
The following mixture models are available for liquid thermal conductivity
Ideal
 
To calculate the mixture liquid thermal conductivity, λl,m, the values of pure liquid thermal conductivity correlations are mixed ideally
(6-151)
The pressure dependence is based on the work of Missenard (Ref. 52) where
(6-152)
where Q is correlated as
(6-153)
and the following mixing rules are used
(6-154)
(6-155)
(6-156)
(6-157)
(6-158)
(6-159)
All liquid thermal conductivity correlations must be specified. All values for critical temperatures, Tc,i, critical volumes, Pc,i and critical compressibility factors, Zc,i must be specified for all species i.
Power Law
 
The values of pure liquid vapor thermal conductivity correlations are mixed according to the following power law
(6-160)
All liquid thermal conductivity correlations must be specified. The model is valid for pure compound thermal conductivity values that are no more apart than a factor of 2 (Ref. 53 and Ref. 54). The pressure dependence is introduced using Equation 6-152 through Equation 6-159. All values for critical temperatures, Tc,i, critical volumes, Pc,i and critical compressibility factors, Zc,i must be specified for all spices i.
Local Composition
 
The local composition model by Rowley (Ref. 54) uses an ideal and excess contribution
(6-161)
The ideal part is based on mass fractions
(6-162)
(6-163)
The excess term is based on NRTL local concentrations
(6-164)
where Gj,i follows from Equation 6-57. The binary interaction terms follow from
(6-165)
which is symmetric, and on the diagonal,
(6-166)
and
(6-167)
with ϖi is the composition in the binary mixture of species i and j and the local composition is equi-molar
(6-168)
All liquid thermal conductivity correlations must be specified. The pressure dependence is introduced using Equation 6-152Equation 6-159. All values for critical temperatures Tc,i, critical volumes Pc,i, critical compressibility factors Zc,i, and molecular weights Mi must be specified for all compounds i. In addition, all NRTL binary interaction parameters Ai,j must be specified. Unspecified values for NRTL interaction parameters Bi,j are set to zero. The randomness parameters αi,j have values of zero on the diagonal and the matrix is symmetric. All off-diagonal values must be specified. NRTL model is presented in Equation 6-53 to Equation 6-64.
Local Composition (Modified)
 
Rowley (Ref. 53) adapted the local composition model by replacing the mixing rule in Equation 6-165 by the following
(6-169)
which he found to produce better model predictions in most cases where both the Local Composition model and Power Law model have trouble. However, the model is not as generally applicable; for instance, systems containing H2O are not well described by this model due to the low molecular weight of H2O.
Water (IAPWS)
 
The International Association of the Properties of Water and Steam recommend an equation Ref. 36 which is valid in the following range:
(6-170)
(6-171)
(6-172)
Viscosity
Vapor
Wilke
 
Wilke, see Ref. 37, based his method for mixture viscosity of the vapor phase on kinetic theory:
(6-173)
(6-174)
The vapor viscosity correlation ηi,v must be available for all species. In addition, all molecular weights Mi must be specified.
Brokaw
 
Brokaw (see Ref. 38) uses the same basic equation as Wilke (Equation 6-173). However, Equation 6-174 is replaced by
(6-175)
and the interaction parameter is defined as
(6-176)
where
(6-177),
The vapor viscosity correlation, ηi,v must be available for all species i. In addition, all molecular weights Mi must be specified. If Lennard-Jones energy εi (see Ref. 39) Stockmayer’s polar parameter δs,i(Ref. 40 and Ref. 41) are specified for both species i and j then
(6-178)
Otherwise,
(6-179)
High Pressure Modification
 
To account for the effect of pressure on vapor viscosity, a pressure correction can be applied. The pressure dependence is based on kinetic gas theory, which adds the following term to the vapor viscosity:
(6-180)
where ξ is calculated from the correlation of Jossi (Ref. 42), which is applicable for ρr < 3.0. It is less accurate for H2, strongly polar gases and gases with a high degree of hydrogen bonding such as H2O and NH3.
The correction factor is due to using pressure, atm, and viscosity, cP, units in Jossi’s correlation. It is expressed as:
(6-181)
(6-182)
(6-183)
The following mixture rules are used
(6-184), ,
(6-185),
The values for critical volumes, Vc,i, critical temperatures, Tc,i, critical compressibility factors, Zc,i and molecular weights, Mi must be specified for all species i.
The high pressure correction is available for the Wilke and Brokaw mixture models. The vapor viscosity follows from
(6-186)
where ηv,Wilke is calculated from Equation 6-173.
Pedersen Corresponding States Model
 
The corresponding states viscosity model of Pedersen (Ref. 43 and Ref. 44) applies to both vapor and liquid phases of hydrocarbon mixtures. The selected reference species is CH4.
The CH4 viscosity is calculated from Ref. 45, modified by Pedersen and Fredunsland (Ref. 46) to avoid issues below 91 K where CH4 becomes solid
(6-187)
where
(6-188)
(6-189)
(6-190)
(6-191)
Here, ρCH4 is used in g/cm3; for the mass-mole conversion of ρCH4, a molecular weight of MCH4 = 16.042568 g/mol is used.
The dilute gas part is given by
(6-192)
The first density correction for the moderately dense gas is given by
(6-193)
The remainder is given by the empirical correlation
(6-194)
The correction term for solid CH4
(6-195)
with the values of the parameters L1 through L25 are listed in Table 6-5 below:
L1
L10
L19
L2
L11
L20
L3
L12
L21
L4
L13
L22
L5
L14
L23
L6
L15
L24
L7
L16
L25
L8
L17
L9
L18
Here, ρCH4 is used in g/cm3; the critical density is given by ρc,CH4 = 0.16284 g/cm3. The following equation by McCarty (Ref. 47) is solved for the density of CH4
(6-196)
where ρCH4 is used in mol/l.
N1
N13
N25
N2
N14
N26
N3
N15
N27
N4
N16
N28
N5
N17
N29
N6
N18
N30
N7
N19
N31
N8
N20
N32
N9
N21
N33
N10
N22
N34
N11
N23
N12
N24
With the viscosity and density of CH4 defined, the viscosity of any mixture, ηm, can be calculated from the corresponding states principle
(6-197)
where CH4 viscosity ρCH4,P0,T0 is calculated at temperature, T0 and pressure P0
(6-198)
(6-199)
The following mixing rules are used for the critical properties, see Ref. 48:
(6-200)
(6-201)
(6-202)
The parameter α is
(6-203)
(6-204)
where
(6-205)
(6-206)
with ρc,CH4 = 0.16284 g/cm3. For CH4
(6-207)
where Equation 6-204 is used. The mixture molecular weight is a function of the weight-averaged molecular weight and the number-averaged molecular weight
(6-208)
(6-209)
(6-210)
where the power in Equation 6-208 is determined by fitting to experimental viscosity data.
Note that pure species vapor viscosity correlations ηi,v are not required. However, for each species i, molecular weight Mi, critical temperature, Tc,i and critical pressure, Pc,i must be specified.
Water (IAPWS)
 
The International Association of the Properties of Water and Steam recommend an equation Ref. 49- Ref. 50 for industrial application which is valid in the following range:
(6-211)
(6-212)
(6-213)
(6-214)
(6-215)
where is the pressure dependent melting temperature and is the triple-point pressure. In accordance with industrial application recommendations, critical region correction is not applied.
Liquid
The following mixture models are available for liquid viscosity.
Molar Logarithmic Mixing
 
The values of pure species log liquid viscosity, ln ηi,l are mixed ideally using mole fractions xi
(6-216)
where ηm,l is the mixture viscosity of liquids.
Mass Logarithmic Mixing
 
The values of pure species log liquid viscosity correlation are mixed ideally using the weight fractions ωi
(6-217)
Pedersen Corresponding States Model
 
The Pedersen Corresponding States Model described above for the gas phase viscosity also applies to the liquid phase. Pure species liquid viscosity correlations are not required. However, for each species i, molecular weight Mi, critical temperature Tc,i, and critical pressure Pc,i, must be specified.
Cubic mixing
 
The mixture viscosity is defined using the cubic root average in terms of the mole fractions xi
(6-218)
The model is noted in Ref. 51 to provide reasonable results for hydrocarbon mixtures of similar components.
The model requires that the log liquid viscosity correlation is available for all species i.
Cubic mass mixing
 
The mass fraction equivalent of the previous model is
(6-219)
The model requires that the log liquid viscosity correlation is available for all species i.
Water (IAPWS)
 
The International Association of the Properties of Water and Steam recommend an equation Ref. 49- Ref. 50 for industrial application which is valid in the following range:
(6-220)
(6-221)
(6-222)
(6-223)
(6-224)
where is the pressure dependent melting temperature and is the triple-point pressure. In accordance with industrial application recommendations, critical region correction is not applied.
Diffusivity
Two types of diffusion coefficients are supported. Diffusion coefficients in infinitely diluted systems, and Maxwell-Stefan diffusion coefficients.
For dilute systems, the binary diffusion coefficient D0i,j represent the diffusivity of species i in a medium consisting of pure species j. This corresponds to the Fickian diffusion coefficient.
For any mixture, the binary Maxwell-Stefan diffusion coefficient , represents the inverse drag coefficient of species i moving past species j (Ref. 55 - Ref. 58). This property is referred to as the Maxwell-Stefan diffusivity. The Maxwell-Stefan diffusivity is symmetric, , and the diagonal elements are not used. T
Gas Phase Diffusion Coefficient at Infinite Dilution
The following models are available for the diffusion coefficients at infinite dilution in the vapor phase:
Automatic
When the Gas diffusivity property model is set to Automatic, the Fuller Schettler Giddings model is used, provided that the Fuller diffusion volume is known for both species (i and j), otherwise the Wilke-Lee model is used.
Fuller Schettler Giddings
 
Fuller et al. (Ref. 59) modified the Chapman-Enskog relation to correlate binary diffusion coefficient for species i and j in the vapor phase according to the Fuller Schettler Giddings (FGS) model:
(6-225)
where T denotes the temperature (K), Mi the molecular weight of species i (g/mol) and P is the pressure (Pa). vi are the atomic diffusion volumes (Fuller diffusion volume), (cm3), which are estimated using group contribution for each species (Ref. 60):
For some simple molecules the below values, determined from regression, are used:
H2
D2
N2
O2
N2O
H2O
Wilke-Lee
 
Wilke and Lee (Ref. 61) also modified the Chapman-Enskog relation to correlate binary diffusion coefficient in vapor phase according to:
(6-226)
where
(6-227)
The length scale for the interaction is taken from the Lennard-Jones diameter parameters of species i and j:
(6-228)
If is not specified in the database, it is instead estimated from:
(6-229)
where Vi,l,b is the molar volume of species at normal boiling point.
The collision integral is evaluated from (Ref. 62):
(6-230)
with
(6-231)
The energy scale for the interaction is taken from the Lennard-Jones energy parameters of species i and j:
(6-232)
If is not specified in the database, it is instead estimated from:
(6-233)
where k is the Boltzmann constant, and is normal boiling point temperature.
Gas Phase Maxwell-Stefan Diffusivity
For gas phase diffusion the Maxwell-Stefan diffusivities are defined from the models for gas phase diffusivity at infinite dilution
(6-234)
All models for gas phase diffusion at infinite dilution are symmetric, the diffusivity of species i in species j equals that of species j in species i. In addition, under the ideal gas assumption, the Maxwell-Stefan diffusion coefficient matches the Fick diffusion coefficient. This implies that the Maxwell-Stefan gas diffusivities provided are independent of composition.
Liquid Phase Diffusion Coefficients at Infinite Dilution
The following models are available for the diffusion coefficients at infinite dilution in liquid phase:
Automatic
When the Liquid diffusivity at infinite dilution property model is set to Automatic, a selection for each solute i in solvent j will be made from the models according to the following rules. Except in the case of a temperature correlation and the Erkey-Rodden-Akgerman model, all of these rules require that the log liquid viscosity correlation, ln ηj is available in the database for species j and that the liquid volume at normal boiling point, Vi,l,b is available for the species i.
If the solvent is water, the Siddiqi-Lucas correlation (for aqueous systems) is used.
If parachors Pi and Pj are both available, and the liquid volume at normal boiling point, Vi,l,b is available for the solvent, the Tyn-Calus correlation is used.
If parachors Pi and Pj are both available, the Hayduk-Minhas correlation is used.
If the molecular weight Mj is available for the solvent, the Wilke-Chang correlation is used.
For the Siddiqi-Lucas method, a molecule is considered organic if it has at least one C atom bound to anything else than O or C atoms. This is determined by the SMILES formula, if available. If, for any pair of species, the required input data for none of the above models is available, the entire property liquid diffusion coefficient at infinite dilution is not available.
Wilke-Chang
 
The correlation by Wilke and Chang (Ref. 63) for liquid phase diffusion coefficients at infinite dilution is:
(6-235)
Molecular weight, Mj and log liquid viscosity correlation, , for species j and liquid molar volume at normal boiling point for species i, , is required. The Wilke-Chang association parameter , if unavailable, is set to:
Species are identified by their CAS number or SMILES formula. The Wilke-Chang correlation is not suitable for diffusion of water. If water is the solute, the correction suggested by Kooijman (Ref. 65) is applied where liquid molar volume of water at normal boiling point, Vwater,l,b is multiplied by 4.5.
Tyn-Calus
 
The Correlation by Tyn and Calus (Ref. 66) for liquid diffusion coefficients at infinite dilution reads:
(6-236)
The log liquid viscosity correlation should be available for species j, The liquid volume at normal boiling point, Vi,l,b and parachor, Pi should be available for both species i and j.
If the solvent is nonpolar (dipole moment is zero), and the solvent is methanol, ethanol or1-butanol or if the solvent is a mono-hydroxy alcohol, both the liquid volume at normal boiling point, Vj,l,b and the parachor, Pj are corrected by a factor of .
If the solute is water, both the liquid volume at normal boiling point, Vi,l,b and parachor, Pi for the solute are corrected by a factor 2. This factor also applies if it is detected from the SMILE formula that the species is an organic acid (a carboxyl group is found), except in the cases where the solvent is water, methanol or n-butanol.
Hayduk-Minhas
 
Hayduk and Minhas (Ref. 67) suggested three different correlations for liquid diffusion coefficients at infinite dilution.
Aqueous Solutions: in case the solvent is water (derived from CAS number or SMILES formula), the correlation reads:
(6-237)
Normal paraffin solutions: In case both the solute and solvent are normal paraffins (derived from SMILES formula), the correlation reads:
(6-238)
For all other systems, the correlation is:
(6-239)
If the solvent is nonpolar (dipole moment is zero), and the solvent is methanol, ethanol or1-butanol or if the solvent is a mono-hydroxy alcohol, both the liquid volume at normal boiling point, Vj,l,b and the parachor, Pj are corrected by a factor of .
If the solute is water, both the liquid volume at normal boiling point, Vi,l,b and parachor, Pi for solute are corrected by a factor 2. This factor also applies if it is detected from the SMILES formula that the species is an organic acid (a carboxyl group is found), except in the cases where the solvent is water, methanol or n-butanol.
Siddiqi-Lucas
 
Siddiqi and Lucas (Ref. 68) suggested correlations for liquid diffusion coefficients at infinite dilution for aqueous system (including gases) and for organic solutions. For normal paraffins systems, they recommended the Hayduk-Minhas model.
For aqueous solutions:
(6-240)
For all other systems, the equation for organic solution is used as:
(6-241)
Where T is temperature, Vi,l,b is liquid volume at normal boiling point and is the viscosity of the solvent.
Erkey-Rodden-Akgerman
 
The correlation by Erkey et al. (Ref. 69) for liquid diffusion coefficients at infinite dilution in normal paraffins is:
(6-242)
where the reference volume is
(6-243)
NAv is Avogadro’s number, and the deviation from closest packing volume is given by
(6-244)
The correlation is fitted to normal paraffins, hydrogen, carbon monoxide and carbon dioxide diffusing in normal paraffins. The molecular weight Mi should be available for both solute and solvent. The Lennard-Jones diameter, σi is estimated from Bondi group contribution method (Ref. 70), and for some species are taken from Ref. 69 and Ref. 71 as:
H2
n-CjH2j+2
To prevent the diffusion coefficients from becoming negative, the minimum difference of (V − Vjref) is considered to be 1012 mol/m3.
Liquid Phase Maxwell-Stefan Diffusivity
The Maxwell-Stefan liquid diffusion coefficients are calculated from the liquid diffusion coefficients at infinite dilution. The diagonal values should be ignored and are set to zero. The Vignes (Ref. 72) model for diffusion in binary solutions can be extended to multi-component systems (Ref. 73):
(6-245)
where denotes the Maxwell-Stefan liquid diffusion coefficient for species i and j in the limited of pure species k. If k=j, it represents a binary system of species i and j where i is infinitely diluted. In addition, at infinite dilution the thermodynamics factor (activity) (Ref. 58) becomes unity and the Maxwell-Stefan diffusivity equals to Fick diffusivity:
(6-246)
Similarly,
(6-247)
Models for the for and remain to be defined in such a way that symmetry is ensured:
(6-248)
and to ensure continuity if both species i and j vanish:
(6-249)
For binary systems, Equation 6-245 reduces to the Vignes interpolation formula.
Wesseling-Krishna
 
Wesseling and Krishna (Ref. 73) proposed:
(6-250)
leading to:
(6-251)
Kooijman-Taylor
 
Kooijman and Taylor (Ref. 74) found on a limited number of systems that this provides better results:
(6-252)
leading to:
(6-253)
Wesselingh and Bollen (Ref. 75) asserted that this is a reasonable estimate.
Krishna-van Baten
 
Krishna and van Baten (Ref. 76), on the basis of data obtained from molecular dynamic simulations, proposed the following Vignes-based (Ref. 72) interpolation:
(6-254)
The model reduces to the Kooijman-Taylor model for xi = xj. The value of is undefined in the limit of both and . From a physical point of view this is inconsequential as the value cancels out in the expressions obtained for the fluxes using the Maxwell-Stefan equations. Nevertheless, in order to obtain well-defined values and composition derivatives of itself, the equation is modified to:
(6-255)
where is taken as a small composition. The limiting case for both and also reduces to the Kooijman-Taylor model. This equation leads to:
(6-256)