AspenPhysPropModelsV7 0-Ref
AspenPhysPropModelsV7 0-Ref
Property System
Physical Property
Models
Version Number: V7.0
July 2008
Copyright (c) 2008 by Aspen Technology, Inc. All rights reserved.
Aspen Physical Property System, the aspen leaf logo and Plantelligence and Enterprise Optimization are trademarks
or registered trademarks of Aspen Technology, Inc., Cambridge, MA.
All other brand and product names are trademarks or registered trademarks of their respective companies.
This document is intended as a guide to using AspenTech's software. This documentation contains AspenTech
proprietary and confidential information and may not be disclosed, used, or copied without the prior consent of
AspenTech or as set forth in the applicable license agreement. Users are solely responsible for the proper use of
the software and the application of the results obtained.
Although AspenTech has tested the software and reviewed the documentation, the sole warranty for the software
may be found in the applicable license agreement between AspenTech and the user. ASPENTECH MAKES NO
WARRANTY OR REPRESENTATION, EITHER EXPRESSED OR IMPLIED, WITH RESPECT TO THIS DOCUMENTATION,
ITS QUALITY, PERFORMANCE, MERCHANTABILITY, OR FITNESS FOR A PARTICULAR PURPOSE.
Contents 1
COSMO-SAC .......................................................................................76
Electrolyte NRTL Activity Coefficient Model ..............................................79
ENRTL-SAC ........................................................................................92
Hansen ..............................................................................................97
Ideal Liquid ........................................................................................98
NRTL (Non-Random Two-Liquid)............................................................98
Using NRTL-SAC................................................................................ 104
Pitzer Activity Coefficient Model ........................................................... 105
Polynomial Activity Coefficient............................................................. 117
Redlich-Kister ................................................................................... 118
Scatchard-Hildebrand ........................................................................ 119
Three-Suffix Margules ........................................................................ 120
UNIFAC Activity Coefficient Model ........................................................ 121
UNIFAC (Dortmund Modified) .............................................................. 123
UNIFAC (Lyngby Modified) .................................................................. 124
UNIQUAC Activity Coefficient Model ..................................................... 125
Van Laar Activity Coefficient Model ...................................................... 126
Wagner Interaction Parameter ............................................................ 127
Wilson Activity Coefficient Model.......................................................... 128
Wilson Model with Liquid Molar Volume................................................. 129
Vapor Pressure and Liquid Fugacity Models ..................................................... 131
Extended Antoine/Wagner/PPDS/IK-CAPE Liquid Vapor Pressure Model ..... 131
API Sour Model ................................................................................. 136
Braun K-10 Model ............................................................................. 136
Chao-Seader Pure Component Liquid Fugacity Model .............................. 137
Grayson-Streed Pure Component Liquid Fugacity Model .......................... 137
Kent-Eisenberg Liquid Fugacity Model................................................... 138
Maxwell-Bonnell Vapor Pressure Model ................................................. 139
Solid Antoine Vapor Pressure Model ..................................................... 139
DIPPR/Watson/PPDS/IK-CAPE Heat of Vaporization Model ................................. 139
DIPPR Heat of Vaporization Equation.................................................... 140
Watson Heat of Vaporization Equation .................................................. 140
PPDS Heat of Vaporization Equation ..................................................... 141
IK-CAPE Heat of Vaporization Equation................................................. 141
NIST TDE Watson Heat of Vaporization Equation .................................... 142
Clausius-Clapeyron Equation ............................................................... 142
Molar Volume and Density Models ................................................................. 143
API Liquid Volume ............................................................................. 143
Brelvi-O'Connell ................................................................................ 144
Clarke Aqueous Electrolyte Volume ...................................................... 145
COSTALD Liquid Volume..................................................................... 147
Debye-Hückel Volume ........................................................................ 149
Liquid Constant Molar Volume Model .................................................... 150
Rackett/DIPPR/PPDS/IK-CAPE Liquid Molar Volume ................................ 150
Rackett/Campbell-Thodos Mixture Liquid Volume ................................... 155
Modified Rackett Liquid Molar Volume................................................... 156
Aspen/DIPPR/IK-CAPE Solid Molar Volume ............................................ 157
Liquid Volume Quadratic Mixing Rule.................................................... 159
Heat Capacity Models ..................................................................................160
Aqueous Infinite Dilution Heat Capacity ................................................ 160
Criss-Cobble Aqueous Infinite Dilution Ionic Heat Capacity ...................... 161
DIPPR/PPDS/IK-CAPE Liquid Heat Capacity ........................................... 161
Aspen/DIPPR/Barin/PPDS/IK-CAPE Ideal Gas Heat Capacity .................... 165
2 Contents
Aspen/DIPPR/Barin/IK-CAPE Solid Heat Capacity ................................... 168
Solubility Correlations .................................................................................. 170
Henry's Constant...............................................................................171
Water Solubility ................................................................................ 172
Hydrocarbon Solubility ....................................................................... 172
Other Thermodynamic Property Models .......................................................... 173
Cavett ............................................................................................. 173
Barin Equations for Gibbs Energy, Enthalpy, Entropy, and Heat Capacity ... 174
Electrolyte NRTL Enthalpy................................................................... 176
Electrolyte NRTL Gibbs Free Energy ..................................................... 178
Liquid Enthalpy from Liquid Heat Capacity Correlation ............................ 179
Enthalpies Based on Different Reference States ..................................... 180
Helgeson Equations of State ............................................................... 184
Quadratic Mixing Rule ........................................................................ 187
Contents 3
Nernst-Hartley .................................................................................. 229
Wilke-Chang (Binary)......................................................................... 230
Wilke-Chang (Mixture) ....................................................................... 231
Surface Tension Models................................................................................ 232
Liquid Mixture Surface Tension ............................................................ 232
API Surface Tension........................................................................... 233
IAPS Surface Tension for Water........................................................... 233
Hakim-Steinberg-Stiel/DIPPR/PPDS/IK-CAPE Liquid Surface Tension......... 233
Onsager-Samaras ............................................................................. 237
Modified MacLeod-Sugden .................................................................. 238
Index ..................................................................................................................264
4 Contents
1 Introduction
This manual describes the property models available in the Aspen Physical
Property System and defines the parameters used in each model. The
description for each model lists the parameter names used to enter values on
the Properties Parameters forms.
This manual also lists the pure component temperature-dependent properties
that the Aspen Physical Property System can calculate from a model that
supports several equations or submodels. See Pure Component Temperature-
Dependent Properties (below).
Many parameters have default values indicated in the Default column. A dash
(–) indicates that the parameter has no default value and you must provide a
value. If a parameter is missing, calculations stop. The lower limit and upper
limit for each parameter, when available, indicate the reasonable bounds for
the parameter. The limits are used to detect grossly erroneous parameter
values.
The property models are divided into the following categories:
• Thermodynamic property models
• Transport property models
• Nonconventional solid property models
The property types for each category are discussed in separate sections. The
following table (below) provides an organizational overview of this manual.
The tables labeled Thermodynamic Property Models, Transport Property
Models, and Nonconventional Solid Property Models present detailed lists of
models. These tables also list the Aspen Physical Property System model
names, and their possible use in different phase types, for pure components
and mixtures.
Electrolyte and conventional solid property models are presented in
Thermodynamic Property Models.
1 Introduction 5
Categories of Models
Category Sections
Thermodynamic Equation-of-State Models
Property Models Activity Coefficient Models (Including Electrolyte Models)
Vapor Pressure and Liquid Fugacity Models
Heat of Vaporization Models
Molar Volume and Density Models
Heat Capacity Models
Solubility Correlations
Other
Transport Property Viscosity Models
Models Thermal Conductivity Models
Diffusivity Models
Surface Tension Models
Nonconventional Solid General Enthalpy and Density Models
Property Models Enthalpy and Density Models for Coal and Char
6 1 Introduction
Pure Component Temperature-
Dependent Properties
The following table lists the pure component temperature-dependent
properties that the Aspen Physical Property System can calculate from a
model that supports several equations or submodels.
For example, the Aspen Physical Property System can calculate heat of
vaporization using these equations:
• Watson
• DIPPR
• PPDS
• IK-CAPE
1 Introduction 7
The selection of submodels is driven by the data hierarchy, and controlled by
the submodel-selection parameters.
The thermodynamic properties use the THRSWT (thermo switch) submodel-
selection parameter, and the transport properties use the TRNSWT (transport
switch) submodel-selection parameter.
As the previous table shows, a property is associated with an element of the
submodel-selection parameter. For example, THRSWT element 1 controls the
submodel for solid volume.
The following table shows the values for THRSWT or TRNSWT, and the
corresponding submodels.
Parameter Values
(Equation Number) Submodel
0 Aspen
1 to 116 DIPPR
200 BARIN
301 to 302 PPDS or property-specific methods
400 PML
401 to 404 IK-CAPE
501 to 515 NIST
8 1 Introduction
Set the fourth element of the THRSWT parameter to 106, and the 1-3 and 5-8
elements to 0. If you want to set the other temperature-dependent properties
to use what is defined for that component in the databank, leave the element
blank.
The following table lists the available DIPPR equations and the corresponding
equation (submodel) number.
101
102
103
104
105
106
107
114
116
Extrapolation Methods
Many temperature dependent property models have upper and lower
temperature limits. The Aspen Physical Property System usually extrapolates
linearly beyond such limits. It calculates the slope of the property-versus-
temperature curve, or the ln(property)-versus-temperature curve for models
expressed in logarithmic form, at the upper or lower temperature limit. For
temperatures beyond the limit, it uses a linear model with this slope which
meets the curve from the equation at the temperature limit. Thus the model
is:
1 Introduction 9
For T beyond the upper or lower limit, where Tlim is that limit and Z is the
property or ln(property) as appropriate. Some liquid molar volume models are
actually molar density models which then return the reciprocal of the density
as the liquid molar volume. In these models, the extrapolation occurs for the
density calculation.
There are certain exceptions, detailed below.
For these models, the extrapolation maintains the slope of ln(property) versus
1/T. This applies to the Extended Antoine vapor pressure equation and the
Andrade and DIPPR liquid viscosity equations. Note that the equation for
Henry's Constant is extrapolated by ln(Henry) versus T.
Exception 3: No Extrapolation
The equations for certain models are used directly at all temperatures, so that
no extrapolation is performed. These models are the Wagner vapor pressure
equation, the Aly and Lee equation for the DIPPR Ideal Gas Heat Capacity
(using the CPIGDP parameter), and the Water Solubility and Hydrocarbon
Solubility models. The equations for temperature-dependent binary
interaction parameters are also used directly at all temperatures with no
extrapolation.
10 1 Introduction
2 Thermodynamic Property
Models
Equation-of-State Models
Property Model Model Name(s) Phase(s) Pure Mixture
ASME Steam Tables ESH2O0,ESH2O VL X —
BWR-Lee-Starling ESBWR0, ESCSTBWR V L X X
Benedict-Webb-Rubin-Starling ESBWRS, ESBWRS0 VL X X
Hayden-O'Connell ESHOC0,ESHOC V X X
HF equation-of-state ESHF0, ESHF V X X
Ideal Gas ESIG0, ESIG V X X
Lee-Kesler ESLK VL — X
Lee-Kesler-Plöcker ESLKP0,ESLKP VL X X
NBS/NRC Steam Tables ESSTEAM0,ESSTEAM VL X —
Nothnagel ESNTH0,ESNTH V X X
Peng-Robinson ESPR0, ESPR VL X X
Standard Peng-Robinson ESPRSTD0,ESPRSTD VL X X
Peng-Robinson-Wong-Sandler ESPRWS0,ESPRWS VL X X
Peng-Robinson-MHV2 ESPRV20,ESPRV2 VL X X
Predictive SRK ESRKSV10, ESRKSV1 V L X X
Other Models
Property Model Model Name Phase(s) Pure Mixture
Cavett Liquid Enthalpy DHL0CVT, DHL2CVT L X X
Departure
BARIN Equations for Gibbs — SLV X —
Energy, Enthalpy, Entropy and
Heat Capacity
Electrolyte NRTL Enthalpy HAQELC, HMXELC L — X
Electrolyte NRTL Gibbs Energy GAQELC, GMXELC L — X
Liquid Enthalpy from Liquid Heat DHL0DIP L X X
Capacity Correlation
Enthalpies Based on Different DHL0HREF LV X X
Reference States
Helgeson Equations of State — L — X
Quadratic Mixing Rule — LV — X
Equation-of-State Models
The Aspen Physical Property System has the following built-in equation-of-
state property models. This section describes the equation-of-state property
models available.
Model Type
ASME Steam Tables Fundamental
BWR-Lee-Starling Virial
Benedict-Webb-Rubin-Starling Virial
Hayden-O'Connell Virial and association
HF Equation-of-State Ideal and association
Huron-Vidal mixing rules Mixing rules
Ideal Gas Ideal
Lee-Kesler Virial
References
ASME Steam Tables, Thermodynamic and Transport Properties of Steam,
(1967).
K. V. Moore, Aerojet Nuclear Company, prepared for the U.S. Atomic Energy
Commision, ASTEM - A Collection of FORTRAN Subroutines to Evaluate the
1967 ASME equations of state for water/steam and derivatives of these
equations.
Where:
Mixing Rules:
Where:
Binary interaction parameters BWRKV and BWRKT are available in the Aspen
Physical Property System for a large number of components. (See Physical
Property Data, Chapter 1).
References
M.R. Brulé, C.T. Lin, L.L. Lee, and K.E. Starling, AIChE J., Vol. 28, (1982) p.
616.
Brulé et al., Chem. Eng., (Nov., 1979) p. 155.
Benedict-Webb-Rubin-Starling
The Benedict-Webb-Rubin-Starling equation-of-state is the basis of the BWRS
property method. It is a modification by Han and Starling of the virial
equation-of-state for pure fluids by Benedict, Webb and Rubin. This equation-
of-state can be used for hydrocarbon systems that include the common light
gases, such as H2S, CO2 and N2.
The form of the equation-of-state is:
Where:
In the mixing rules given above, A0i, B0i, C0i, D0i, E0i, ai, bi, ci, di, αi, γi are pure
component constants which can be input by the user. If the values of these
parameters are not given, the Aspen Physical Property System will calculate
them using the critical temperature, the critical volume (or critical density),
the acentric factor and generalized correlations given by Han and Starling.
When water is present, by default Benedict-Webb-Rubin-Starling uses the
steam table to calculate the enthalpy, entropy, Gibbs energy, and molar
volume of water. The total properties are mole-fraction averages of these
values with the properties calculated by the equation of state for other
components. Fugacity coefficient is not affected. An option code can disable
this use of the steam table.
For best results, the binary parameter kij must be regressed using phase-
equilibrium data such as VLE data.
Parameter SymbolDefault MDS Lower Upper Units
Name/ Limit Limit
Element
BWRSTC Tci TC x 5.0 2000.0 TEMPERATURE
BWRSVC Vci VC x 0.001 3.5 MOLE-VOLUME
BWRSOM ωi OMEGA x –0.5 2.0 –
BWRSA/1 B0i fcn(ωi ,Vci , Tci) x – – MOLE-VOLUME
BWRSA/2 A0i fcn(ωi ,Vci , Tci) x – – PRESSURE * MOLE-
VOL^2
BWRSA/3 C0i fcn(ωi ,Vci , Tci) x – – PRESSURE *
TEMPERATURE^2 *
MOLE-VOLUME^2
BWRSA/4 γi fcn(ωi ,Vci , Tci) x – – MOLE-VOLUME^2
References
M. Benedict, G. B. Webb, and L. C. Rubin, J. Chem. Phys., Vol. 8, (1940), p.
334.
M. S. Han, and K. E. Starling , "Thermo Data Refined for LPG. Part 14:
Mixtures", Hydrocarbon Processing, Vol. 51, No. 5, (1972), p.129.
K. E. Starling, "Fluid Themodynamic Properties for Light Petroleum Systems",
Gulf Publishing Co., Houston, Texas (1973).
Where:
, with
, where
, with
, where
, and
Cross-Interactions
The previous equations are valid for dimerization and cross-dimerization if
these mixing rules are applied:
Chemical Theory
with
In this case, molar volume is equal to V/nt.
This represents true total volume over the true number of species nt.
However, the reported molar volume is V/na.
This represents the true total volume over the apparent number of species na.
If dimerization does not occur, na is defined as the number of species. V/na
reflects the apparently lower molar volume of an associating gas mixture.
The chemical equilibrium constant for the dimerization reaction on pressure
basis Kp, is related to the true mole fractions and fugacity coefficients:
Where:
yi and yj = True mole fractions of monomers
yij = True mole fraction of dimer
Apparent mole fractions yia are reported, but in the calculation real mole
fractions yi, yj, and yij are used.
The heat of reaction due to each dimerization is calculated according to:
The sum of the contributions of all dimerization reactions, corrected for the
ratio of apparent and true number of moles is added to the molar enthalpy
departure .
Parameter Name/ Symbol Default MDS Lower Upper Units
Element Limit Limit
TC Tci — — 5.0 2000.0 TEMPERATURE
PC pci — — 105 108 PRESSURE
RGYR rigyr — — 10 -11
5x10 -9
LENGTH
-24
MUP pi — — 0.0 5x10 DIPOLEMOMENT
HOCETA η 0.0 X — — —
The binary parameters HOCETA for many component pairs are available in the
Aspen Physical Property System. These parameters are retrieved
automatically when you specify any of the following property methods: NRTL-
HOC, UNIF-HOC, UNIQ-HOC, VANL-HOC, and WILS-HOC.
References
J.G. Hayden and J.P. O'Connell, "A Generalized Method for Predicting Second
Virial Coefficients," Ind. Eng. Chem., Process Des. Dev., Vol. 14,No. 3,
(1974), pp. 209 – 216.
HF Equation-of-State
HF forms oligomers in the vapor phase. The non-ideality in the vapor phase is
found in important deviations from ideality in all thermodynamic properties.
The HF equation accounts for the vapor phase nonidealities. The model is
based on chemical theory and assumes the formation of hexamers.
Species like HF that associate linearly behave as single species. For example,
they have a vapor pressure curve, like pure components. The component on
which a hypothetical unreacted system is based is often called the apparent
(or parent) component. Apparent components react to the true species.
Electrolyte Calculation in Physical Property Methods discusses apparent and
The constants C0 and C1 are from Long et al. (1943), and C2 and C3 are set to
0. The correlation is valid between 270 and 330 K, and can be extrapolated to
about 370 K (cf. sec. 4). Different sets of constants can be determined by
experimental data regression.
Where the true number of species is given by 1/f. The association factor is
easily determined from (p, V, T) experiments. For a critical evaluation of data
refer to Vanderzee and Rodenburg (1970).
If only one reaction is assumed for a mixture of HF and its associated species,
(refer to Long et al., 1943), then:
(3)
K is known from Long et al., or can be regressed from (p,V,T) data. The
apparent mole fraction of HF, ya, is known to the user and the simulator, but
p1, or y = p1/p must also be known in order to calculate the thermodynamic
properties of the mixture. Equation 9 must be solved for p1.
Equation 9 can be written as a polynomial in p1 of degree 6:
K(6 - 5ya)(p1)6 + p1 - pya = 0 (9a)
A second order Newton-Raphson technique is used to determine p1. Then p6
can be calculated by equation 5, and f is known (equation 7).
From the fugacity coefficient, the Gibbs energy departure of the mixture or
pure apparent components can be calculated:
(12)
(12a)
(14)
(17)
From the Gibbs energy departure and enthalpy departure, the entropy
departure can be calculated:
(18)
Usage
The HF equation-of-state should only be used for vapor phase calculations. It
is not suited for liquid phase calculations.
References
M. M. Abbott and H. C. van Ness, "Thermodynamics of Solutions Containing
Reactive Species, a Guide to Fundamentals and Applications," Fluid Phase Eq.,
Vol. 77, (1992) pp. 53 – 119.
V. V. De Leeuw and S. Watanasiri, "Modelling Phase Equilibria and Enthalpies
of the System Water and Hydroflouric Acid Using an HF Equation-of-state in
Conjunction with the Electrolyte NRTL Activity Coefficient Model," Paper
presented at the 13th European Seminar on Applied Thermodynamics, June 9
– 12, Carry-le-Rouet, France, 1993.
R. W. Long, J. H. Hildebrand, and W. E. Morrell, "The Polymerization of
Gaseous Hydrogen and Deuterium Flourides," J. Am. Chem. Soc., Vol. 65,
(1943), pp. 182 – 187.
C. E. Vanderzee and W. WM. Rodenburg, "Gas Imperfections and
Thermodynamic Excess Properties of Gaseous Hydrogen Fluoride," J. Chem.
Thermodynamics, Vol. 2, (1970), pp. 461 – 478.
Ideal Gas
The ideal gas law (ideal gas equation-of-state) combines the laws of Boyle
and Gay-Lussac. It models a vapor as if it consisted of point masses without
any interactions. The ideal gas law is used as a reference state for equation-
of-state calculations, and can be used to model gas mixtures at low pressures
(without specific gas phase interactions).
The equation is:
p = RT / Vm
Lee-Kesler
This equation-of-state model is based on the work of Lee and Kesler (1975).
In this equation, the volumetric and thermodynamic properties of fluids based
on the Curl and Pitzer approach (1958) have been analytically represented by
a modified Benedict-Webb-Rubin equation-of-state (1940). The model
calculates the molar volume, enthalpy departure, Gibbs free energy
departure, and entropy departure of a mixture at a given temperature,
pressure, and composition for either a vapor or a liquid phase. Partial
In the above equation both Z(0) and Z(1) are represented as generalized
equations of the BWR form in terms of reduced temperature and pressure.
Thus,
Equations for the enthalpy departure, Gibbs free energy departure, and
entropy departure are obtained from the compressibility factor using standard
thermodynamic relationships, thus ensuring thermodynamic consistency.
In the case of mixtures, mixing rules (without any binary parameters) are
used to obtain the mixture values of the critical temperature and pressure,
and the acentric factor.
This equation has been found to provide a good description of the volumetric
and thermodynamic properties of mixtures containing nonpolar and slightly
polar components.
Symbol Parameter Name Default Definition
Tc TCLK TC Critical temperature
Pc PCLK PC Critical pressure
References
M. Benedict, G. B. Webb, and L. C. Rubin, J. Chem. Phys., Vol. 8, (1940), p.
334.
R. F. Curl and K.S. Pitzer, Ind. Eng. Chem., Vol. 50, (1958), p. 265.
B. I. Lee and M.G. Kesler, AIChE J., Vol. 21, (1975), p. 510.
Lee-Kesler-Plöcker
The Lee-Kesler-Plöcker equation-of-state is the basis for the LK-PLOCK
property method. This equation-of-state applies to hydrocarbon systems that
include the common light gases, such as H2S and CO2. It can be used in gas-
processing, refinery, and petrochemical applications.
The general form of the equation is:
Where:
The fo and fR parameters are functions of the BWR form. The fo parameter is
for a simple fluid, and fR is for reference fluid n-octane.
ω =
Zm =
Where:
Vcij =
Tcij =
kij = kji
The binary parameter kij is determined from phase-equilibrium data
regression, such as VLE data. The Aspen Physical Property System stores the
binary parameters for a large number of component pairs. These binary
parameters are used automatically with the LK-PLOCK property method. If
binary parameters for certain component pairs are not available, they can be
estimated using built-in correlations. The correlations are designed for binary
interactions among the components CO, CO2, N2, H2, CH4 alcohols and
hydrocarbons. If a component is not CO, CO2, N2, H2, CH4 or an alcohol, it is
assumed to be a hydrocarbon.
Parameter Symbol Default MDS Lower Upper Units
Name/ Limit Limit
Element
TCLKP Tci TC x 5.0 2000.0 TEMPERATURE
PCLKP pci PC x PRESSURE
Method 1 is the default for LKPZC; Method 2 can be invoked by setting the
value of LKPZC equal to zero.
Binary interaction parameters LKPKIJ are available for a large number of
components in the Aspen Physical Property System.
References
B.I. Lee and M.G. Kesler, AIChE J., Vol. 21, (1975) p. 510; errata: AIChE J.,
Vol. 21, (1975) p. 1040.
V. Plöcker, H. Knapp, and J.M. Prausnitz, Ind. Eng. Chem., Process Des. Dev.,
Vol. 17, (1978), p. 324.
References
L. Haar, J.S. Gallagher, and J.H. Kell, "NBS/NRC Steam Tables," (Washington:
Hemisphere Publishing Corporation, 1984).
Nothnagel
The Nothnagel equation-of-state calculates thermodynamic properties for the
vapor phase. It is used in property methods NRTL-NTH, UNIQ-NTH, VANL-
NTH, and WILS-NTH. It is recommended for systems that exhibit strong vapor
phase association. The model incorporates the chemical theory of
dimerization to account for strong association and solvation effects, such as
those found in organic acids, like acetic acid. The equation-of-state is:
Where:
b =
bij =
Where:
yi and yj = True mole fractions of monomers
yij = True mole fraction of dimer
ϕi = True fugacity coefficient of component i
Kij = Equilibrium constant for the dimerization of i and j, on a
pressure basis
The sum of the contributions of all dimerization reactions, corrected for the
ratio of apparent and true number of moles, is added to the molar enthalpy
departure:
In this method, you enter Ai, Bi, Ci, and Di on the Properties Parameters
Unary.T-Dependent form. The units for Kii is pressure-1; use absolute units for
temperature. If you enter Kii and Kjj, then Kij is computed from
If you enter Ai, Bi, Ci, and Di, the equilibrium constants are computed using
the parameters you entered. Otherwise the equilibrium constants are
computed using built-in correlations.
Parameter Symbol Default Lower Upper Limit Units
Name/Element Limit
TC Tci — 5.0 2000.0 TEMPERATURE
TB Tbi — 4.0 2000.0 TEMPERATURE
5 8
PC pci — 10 10 PRESSURE
NTHA/1 bi 0.199 RTci / pci 0.01 1.0 MOLE-VOLUME
NTHA/2 di 0.33 0.01 3.0 —
† For the following systems, the values given in Nothnagel et al., 1973 are
used by default:
• Methyl chloride/acetone
• Acetonitrile/acetaldehyde
• Acetone/chloroform
• Chloroform/diethyl amine
• Acetone/benzene
• Benzene/chloroform
• Chloroform/diethyl ether
• Chloroform/propyl formate
• Chloroform/ethyl acetate
• Chloroform/methyl acetate
• Chloroform/methyl formate
• Acetone/dichloro methane
• n-Butane/n-perfluorobutane
• n-Pentane/n-perfluoropentane
• n-Pentane/n-perfluorohexane
References
K.-H. Nothnagel, D. S. Abrams, and J.M. Prausnitz, "Generalized Correlation
for Fugacity Coefficients in Mixtures at Moderate Pressures," Ind. Eng. Chem.,
Process Des. Dev., Vol. 12, No. 1 (1973), pp. 25 – 35.
Peng-Robinson
The Peng-Robinson equation-of-state is the basis for the PENG-ROB and PR-
BM property methods. The model has been implemented with choices of
different alpha functions (see Peng-Robinson Alpha Functions) and has been
extended to include advanced asymmetric mixing rules.
Note: You can choose any of the available alpha functions, but you cannot
define multiple property methods based on this model using different alpha
functions within the same run.
By default, the PENG-ROB property method uses the literature version of the
alpha function and mixing rules. The PR-BM property method uses the
Boston-Mathias alpha function and standard mixing rules. These default
Where:
b =
c =
a = a0+a1
a0 =
kij = kji
a1
In general, .
ai =
bi =
ci =
For best results, the binary parameters kij and lij must be determined from
regression of phase equilibrium data such as VLE data. The Aspen Physical
Property System also has built-in kij and lij for a large number of component
pairs in the EOS-LIT databank. These parameters are used automatically with
the PENG-ROB property method. Values in the databank can be different than
References
D.-Y. Peng and D. B. Robinson, "A New Two-Constant Equation-of-state," Ind.
Eng. Chem. Fundam., Vol. 15, (1976), pp. 59–64.
P.M. Mathias, H.C. Klotz, and J.M. Prausnitz, "Equation of state mixing rules
for multicomponent mixtures: the problem of invariance," Fluid Phase
Equilibria, Vol 67, (1991), pp. 31-44.
Standard Peng-Robinson
The Standard Peng-Robinson equation-of-state is the original formulation of
the Peng-Robinson equation of state with the standard alpha function (see
Peng-Robinson Alpha Functions). It is recommended for hydrocarbon
processing applications such as gas processing, refinery, and petrochemical
processes. Its results are comparable to those of the standard Redlich-
Kwong-Soave equation of state.
Note: You can choose any of the available alpha functions, but you cannot
define multiple property methods based on this model using different alpha
functions within the same run.
The equation for this model is:
a =
ai =
bi =
kij =
The model has option codes which can be used to customize the model, by
selecting a different alpha function and other model options. See Peng-
Robinson Alpha Functions for a description of the alpha functions. See Option
Codes for Equation of State Models (under ESPRSTD) for a list of the option
codes.
For best results, the binary parameter kij must be determined from regression
of phase equilibrium data such as VLE data. The Aspen Physical Property
System also has built-in kij for a large number of component pairs in the EOS-
LIT databank. These parameters are used automatically with the PENG-ROB
property method. Values in the databank can be different than those used
with other models such as Soave-Redlich-Kwong or Redlich-Kwong-Soave,
and this can produce different results.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
TCPRS Tci TC x 5.0 2000.0 TEMPERATURE
5 8
PCPRS pci PC x 10 10 PRESSURE
OMGPRS ωi OMEGA x -0.5 2.0 —
PRKBV/1 kij(1) 0 x - - -
PRKBV/2 kij(2) 0 x - - TEMPERATURE
PRKBV/3 kij(3) 0 x - - TEMPERATURE
PRKBV/4 Tlower 0 x - - TEMPERATURE
PRKBV/5 Tupper 1000 x - - TEMPERATURE
References
D.-Y. Peng and D. B. Robinson, "A New Two-Constant Equation-of-state," Ind.
Eng. Chem. Fundam., Vol. 15, (1976), pp. 59–64.
Note: You can choose any of the available alpha functions, but you cannot
define multiple property methods based on this model using different alpha
functions within the same run.
The PSRK method is an example of modified Huron-Vidal mixing rules. A brief
introduction is provided in Huron-Vidal Mixing Rules. For more details, see
Predictive Soave-Redlich-Kwong-Gmehling Mixing Rules.
Peng-Robinson-Wong-Sandler
This model uses the Peng-Robinson equation-of-state for pure compounds.
The mixing rules are the predictive Wong-Sandler rules. Several alpha
functions can be used in the Peng-Robinson-Wong-Sandler equation-of-state
model for a more accurate description of the pure component behavior. The
pure component behavior and parameter requirements are described in Peng-
Robinson, and in Peng-Robinson Alpha Functions.
Note: You can choose any of the available alpha functions, but you cannot
define multiple property methods based on this model using different alpha
functions within the same run.
The Wong-Sandler mixing rules are an example of modified Huron-Vidal
mixing rules. A brief introduction is provided in Huron-Vidal Mixing Rules. For
more details see Wong-Sandler Mixing Rules., this chapter.
Where:
=
b =
ai =
bi =
References
O. Redlich and J.N.S. Kwong, "On the Thermodynamics of Solutions V. An
Equation-of-state. Fugacities of Gaseous Solutions," Chem. Rev., Vol. 44,
(1979), pp. 223 – 244.
Redlich-Kwong-Aspen
The Redlich-Kwong-Aspen equation-of-state is the basis for the RK-ASPEN
property method. It can be used for hydrocarbon processing applications. It is
also used for more polar components and mixtures of hydrocarbons, and for
light gases at medium to high pressures.
The equation is the same as Redlich-Kwong-Soave:
bi =
kb,ij =
For best results, binary parameters kij must be determined from phase-
equilibrium data regression, such as VLE data.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
TCRKA Tci TC x 5.0 2000.0 TEMPERATURE
PCRKA pci PC x 105 108 PRESSURE
OMGRKA ωi OMEGA x -0.5 2.0 —
References
Mathias, P.M., "A Versatile Phase Equilibrium Equation-of-state", Ind. Eng.
Chem. Process Des. Dev., Vol. 22, (1983), pp. 385 – 391.
Redlich-Kwong-Soave
This is the standard Redlich-Kwong-Soave equation-of-state, and is the basis
for the RK-SOAVE property method. It is recommended for hydrocarbon
processing applications, such as gas-processing, refinery, and petrochemical
Where:
b =
ai =
bi =
kij = kji
; ;
The parameter ai is calculated according to the standard Soave formulation
(see Soave Alpha Functions, equations 1, 2, 3, 5, and 6).
Note: You can choose any of the available alpha functions, but you cannot
define multiple property methods based on this model using different alpha
functions within the same run.
The model uses option codes which are described in Soave-Redlich-Kwong
Option Codes.
For best results, binary parameters kij must be determined from phase-
equilibrium data regression (for example, VLE data). The Aspen Physical
Property System also has built-in kij for a large number of component pairs in
the EOS-LIT databank. These binary parameters are used automatically with
the RK-SOAVE property method. Values of kij in the databank can be different
than those used with other models such as Soave-Redlich-Kwong or Standard
Peng-Robinson, and this can produce different results.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
TCRKSS Tci TC x 5.0 2000.0 TEMPERATURE
References
G. Soave, "Equilibrium Constants for Modified Redlich-Kwong Equation-of-
state," Chem. Eng. Sci., Vol. 27, (1972), pp. 1196 – 1203.
J. Schwartzentruber and H. Renon, "Extension of UNIFAC to High Pressures
and Temperatures by the Use of a Cubic Equation-of-state," Ind. Eng. Chem.
Res., Vol. 28, (1989), pp. 1049 – 1955.
A. Peneloux, E. Rauzy, and R. Freze, "A Consistent Correction For Redlich-
Kwong-Soave Volumes", Fluid Phase Eq., Vol. 8, (1982), pp. 7–23.
Redlich-Kwong-Soave-Boston-Mathias
The Redlich-Kwong-Soave-Boston-Mathias equation-of-state is the basis for
the RKS-BM property method. It is the Redlich-Kwong-Soave equation-of-
state with the Boston-Mathias alpha function (see Soave Alpha Functions). It
is recommended for hydrocarbon processing applications, such as gas-
processing, refinery, and petrochemical processes. Its results are comparable
to those of the Peng-Robinson-Boston-Mathias equation-of-state.
Note: You can choose any of the available alpha functions, but you cannot
define multiple property methods based on this model using different alpha
functions within the same run.
The equation is:
p =
Where:
b =
ai =
bi =
kij = kji
; ;
The parameter ai is calculated by the standard Soave formulation at
supercritical temperatures. If the component is supercritical, the Boston-
Mathias extrapolation is used (see Soave Alpha Functions).
The model uses option codes which are described in Soave-Redlich-Kwong
Option Codes.
For best results, binary parameters kij must be determined from phase-
equilibrium data regression (for example, VLE data).
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
TCRKS Tci TC x 5.0 2000.0 TEMPERATURE
PCRKS pci PC x 105 108 PRESSURE
OMGRKS ωi OMEGA x -0.5 2.0 —
References
G. Soave, "Equilibrium Constants for Modified Redlich-Kwong Equation-of-
state," Chem. Eng. Sci., Vol. 27, (1972), pp. 1196 – 1203.
Note: You can choose any of the available alpha functions, but you cannot
define multiple property methods based on this model using different alpha
functions within the same run.
The Wong-Sandler mixing rules are an example of modified Huron-Vidal
mixing rules. A brief introduction is provided in Huron-Vidal Mixing Rules. For
more details, see Wong-Sandler Mixing Rules.
Redlich-Kwong-Soave-MHV2
This equation-of-state model uses the Redlich-Kwong-Soave equation-of-state
for pure compounds. The predictive MHV2 mixing rules are used. Several
alpha functions can be used in the RK-Soave-MHV2 equation-of-state model
for a more accurate description of the pure component behavior. The pure
component behavior and its parameter requirements are described in
Standard Redlich-Kwong-Soave, and in Soave Alpha Functions.
Note: You can choose any of the available alpha functions, but you cannot
define multiple property methods based on this model using different alpha
functions within the same run.
The MHV2 mixing rules are an example of modified Huron-Vidal mixing rules.
A brief introduction is provided in Huron-Vidal Mixing Rules. For more details,
see MHV2 Mixing Rules.
Schwartzentruber-Renon
The Schwartzentruber-Renon equation-of-state is the basis for the SR-POLAR
property method. It can be used to model chemically nonideal systems with
the same accuracy as activity coefficient property methods, such as the
WILSON property method. This equation-of-state is recommended for highly
non-ideal systems at high temperatures and pressures, such as in methanol
synthesis and supercritical extraction applications.
The equation for the model is:
p =
Where:
b =
c =
ai =
bi =
ci =
ka,ij =
lij =
kb,ij =
ka,ij = ka,ji
lij = -lji
kb,ij = kb,ji
The binary parameters ka,ij, kb,ij, and lij are temperature-dependent. In most
cases, ka,ij0 and lij0 are sufficient to represent the system of interest.
VLE calculations are independent of c. However, c does influence the fugacity
values and can be adjusted to (liquid) molar volumes. For a wide temperature
range, adjust cio to the molar volume at 298.15K or at boiling temperature.
Warning: Using c1i and c2i can cause anomalous behavior in enthalpy and
heat capacity. Their use is strongly discouraged.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
TCRKU Tci TC x 5.0 2000.0 TEMPERATURE
PCRKU pci PC x 105 108 PRESSURE
OMGRKU ωi OMEGA x -0.5 2.0 —
RKUPP0 †† q0i — x — — —
RKUPP1 †† q1i 0 x — — —
RKUPP2 †† q2i 0 x — — —
RKUC0 c0i 0 x — — —
RKUC1 c1i 0 x — — —
RKUC2 c2i 0 x — — —
RKUKA0 ††† ka,ij0 0 x — — —
RKUKA1 ††† ka,ij1 0 x — — TEMPERATURE
RKUKA2 ††† ka,ij2 0 x — — TEMPERATURE †
† Absolute temperature units are assumed for ka,ij2, lij2, and kb,ij2.
†† For polar components (dipole moment >> 0), if you do not enter q0i, the
system estimates q0i, q1i, q2i from vapor pressures using the Antoine vapor
pressure model.
††† If you do not enter at least one of the binary parameters ka,ij0, ka,ij2, lij0,
lij2, kb,ij0, or kb,ij2 the system estimates ka,ij0, ka,ij2, lij0, and lij2 from the UNIFAC
or Hayden O'Connell models.
References
G. Soave, "Equilibrium Constants for Modified Redlich-Kwong Equation-of-
state," Chem. Eng. Sci., Vol. 27, (1972), pp. 1196 - 1203.
J. Schwartzentruber and H. Renon, "Extension of UNIFAC to High Pressures
and Temperatures by the Use of a Cubic Equation-of-State," Ind. Eng. Chem.
Res., Vol. 28, (1989), pp. 1049 – 1955.
A. Peneloux, E. Rauzy, and R. Freze, "A Consistent Correction For Redlich-
Kwong-Soave Volumes", Fluid Phase Eq., Vol. 8, (1982), pp. 7–23.
Soave-Redlich-Kwong
The Soave-Redlich-Kwong equation-of-state is the basis of the SRK property
method. This model is based on the same equation of state as the Redlich-
Kwong-Soave model. However, this model has several important differences.
• A volume translation concept introduced by Peneloux and Rauzy is used to
improve molar liquid volume calculated from the cubic equation of state.
• Improvement in water properties is achieved by using the NBS Steam
Table.
• Improvement in speed of computation for equation based calculation is
achieved by using composition independent fugacity.
• Optional Kabadi-Danner mixing rules for improved phase equilibrium
calculations in water-hydrocarbon systems (see SRK-Kabadi-Danner)
• Optional Mathias alpha function
Note: You can choose any of the available alpha functions, but you cannot
define multiple property methods based on this model using different alpha
functions within the same run.
This equation-of-state can be used for hydrocarbon systems that include the
common light gases, such as H2S, CO2 and N2.
Where:
Where:
;
a1 is an additional, asymmetric (polar) term:
Where:
; ;
The enthalpy, entropy, Gibbs energy, and molar volume of water are
calculated from the steam tables. The total properties are mole-fraction
averages of these values with the properties calculated by the equation of
state for other components. Fugacity coefficient is not affected.
For best results, the binary parameter kij must be determined from phase
equilibrium data regression (for example, VLE data). The Aspen Physical
Property System also has built-in kij for a large number of component pairs in
the SRK-ASPEN databank. These parameters are used automatically with the
SRK property method. Values of kij in the databank can be different than
those used with other models such as Standard Redlich-Kwong-Soave or
Standard Redlich-Kwong-Soave, and this can produce different results.
The model uses option codes which are described in Soave-Redlich-Kwong
Option Codes.
References
G. Soave, "Equilibrium Constants for Modified Redlich-Kwong Equation-of-
state," Chem. Eng. Sci., Vol. 27, (1972), pp. 1196 - 1203.
A. Peneloux, E. Rauzy, and R. Freze, "A Consistent Correction For Redlich-
Kwong-Soave Volumes", Fluid Phase Eq., Vol. 8, (1982), pp. 7–23.
P.M. Mathias, H.C. Klotz, and J.M. Prausnitz, "Equation of state mixing rules
for multicomponent mixtures: the problem of invariance," Fluid Phase
Equilibria, Vol 67, (1991), pp. 31-44.
SRK-Kabadi-Danner
The SRK-Kabadi-Danner property model uses the SRK equation-of-state with
improved phase equilibrium calculations for mixtures containing water and
hydrocarbons. These improvements are achieved by using the Kabadi-Danner
mixing rules.
The form of the equation-of-state is:
Where:
;
The best values of kwj (w = water) were obtained from experimental data.
Results are given for seven homologous series.
Where:
† This value is obtained from very little data. Might not be reliable.
The model uses option codes which are described in Soave-Redlich-Kwong
Option Codes.
SRK-Kabadi-Danner uses the same parameters as SRK, with added
interaction parameters. Do not specify values for the SRKLIJ parameters
when using SRK-KD.
Parameter Symbol Default MDS Lower Upper Units
Name/ Limit Limit
Element
SRKTC Tci TC x 5.0 2000.0 TEMPERATURE
References
V. Kabadi and R. P. Danner, "A Modified Soave-Redlich-Kwong Equation of
State for Water-Hydrocarbon Phase Equilibria", Ind. Eng. Chem. Process Des.
Dev., Vol. 24, No. 3, (1985), pp. 537-541.
SRK-ML
The SRK-ML property model is the same as the standard Soave-Redlich-
Kwong model with two exceptions:
• kij may not equal kji; they are unsymmetric, and a different set of
parameters are used, as shown below.
• The lij are calculated from the equation lij = kji - kij
VPA/IK-CAPE Equation-of-State
The VPA/IK-CAPE equation of state is similar to the HF equation of state but
allows dimerization, tetramerization and hexamerization to occur
simultaneously. The main assumption of the model is that only molecular
association causes the gas phase nonideality. Attractive forces between the
molecules and the complexes are neglected.
There are three kinds of associations, which can be modeled:
• Dimerization (examples: formic acid, acetic acid)
• Tetramerization (example: acetic acid)
• Hexamerization (example: hydrogen fluoride)
To get the largest possible flexibility of the model all these kinds of
association can occur simultaneously, for example, in a mixture containing
acetic acid and HF. Up to five components can associate, and any number of
inert components are allowed. This is the only difference between this model
and the HF equation of state, which account only the hexamerization of HF.
Symbols
In the following description, these symbols are used:
yi = Apparent concentration
zin = True concentration, for component i and degree of
association n=1, 2, 4, 6
(2)
(4)
By setting
(5)
(6)
(7)
(8)
(9)
Fugacity Coefficient
As is well-known from thermodynamics, the fugacity coefficient can be
calculated by
(11)
with
(13)
and
(16)
(2)
(18)
(20)
and
(23)
(24)
Using the association model, more different species occur than can be
distinguished. Thus, the equivalent expression for the entropy of mixing
should be written with the true concentrations. As eq. 24 refers to 1 mole
monomers, the expression should be weighted by the compressibility factor
representing the true number of moles. The new expression is
(25)
and, analogously,
(27)
References
M. M. Abbott and H. C. van Ness, "Thermodynamics of Solutions Containing
Reactive Species, a Guide to Fundamentals and Applications," Fluid Phase Eq.,
Vol. 77, (1992) pp. 53–119.
V. V. De Leeuw and S. Watanasiri, "Modeling Phase Equilibria and Enthalpies
of the System Water and Hydrofluoric Acid Using an HF Equation-of-state in
Conjunction with the Electrolyte NRTL Activity Coefficient Model," Paper
Presented at the 13th European Seminar on Applied Thermodynamics, June
9–12, Carry-le-Rouet, France, 1993.
(2)
Note: You can choose any of the alpha functions described here, but you
cannot define multiple property methods based on this model using different
alpha functions within the same run.
This approach was also adopted by Peng and Robinson:
(4)
With
=
For c2,i = 0 and c3,i = 0, this expression reduces to the standard Peng-
Robinson formulation if c2,i = mi. You can use vapor pressure data if the
temperature is subcritical to regress the constants. If the temperature is
supercritical, c2,i and c3,i are set to 0.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
TCPR Tci TC X 5.0 2000.0 TEMPERATURE
PCPR pci PC X 105 108 PRESSURE
Where mi is computed by equation 5. The polar parameters p1,i, p2,i and p3,i
are comparable with the c parameters of the Mathias-Copeman expression.
Equation 8 reduces to the standard Peng-Robinson formulation if the polar
parameters are zero. Equation 8 is used only for below-critical temperatures.
For above-critical temperatures, the Boston-Mathias extrapolation is used.
Use equation 6 with:
(9)
(10)
PRSRP/1 — X — — —
PRSRP/2 0 X — — —
PRSRP/3 0 X — — —
References
J. F. Boston and P.M. Mathias, "Phase Equilibria in a Third-Generation Process
Simulator" in Proceedings of the 2nd International Conference on Phase
Equilibria and Fluid Properties in the Chemical Process Industries, West Berlin,
(17-21 March 1980) pp. 823-849.
D.-Y. Peng and D.B. Robinson, "A New Two-Constant Equation-of-state," Ind.
Eng. Chem. Fundam., Vol. 15, (1976), pp. 59-64.
P.M. Mathias and T.W. Copeman, "Extension of the Peng-Robinson Equation-
of-state To Complex Mixtures: Evaluation of the Various Forms of the Local
Composition Concept",Fluid Phase Eq., Vol. 13, (1983), p. 91.
J. Schwartzentruber, H. Renon, and S. Watanasiri, "K-values for Non-Ideal
Systems:An Easier Way," Chem. Eng., March 1990, pp. 118-124.
G. Soave, "Equilibrium Constants for a Modified Redlich-Kwong Equation-of-
state," Chem Eng. Sci., Vol. 27, (1972), pp. 1196-1203.
C.H. Twu, J. E. Coon, and J.R. Cunningham, "A New Cubic Equation of State,"
Fluid Phase Equilib., Vol. 75, (1992), pp. 65-79.
C.H. Twu, D. Bluck, J.R. Cunningham, and J.E. Coon, "A Cubic Equation of
State with a New Alpha Function and a New Mixing Rule," Fluid Phase Equilib.,
Vol. 69, (1991), pp. 33-50.
(2)
Note: You can choose any of the alpha functions described here, but you
cannot define multiple property methods based on this model using different
alpha functions within the same run.
In the Redlich-Kwong equation-of-state, alpha is:
Soave Modification
The parameter αi is a temperature function introduced by Soave in the
Redlich-Kwong equation-of-state to improve the correlation of the pure
component vapor pressure:
(5)
Equation 3 still holds. The parameter mi can be correlated with the acentric
factor:
(6)
Boston-Mathias Extrapolation
For light gases at high reduced temperatures (> 5), equation 5 gives
unrealistic results. The boundary conditions are that attraction between
molecules should vanish for extremely high temperatures, and reduces
asymptotically to zero. Boston and Mathias derived an alternative function for
temperatures higher than critical:
(7)
With
di =
ci =
Where:
(10)
† RKAPOL for Redlich-Kwong-Aspen, SRKPOL for SRK and SRKKD, SRKMLP for
SRK-ML, RKSPOL for RKS-BM, or RKSSPO for RK-SOAVE.
Where mi is computed by equation 6. If the polar parameters p1,i, p2,i and p3,i
are zero, equation 11 reduces to the standard Redlich-Kwong-Soave
formulation. You can use vapor pressure data to regress the constants if the
temperature is subcritical. Equation 11 is used only for temperatures below
critical.
The Boston-Mathias extrapolation is used for temperatures above critical, that
is, with:
(12)
For c2,i=0 and c3,i=0 this expression reduces to the standard Redlich-Kwong-
Soave formulation if c1,i=mi. If the temperature is subcritical, use vapor
pressure data to regress the constants. If the temperature is supercritical, set
c2,i and c3,i to 0.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
TCRKS Tci TC X 5.0 2000.0 TEMPERATURE
5 8
PCRKS pci PC X 10 10 PRESSURE
RKSMCP/1 c1,i — X — — —
RKSMCP/2 c2,i 0 X — — —
RKSMCP/3 c3,i 0 X — — —
Where mi is computed by equation 6 and the polar parameters p1,i, p2,i and p3,i
are comparable with the c parameters of the Mathias-Copeman expression.
Equation 15 reduces to the standard Redlich-Kwong-Soave formulation if the
polar parameters are zero. Equation 15 is very similar to the extended
Mathias equation, but it is easier to use in data regression. It is used only for
temperatures below critical. The Boston-Mathias extrapolation is used for
temperatures above critical, that is, use equation 7 with:
(16)
References
J. F. Boston and P.M. Mathias, "Phase Equilibria in a Third-Generation Process
Simulator" in Proceedings of the 2nd International Conference on Phase
Equilibria and Fluid Properties in the Chemical Process Industries, West Berlin,
(17-21 March 1980), pp. 823-849.
P. M. Mathias, "A Versatile Phase Equilibrium Equation-of-state", Ind. Eng.
Chem. Process Des. Dev., Vol. 22, (1983), pp. 385–391.
P.M. Mathias and T.W. Copeman, "Extension of the Peng-Robinson Equation-
of-state To Complex Mixtures: Evaluation of the Various Forms of the Local
Composition Concept", Fluid Phase Eq., Vol. 13, (1983), p. 91.
O. Redlich and J. N. S. Kwong, "On the Thermodynamics of Solutions V. An
Equation-of-state. Fugacities of Gaseous Solutions," Chem. Rev., Vol. 44,
(1949), pp. 223–244.
J. Schwartzentruber and H. Renon, "Extension of UNIFAC to High Pressures
and Temperatures by the Use of a Cubic Equation-of-state," Ind. Eng. Chem.
Res., Vol. 28, (1989), pp. 1049–1055.
J. Schwartzentruber, H. Renon, and S. Watanasiri, "K-values for Non-Ideal
Systems:An Easier Way," Chem. Eng., March 1990, pp. 118-124.
G. Soave, "Equilibrium Constants for a Modified Redlich-Kwong Equation-of-
state," Chem Eng. Sci., Vol. 27, (1972), pp. 1196-1203.
C.H. Twu, W.D. Sim, and V. Tassone, "Getting a Handle on Advanced Cubic
Equations of State", Chemical Engineering Progress, Vol. 98 #11 (November
2002) pp. 58-65.
(3)
Where:
(5)
Redlich-Kwong-Soave 1 0
This expression can be used at any pressure as a mixing rule for the
parameter. The mixing rule for b is fixed by equation 3. Even when used at
other pressures, this expression contains the excess Gibbs energy at infinite
pressure. You can use any activity coeffecient model to evaluate the excess
Gibbs energy at infinite pressure. Binary interaction coefficients must be
regressed. The mixing rule used contains as many binary parameters as the
activity coefficient model chosen.
This mixing rule has been used successfully for polar mixtures at high
pressures, such as systems containing light gases. In theory, any activity
coefficient model can be used. But the NRTL equation (as modified by Huron
and Vidal) has demonstrated better performance.
The Huron-Vidal mixing rules combine extreme flexibility with thermodynamic
consistency, unlike many other mole-fraction-dependent equation-of-state
mixing rules. The Huron-Vidal mixing rules do not allow flexibility in the
description of the excess molar volume, but always predict reasonable excess
volumes.
The Huron-Vidal mixing rules are theoretically incorrect for low pressure,
because quadratic mole fraction dependence of the second virial coefficient (if
derived from the equation-of-state) is not preserved. Since equations of state
are primarily used at high pressure, the practical consequences of this
drawback are minimal.
The Gibbs energy at infinite pressure and the Gibbs energy at an arbitrary
high pressure are similar. But the correspondence is not close enough to
make the mixing rule predictive.
There are several methods for modifying the Huron-Vidal mixing rule to make
it more predictive. The following three methods are used in Aspen Physical
Property System equation-of-state models:
References
M.- J. Huron and J. Vidal, "New Mixing Rules in Simple Equations of State for
representing Vapour-liquid equilibria of strongly non-ideal mixtures," Fluid
Phase Eq., Vol. 3, (1979), pp. 255-271.
The advantage is that the expressions for mixture and pure component
fugacities do not contain the pressure. They are functions of compacity V/b
and α:
(2)
Where:
(3)
and
(4)
With:
(5)
The constants λ1 and λ2, which depend only on the equation-of-state (see
Huron-Vidal Mixing Rules) occur in equations 2 and 4.
Instead of using infinite pressure for simplification of equation 1, the condition
of zero pressure is used. At p = 0 an exact relationship between the
Predictive Soave-Redlich-Kwong-Gmehling
Mixing Rules
These mixing rules by Holderbaum and Gmehling (1991) use a relationship
between the excess Helmholtz energy and equation-of-state. They do not use
a relationship between equation-of-state properties and excess Gibbs energy,
as in the Huron-Vidal mixing rules. The pressure-explicit expression for the
equation-of-state is substituted in the thermodynamic equation:
Where both Ai* and Am are calculated by using equation 1. Ai* and Am are
written in terms of equation-of-state parameters.
The simplification of constant packing fraction (Vm / b) is used:
(3)
With:
(4)
Therefore:
(5)
Where Λ' is slightly different from Λ for the Huron-Vidal mixing rule:
(7)
The result is a predictive mixing rule for cubic equations of state. But the
original UNIFAC formulation gives the best performance for any binary pair
with interactions available from UNIFAC. Gas-solvent interactions are
unavailable.
At the University of Oldenburg in Germany, the UNIFAC groups were
extended with often-occurring gases. New group interactions were
determined from gas-solvent data, specific to the Redlich-Kwong-Soave
equation-of-state. The new built-in parameters to the Aspen Physical Property
System are activated when using the PSRK equation-of-state model.
The PSRK method has a lot in common with the Huron-Vidal mixing rules. The
mole fraction is dependent on the second virial coefficient and excess volume
is predicted. These are minor disadvantages.
References
K. Fischer, "Die PSRK-Methode: Eine Zustandsgleichung unter Verwendung
des UNIFAC-Gruppenbeitragsmodells," (Düsseldorf: VDI Fortschrittberichte,
Reihe 3: Verfahrenstechnik, Nr. 324, VDI Verlag GmbH, 1993).
T. Holderbaum and J. Gmehling, "PSRK: A Group Contribution Equation-of-
state based on UNIFAC," Fluid Phase Eq., Vol. 70, (1991), pp. 251-265.
Where both Ai* and Am are calculated by using equation 1. Ai* and Am are
written in terms of equation-of-state parameters.
Like Huron and Vidal, the limiting case of infinite pressure is used. This
simplifies the expressions for Ai* and Am. Equation 2 becomes:
(3)
With:
(5)
The relationship between the equation-of-state at low pressure and the virial
coefficient is:
(6)
(7)
Wong and Sandler discovered the following mixing rule to satisfy equation 4
(using equations 6 and 7):
References
D. S. Wong and S. I. Sandler, "A Theoretically Correct New Mixing Rule for
Cubic Equations of State for Both Highly and Slightly Non-ideal Mixtures,"
AIChE J., Vol. 38, (1992), pp. 671 – 680.
Parameters β(0), β(1), β(2), β(3), and θ have five elements to account for
temperature dependencies. Elements P1 through P5 follow the temperature
dependency relation:
Where:
Tref = 298.15K
The user must:
• Supply these elements using a Properties Parameters Binary T-Dependent
form.
• Specify Comp ID i and Comp ID j on this form, using the same order that
appears on the Components Specifications Selection sheet form.
Parameter Name Symbol No. of Elements Default Units
Ionic Unary Parameters
GMBPB βion 1 0 —
GMBPD δion 1 0 —
Cation-Anion Parameters
GMBPB0 β(0) 5 0 —
GMBPB1 β(1) 5 0 —
GMBPB2 β(2) 5 0 —
GMBPB3 β(3) 5 0 —
Cation-Cation Parameters
GMBPTH θcc' 5 0 —
Anion-Anion Parameters
GMBPTH θaa' 5 0 —
Molecule-Ion and Molecule-Molecule Parameters
GMBPB0 β(0) 5 0 —
GMBPB1 β(1) 5 0 —
Working Equations
The complete Pitzer equation (Fürst and Renon, 1982) for the excess Gibbs
energy is (see also equation 4):
Where:
xi = Mole fraction of ion i
xw = Mole fraction of water
Mw = Molecular weight of water (g/mol)
ni = Moles of ion i
The C term and the Ψ term are dropped from equation 1 to give the
simplified Pitzer equation.
(2)
Where:
Bij = f(βij(0),βij(1),βij(2),βij(3))
Therefore, the simplified Pitzer equation has two types of binary interaction
parameters, β 's and θ''s. There are no ternary interaction parameters with
the simplified Pitzer equation.
Note that the Pitzer model parameter databank described in Physical Property
Data, Chapter 1, is not applicable to the simplified Pitzer equation.
A built-in empirical correlation estimates the β(0) and β(1) parameters for
cation-anion pairs from the Bromley ionic parameters, βion and δion (Bromley,
1973). The estimated values of β(0)'s and β(1)'s are overridden by the user's
References
L.A. Bromley, "Thermodynamic Properties of Strong Electrolytes in Aqueous
Solution, " AIChE J., Vol. 19, No. 2, (1973), pp. 313 – 320.
W. Fürst and H. Renon, "Effects of the Various Parameters in the Application
of Pitzer's Model to Solid-Liquid Equilibrium. Preliminary Study for Strong 1-1
Electrolytes," Ind. Eng. Chem. Process Des. Dev., Vol. 21, No. 3, (1982),
pp. 396-400.
Parameter Conversion
For n-m electrolytes, n and m>1 (2-2, 2-3, 3-4, and so on), the parameter
β(3) corresponds to Pitzer's β(1); β(2) is the same in both Aspen Physical
Property System and original Pitzer models. Pitzer refers to the n-m
electrolyte parameters as β(1), β(2), β(0). β(0) and β(2) retain their meanings in
both models, but Pitzer's β(1) is Aspen Physical Property System β(3). Be
careful to make this distinction when entering n-m electrolyte parameters.
Chien-Null
The Chien-Null model calculates liquid activity coefficients and it can be used
for highly non-ideal systems. The generalized expression used in its derivation
can be adapted to represent other well known formalisms for the activity
coefficient by properly defining its binary terms. This characteristic allows the
model the use of already available binary parameters regressed for those
other liquid activity models with thermodynamic consistency.
The equation for the Chien-Null liquid activity coeficient is:
Where:
Rji = Aji / Aij
Aii = 0
Aji = 2τjiGji
Vji = Gji
The following are defined for the Non-Random Two-Liquid activity coefficient
model, where:
γi = ai
Parameter Symbol Default MDS Upper Lower Units
Name/Element Limit Limit
GMCONS ai 1.0 x — — —
COSMO-SAC
Cosmo-SAC is a solvation model that describes the electric fields on the
molecular surface of species that are polarizable. It requires a fairly
complicated quantum mechanical calculation, but this calculation must be
done only once for a particular molecule; then the results can be stored. In its
final form, it uses individual atoms as the building blocks for predicting phase
equilibria instead of functional groups. This model formulation provides a
considerably larger range of applicability than group-contribution methods.
The calculation for liquid nonideality is only slightly more computationally
intensive than activity-coefficient models such as NRTL or UNIQUAC. Cosmo-
SAC complements the UNIFAC group-contribution method, because it is
applicable to virtually any mixture.
The Cosmo-SAC model calculates liquid activity coefficients. The equation for
the Cosmo-SAC model is:
With
The Cosmo-SAC model does not require binary parameters. For each
component, it has six input parameters. CSACVL is the component volume
parameter which is always defined in cubic angstroms, regardless of chosen
units sets. SGPRF1 to SGPRF5 are five component sigma profile parameters;
each can store up to 12 points of sigma profile values. All six input
parameters are obtained from COSMO calculation. The Aspen Physical
Property System includes a database of sigma profiles for over 1400
compounds from Mullins et al. (2006). The parameters were obtained by
permission from the Virginia Tech Sigma Profile Database website
(http://www.design.che.vt.edu/VT-2004.htm). Aspen Technology, Inc. does
not claim proprietary rights to these parameters.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
CSACVL Vi — x — — VOLUME (Å3)
SGPRF1 Ai pi(1-12) — x — — —
SGPRF2 Ai pi(13-24) — x — — —
SGPRF3 Ai pi(25-36) — x — — —
SGPRF4 Ai pi(37-48) — x — — —
SGPRF5 Ai pi(49-51) — x — — —
Option Codes
The primary version of COSMO-SAC is the model by Lin and Sandler (2002).
Two other versions are available using an option code, as detailed in the table
below:
References
A. Klamt and F. Eckert, "COSMO-RS: A Novel and Efficient Method for the a
Priori Prediction of Thermophysical Data of Liquids," Fluid Phase Equilibria 43,
172 (2000).
S.-T. Lin, P. M. Mathias, Y. Song, C.-C. Chen, and S. I. Sandler,
"Improvements of Phase-Equilibrium Predictions for Hydrogen-Bonding
Systems from a New Expression for COSMO Solvation Models," presented at
the AIChE Annual Meeting, Indianapolis, IN, 3-8 November (2002).
S.-T. Lin and S. I. Sandler, "A Priori Phase Equilibrium Prediction from a
Segment Contribution Solvation Model," Ind. Eng. Chem. Res. 41, 899
(2002).
E. Mullins, et al. "Sigma-Profile Database for Using COSMO-Based
Thermodynamic Methods," Ind. Eng. Chem. Res. 45, 4389 (2006).
Where:
Tref = 298.15K
Many parameter pairs are included in the electrolyte NRTL model parameter
databank (see Physical Property Data, Chapter 1).
Option codes can affect the performance of this model. See Option Codes for
Activity Coefficient Models for details.
Parameter Symbol No. of Default MDS Units
Name Elements
Dielectric Constant Unary Parameters
CPDIEC AB 1 — — —
BB 1 0 — —
CB 1 298.15 — TEMPERATURE†
NRTL/4 — — 0 x TEMPERATURE
NRTL/5 FBB' — 0 x TEMPERATURE
FB'B — 0 x TEMPERATURE
NRTL/6 GBB' — 0 x TEMPERATURE
GB'B — 0 x TEMPERATURE
Electrolyte-Molecule Pair Parameters
GMELCC Cca,B 1 0 x —
CB,ca 1 0 x —
GMELCD Dca,B 1 0 x TEMPERATURE†
DB,ca 1 0 x TEMPERATURE†
GMELCE Eca,B 1 0 x —
EB,ca 1 0 x —
GMELCN αca,B = αB,ca 1 .2 x —
Electrolyte-Electrolyte Pair Parameters
GMELCC Cca',ca'' 1 0 x —
Cca'',ca' 1 0 x —
Cc'a,c''a 1 0 x —
Cc''a,c'a 1 0 x —
†
GMELCD Dca',ca'' 1 0 x TEMPERATURE
†
Dca'',ca' 1 0 x TEMPERATURE
Dc'a,c''a 1 0 x TEMPERATURE†
Dc''a,c'a 1 0 x TEMPERATURE†
GMELCE Eca',ca'' 1 0 x —
Eca'',ca' 1 0 x —
Ec'a,c''a 1 0 x —
Ec''a,c'a 1 0 x —
GMELCN αca',ca'' = αca'',ca' 1 .2 x —
αc'a,c''a = αc''a,c'a 1 .2 x —
†
Certain Electrolyte NRTL activity coefficient model parameters are used with
reciprocal temperature terms:
• CPDIEC
• NRTL/2
• GMELCD for electrolyte-electrolyte or electrolyte-molecule pairs
This leads to
(2)
References
C.-C. Chen, H.I. Britt, J.F. Boston, and L.B. Evans, "Local Compositions Model
for Excess Gibbs Energy of Electrolyte Systems: Part I: Single Solvent, Single
Completely Dissociated Electrolyte Systems:, AIChE J., Vol. 28, No. 4, (1982),
p. 588-596.
C.-C. Chen, and L.B. Evans, "A Local Composition Model for the Excess Gibbs
Energy of Aqueous Electrolyte Systems," AIChE J., Vol. 32, No. 3, (1986), p.
444-459.
B. Mock, L.B. Evans, and C.-C. Chen, "Phase Equilibria in Multiple-Solvent
Electrolyte Systems: A New Thermodynamic Model," Proceedings of the 1984
Summer Computer Simulation Conference, p. 558.
B. Mock, L.B. Evans, and C.-C. Chen, "Thermodynamic Representation of
Phase Equilibria of Mixed-Solvent Electrolyte Systems," AIChE J., Vol. 32, No.
10, (1986), p. 1655-1664.
with
(4)
(5)
Where:
Aϕ = Debye-Hückel parameter
NA = Avogadro's number
ds = Mass density of solvent
Qe = Electron charge
T = Temperature
k = Boltzmann constant
Ix = Ionic strength (mole fraction scale)
xi = Mole fraction of component i
zi = Charge number of ion i
The Born equation is used to account for the Gibbs energy of transfer of ionic
species from the infinite dilution state in a mixed-solvent to the infinite
dilution state in aqueous phase.
(7)
Where:
The Debye-Hückel theory is based on the infinite dilution reference state for
ionic species in the actual solvent media. For systems with water as the only
solvent, the reference state is the infinite dilution aqueous solution. For
mixed-solvent systems, the reference state for which the Debye-Hückel
theory remains valid is the infinite dilution solution with the corresponding
mixed-solvent composition. However, the molecular weight Ms, the mass
density ds, and the dielectric constant εs for the single solvent need to be
(7a)
(8a)
(8b)
(8c)
(8d)
Where:
xm = Mole fraction of the solvent m in the solution
Mm = Molecular weight of the solvent m
Vml = Molar volume of the solvent mixture
Where:
Xj = xjCj
(Cj = zj for ions and Cj = unity for molecules)
Gji =
τji =
gji and gii are energies of interaction between species j and i, and i and i,
respectively. Both gij and αij are inherently symmetric (gij = gji and αij = αji).
Similarly,
(10)
Where:
Gji,ki =
τji,ki =
i = c, a, or B
(15)
(16)
To obtain an expression for the excess Gibbs energy, let the residual Gibbs
energies, per mole of cells of central cation, anion, or solvent, respectively, be
, , and . These are then related to the
effective local mole fractions:
(17)
(18)
(19)
(21)
(22)
Where:
Using the previous relation for the excess Gibbs energy and the expressions
for the residual and reference Gibbs energy (equations 17 to 19 and 20 to
22), the following expression for the excess Gibbs energy is obtained:
(24)
Combining this expression with the expression for the effective local mole
fractions given in equations 9 and 10, the following equality is obtained:
(26)
(28)
and,
(29)
(31)
The binary parameters αca,B , τca,B and τB,ca are now the adjustable parameters
for an apparent binary system of a single electrolyte and a single solvent.
The excess Gibbs energy expression (equation 24) must now be normalized
to the infinite dilution reference state for ions:
(32)
(35)
(36)
Multicomponent Systems
The Electrolyte NRTL model can be extended to handle multicomponent
systems.
The excess Gibbs energy expression is:
(37)
Where:
j and k can be any species (a, C, or B)
The activity coefficient equation for molecular components is given by:
Where:
(41)
(42)
(43)
(45)
(46)
(47)
(48)
Parameters
The model adjustable parameters include:
• Pure component dielectric constant coefficient of nonaqueous solvents
• Born radius of ionic species
• NRTL interaction parameters for molecule-molecule, molecule-electrolyte,
and electrolyte-electrolyte pairs
Note that for the electrolyte-electrolyte pair parameters, the two electrolytes
must share either one common cation or one common anion.
Each type of the electrolyte NRTL parameter consists of both the
nonrandomness factor, α, and energy parameters, τ.
The pure component dielectric constant coefficients of nonaqueous solvents
and Born radius of ionic species are required only for mixed-solvent
electrolyte systems.
The temperature dependency relations of these parameters are given in
Electrolyte NRTL Activity Coefficient Model.
Heat of mixing is calculated from temperature derivatives of activity
coefficients. Heat capacity is calculated from secondary temperature
derivative of the activity coefficient. As a result, the temperature dependent
parameters are critical for modeling enthalpy correctly. It is recommended
that enthalpy data and heat capacity data be used to obtain these
temperature dependency parameters. See also Electrolyte NRTL Enthalpy and
Electrolyte NRTL Gibbs Energy.
ENRTL-SAC
eNRTL-SAC (ENRTLSAC, patent pending) is an extension of the nonrandom
two-liquid segment activity coefficient model (NRTL-SAC, patent pending) by
Chen and Song (Ind. Eng. Chem. Res., 2004, 43, 8354) to include
electrolytes in the solution. It can be used in usable in Aspen Properties and
Aspen Polymers. It is intended for the computation of ionic activity
coefficients and solubilities of electrolytes, organic and inorganic, in common
solvents and solvent mixtures. In addition to the three types of molecular
parameters defined for organic nonelectrolytes in NRTL-SAC (hydrophobicity
X, hydrophilicity Z, and polarity Y- and Y+), an electrolyte parameter, E, is
introduced to characterize both local and long-range ion-ion and ion-molecule
interactions attributed to ionized segments of electrolytes.
In applying the segment contribution concept to electrolytes, a new
conceptual electrolyte segment e corresponding to the electrolyte parameter
E, is introduced. This conceptual segment e would completely dissociate to a
cationic segment (c) and an anionic segment (a), both of unity charge. All
electrolytes, organic or inorganic, symmetric or unsymmetric, univalent or
multivalent, are to be represented with this conceptual uni-univalent
With
Aϕ = Debye-Hückel parameter
τ21 5 -3 -3 -4.549
Option codes
There are three option codes in ENRTLSAC. The first is used to enable or
disable the Flory-Huggins term. The other two are only used internally and
you should not change their values. The Flory-Huggins term is included by
default in eNRTL-SAC model. You can remove this term using the first option
code. The table below lists the values for the first option code.
0 Flory-Huggins term included (default)
References
C.-C. Chen and Y. Song, "Solubility Modeling with a Nonrandom Two-Liquid
Segment Activity Coefficient Model," Ind. Eng. Chem. Res. 43, 8354 (2004).
C.-C. Chen and Y. Song, "Extension of Nonrandom Two-Liquid Segment
Activity Coefficient Model for Electrolytes," Ind. Eng. Chem. Res. 44, 8909
(2005).
Hansen
Hansen is a solubility parameter model and is commonly used in the solvent
selection process. It is based on the regular solution theory and Hansen
solubility parameters. This model has no binary parameters and its application
merely follows the empirical guide like dissolves like.
The Hansen model calculates liquid activity coefficients. The equation for the
Hansen model is:
with
Where:
HANVOL Vi — x — — VOLUME
Option codes
The Hansen volume is implemented as an input parameter. If the Hansen
volume is not input by the user it will be calculated by an Aspen Plus internal
method. You can also request the Aspen Plus method using Option Codes in
Aspen Plus Interface. The table below lists the option codes.
Reference
Frank, T. C.; Downey, J. R.; Gupta, S. K. "Quickly Screen Solvents for
Organic Solids," Chemical Engineering Progress 1999, December, 41.
Hansen, C. M. Hansen Solubility Parameters: A User’s Handbook; CRC Press,
2000.
Ideal Liquid
This model is used in Raoult's law. It represents ideality of the liquid phase.
This model can be used for mixtures of hydrocarbons of similar carbon
number. It can be used as a reference to compare the results of other activity
coefficient models.
The equation is:
ln γi = 0
Where:
Gij =
τij =
αij =
τii = 0
Gii = 1
aij, bij, eij, and fij are unsymmetrical. That is, aij may not be equal to aji, etc.
The binary parameters aij, bij, cij, dij, eij, and fij can be determined from VLE
and/or LLE data regression. The Aspen Physical Property System has a large
number of built-in binary parameters for the NRTL model. The binary
parameters have been regressed using VLE and LLE data from the Dortmund
Databank. The binary parameters for the VLE applications were regressed
using the ideal gas, Redlich-Kwong, and Hayden O'Connell equations of state.
See Physical Property Data, Chapter 1, for details.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
NRTL/1 aij 0 x -100.0 100.0 —
NRTL/2 bij 0 x -30000 30000.0 TEMPERATURE
NRTL/3 cij 0.30 x 0.0 1.0 —
NRTL/4 dij 0 x -0.02 0.02 TEMPERATURE
NRTL/5 eij 0 x — — TEMPERATURE
NRTL/6 fij 0 x — — TEMPERATURE
The NRTL-2 property method uses data set 2 for NRTL. All other NRTL
methods use data set 1.
NRTL-SAC
NRTL-SAC (patent pending) is a segment contribution activity coefficient
model, derived from the Polymer NRTL model, usable in Aspen Properties and
Aspen Polymers. NRTL-SAC can be used for fast, qualitative estimation of the
solubility of complex organic compounds in common solvents. Conceptually,
the model treats the liquid non-ideality of mixtures containing complex
organic molecules (solute) and small molecules (solvent) in terms of
interactions between three pairwise interacting conceptual segments:
hydrophobic segment (x), hydrophilic segment (z), and polar segments (y-
and y+). In practice, these conceptual segments become the molecular
descriptors used to represent the molecular surface characteristics of each
solute or solvent molecule. Hexane, water, and acetonitrile are selected as
the reference molecules for the hydrophobic, hydrophilic, and polar segments,
respectively. The molecular parameters for all other solvents can be
determined by regression of available VLE or LLE data for binary systems of
solvent and the reference molecules or their substitutes. The treatment
results in four component-specific molecular parameters: hydrophobicity X,
hydrophilicity Z, and polarity Y- and Y+. The two types of polar segments, Y-
and Y+, are used to reflect the wide variations of interactions between polar
molecules and water.
The conceptual segment contribution approach in NRTL-SAC represents a
practical alternative to the UNIFAC functional group contribution approach.
This approach is suitable for use in the industrial practice of carrying out
measurements for a few selected solvents and then using NRTL-SAC to
quickly predict other solvents or solvent mixtures and to generate a list of
suitable solvent systems.
The NRTL-SAC model calculates liquid activity coefficients. The equation for
the NRTL-SAC model is:
with
Option codes
The Flory-Huggins term is included by default in the NRTL-SAC model. You
can remove this term using the first option code. The table below lists the
values for this option code.
0 Flory-Huggins term included (default)
Others Flory-Huggins term removed
Reference
C.-C. Chen and Y. Song, "Solubility Modeling with a Nonrandom Two-Liquid
Segment Activity Coefficient Model," Ind. Eng. Chem. Res. 43, 8354 (2004).
Using NRTL-SAC
NRTL-SAC (patent pending) is a segment contribution activity coefficient
model, derived from the Polymer NRTL model, usable in Aspen Properties and
Aspen Polymers. NRTL-SAC can be used for fast, qualitative estimation of the
solubility of complex organic compounds in common solvents. For more
information about the model, see NRTL-SAC.
Where:
Tref = 298.15 K
The user must:
• Supply these elements for the binary parameters using a Properties |
Parameters | Binary | T-Dependent form.
• Supply these elements for Ψ on the Properties | Parameters | Electrolyte
Ternary form.
• Specify Comp ID i and Comp ID j (and Comp ID k for Ψ) on these forms,
using the same order that appears on the Components Specifications
Selection sheet.
The parameters are summarized in the following table. There is a Pitzer
parameter databank in the Aspen Physical Property System (see Physical
Property Data).
Parameter Provides No. of Default MDS Units
Name P1 - P5 for Elements
Cation-Anion Parameters
GMPTB0 β(0) 5 0 x —
GMPTB1 β(1) 5 0 x —
GMPTB2 β(2) 5 0 x —
GMPTB3 β(3) 5 0 x —
GMPTC ϕ 5 0 x —
C
Cation-Cation Parameters
GMPTTH θcc' 5 0 x —
Anion-Anion Parameters
GMPTTH θaa' 5 0 x —
Ternary Parameters
GMPTPS, Ψijk 1 (in each 0 x —
GMPTP1, parameter)
GMPTP2,
GMPTP3,
GMPTP4
Molecule-Ion and Molecule-Molecule Parameters
GMPTB0 β(0) 5 0 x —
GMPTB1 β(1) 5 0 x —
Model Development
The Pitzer model analyzes "hard-core" effects in the Debye-Hückel theory. It
uses the following expansion as a radial distribution function:
(1)
Where:
gij = Distribution function
r = Radius
qij =
k = Boltzmann's constant
T = Temperature
This radial distribution function is used in the so-called pressure equation that
relates this function and the intermolecular potential to thermodynamic
properties. From this relation you can obtain an expression for the osmotic
coefficient.
Pitzer proposes a general equation for the excess Gibbs energy. The basic
equation is:
(2)
Where:
GE = Excess Gibbs energy
R = Gas constant
T = Temperature
nw = Kilograms of water
mi =
(molality of ion i)
Subscripts c, c', and a, a' denote cations and anions of the solution. B, C, θ,
and Ψ are interaction parameters. f(I) is an electrostatic term as a function of
ionic strength. The cation-anion parameters B and C are characteristic for an
aqueous single-electrolyte system. These parameters can be determined by
the properties of pure (apparent) electrolytes. B is expressed as a function of
β(0) and β(1), or of β(0), β(2), and β(3) (see equations 11 through 15).
The parameters θ and Ψ are for the difference of interaction of unlike ions of
the same sign from the mean of like ions. These parameters can be measured
from common-ion mixtures. Examples are NaCl + KCl + H2O or NaCl + NaNO3
+ H2O (sic, Pitzer, 1989). These terms are discussed in detail later in this
section.
Fürst and Renon (1982) propose the following expression as the Pitzer
equation for the excess Gibbs energy:
(4)
The difference between equations 3 and 4 is that Pitzer orders cation before
anions. Fürst and Renon do not. All summations are taken over all ions i and j
(both cations and anions). This involves making the parameter matrices Bij,
Cij, θij, and Ψijk symmetric, as follows:
Second-order parameters are written Bij if i and j are ions of different sign. Bij
= 0 if the sign of zi = sign of zj, and Bii = 0. Since cations are not ordered
before anions, Bij = Bji. This eliminates the 2 in the second term in brackets in
Pitzer's original expression (equation 3). Second-order parameters are written
θij if i and j are ions of the same sign. Thus θij = 0 if the sign of zi is different
from the sign of zj, and θii = 0 with θij = θji.
Third-order parameters are written Cij if i and j are ions with different signs.
Cij = 0 if the sign of zi = sign of zj, and Cii = 0 with Cij = Cji. The factor of 2 in
the fifth bracketed term in Pitzer's original expression (equation 3) becomes
Ψijk is written for three different ions Ψijk = Ψkij = Ψjki , and Ψikk = 0. Ψijk = 0
if the sign of zi =sign of zj =sign of zk. The factor of 1/6 is different from 1/2
in the last term in brackets in Pitzer's original expression. Pitzer distinguishes
between cations and anions. In Pitzer's original model this parameter appears
twice, as Ψcc'a and Ψc'ca. In this modified model, it appears six times, as Ψcc'a;
Ψc'ca; Ψacc'; Ψac'c; Ψcac'; and Ψc'ac. Fürst and Renon's expression, equation 4,
calculates the expressions for activity coefficients and osmotic coefficients.
Pitzer (1975) modified his model by adding the electrostatic unsymmetrical
mixing effects, producing this modified Pitzer equation for the excess Gibbs
energy:
Where θ is neglected and Φij and Φ'ij are the electrostatic unsymmetric mixing
effects:
So that:
(10)
Where:
NA = Avogadro's constant
dw = Water density
with α1=2.0.
For n-m electrolytes, n and m>1 (2-2, 2-3, 3-4, and so on), B is determined
by the following expression:
(13)
The parameters β(0), β(1), β(2), β(3) and also C, θ, and Ψ can be found in
Pitzer's articles .
After the activity coefficients are calculated, they can be converted to the
mole fraction scale from the molality scale by the following relations:
For solutes:
(16)
Where:
Ψiii = 6μiii
However, molecule-molecule interactions were not taken into account by
Pitzer and coworkers. So μiii is an artificially introduced quantity.
Parameters
The Pitzer model in the Aspen Physical Property System involves user-
supplied parameters. These parameters are used in the calculation of binary
and ternary parameters for the electrolyte system. These parameters include
ϕ
the cation-anion parameters β(0), β(1), β(2), β(3) and C , cation-cation
parameter θcc', anion-anion parameter θaa', cation1-cation2-common anion
parameter Ψcc'a, anion1-anion2-common cation parameter Ψcaa', and the
ϕ
molecule-ion and molecule-molecule parameters β(0), β(1), and, C . The
parameter names in the Aspen Physical Property System and their
requirements are discussed in Pitzer Activity Coefficient Model.
Parameter Conversion
For n-m electrolytes, n and m>1 (2-2, 2-3, 3-4, and so on), the parameter
β(3) corresponds to Pitzer's β(1). β(2) is the same in both the Aspen Physical
Property System and original Pitzer models. Pitzer refers to the n-m
electrolyte parameters as β(1), β(2), β(0). β(0) and β(2) retain their meanings in
both models, but Pitzer's β(1) is β(3) in the Aspen Physical Property System. Be
careful to make this distinction when entering n-m electrolyte parameters.
ϕ
Pitzer often gives values of β(0), β(1), β(2), β(3), and C that are corrected by
some factors (see Pitzer and Mayorga (1973) for examples). These factors
originate from one of Pitzer's earlier expressions for the excess Gibbs energy:
(18)
Where:
=
Here β(0), β(1), β(2), and β(3) are multiplied by a factor of 2ncna. C is multiplied
by a factor of 2(ncna)3/2.
Aspen Physical Property System accounts for these correcting factors. Enter
the parameters without their correcting factors.
For example, Pitzer gives the values of parameters for MgCl2 as:
4/3β(0) = 0.4698
4/3β(1) = 2.242
= 1.6815
= 0.00520
Parameter Sources
Binary and ternary parameters for the Pitzer model for various electrolyte
systems are available from Pitzer's series on the thermodynamics of
electrolytes. These papers and the electrolyte parameters they give are:
Reference Parameters available
ϕ
(Pitzer, 1973) Binary parameters (β(0), β(1), C ) for 13
dilute aqueous electrolytes
(Pitzer and Mayorga, 1973) Binary parameters for 1-1 inorganic
electrolytes, salts of carboxylic acids (1-1),
tetraalkylammonium halids, sulfonic acids
and salts, additional 1-1 organic salts, 2-1
inorganic compounds, 2-1 organic
electrolytes, 3-1 electrolytes, 4-1 and 5-1
electrolytes
(Pitzer and Mayorga, 1974) Binary parameters for 2-2 electrolytes in
water at 25°C
(Pitzer and Kim, 1974) Binary and ternary parameters for mixed
electrolytes, binary mixtures without a
common ion, mixed electrolytes with three
or more solutes
(Pitzer, 1975) Ternary parameters for systems mixing
doubly and singly charged ions
(Pitzer and Silvester, 1976) Parameters for phosphoric acid and its buffer
solutions
(Pitzer, Roy and Silvester, 1977) Parameters and thermodynamic properties
for sulfuric acid
(Silvester and Pitzer, 1977) Data for NaCl and aqueous NaCl solutions
(Pitzer, Silvester, and Peterson, 1978) Rare earth chlorides, nitrates, and
perchlorates
(Peiper and Pitzer, 1982) Aqueous carbonate solutions, including
mixtures of sodium carbonate, bicarbonate,
and chloride
(Phutela and Pitzer, 1983) Aqueous calcium chloride
(Conceicao, de Lima, and Pitzer, 1983) Saturated aqueous solutions, including
mixtures of sodium chloride, potassium
chloride, and cesium chloride
Pitzer References
Conceicao, M., P. de Lima, and K.S. Pitzer, "Thermodynamics of Saturated
Aqueous Solutions Including Mixtures of NaCl, KCl, and CsCl, "J. Solution
Chem, Vol. 12, No. 3, (1983), pp. 171-185.
Fürst, W. and H. Renon, "Effects of the Various Parameters in the Application
of Pitzer's Model to Solid-Liquid Equilibrium. Preliminary Study for Strong 1-1
Electrolytes," Ind. Eng. Chem. Process Des. Dev., Vol. 21, No. 3, (1982),
pp. 396-400.
Guggenheim, E.A., Phil. Mag., Vol. 7, No. 19, (1935), p. 588.
Guggenheim, E.A. and J.C. Turgeon, Trans. Faraday Soc., Vol. 51, (1955), p.
747.
Kim, H. and W.J. Frederick, "Evaluation of Pitzer Ion Interaction Parameters
of Aqueous Mixed Electrolyte Solution at 25C, Part 2: Ternary Mixing
Parameters," J. Chem. Eng. Data, 33, (1988), pp. 278-283.
Pabalan, R.T. and K.S. Pitzer, "Thermodynamics of Concentrated Electrolyte
Mixtures and the Prediction of Mineral Solubilities to High Temperatures for
Mixtures in the system Na-K-Mg-Cl-SO4-OH-H2O," Geochimica Acta, 51,
(1987), pp. 2429-2443.
Peiper, J.C. and K.S. Pitzer, "Thermodynamics of Aqueous Carbonate
Solutions Including Mixtures of Sodium Carbonate, Bicarbonate, and
Chloride," J. Chem. Thermodynamics, Vol. 14, (1982), pp. 613-638.
Phutela, R.C. and K.S. Pitzer, "Thermodynamics of Aqueous Calcium
Chloride," J. Solution Chem., Vol. 12, No. 3, (1983), pp. 201-207.
Pitzer, K.S., "Thermodynamics of Electrolytes. I. Theoretical Basis and
General Equations, " J. Phys. Chem., Vol. 77, No. 2, (1973), pp. 268-277.
Pitzer, K.S., J. Solution Chem., Vol. 4, (1975), p. 249.
Pitzer, K.S., "Fluids, Both Ionic and Non-Ionic, over Wide Ranges of
Temperature and Composition," J. Chen. Thermodynamics, Vol. 21, (1989),
pp. 1-17. (Seventh Rossini lecture of the commission on Thermodynamics of
the IUPAC, Aug. 29, 1988, Prague, ex-Czechoslovakia).
Pitzer, K.S. and J.J. Kim, "Thermodynamics of Electrolytes IV; Activity and
Osmotic Coefficients for Mixed Electrolytes," J.Am. Chem. Soc., Vol. 96
(1974), p. 5701.
Pitzer, K.S. and G. Mayorga, "Thermodynamics of Electrolytes II; Activity and
Osmotic Coefficients for Strong Electrolytes with One or Both Ions Univalent,"
J. Phys. Chem., Vol. 77, No. 19, (1973), pp. 2300-2308.
Pitzer, K.S. and G. Mayorga, J. Solution Chem., Vol. 3, (1974), p. 539.
Where:
Ai =
Bi =
Ci =
Di =
Ei =
For any component i, the value of the activity coefficient can be fixed:
γi = fi
This model is not part of any property method. To use it:
1. On the Properties | Specifications sheet, specify an activity coefficient
model, such as NRTL.
2. Click the Properties | Property Methods folder.
3. In the Object Manager, click New.
4. In the Create New ID dialog box, enter a name for the new method.
5. In the Base Property Method field, select NRTL.
6. Click the Models tab.
7. Change the Model Name for GAMMA from GMRENON to GMPOLY.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
GMPLYP/1 ai1 0 x — — —
GMPLYP/2 ai2 0 x — — —
GMPLYP/3 ai3 0 x — — —
Redlich-Kister
This model calculates activity coefficients. It is a polynomial in the difference
between mole fractions in the mixture. It can be used for liquid and solid
mixtures (mixed crystals).
The equation is:
Where:
nc = Number of components
A1,ij = aij / T + bij
A2,ij = cij / T + dij
A3,ij = eij / T + fij
γi = vi
Parameter Name/ Symbol Default MDS Lower Upper Units
Element Limit Limit
GMRKTB/1 aij 0 x — — —
GMRKTB/2 bij 0 x — — —
GMRKTB/3 cij 0 x — — —
GMRKTB/4 dij 0 x — — —
GMRKTB/5 eij 0 x — — —
GMRKTB/6 fij 0 x — — —
GMRKTB/7 gij 0 x — — —
GMRKTB/8 hij 0 x — — —
GMRKTB/9 mij 0 x — — —
GMRKTB/10 nij 0 x — — —
GMRKTO vi — x — — —
Scatchard-Hildebrand
The Scatchard-Hildebrand model calculates liquid activity coefficients. It is
used in the CHAO-SEA property method and the GRAYSON property method.
The equation for the Scatchard-Hildebrand model is:
Where:
Aij =
ϕi =
Vm*,l =
Three-Suffix Margules
This model can be used to describe the excess properties of liquid and solid
solutions. It does not find much use in chemical engineering applications, but
is still widely used in metallurgical applications. Note that the binary
parameters for this model do not have physical significance.
The equation is:
For any component i, the value of the activity coefficient can be fixed:
γi = di
Parameter Name/ Symbol Default MDS Lower Upper Units
Element Limit Limit
GMMRGB/1 aij 0 x — — TEMPERATURE
GMMRGB/2 bij 0 x — — —
GMMRGB/3 cij 0 x — — —
GMMRGO di — x — — —
References
M. Margules, "Über die Zusammensetzung der gesättigten Dämpfe von
Mischungen," Sitzungsber. Akad. Wiss. Vienna, Vol. 104, (1895), p. 1293.
D.A. Gaskell, Introduction to Metallurgical Thermodyanics, 2nd ed., (New
York: Hemisphere Publishing Corp., 1981), p. 360.
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987).
ln γ = ln γic + ln γir
ln γic =
and
Where nc is the number of components in the mixture. The coordination
number z is set to 10. The parameters ri and qi are calculated from the group
volume and area parameters:
and
With:
GMUFQ Qk — — — — —
GMUFR Rk — — — — —
GMUFB bkn — — — — TEMPERATURE
The parameter UFGRP stores the UNIFAC functional group number and
number of occurrences of each group. UFGRP is stored in the Aspen Physical
Property System pure component databank for most components. For
nondatabank components, enter UFGRP on the Properties Molecular Structure
Functional Group sheet. See Physical Property Data, Chapter 3, for a list of
the UNIFAC functional groups.
UNIFAC-PSRK
The PSRK property method uses GMUFPSRK, the UNIFAC-PSRK model, which
is a variation on the standard UNIFAC model. UNIFAC-PSRK has special
groups defined for the light gases CO2, H2, NH3, N2, O2, CO, H2S, and argon,
and the group binary interaction parameters are temperature-dependent,
using the values in parameter UNIFPS, instead of the constant value from
GMUFB used above, so that:
References
Aa. Fredenslund, J. Gmehling and P. Rasmussen, "Vapor-Liquid Equilibria
using UNIFAC," (Amsterdam: Elsevier, 1977).
Aa. Fredenslund, R.L. Jones and J.M. Prausnitz, AIChE J., Vol. 21, (1975), p.
1086.
With:
GMUFDQ Qk — — — — —
GMUFDR Rk — — — — —
UNIFDM/1 amn,1 0 — — — TEMPERATURE
UNIFDM/2 amn,2 0 — — — TEMPERATURE
UNIFDM/3 amn,3 0 — — — TEMPERATURE
The parameter UFGRPD stores the group number and the number of
occurrences of each group. UFGRPD is stored in the Aspen Physical Property
System pure component databank. For nondatabank components, enter
UFGRPD on the Properties Molecular Structure Functional Group sheet. See
Physical Property Data, Chapter 3, for a list of the Dortmund modified UNIFAC
functional groups.
References
U. Weidlich and J. Gmehling, "A Modified UNIFAC Model 1. Prediction of VLE,
hE and ," Ind. Eng. Chem. Res., Vol. 26, (1987), pp. 1372–1381.
J. Gmehling, J. Li, and M. Schiller, "A Modified UNIFAC Model. 2. Present
Parameter Matrix and Results for Different Thermodynamic Properties," Ind.
Eng. Chem. Res., Vol. 32, (1993), pp. 178–193.
With:
Where Γk and Γki have the same meaning as in the original UNIFAC, but
defined as:
With:
GMUFLQ Qk — — — — —
The parameter UFGRPL stores the modified UNIFAC functional group number
and the number of occurrences of each group. UFGRPL is stored in the Aspen
Physical Property System pure component databank. For nondatabank
components, enter UFGRP on the Properties | Molecular Structure |
Functional Group sheet. See Physical Property Data, Chapter 3, for a list of
the Larsen modified UNIFAC functional groups.
Reference: B. Larsen, P. Rasmussen, and Aa. Fredenslund, "A Modified
UNIFAC Group-Contribution Model for Prediction of Phase Equilibria and Heats
of Mixing," Ind. Eng. Chem. Res., Vol. 26, (1987), pp. 2274 – 2286.
Where:
θi =
θi ' =
Φi =
li =
ti ' =
τij =
z = 10
References
D.S. Abrams and J.M. Prausnitz, "Statistical Thermodynamics of liquid
mixtures: A new expression for the Excess Gibbs Energy of Partly or
Completely Miscible Systems," AIChE J., Vol. 21, (1975), p. 116.
A. Bondi, "Physical Properties of Molecular Crystals, Liquids and Gases," (New
York: Wiley, 1960).
Simonetty, Yee and Tassios, "Prediction and Correlation of LLE," Ind. Eng.
Chem. Process Des. Dev., Vol. 21, (1982), p. 174.
Where:
Ai =
Bi =
Ci =
Aij =
Cij =
Cij = Cji
Aii = Bii = Cii = 0
aij and bij are unsymmetrical. That is, aij may not be equal to aji, and bij may
not be equal to bji.
Parameters Symbol Default MDS Lower Limit Upper Units
Name/Element Limit
VANL/1 aij 0 x -50.0 50.0 —
VANL/2 bij 0 x -15000.0 15000.0 TEMPERATURE
VANL/3 cij 0 x -50.0 50.0 —
VANL/4 dij 0 x -15000.0 15000.0 TEMPERATURE
The VANL-2 property method uses data set 2 for VANL. All other Van Laar
methods use data set 1.
References
J.J. Van Laar, "The Vapor Pressure of Binary Mixtures," Z. Phys. Chem., Vol.
72, (1910), p. 723.
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed. (New York: McGraw-Hill, 1987).
Where:
For any component i, the value of the activity coefficient can be fixed:
γi = gi
This model is recommended for dilute solutions.
Parameter Name/ Symbol Default MDS Lower Upper Units
Element Limit Limit
GMWIPR/1 ai 0 x — — TEMPERATURE
GMWIPR/2 bi 0 x — — —
GMWIPR/3 ci 0 x — — —
GMWIPB/1 dij 0 x — — TEMPERATURE
GMWIPB/2 eij 0 x — — —
GMWIPB/3 fij 0 x — — —
GMWIPO gi — x — — —
GMWIPS — 0 x — — —
GMWIPS is used to identify the solvent component. You must set GMWIPS to
1.0 for the solvent component. This model allows only one solvent.
References
A.D. Pelton and C. W. Bale, "A Modified Interaction Parameter Formalism for
Non-Dilute Solutions," Metallurgical Transactions A, Vol. 17A, (July 1986),
p. 1211.
Where:
The WILS-2 property method uses data set 2 for WILSON. All other Wilson
methods use data set 1.
† In the original formulation of the Wilson model, aij = ln Vj/Vi, cij = dij = eij =
0, and
bij = -(λij - λii)/R, where Vj and Vi are pure component liquid molar volume at
25°C.
†† If any of biA, ciA, and eiA are non-zero, absolute temperature units are
assumed for all coefficients. If biA, ciA, and eiA are all zero, the others are
interpreted in input units. The temperature limits are always interpreted in
input units.
References
G.M. Wilson, J. Am. Chem. Soc., Vol. 86, (1964), p. 127.
Where:
†
ln Aij =
Vj and Vi are pure component liquid molar volume at the system temperature
calculated using the Rackett/DIPPR/IK-CAPE model. The extended form of ln
Aij provides more flexibility in fitting phase equilibrium and enthalpy data. aij,
bij, cij, dij, and eij are unsymmetrical. That is, aij may not be equal to aji, etc.
The binary parameters aij, bij, cij, dij, and eij must be determined from data
regression of VLE and/or heat-of-mixing data. There are no built in binary
parameters for this model.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
WSNVOL/1 aij 0 x -50.0 50.0 —
WSNVOL/2 bij 0 x -15000.0 15000.0 TEMPERATURE
††
WSNVOL/3 cij 0 x -— — TEMPERATURE
††
WSNVOL/4 dij 0 x — — TEMPERATURE
††
WSNVOL/5 eij 0 x — — TEMPERATURE
††
WSNVOL/6 Tlower 0K x — — TEMPERATURE
WSNVOL/7 Tupper 1000 K x — — TEMPERATURE
Extended Antoine/Wagner/PPDS/IK-CAPE
Liquid Vapor Pressure Model
The Aspen Physical Property System has several submodels for calculating
vapor pressure of a liquid. It uses parameter THRSWT/3 to determine which
submodel is used. See Pure Component Temperature-Dependent Properties
for details.
If THRSWT/3 is Then this equation is used
0 Extended Antoine
200 BARIN
301 Wagner
302 PPDS Modified Wagner
400 PML
401 IK-CAPE
501 NIST TDE Polynomial
502 NIST Wagner 25
If C5i, C6i, or C7i is non-zero, absolute temperature units are assumed for all
coefficients C1i through C7i. The temperature limits are always in user input
units.
Where:
Tri = T / Tci
pri*,l = pi*,l / pci
Linear extrapolation of ln pi*,l versus T occurs outside of temperature bounds.
Where:
Tri = T / Tci
pri*,l = pi*,l / pci
Linear extrapolation of ln pi*,l versus T occurs outside of temperature bounds.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
WAGNER/1 C1i x
WAGNER/2 C2i 0 x
WAGNER/3 C3i 0 x
WAGNER/4 C4i 0 x
TC Tci TEMPERATURE
PC pci PRESSURE
Where:
Tri = T / Tci
Linear extrapolation of ln pi*,l versus T occurs outside of temperature bounds.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
WAGNER25/1 C1i x
WAGNER25/2 C2i 0 x
WAGNER25/3 C3i 0 x
WAGNER25/4 C4i 0 x
WAGNER25/5 ln pci 0 x
Where:
Tri = T / Tci
pri*,l = pi*,l / pci
LNVPEQ/2 and LNVPEQ/3 are the lower and upper temperature limits,
respectively.
In equations 1-4, if elements 4, 7, or 8 are non-zero, absolute temperature
units are assumed for all elements.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
LNVP1/1, ...,8 C1i, ..., C8i x PRESSURE
TEMPERATURE
LOGVP1/1, ..., 8 C1i, ..., C8i x PRESSURE
TEMPERATURE
LNPR1/1, ..., 8 C1i, ..., C8i x PRESSURE
TEMPERATURE
LOGPR1/1, ..., 8 C1i, ..., C8i x PRESSURE
TEMPERATURE
LNPR2/1,2 C1i, C2i x
TC Tci TEMPERATURE
PC pci PRESSURE
PLTDEPOL/3 C3i 0 X
References
Reid, Prausnitz, and Poling, The Properties of Gases and Liquids, 4th ed.,
(New York: McGraw-Hill, 1987).
Harlacher and Braun, "A Four-Parameter Extension of the Theorem of
Corresponding States," Ind. Eng. Chem. Process Des. Dev., Vol. 9, (1970), p.
479.
W. Wagner, Cryogenics, Vol. 13, (1973), pp. 470-482.
D. Ambrose, M. B. Ewing, N. B. Ghiassee, and J. C. Sanchez Ochoa "The
ebulliometric method of vapor-pressure measurement: vapor pressures of
benzene, hexafluorobenzene, and naphthalene," J. Chem. Thermodyn. 22
(1990), p. 589.
Detail of the model is described in the reference below and is too involved to
be reproduced here.
Reference
New Correlation of NH3, CO2, and H2S Volatility Data from Aqueous Sour
Water Systems, API Publication 955, March 1978 (American Petroleum
Institute).
References
B.C. Cajander, H.G. Hipkin, and J.M. Lenoir, "Prediction of Equilibrium Ratios
from Nomograms of Improved Accuracy," Journal of Chemical Engineering
Data, vol. 5, No. 3, July 1960, p. 251-259.
J.M. Lenoir, "Predict K Values at Low Temperatures, part 1," Hydrocarbon
Processing, p. 167, September 1969.
J.M. Lenoir, "Predict K Values at Low Temperatures, part 2," Hydrocarbon
Processing, p. 121, October 1969.
Where:
=
References
K.C. Chao and J.D. Seader, "A General Correlation of Vapor-Liquid Equilibria
in Hydrocarbon Mixtures," AIChE J., Vol. 7, (1961), p. 598.
Where:
References
H.G. Grayson and C.W. Streed, Paper 20-PO7, Sixth World Petroleum
Conference, Frankfurt, June 1963.
Where:
R' = Alcohol substituted alkyl groups
The equilibrium constants are given by:
The apparent fugacities and partial molar enthalpies, Gibbs energies and
entropies of H2S and CO2 are calculated by standard thermodynamic
relationships. The chemical reactions are always considered.
The values of the coefficients for the seven equilibrium constants (A1i, ... A5i)
and for the two Henry's constants B1i and B2i are built into the Aspen Physical
Property System. The coefficients for the equilibrium constants were
determined by regression. All available data for the four amines were used:
monoethanolamine, diethanolamine, disopropanolamine and diglycolamine.
You are not required to enter any parameters for this model.
References
API procedure 5A1.15 and 5A1.13.
DIPPR/Watson/PPDS/IK-CAPE
Heat of Vaporization Model
The Aspen Physical Property System has several submodels for calculating
pure component heat of vaporization. It uses parameter THRSWT/4 to
Where:
Tri = T / Tci
.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
DHVLDP/1 C1i — x — — MOLE-ENTHALPY
DHVLDP/2, ..., 5 C2i, ..., C5i 0 x — — —
DHVLDP/6 C6i 0 x — — TEMPERATURE
DHVLDP/7 C7i 1000 x — — TEMPERATURE
TC Tci — — 5.0 2000.0 TEMPERATURE
Clausius-Clapeyron Equation
The Aspen Physical Property System can calculate heat of vaporization using
the Clausius Clapeyron equation:
Where:
Where:
fcn = A correlation based on API Figure 6A3.5 (API Technical Data Book,
Petroleum Refining, 4th edition)
At high density, the Ritter equation is used (adapted from Ritter, Lenoir, and
Schweppe, Petrol. Refiner 37 [11] 225 (1958)):
Brelvi-O'Connell
The Brelvi-O'Connell model calculates partial molar volume of a supercritical
component i at infinite dilution in pure solvent A. Partial molar volume at
infinite dilution is required to compute the effect of pressure on Henry's
constant. (See Henry's Constant.)
The general form of the Brelvi-O'Connell model is:
Where:
References
S.W. Brelvi and J.P. O'Connell, AIChE J., Vol. 18, (1972), p. 1239.
S.W. Brelvi and J.P. O'Connell, AIChE J., Vol. 21, (1975), p. 157.
Where:
Vml = Liquid molar volume for electrolyte solutions.
Where:
xw = Mole fraction of water.
For electrolytes:
(3)
(4)
(5)
Where:
xca = Apparent mole fraction of electrolyte ca
Liquid molar volume for electrolyte ca
The mole fractions xca are reconstituted arbitrarily from the true ionic
concentrations, even if you use the apparent component approach. This
technique is explained in Electrolyte Calculation in Physical Property Methods.
The result is that electrolytes are generated from all possible combinations of
ions in solution. The following equation is consistently applied to determine
the amounts of each possible apparent electrolyte nca:
(6)
Where:
nca = Number of moles of apparent electrolyte ca
zc = Charge of c
zfactor = zc if c and a have the same number of
charges; otherwise 1.
nc = Number of moles of cation c
na = Number of moles of anion a
For example: given an aqueous solution of Ca2+, Na+, SO42-, Cl- four
electrolytes are found: CaCl2, Na2SO4, CaSO4, and NaCl. The Clarke
parameters of all four electrolytes are used. You can rely on the default,
which calculates the Clarke parameters from ionic parameters. Otherwise, you
Where:
Vml,t = Liquid volume per number of true species
Vml,a = Liquid volume per number of apparent species, Vml
of equation 1
na = Number of apparent species
t
n = Number of true species
The apparent molar volume is calculated as explained in the preceding
subsection.
Temperature Dependence
The temperature dependence of the molar volume of the solution is
approximately equal to the temperature dependence of the molar volume of
the solvent mixture:
(8)
Reference
C.C. Chen, private communication.
Where B and C are functions of T, ω, Tc, Pc and Psat is the saturated pressure
at T.
Where:
References
R.W. Hankinson and G.H. Thomson, AIChE J., Vol. 25, (1979), p. 653.
G.H. Thomson, K.R. Brobst, and R.W. Hankinson, AIChE J., Vol. 28, (1982),
p. 4, p. 671.
Where:
∞
Vk is the molar volume for water and is calculated from the ASME steam
table.
Vk is calculated from the Debye-Hückel limiting law for ionic species. It is
assumed to be the infinite dilution partial volume for molecular solutes.
Where:
∞ = Partial molar ionic volume at infinite dilution
Vk
zk = Charge number of ion k
Av = Debye-Hückel constant for volume
b = 1.2
I =
the ionic strength, with
mk = Molarity of ion k
Av is computed as follows:
Where:
Rackett/DIPPR/PPDS/IK-CAPE Liquid
Molar Volume
The Aspen Physical Property System has several submodels for calculating
liquid molar volume. It uses parameter THRSWT/2 to determine which
submodel is used. See Pure Component Temperature-Dependent Properties
for details.
If THRSWT/2 is This equation is used And this parameter is used
0 Rackett RKTZRA
100-116 DIPPR DNLDIP
301 PPDS DNLPDS
401 IK-CAPE VLPO
503 NIST ThermoML DNLTMLPO
Polynomial
504 NIST TDE expansion DNLEXSAT
514 NIST TDE Rackett DNLRACK
515 NIST COSTALD DNLCOSTD
For liquid molar volume of mixtures, the Rackett mixture equation is always
used by default. This is not necessarily consistent with the pure component
molar volume or density. If you need consistency, select route VLMX26 on the
Properties | Property Methods form. This route calculates mixture molar
volume from the mole-fraction average of pure component molar volumes.
Many of these equations calculate density first, and return calculate liquid
molar volume based on that density:
τ = 1 - T / Tc
Other DIPPR equations, such as equation 100, may be used for some
substances. Check the value of THRSWT/2 to determine the equation used.
See Pure Component Temperature-Dependent Properties for details about
DIPPR equations.
† For equation 116, the units are MOLE-DENSITY for parameters DNLDIP/1
through DNLDIP/5 and the default for DNLDIP/3 is 0. For equation 105,
DNLDIP/5 is not used, and absolute temperature units are assumed for
DNLDIP/3.
PPDS
The PPDS equation is:
IK-CAPE
The IK-CAPE equation is:
Where:
Tr = T / Tci
References
H.G. Rackett, J.Chem. Eng. Data., Vol. 15, (1970), p. 514.
C.F. Spencer and R.P. Danner, J. Chem. Eng. Data, Vol. 17, (1972), p. 236.
Rackett
The equation for the Rackett model is:
Where:
Tc =
ZmRA =
Vcm =
Tr = T / Tc
References
H.G. Rackett, J.Chem, Eng. Data., Vol. 15, (1970), p. 514.
C.F. Spencer and R.P. Danner, J. Chem. Eng. Data, Vol. 17, (1972), p. 236.
Where:
Tc =
kij =
ZmRA =
Zi*,RA =
Vcm =
Tr = T / Tc
MRKKIJ/2 Bij 0 x — — —
MRKKIJ/3 Cij 0 x — — —
References
H.G. Rackett, J.Chem, Eng. Data., Vol. 15, (1970), p. 514.
C.F. Spencer and R.P. Danner, J. Chem. Eng. Data, Vol. 17, (1972), p. 236.
IK-CAPE Equation
The IK-CAPE equation is:
DIPPR
The DIPPR equation is:
Option Codes
Option Code Value Descriptions
1 0 Use normal pure component liquid volume model for all
components (default)
1 Use steam tables for water
2 0 Use mole basis composition (default)
1 Use mass basis composition
Parameter
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
VLQKIJ Kij - x - - —
versus T is linearly extrapolated using the slope at C7i for T < C7i
versus T is linearly extrapolated using the slope at C8i for T < C8i
Parameter Applicable Symbol Default Units
Name/Element Components
CPAQ0/1 Ions, molecular solutes C1i — TEMPERATURE and
HEAT CAPACITY
CPAQ0/2,…,6 Ions, molecular solutes C2i, ..., C6i 0 TEMPERATURE and
HEAT CAPACITY
CPAQ0/7 Ions, molecular solutes C7i 0 TEMPERATURE
CPAQ0/8 Ions, molecular solutes C8i 1000 TEMPERATURE
Cations MOLE-ENTROPY
= Enthalpy of vaporization
To specify that the model is used for a component, enter a value of 1.0 for
COMPHL.
(Other DIPPR equations may sometimes be used. See Pure Component
Temperature-Dependent Properties for details.)
† If C5i is non-zero, absolute temperature units are assumed for C2i through
C5i. Otherwise, user input units for temperature are used. The temperature
limits are always interpreted in user input units.
When THRSWT/6 is 401, the IK-CAPE heat capacity polynomial is used.The
equation is:
When THRSWT/6 is 506, the TDE equation is used to calculate liquid heat
capacity:
Aspen/DIPPR/Barin/PPDS/IK-CAPE Ideal
Gas Heat Capacity
The Aspen Physical Property System has several submodels for calculating
ideal gas heat capacity. It uses parameter THRSWT/7 to determine which
submodel is used. See Pure Component Temperature-Dependent Properties
for details.
If THRSWT/7 is This equation is used And this parameter is
used
0 Ideal gas heat capacity CPIG
polynomial
107 DIPPR CPIGDP
200 Barin CPIXP1, CPIXP2, CPIXP3
301 PPDS CPIGDS
401 IK-CAPE heat capacity CPIGPO
polynomial
503 NIST ThermoML CPITMLPO
polynomial
513 NIST Aly-Lee CPIALEE
These equations are also used to calculate ideal gas enthalpies, entropies,
and Gibbs energies.
† If C10i or C11i is non-zero, then absolute temperature units are assumed for
C9i through C11i. Otherwise, user input temperature units are used for all
parameters. User input temperature units are always used for C1i through C8i.
DIPPR
The DIPPR ideal gas heat capacity equation by Aly and Lee 1981 is:
†† Absolute temperature units are assumed for C3i and C5i. The temperature
limits are always interpreted in user input units.
NIST Aly-Lee
This equation is the same as the DIPPR Aly and Lee equation above, but it
uses a different parameter set. Note that elements 6 and 7 of the CPIALEE
parameter are not used in the equation.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
CPIALEE/1 C1i — x — — MOLE-HEAT-
CAPACITY
CPIALEE/2 C2i 0 x — — MOLE-HEAT-
CAPACITY
CPIALEE/3 C3i 0 x — — TEMPERATURE ††
CPIALEE/4 C4i 0 x — — MOLE-HEAT-
CAPACITY
CPIALEE/5 C5i 0 x — — TEMPERATURE ††
CPIALEE/8 C6i 0 x — — TEMPERATURE
CPIALEE/9 C7i 1000 x — — TEMPERATURE
†† Absolute temperature units are assumed for C3i and C5i. The temperature
limits are always interpreted in user input units.
PPDS
The PPDS equation is:
References
Data for the Ideal Gas Heat Capacity Polynomial: Reid, Prausnitz and Poling,
The Properties of Gases and Liquids, 4th ed., (New York: McGraw-Hill, 1987).
The Aspen Physical Property System combustion data bank, JANAF
Thermochemical Data, Compiled and calculated by the Thermal Research
Laboratory of Dow Chemical Company.
F. A. Aly and L. L. Lee, "Self-Consistent Equations for Calculating the Ideal
Gas Heat Capacity, Enthalpy, and Entropy, Fluid Phase Eq., Vol. 6, (1981), p.
169.
The enthalpy, entropy, and Gibbs energy of solids are also calculated from
these equations:
† The units are TEMPERATURE and HEAT-CAPACITY for all elements. If any of
C4i through C6i are non-zero, absolute temperature units are assumed for
elements C1i through C6i. Otherwise, user input temperature units are
assumed for all elements. The temperature limits are always interpreted in
user input units.
DIPPR
The DIPPR equation is:
Solubility Correlations
The Aspen Physical Property System has three built-in solubility correlation
models. This section describes the solubility correlation models available.
Model Type
Henry's constant Gas solubility in liquid
Water solubility Water solubility in organic liquid
Hydrocarbon solubility Hydrocarbon solubility in water-rich
liquid
Where:
wA =
ln HiA(T, pA*,l) =
HiA(T, P) =
† If any of biA, ciA, and eiA are non-zero, absolute temperature units are
assumed for all coefficients. If biA, ciA, and eiA are all zero, the others are
interpreted in input units. The temperature limits are always interpreted in
input units.
Water Solubility
This model calculates solubility of water in a hydrocarbon-rich liquid phase.
The model is used automatically when you model a hydrocarbon-water
system with the free-water option. See Free-Water Immiscibility Simplification
in Free-Water and Three-Phase Calculations for details.
The expression for the liquid mole fraction of water in the ith hydrocarbon
species is:
Hydrocarbon Solubility
This model calculates solubility of hydrocarbon in a water-rich liquid phase.
The model is used automatically when you model a hydrocarbon-water
system with the dirty-water option. See Free-Water Immiscibility
Simplification in Free-Water and Rigorous Three-Phase Calculations for
details.
The expression for the liquid mole fraction of the ith hydrocarbon species in
water is:
Reference
C. Tsonopoulos, Fluid Phase Equilibria, 186 (2001), 185-206.
Cavett
The general form for the Cavett model is:
• Enthalpy:
(2)
• Entropy:
(3)
• Heat capacity:
(4)
α refers to an arbitrary phase which can be solid, liquid, or ideal gas. For each
phase, multiple sets of parameters from 1 to n are present to cover multiple
temperature ranges. The value of the parameter n depends on the phase.
(See tables that follow.) When the temperature is outside all these
temperature ranges, linear extrapolation of properties versus T is performed
using the slope at the end of the nearest temperature bound.
The four properties Cp, H, S, and G are interrelated as a result of the
thermodynamic relationships:
Solid Phase
The parameters in range n are valid for temperature: Tn,ls < T < Tn,hs
When you specify this parameter, be sure to specify at least elements 1
through 3.
Parameter Name † Symbol Default MDS Lower Upper Units
/Element Limit Limit
CPSXPn/1 Tn,ls — x — — TEMPERATURE
CPSXPn/2 Tn,hs — x — — TEMPERATURE
CPSXPn/3 an,is — x — — ††
CPSXPn/4 bn,is 0 x — — ††
CPSXPn/5 cn,is 0 x — — ††
CPSXPn/6 dn,is 0 x — — ††
s
CPSXPn/7 en,i 0 x — — ††
CPSXPn/8 fn,is 0 x — — ††
CPSXPn/9 gn,is 0 x — — ††
CPSXPn/10 hn,is 0 x — — ††
Liquid Phase
The parameters in range n are valid for temperature: Tn,ll < T < Tn,hl
When you specify this parameter, be sure to specify at least elements 1
through 3.
Parameter Name † Symbol Default MDS Lower Upper Units
/Element Limit Limit
CPLXPn/1 Tn,ll — x — — TEMPERATURE
CPLXPn/2 Tn,hl — x — — TEMPERATURE
l
CPLXPn/3 an,i — x — — ††
CPLXPn/4 bn,il 0 x — — ††
CPLXPn/5 cn,il 0 x — — ††
CPLXPn/6 dn,il 0 x — — ††
l
CPLXPn/7 en,i 0 x — — ††
CPLXPn/8 fn,il 0 x — — ††
† n is 1 through 3. CPIXP1 vector stores ideal gas parameters for the first
temperature range. CPIXP2 vector stores ideal gas parameters for the second
temperature range, and so on.
†† TEMPERATURE, ENTHALPY, ENTROPY
The molar enthalpy Hm* and the molar excess enthalpy Hm*E are defined with
the asymmetrical reference state: the pure solvent water and infinite dilution
of molecular solutes and ions. (here * refers to the asymmetrical reference
state.)
Hw* is the pure water molar enthalpy, calculated from the Ideal Gas model
and the ASME Steam Table equation-of-state. (here * refers to pure
component.)
Hm*E is excess enthalpy and is calculated from the electrolyte NRTL activity
coefficient model.
See Criss-Cobble model and Henry's law model for more information.
Option codes can improve the performance of this model when using mixed
solvents and Henry components. See Option Codes for Enthalpy Models for
details.
Parameter Applicable Symbol Default Units
Name Components
IONTYP Ions † Ion 0 —
SO25C Cations ∞,aq — MOLE-ENTROPY
Sc (T=298)
Anions ∞,aq — MOLE-ENTROPY
Sa (T=298)
DHAQFM Ions, Molecular Solutes Δf Hk∞,aq — MOLE-ENTHALPY
††
CPAQ0 Ions, Molecular Solutes ∞,aq — HEAT-CAPACITY
Cp,k
††
DHFORM Molecular Solutes †† Δf Hi*,ig — MOLE-ENTHALPY
The molar Gibbs free energy and the molar excess Gibbs free energy Gm* and
Gm*E are defined with the asymmetrical reference state: as pure water and
infinite dilution of molecular solutes and ions. (* refers to the asymmetrical
reference state.) The ideal mixing term is calculated normally, where j refers
to any component. The molar Gibbs free energy of pure water (or
thermodynamic potential) μw* is calculated from the ideal gas contribution.
This is a function of the ideal gas heat capacity and the departure function.
(here * refers to the pure component.)
∞ ∞
Because Δf Hk ,aq and Δf Gk ,aq are based on a molality scale, and μk∞ is based
on mole fraction scale, the term RT ln(1000/Mw) is added.
† IONTYP and SO25C are not needed if CPAQ0 is given for ions.
†† DGFORM is not needed if DHAQFM and CPAQ0 are given for molecular
solutes.
††† The unit keywords for CPIG are TEMPERATURE and HEAT-CAPACITY. If
CPIG/10 or CPIG/11 is non-zero, then absolute temperature units are
assumed for CPIG/9 through CPIG/11. Otherwise, user input temperature
units are used for all elements of CPIG. User input temperature units are
always used for other elements of CPIG.
Where:
Hi*,ig = Ideal gas enthalpy
Hi*,v - Hi*,ig = Vapor enthalpy departure from equation of state
Where:
Hmig = Enthalpy of ideal gas mixture
=
HmE,l =
Where:
Hiref,l = Reference enthalpy for liquid state at Tiref,l
Cp,i*,l = Liquid heat capacity of component i
The saturated liquid Gibbs free energy is calculated as follows:
Where:
Giref,l = Reference Gibbs free energy for liquid state at Tiref,l
p = System pressure
ref
p = Reference pressure (=101325 N/m2)
For the WILS-LR property method, Hiref,l and Giref,l default to zero (0). The
reference temperature Tiref,l defaults to 273.15K.
For the WILS-GLR property method, the reference temperature Tiref,l defaults
to 298.15K and Hiref,l defaults to:
Where:
Hiref,ig = Ideal gas enthalpy of formation for liquid state at
298.15K
Giref,ig = Ideal gas Gibbs free energy of formation for liquid
state at 298.15K
Note that we cannot default the liquid reference enthalpy and Gibbs free
energy to zero, as is the case for WILS-LR, because it will cause inconsistency
with the enthalpy of components that use ideal-gas reference state. The
default values used result in the enthalpies of all components being on the
same basis. In fact, if you enter values for Hiref,l and Giref,l for a liquid-
reference state component you must make sure that they are consistent with
each other and are consistent with the enthalpy basis of the remaining
components in the mixture. If you enter a value for Hiref,l, you should also
enter a value for Giref,l to ensure consistency.
When the liquid-reference state is used, the ideal gas enthalpy at
temperature T is not calculated from the integration of the ideal gas heat
capacity equation (see Ideal Gas as Reference State section below). For
consistency, it is calculated from liquid enthalpy as follows:
Where:
Ticon,l = Temperature of conversion from liquid to vapor
enthalpy for component i
Ticon,l is the temperature at which one crosses from liquid state to the vapor
state. This is a user defined temperature that defaults to the system
temperature T. Ticon,l may be selected such that heat of vaporization for
component i at the temperature is most accurate.
The vapor enthalpy is calculated from ideal gas enthalpy as follows:
Where:
Where:
Hiref,ig = Reference state enthalpy for ideal gas at Tiref,ig
= Heat of formation of ideal gas at 298.15 K by default
Tiref,ig = Reference temperature corresponding to Hiref,ig.
Defaults to 298.15 K
Ticon,ig = The temperature at which one crosses from vapor
state to liquid state. This is a user defined
temperature that defaults to the system temperature
T. Ticon,igmay be selected such that heat of
vaporization of component i at the temperature is
most accurate.
The ideal gas enthalpy is calculated as follows:
The liquid heat capacity and the ideal gas heat capacity can be calculated
from the ASPEN, DIPPR, IK-CAPE, BARIN, NIST/TDE, or other available
models. The heat of vaporization can be calculated from the Watson, DIPPR,
IK-CAPE, NIST/TDE, or other available models. The enthalpy departure is
ω = Born coefficient
References
Tanger, J.C. IV and H.C. Helgeson, "Calculation of the thermodynamic and
transport properties of aqueous species at high pressures and temperatures:
Revised equation of state for the standard partial properties of ions and
electrolytes," American Journal of Science, Vol. 288, (1988), p. 19-98.
Shock, E.L. and H.C. Helgeson, "Calculation of the thermodynamic and
transport properties of aqueous species at high pressures and temperatures:
The pure component properties Qi and Qj are calculated by the default model
for that property, unless modified by option codes. Composition xi and xj is in
mole fraction unless modified by option codes. Kij is a binary parameter
specific to the mixing rule for each property. A variation on the equation using
the logarithm of the property is used for viscosity.
This section describes the transport property models available in the Aspen
Physical Property System. The following table provides an overview of the
available models. This table lists the Aspen Physical Property System model
names, and their possible use in different phase types, for pure components
and mixtures.
Viscosity Models
The Aspen Physical Property System has the following built-in viscosity
models:
Model Type
Andrade Liquid Mixture Viscosity Liquid
Andrade/DIPPR/PPDS/IK-CAPE Pure component liquid
API liquid viscosity Liquid
API 1997 liquid viscosity Liquid
Chapman-Enskog-Brokaw/DIPPR/IK-CAPE Low pressure vapor,
pure components
Chapman-Enskog-Brokaw-Wilke Mixing Rule Low pressure vapor,
mixture
Chung-Lee-Starling Low Pressure Low pressure vapor
Chung-Lee-Starling Liquid or vapor
Dean-Stiel Pressure correction Vapor
IAPS viscosity Water or steam
Jones-Dole Electrolyte Correction Electrolyte
Letsou-Stiel High temperature liquid
Lucas Vapor
TRAPP viscosity Vapor or liquid
Aspen Liquid Mixture Viscosity Liquid
ASTM Liquid Mixture Viscosity Liquid
Twu liquid viscosity Liquid
Viscosity quadratic mixing rule Liquid
Where:
kij =
mij =
Note that the Andrade liquid mixture viscosity model is called from other
models. The first option codes of these models cause fi to be mole or mass
fraction when Andrade is used in the respective models. To maintain
consistency across models, if you set the first option code for MUL2ANDR to
1, you should also the set the first option code of the other models to 1, if
they are used in your simulation.
Model Model Name
MUL2JONS Jones-Dole Electrolyte Viscosity model
DL0WCA Wilke-Chang Diffusivity model (binary)
DL1WCA Wilke-Chang Diffusivity model (mixture)
DL0NST Nernst-Hartley Electrolyte Diffusivity model (binary)
DL1NST Nernst-Hartley Electrolyte Diffusivity model (mixture)
The binary parameters kij and mij allow accurate representation of complex
liquid mixture viscosity. Both binary parameters default to zero.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
ANDKIJ/1 aij 0 — — — —
ANDKIJ/2 bij 0 — — — —
ANDMIJ/1 cij 0 — — — —
ANDMIJ/2 dij 0 — — — —
Andrade/DIPPR/PPDS/IK-CAPE Pure
Component Liquid Viscosity
The Aspen Physical Property System has several submodels for calculating
pure component liquid viscosity. It uses parameter TRNSWT/1 to determine
which submodel is used. See Pure Component Temperature-Dependent
Properties for details.
If TRNSWT/1 is This equation is used And this parameter is used
0 Andrade MULAND
101 DIPPR MULDIP
301 PPDS MULPDS
401 IK-CAPE polynomial equation MULPO
†† If any of C3i through C5i are non-zero, absolute temperature units are
assumed for C3i through C5i. Otherwise, all coefficients are interpreted in user
input temperature units. The temperature limits are always interpreted in
user input units.
Exponential
††† Absolute temperature units are assumed for C2i. The temperature limits
are always interpreted in user input units.
Polynomial
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 439.
Where:
fcn = A correlation based on API Procedures and Figures 11A4.1, 11A4.2,
and 11A4.3 (API Technical Data Book, Petroleum Refining, 4th
edition)
Vml is obtained from the API liquid volume model.
For real components, the Andrade/DIPPR model is used.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
TB Tbi — — 4.0 2000.0 TEMPERATURE
API APIi — — -60.0 500.0 —
Where:
fcn = A correlation based on API Procedures and Figures 11A4.2, 11A4.3,
and 11A4.4 (API Technical Data Book, Petroleum Refining, 1997
edition)
Vml is obtained from the API liquid volume model.
For real components, the Andrade/DIPPR model is used.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
TB Tbi — — 4.0 2000.0 TEMPERATURE
API APIi — — -60.0 500.0 —
Where:
Xi = Mole fraction or weight fraction of component i
kij = Symmetric binary parameter (kij = kji)
lij = Antisymmetric binary parameter (lij = -lji)
Where:
wi = Weight fraction of component i
Chapman-Enskog-
Brokaw/DIPPR/PPDS/IK-CAPE Vapor
Viscosity
The Aspen Physical Property System has several submodels for calculating
pure component low pressure vapor viscosity. It uses parameter TRNSWT/2
to determine which submodel is used. See Pure Component Temperature-
Dependent Properties for details.
If TRNSWT/2 is This equation is used And this parameter is
used
0 Chapman-Enskog- —
Brokaw
102 DIPPR MUVDIP
301 PPDS MUVPDS
302 PPDS kinetic theory MUVCEB
401 IK-CAPE polynomial MUVPO
equation
402 IK-CAPE Sutherland MUVSUT
equation
503 NIST ThermoML MUVTMLPO
polynomial
Chapman-Enskog-Brokaw
The equation for the Chapman-Enskog model is:
Where:
Ωη =
† If any of C2i through C4i are non-zero, absolute temperature units are
assumed for C1i through C4i. Otherwise, all coefficients are interpreted in user
input temperature units. The temperature limits are always interpreted in
user input units.
†† MUVDIP/5 is not used in this equation. It is normally set to zero. The
parameter is provided for consistency with other DIPPR equations.
PPDS
The PPDS submodel includes both the basic PPDS vapor viscosity equation
and the PPDS kinetic theory vapor viscosity equation. For either equation,
linear extrapolation of viscosity versus T occurs for temperatures outside
bounds.
† Absolute temperature units are assumed for C2i . The temperature limits are
always interpreted in user input units.
Sutherland Equation
†† If C2i is non-zero, absolute temperature units are assumed for C1i and C2i.
Otherwise, all coefficients are interpreted in user input temperature units. The
temperature limits are always interpreted in user input units.
Polynomial
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling. The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 392.
Chapman-Enskog-Brokaw-Wilke Mixing
Rule
The low pressure vapor mixture viscosity is calculated by the Wilke
approximation of the Chapman-Enskog equation:
Where:
Where p is in debye.
The pure component vapor viscosity ηi*,v (p = 0) can be calculated using the
Chapman- Enskog-Brokaw/DIPPR (or another) low pressure vapor viscosity
model.
References
R.C. Reid, J.M. Prausnitz, and T.K. Sherwood, The Properties of Gases and
Liquids, 3rd ed., (New York: McGraw-Hill, 1977), pp. 410–416.
C1 is a constant of correlation.
The polar parameter κ is tabulated for certain alcohols and carboxylic acids.
The previous equations can be used for mixtures when applying these mixing
rules:
Tcij =
ωij =
Mij =
κij =
Chung-Lee-Starling Viscosity
The Chung-Lee-Starling viscosity equation for vapor and liquid, high and low
pressure is:
With:
f1 =
f2 =
FC =
The polar parameter κ is tabulated for certain alcohols and carboxylic acids.
For low pressures, f1 is reduced to 1.0 and f2 becomes negligible. The
equation reduces to the low pressure vapor viscosity model by Chung-Lee and
Starling.
The previous equations can be used for mixtures when applying these mixing
rules:
Tcij =
ωij =
Mij =
κij =
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 427.
Tc =
M =
pc =
Vcm =
Zcm =
ηw=fcn(T, p)
Where:
fcn = Correlation developed by IAPS
The models are only applicable to water. There are no parameters required
for the models.
Where:
Jones-Dole
The Jones-Dole equation is:
(1)
Where:
= Concentration of apparent electrolyte ca (2)
La = (5)
Bca = (7)
Breslau-Miller
The Breslau-Miller equation is:
(8)
Carbonell
The Carbonell equation is:
(10)
Where:
Mk = Molecular weight of an apparent electrolyte
component k
You must provide parameters for the model used for the calculation of the
liquid mixture viscosity of the solvent mixture.
Parameter Symbol Default Lower Upper Units
Name/Element Limit Limit
CHARGE z 0.0 — — —
MW M — 1.0 5000.0 —
IONMOB/1 l1 — — — AREA, MOLES
IONMOB/2 l2 0.0 — — AREA, MOLES,
TEMPERATURE
IONMUB/1 b1 — — — MOLE-VOLUME
IONMUB/2 b2 0,0 — — MOLE-VOLUME,
TEMPERATURE
References
A. L. Horvath, Handbook of Aqueous Electrolyte Solutions, (Chichester: Ellis
Horwood, 1985).
(ηlε)0 =
(ηlε)1 =
ε =
ω =
References
R.C. Reid, J.M. Pransnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 471.
The polar and quantum correction factors at high and low pressure are:
FP =
FQ =
FPi (p = 0) =
pc =
M =
FP (p = 0) =
FQ (p = 0) =
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 421, 431.
Where:
The parameter is the mole fraction vector; fcn is a corresponding states
correlation based on the model for vapor and liquid viscosity TRAPP, by the
National Bureau of Standards (NBS, currently NIST) . The model can be used
for both pure components and mixtures. The model should be used for
nonpolar components only.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
MW Mi — — 1.0 5000.0 —
TCTRAP Tci TC x 5.0 2000.0 TEMPERATURE
5 8
PCTRAP pci PC x 10 10 PRESSURE
VCTRAP Vci VC x 0.001 3.5 MOLE-VOLUME
ZCTRAP Zci ZC x 0.1 1.0 —
OMGRAP ωi OMEGA x -0.5 3.0 —
References
J.F. Ely and H.J.M. Hanley, "Prediction of Transport Properties. 1. Viscosities
of Fluids and Mixtures," Ind. Eng. Chem. Fundam., Vol. 20, (1981), pp. 323–
332.
Option Codes
Option Code Value Descriptions
1 0 Use mole basis composition (default)
1 Use mass basis composition
Parameter
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
MLQKIJ Kij - x - - —
Where:
f1 =
f2 =
Ψ =
The polar parameter κ is tabulated for certain alcohols and carboxylic acids.
For low pressures, f1 is reduced to 1.0 and f2 is reduced to zero. This gives
the Chung-Lee-Starling expression for thermal conductivity of low pressure
gases.
The same expressions are used for mixtures. The mixture expression for η (p
= 0) must be used. (See Chung-Lee-Starling Low-Pressure Vapor Viscosity.)
Where:
Tcij =
ωij =
Mij =
κij =
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 505, 523.
λw=fcn(T, p)
Li Mixing Rule
Liquid mixture thermal conductivity is calculated using Li equation (Reid
et.al., 1987):
Where:
The pure component liquid molar volume Vi*,l is calculated from the Rackett
model.
Where:
Sato-Riedel/DIPPR/IK-CAPE Liquid
Thermal Conductivity
The Aspen Physical Property System has several submodels for calculating
pure component liquid thermal conductivity. It uses parameter TRNSWT/3 to
determine which submodel is used. See Pure Component Temperature-
Dependent Properties for details.
If TRNSWT/3 is This equation is used And this parameter is used
0 Sato-Riedel —
100 DIPPR KLDIP
301 PPDS KLPDS
401 IK-CAPE KLPO
503 NIST ThermoML KLTMLPO
polynomial
510 NIST PPDS8 equation KLPPDS8
Sato-Riedel
The Sato-Riedel equation is (Reid et al., 1987):
PPDS
The equation is:
IK-CAPE
The IK-CAPE equation is:
References
R.C. Reid, J.M. Prausnitz, and T.K. Sherwood, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1977), p. 533.
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 550.
For mixtures:
Stiel-Thodos
The Stiel-Thodos equation is:
Where:
† If any of C2i through C4i are non-zero, absolute temperature units are
assumed for C1i through C4i. Otherwise, all coefficients are interpreted in user
input temperature units. The temperature limits are always interpreted in
user input units.
PPDS
The equation is:
IK-CAPE Polynomial
References
R.C. Reid, J.M. Praunitz, and B.E. Poling, The Properties of Gases and Liquid,
4th ed., (New York: McGraw-Hill, 1987), p. 494.
Where:
ρrm =
Where:
wi = Liquid phase weight fraction of component i
For most systems, the ratio of maximum to minimum pure component liquid
thermal conductivity is between 1 and 2, where the exponent -2 is
recommended, and is the default value used.
References
J.F. Ely and H.J. M. Hanley, "Prediction of Transport Properties. 2. Thermal
Conductivity of Pure Fluids and Mixtures," Ind. Eng. Chem. Fundam., Vol. 22,
(1983), pp. 90–97.
Where:
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), pp. 530–531.
Diffusivity Models
The Aspen Physical Property System has seven built-in diffusivity models. This
section describes the diffusivity models available.
Model Type
Chapman-Enskog-Wilke-Lee (Binary) Low pressure vapor
Chapman-Enskog-Wilke-Lee (Mixture) Low pressure vapor
Dawson-Khoury-Kobayashi (Binary) Vapor
Dawson-Khoury-Kobayashi (Mixture) Vapor
Nernst-Hartley Electrolyte
Wilke-Chang (Binary) Liquid
Wilke-Chang (Mixture) Liquid
Chapman-Enskog-Wilke-Lee (Binary)
The binary diffusion coefficient at low pressures is calculated using
the Chapman-Enskog-Wilke-Lee model:
Where:
ΩD =
σij =
References
R.C. Reid, J.M. Praunitz, and B.E. Poling, The Properties of Gases and Liquids,
4th ed., (New York: McGraw-Hill, 1987), p. 587.
Chapman-Enskog-Wilke-Lee (Mixture)
The diffusion coefficient of a gas into a gas mixture at low pressures is
calculated using an equation of Bird, Stewart, and Lightfoot by default (option
code 0):
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 597.
Dawson-Khoury-Kobayashi (Binary)
The binary diffusion coefficient Dijv at high pressures is calculated from the
Dawson-Khoury-Kobayashi model:
The parameters ρmv and Vmv are obtained from the Redlich-Kwong equation-
of-state model.
You must supply parameters for these two models.
Subscript i denotes a diffusing component. j denotes a solvent.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
VC Vci — x 0.001 3.5 MOLE-VOLUME
References
R.C. Reid, J.M. Prausnitz, and T.K. Sherwood. The Properties of Gases and
Liquids, 3rd ed., (New York: McGraw-Hill, 1977), pp. 560-565.
Dawson-Khoury-Kobayashi (Mixture)
The diffusion coefficient of a gas into a gas mixture at high pressure is
calculated using an equation of Bird, Stewart, and Lightfoot by default (option
code 0):
The binary diffusion coefficient Dijv at high pressures is calculated from the
Dawson-Khoury-Kobayashi model. (See Dawson-Khoury-Kobayashi (Binary).)
At low pressures (up to 1 atm) the binary diffusion coefficient is instead
calculated by the Chapman-Enskog-Wilke-Lee (Binary) model.
You must provide parameters for this model.
Parameter Symbol Default MDS Lower Upper Units
Name/Element Limit Limit
DVBLNC — 1 — — — —
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 597.
Nernst-Hartley
The effective diffusivity of an ion i in a liquid mixture with electrolytes can be
calculated using the Nernst-Hartley model:
(1)
Where:
F = 9.65x107C/kmole (Faraday's number)
xk = Mole fraction of any molecular species k
zi = Charge number of species i
The binary diffusion coefficient of the ion with respect to a molecular species
is set equal to the effective diffusivity of the ion in the liquid mixture:
(2)
The binary diffusion coefficient of an ion i with respect to an ion j is set to the
mean of the effective diffusivities of the two ions:
References
A. L. Horvath, Handbook of Aqueous Electrolyte Solutions, (Chichester: Ellis
Horwood, Ltd, 1985).
Wilke-Chang (Binary)
The Wilke-Chang model calculates the liquid diffusion coefficient of
component i in a mixture at finite concentrations:
Dijl = Djil
The equation for the Wilke-Chang model at infinite dilution is:
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 598–600.
Wilke-Chang (Mixture)
The Wilke-Chang model calculates the liquid diffusion coefficient of
component i in a mixture.
The equation for the Wilke-Chang model is:
With:
Where:
References
R.C. Reid, J.M. Praunsnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th ed., (New York: McGraw-Hill, 1987), p. 618.
Where:
x = Mole fraction
r = Exponent (specified by the option code)
Hadden (S. T. Hadden, Hydrocarbon Process Petrol Refiner, 45(10), 1966, p.
161) suggested that the exponent value r =1 should be used for most
hydrocarbon mixtures. However, Reid recommended the value of r in the
range of -1 to -3. The exponent value r can be specified using the model’s
Option Code (option code = 1, -1, -2, ..., -9 corresponding to the value of r).
The default value of r for this model is 1.
The pure component liquid surface tension σi*,l can be calculated by these
models:
• Hakim-Steinberg-Stiel
• DIPPR
Where:
fcn = A correlation based on API Procedure 10A32 (API Technical Data
Book, Petroleum Refining, 4th edition)
σw=fcn(T, p)
Where:
fcn = Correlation developed by IAPS
The model is only applicable to water. No parameters are required.
Hakim-Steinberg-Stiel/DIPPR/PPDS/IK-
CAPE Liquid Surface Tension
The Aspen Physical Property System has several submodels for calculating
liquid surface tension. It uses parameter TRNSWT/5 to determine which
submodel is used. See Pure Component Temperature-Dependent Properties
for details.
If TRNSWT/5 is This equation is used And this parameter is
used
0 Hakim-Steinberg-Stiel —
Hakim-Steinberg-Stiel
The Hakim-Steinberg-Stiel equation is:
Where:
Qpi =
mi =
CHI χi 0 — — —
Where:
Tri = T / Tci
PPDS
The equation is:
IK-CAPE Polynomial
The IK-CAPE equation is:
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 4th. ed., (New York: McGraw-Hill, 1987), p. 638.
Where:
Where:
ccaa =
References
A. L. Horvath, Handbook of Aqueous Electrolyte Solutions, (Chichester: Ellis,
Ltd. 1985).
Where:
References
R.C. Reid, J.M. Prausnitz, and B.E. Poling, The Properties of Gases and
Liquids, 3rd. ed., (New York: McGraw-Hill, 1977).
This section describes the nonconventional solid density and enthalpy models
available in the Aspen Physical Property System. The following table lists the
available models and their model names. Nonconventional components are
solid components that cannot be characterized by a molecular formula. These
components are treated as pure components, though they are complex
mixtures.
Where:
wij = Mass fraction of the jth constituent in component i
The elements of GENANAL are used to input the mass percentages of the
constituents. The structure for HCGEN is: Elements 1 to 4 are the four
coefficients for the first constituent, elements 5 to 8 are the coefficients for
the second constitutent, and so on, for up to 20 constituents.
Notation
Most correlations for the calculation of coal properties require proximate,
ultimate, and other analyses. These are converted to a dry, mineral-matter-
free basis. Only the organic portion of the coal is considered.
Moisture corrections are made for all analyses except hydrogen, according to
the formula:
Where:
w = The value determined for weight fraction
d
w = The value on a dry basis
= The moisture weight fraction
The mineral matter content is calculated using the modified Parr formula:
The ash term corrects for water lost by decomposition of clays in the ash
determination. The average water constitution of clays is assumed to be 11.2
percent. The sulfur term allows for loss in weight of pyritic sulfur when pyrite
is burned to ferric oxide. The original Parr formula assumed that all sulfur is
pyritic sulfur. This formula included sulfatic and organic sulfur in the mineral-
matter calculation. When information regarding the forms of sulfur is
available, use the modified Parr formula to give a better approximation of the
percent of inorganic material present. Because chlorine is usually small for
United States coals, you can omit chlorine from the calculation.
Correct analyses from a dry basis to a dry, mineral-matter-free basis, using
the formula:
Where:
The oxygen and organic sulfur contents are usually calculated by difference
as:
Where:
Cp = Heat capacity / (J/kgK)
cp = Heat capacity / (cal/gC)
h = Specific enthalpy
ρ = Specific density
Subscripts:
A = Ash
C = Carbon
Cl = Chlorine
FC = Fixed carbon
H = Hydrogen
H2O = Moisture
MM = Mineral matter
N = Nitrogen
O = Oxygen
So = Organic sulfur
Sp = Pyritic sulfur
St = Total sulfur
Dulong Correlation:
IGT Correlation:
Where:
Where:
i = Component index
j = Constituent index j = 1, 2 , ... , ncn
Where the values of j represent:
1 Moisture
2 Fixed carbon
3 Primary volatile matter
4 Secondary volatile matter
5 Ash
wj = Mass fraction of jth constituent on dry basis
This correlation calculates heat capacity in cal/gram-C using temperature in
C. The parameters must be specified in appropriate units for this conversion.
CP1C/7 ai,23
CP1C/8 ai,24 0
CP1C/9 ai,31 0.395
CP1C/10 ai,32
CP1C/11 ai,33 0
CP1C/12 ai,34 0
CP1C/13 ai,41 0.71
CP1C/14 ai,42
CP1C/15 ai,43 0
CP1C/16 ai,44 0
CP1C/17 ai,51 0.18
CP1C/18 ai,52
CP1C/19 ai,53 0
CP1C/20 ai,54 0
CP2C/3 a3i
CP2C/4 a4i
The default values of the parameters were developed by Gomez, Gayle, and
Taylor (1965). They used selected data from three lignites and a
subbituminous B coal, over a temperature range from 32.7 to 176.8°C.
References
Gomez, M., J.B. Gayle, and A.R. Taylor, Jr., Heat Content and Specific Heat of
Coals and Related Products, U.S. Bureau of Mines, R.I. 6607, 1965.
IGT (Institute of Gas Technology), Coal Conversion Systems Technical Data
Book, Section PMa. 44.1, 1976.
Kirov, N.Y., "Specific Heats and Total Heat Contents of Coals and Related
Materials are Elevated Temperatures," BCURA Monthly Bulletin, (1965),
pp. 29, 33.
Swanson, V.E. et al., Collection, Chemical Analysis and Evaluation of Coal
Samples in 1975, U.S. Geological Survey, Open-File Report (1976), pp. 76–
468.
The equation for ρidm is good for a wide range of hydrogen contents, including
anthracities and high temperature cokes. The standard deviation of this
correlation for a set of 190 points collected by IGT from the literature was
12x10-6 m3/kg. The points are essentially uniform over the whole range. This
is equivalent to a standard deviation of about 1.6% for a coal having a
hydrogen content of 5%. It increases to about 2.2% for a coke or anthracite
having a hydrogen content of 1%.
Parameter Name/Element Symbol Default
DENIGT/1 a1i 0.4397
DENIGT/2 a2i 0.1223
DENIGT/3 a3i -0.01715
DENIGT/4 a4i 0.001077
Reference
IGT (Institute of Gas Technology), Coal Conversion Systems Technical Data
Book, Section PMa. 44.1, 1976.
References
I.M. Chang, B.S. Thesis, Massachusetts Institute of Technology, 1979.
M. Gomez, J.B. Gayle, and A.R. Taylor, Jr., Heat Content and Specific Heat of
Coals and Related Products, U.S. Bureau of Mines, R.I. 6607, 1965.
IGT (Institute of Gas Technology), Coal Conversion Systems Technical Data
Book, Section PMa. 44.1, 1976.
N.Y. Kirov, "Specific Heats and Total Heat Contents of Coals and Related
Materials are Elevated Temperatures," BCURA Monthly Bulletin, (1965),
pp. 29, 33.
V.E. Swanson et al., Collection, Chemical Analysis and Evaluation of Coal
Samples in 1975, U.S. Geological Survey, Open-File Report (1976), pp. 76–
468.
Note: The enthalpy, entropy, Gibbs energy, and molar volume of water are
calculated from the steam tables when the relevant option is enabled. The
total properties are mole-fraction averages of these values with the properties
calculated by the equation of state for other components. Fugacity coefficient
is not affected.
Soave-Redlich-Kwong Option
Codes
There are five related models all based on the Soave-Redlich-Kwong equation
of state which are very flexible and have many options. These models are:
• Standard Redlich-Kwong-Soave (ESRKSTD0, ESRKSTD)
• Redlich-Kwong-Soave-Boston-Mathias (ESRKS0, ESRKS)
• Soave-Redlich-Kwong (ESSRK, ESSRK0)
• SRK-Kabadi-Danner (ESSRK, ESSRK0)
• SRK-ML (ESRKSML, ESRKSML0)
The options for these models can be selected using the option codes
described in the following table:
Option Value Description
Code
1 0 Standard SRK alpha function for Tr < 1, Boston-Mathias alpha
function for Tr > 1
1 Standard SRK alpha function for all
2 Grabovsky – Daubert alpha function for H2 and standard SRK
alpha function for others (default)
3 Extended Gibbons-Laughton alpha function for all components
(see notes 1, 2, 3)
4 Mathias alpha function
5 Twu generalized alpha function
2 0 Standard SRK mixing rules (default except for SRK-Kabadi-
Danner)
1 Kabadi – Danner mixing rules (default for SRK-Kabadi-Danner)
(see notes 3, 4, 5)
3 0 Do not calculate water properties from steam table (default for
Redlich-Kwong-Soave models)
Notes
1. The standard alpha function is always used for Helium.
2. If extended Gibbons-Laughton alpha function parameters are missing, the
Boston-Mathias extrapolation will be used if T > Tc, and the standard
alpha function will be used if T < Tc.
3. The extended Gibbons-Laughton alpha function should not be used with
the Kabadi-Danner mixing rules.
4. The Kabadi-Danner mixing rules should not be used if Lij parameters are
provided for water with any other components.
5. It is recommended that you use the SRK-KD property method rather than
change this option code.
6. The enthalpy, entropy, Gibbs energy, and molar volume of water are
calculated from the steam tables when this option is enabled. The total
properties are mole-fraction averages of these values with the properties
calculated by the equation of state for other components. Fugacity
coefficient is not affected.
264 Index
DIPPR model 140, 150, 165, 191, G
197, 217, 221, 233
General models for
heat of vaporization 140
nonconventional components
ideal gas heat capacity 165
239, 240, 244
liquid surface tension 233
coal model for enthalpy 244
liquid thermal conductivity 217
density polynomial model 240
liquid volume 150
enthalpy and density models list
surface tension 233
239
vapor thermal conductivity 221
heat capacity polynomial model
vapor viscosity 197
240
viscosity 191, 197
Grayson-Streed fugacity model 137
DIPPR/IK-CAPE model 161
Group contribution activity
liquid heat capacity 161
coefficient models 121, 123,
DNSTYGEN nonconventional
124
component density model 240
Dortmund-modified UNIFAC 123
Lyngby-modified UNIFAC 124
E UNIFAC 121
Electrolyte models 145, 176, 178,
207, 216, 229, 237 H
Clarke liquid volume 145
Hakim-Steinberg-Stiel/DIPPR
electrolyte NRTL enthalpy 176
surface tension 233
Gibbs energy 178
Hansen solubility parameter model
Jones-Dole viscosity 207
97
Nernst-Hartley diffusion 229
Hayden-O'Connell 20
Onsager-Samaras surface
property model 20
tension 237
HCOALGEN general coal model for
Riedel thermal conductivity 216
enthalpy 244
Electrolyte NRTL 79, 176, 178
Heat capacity models list 160
activity coefficient model 79
Heat of vaporization 139
enthalpy thermodynamic
model 139
property model 176
Helgeson thermodynamic property
Gibbs energy thermodynamic
model 184
property model 178
Henry's constant solubility
eNRTL-SAC activity coefficient
correlation model 171
model 92
HF equation of state 22
Enthalpies model based on
property model 22
different reference states 180
Huron-Vidal mixing rules 64
Enthalpy and density models for
coal and char 241
ENTHGEN nonconventional I
component heat capacity IAPS models for water 206, 215,
model 240 233
Equation-of-state 14 surface tension 233
property models 14 thermal conductivity 215
extrapolation 9 viscosity 206
Ideal gas law 26
F property model 26
Ideal gas/DIPPR heat capacity
Fugacity models list 131
model 165
Ideal liquid activity coefficient
model 98
IGT density model for 250
Index 265
char 250 property model 79
coal 250 NRTL-SAC 100
Using 104
J
O
Jones-Dole electrolyte correction
viscosity model 207 Onsager-Samaras electrolyte
surface tension model 237
K option codes
activity coefficient models 253
Kent-Eisenberg fugacity model 138 enthalpy models 260
equation of state models 255
L gibbs energy models 261
soave-redlich-kwong models 258
Lee-Kesler Plöcker property model
transport property models 252
28
Lee-Kesler property model 26
Letsou-Stiel viscosity model 209 P
Li mixing rule thermal conductivity Peng-Robinson 32, 34, 36, 54
model 216 alpha functions 54
Liquid Constant Molar Volume MHV2 property model 36
Model 150 property model 34
Liquid enthalpy 179 Wong-Sandler property model 36
thermodynamic property model Physical properties 11, 14
179 models 11, 14
Lucas vapor viscosity model 209 Pitzer activity coefficient model 105
Polynomial activity coefficient
M model 117
PPDS Heat of Vaporization
Mathias alpha function 58
Equation 141
Mathias-Copeman alpha function
PPDS Liquid Heat Capacity 163
54, 58
Predictive Soave-Redlich-Kwong-
Maxwell-Bonnell vapor pressure
Gmehling mixing rules 67
model 139
Predictive SRK property model
MHV2 mixing rules 66
(PSRK) 36
Modified MacLeod-Sugden surface
Property models 5, 11, 14
tension model 238
equation-of-state list 14
Modified Rackett model for molar
list of 5
volume 156
thermodynamic list 11
Molar volume and density models
PSRK 36
list 143
property model 36
N
R
Nernst-Hartley electrolyte diffusion
Rackett mixture liquid volume
model 229
model 155
Nonconventional solid property
Rackett/DIPPR pure component
models 239
liquid volume model 150
density 239
Redlich-Kister activity coefficient
enthalpy 239
model 118
list of 239
Redlich-Kwong 37, 58
Nothnagel 30
alpha function 58
property model 30
property model 37
NRTL 79
266 Index
Redlich-Kwong-Aspen property mixtures 232
model 37 models list 232
Redlich-Kwong-Soave 40, 42, 44, Modified MacLeod-Sugden 238
58
alpha function equations 58 T
alpha function list 58
Boston-Mathias property model Thermal conductivity 213, 220
40 models list 213
MHV2 property model 42 solids 220
Soave-Redlich-Kwong property thermo switch 7
model 44 Thermodynamic property 11, 173
Wong-Sandler property model 42 list of additional models 173
Redlich-Kwong-Soave property models list 11
model 38 Three-suffix Margules activity
Riedel electrolyte correction coefficient model 120
thermal conductivity model THRSWT 7
216 Transport property 188
models list 188
transport switch 7
S
TRAPP 211, 224
Sato-Riedel/DIPPR thermal thermal conductivity model 224
conductivity model 217 viscosity model 211
Scatchard-Hildebrand activity TRNSWT 7
coefficient model 119 Twu Liquid Viscosity 211
Schwartzentruber-Renon property
model 42 U
Soave-Redlich-Kwong property
model 44 UNIFAC 121, 123, 124
Solid Antoine vapor pressure activity coefficient model 121
models 139 Dortmund modified activity
Solid Thermal Conductivity coefficient model 123
Polynomial 220 Lyngby modified activity
Solids polynomial heat capacity coefficient model 124
model 168 UNIQUAC 125
Solubility correlation models 170 activity coefficient model 125
list 170 Units for Temperature-Dependent
SRK-Kabadi-Danner property Properties 6
model 46
SRK-ML property model 48 V
Standard Peng-Robinson property
Van Laar activity coefficient model
model 34
126
Standard Redlich-Kwong-Soave
Vapor pressure model list 131
property model 38
Viscosity 190
Steam tables 15, 29
models 190
NBS/NRC 29
property models 15
STEAMNBS property method 29 W
Stiel-Thodos pressure correction Wagner Interaction Parameter
thermal conductivity model activity coefficient model 127
223 Wagner vapor pressure model 131
Stiel-Thodos/DIPPR thermal Wassiljewa-Mason-Saxena mixing
conductivity model 221 rule for thermal conductivity
Surface tension 232, 238 225
Index 267
Water solubility model 172
Watson equation for heat of
vaporization 140
Wilke-Chang diffusion model 230,
231
binary 230
mixture 231
WILS-GLR property method 180
WILS-LR property method 180
Wilson (liquid molar volume)
activity coefficient model 129
Wilson activity coefficient model
128
Wong-Sandler mixing rules 69
268 Index