Water and Steam Properties: Difference between revisions
Merry.Huang (talk | contribs) |
No edit summary |
||
Line 7: | Line 7: | ||
== Introduction == | == Introduction == | ||
As H<sub>2</sub>O (steam | As H<sub>2</sub>O (water / steam) is a common and important component of most models, the property equations for both water and steam are built into SysCAD. These properties are hardcoded and cannot be changed. | ||
However, the user may select the [[Plant_Model_-_Species#Species_Property_Overrides|properties]] of H<sub>2</sub>O to be calculated using either [[#High Fidelity (IF97)|High Fidelity (IF97)]], [[#LowFid|Fast IF97 Sat]] or [[#Simple|Simple]] equations | However, the user may select the [[Plant_Model_-_Species#Species_Property_Overrides|properties]] of H<sub>2</sub>O to be calculated using either [[#High Fidelity (IF97)|'''High Fidelity (IF97)''']], [[#LowFid|'''Fast IF97 Sat''']] or [[#Simple|'''Simple''']] equations (discussed in more detail below). High Fidelity equations will return the most accurate results for T and P, but the calculations do require more computing power, slowing the solver speed. | ||
'''Notes:''' | '''Notes:''' | ||
# Properties for solid | # Properties for solid H<sub>2</sub>O (ice, or water trapped in a solid crystalline structure) are not included. To include these, add a species to the Species Database in the normal way with Phase = 'Solid' but the compound must NOT be 'H2O' and the Definition must NOT be H2O1, rather use '''O1H2'''. | ||
# If alternate user defined properties for water or steam are required, | # If alternate user-defined properties for water or steam are required, it is possible to "trick" SysCAD using a species with Definition 'O1H2' and a Component that is not "H2O". | ||
#* IMPORTANT: Caution should be used when including alternate properties for steam/water as the hardwired steam/water will also exist in the project! | #* IMPORTANT: Caution should be used when including alternate properties for steam/water as the hardwired steam/water will also exist in the project! (H2O(l) and H2O(g) cannot be removed from the species database.) | ||
# SysCAD does a range check on streams and many units to determine if liquid water, H2O(l), is present above the critical temperature of water. If it detects water above Tc it will generate a warning message. Please also see [[Plant Model - RC]] for more information. | # SysCAD does a range check on streams and many units to determine if liquid water, H2O(l), is present above the critical temperature (Tc) of water. If it detects water above Tc it will generate a warning message. Please also see [[Plant Model - RC]] for more information. | ||
== High Fidelity (IF97) == | == High Fidelity (IF97) == | ||
Line 21: | Line 21: | ||
=== IF97 === | === IF97 === | ||
* | *Full details of the algorithms used to calculate the thermodynamic properties of steam/water, such as heat capacity and enthalpy, saturated temperature/pressure, etc., are not presented here, as they are available in the listed references. | ||
*The implementation in SysCAD has been optimised to maximise performance | *The implementation in SysCAD has been optimised to maximise performance and solver speed. | ||
* | *Calculations are based on the IAPWS-IF97 standard, described in a paper by W. Wagner et al. "The IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam", published in ASME J. Eng. Gas Turbines and Power, Vol. 122 (2000). | ||
The valid temperature range is: | |||
:Steam: 273.16 to 3000 K | |||
:Water: 273.16 to 647 K | |||
Some good internet references for further information are: | Some good internet references for further information are: | ||
:http://www.iapws.org/ | :[http://www.iapws.org/ The International Association for the Properties of Water and Steam (http://www.iapws.org/)] | ||
:http://www.cheresources.com/iapwsif97.shtml | :[http://www.cheresources.com/iapwsif97.shtml Thermodynamic and Transport Properties of Water and Steam (http://www.cheresources.com/iapwsif97.shtml)] | ||
===IF97 Equation Names=== | ===IF97 Equation Names=== | ||
In the [[Species Table - Density|SysCAD species database]], equation names for full IF97 correlations are: | In the [[Species Table - Density|SysCAD species database]], equation names for full IF97 correlations are: | ||
<syntaxhighlight | <syntaxhighlight inline>IF97WaterRho(), IF97WaterCp(), IF97SteamRho() and IF97SteamCp() </syntaxhighlight> | ||
'''Note:''' Due to differences in reference points, numerical methods and implementation techniques there can be very small (generally much less than 0.1%) differences between IAPWS results and SysCAD results. | |||
=== Water Properties below 0°C === | |||
A set of equations for water Cp, H, S and density between 235K and 273.15K are {{Available139||y}}. | |||
To enable these equations, please check the tick boxes for '''H2O(l).UseNewLowT.Cp''' and or '''H2O(l).UseNewLowT.Rho''' on the View - [[Plant_Model_-_Species#Water_properties_below_0C|Plant Model - Species]] Tab. | To enable these equations, please check the tick boxes for '''H2O(l).UseNewLowT.Cp''' and or '''H2O(l).UseNewLowT.Rho''' on the View - [[Plant_Model_-_Species#Water_properties_below_0C|Plant Model - Species]] Tab. | ||
*If these options are selected, then water Cp, H, S and Density between 235K and 273.15K are calculated based on | *If these options are selected, then water Cp, H, S and Density between 235K and 273.15K are calculated based on references below. For temperatures below 235K (-38.15°C), the properties are extrapolated. | ||
*If UseNewLowT options are not selected, then water Cp and Density are extrapolated below 0°C. (This is original method used in | *If UseNewLowT options are not selected, then water Cp and Density are extrapolated below 0°C. (This is original method used in {{Removed138|||y}} and earlier, retained here for backward compatibility purposes.) | ||
*UseNewLowT options are auto selected for any new project or existing project with a minimum temperature >0°C. | *UseNewLowT options are auto selected for any new project or existing project with a minimum temperature >0°C. | ||
*For projects with below 0°C temperatures, user can select this option to improve the accuracy of the enthalpy and density values between 235K and 273.15K. | *For projects with below 0°C temperatures, user can select this option to improve the accuracy of the enthalpy and density values between 235K and 273.15K. | ||
Line 51: | Line 52: | ||
References: | References: | ||
# V. Holten, C. E. Bertrand, M. A. Anisimov, and J. V. Sengers | # V. Holten, C. E. Bertrand, M. A. Anisimov, and J. V. Sengers (2012) "Thermodynamics of Supercooled Water", J. Chem. Phys. 136, 094507 | ||
# | # Smith, Van Ness, Abbott (1996) “Introduction to Chemical Engineering Thermodynamics", 5th Ed | ||
===IF97 Equation Version Change Notes=== | ===IF97 Equation Version Change Notes=== | ||
From {{Available139|27886||y}}, a single equation, "IF97" replaces the old equations ("IF97" and "IF97_2") that were available in previous versions of SysCAD. Projects created using older builds will automatically upgrade to always use the updated "IF97" form of the equation. However, it is important to note that this change may affect results if the project was previously using the old equation with steam at high temperatures. | |||
{{collapse | title = IF97 Equation Change - | {{collapse | title = IF97 Equation Change - Detailed notes on how it may affect projects upgraded from Build 138 and earlier. | content = | ||
Prior to {{Removed139|||y}}, there were two equations for Steam Density and Steam Cp: "IF97" and "IF97_2". | |||
* | * The original "IF97" (from SysCAD 9.2) was found to have errors in steam Enthalpy and Density equations at very high temperatures (Region 5). | ||
* A corrected version | * A corrected version, "IF97_2", was implemented and the original "IF97" was retained for backwards compatibility. | ||
* In the SysCAD species database, equation names for | * In the SysCAD species database, equation names for the "IF97" and "IF97_2" correlations were: IF97WaterRho(), IF97WaterCp(), IF97SteamRho(), IF97SteamRho2(), IF97SteamCp() and IF97SteamCp2() | ||
From {{Available139|27886||y}}, the "IF97" and "IF97_2" equations have been merged into a single equation "IF97" (equivalent to the old "IF97_2"). | |||
* | * "IF97" is now the only choice. This uses the correct steam Enthalpy and Density equations for Region 5. | ||
* Projects and | * Projects and CFG files are auto-upgraded to always use "IF97". | ||
* | * "IF97_2" is removed as an option from various dropdown lists (e.g. Plant Model - Species Tab). | ||
Notes for project upgrade: | Notes for project upgrade: | ||
* For project upgrade from Build 138 to Build 139, the CFG files are auto-upgraded to use the new "IF97" form. The updated CFG files are '''not''' backward compatible, so if they are reused in an earlier build there will be load errors/warnings. | |||
* For project upgrade from Build 138 to Build 139, the | |||
* The change of equation may affect results if the project was previously using the old equation with steam at high temperatures. | * The change of equation may affect results if the project was previously using the old equation with steam at high temperatures. | ||
}} | |||
== Fast IF97 Sat == | == Fast IF97 Sat == | ||
Line 87: | Line 85: | ||
|- | |- | ||
| ''Water Density'' | | ''Water Density'' | ||
| | | Special density function '''<font color="green">LiqH2ODensity()</font>'''<br> See [[Fast IF97 Sat Properties]]. | ||
|The equation for water density is an equation fitted to IF97 water density along the saturation line. The valid temperature range is 260K to 647.096K. <br>If the temperature is outside of the valid range then the minimum or maximum temperatures are used to calculate the density. <br>In the [[Species Table - Density|SysCAD species database]], any species can be assigned this density by use of a special density function '''<font color="green">LiqH2ODensity()</font>. | |The equation for water density is an equation fitted to IF97 water density along the saturation line. The valid temperature range is 260K to 647.096K. <br>If the temperature is outside of the valid range then the minimum or maximum temperatures are used to calculate the density. <br>In the [[Species Table - Density|SysCAD species database]], any species can be assigned this density by use of a special density function '''<font color="green">LiqH2ODensity()</font>. | ||
|- | |- | ||
| ''Water Cp'' | | ''Water Cp'' | ||
| | |Special Cp function '''<font color="green">LiqH2OCp()</font>'''. <br> See [[Fast IF97 Sat Properties]]. | ||
|The equation for water Cp is an equation fitted to IF97 water heat capacity along the saturation line. The valid temperature range is 260K to 573.15K. <br>In the [[Species Table - Density|SysCAD species database]], any species can be assigned this Cp by use of a special Cp function '''<font color="green">LiqH2OCp()</font>'''. | |The equation for water Cp is an equation fitted to IF97 water heat capacity along the saturation line. The valid temperature range is 260K to 573.15K. <br>In the [[Species Table - Density|SysCAD species database]], any species can be assigned this Cp by use of a special Cp function '''<font color="green">LiqH2OCp()</font>'''. | ||
|- | |- | ||
| ''Steam Density'' | | ''Steam Density'' | ||
| | | Special Density function '''<font color="green">VapH2ODensity()</font>'''.<br> See [[Fast IF97 Sat Properties]]. | ||
|The density is calculated using the IdealGasDensity equation with the temperature ranged between 260K and 4000K. <br>In the [[Species Table - Density|SysCAD species database]], any species can be assigned this Density by use of a special Density function '''<font color="green">VapH2ODensity()</font>'''. | |The density is calculated using the IdealGasDensity equation with the temperature ranged between 260K and 4000K. <br>In the [[Species Table - Density|SysCAD species database]], any species can be assigned this Density by use of a special Density function '''<font color="green">VapH2ODensity()</font>'''. | ||
|- | |- | ||
| ''Steam Cp'' | | ''Steam Cp'' | ||
| | | Special Cp function '''<font color="green">VapH2OCp()</font>'''.<br> See [[Fast IF97 Sat Properties]]. | ||
|The equation for steam Cp is an equation fitted to IF97 steam heat capacity along the saturation line. <br>In the [[Species Table - Density|SysCAD species database]], any species can be assigned this Cp by use of a special Cp function '''<font color="green">VapH2OCp()</font>'''. | |The equation for steam Cp is an equation fitted to IF97 steam heat capacity along the saturation line. <br>In the [[Species Table - Density|SysCAD species database]], any species can be assigned this Cp by use of a special Cp function '''<font color="green">VapH2OCp()</font>'''. | ||
|} | |} | ||
==Simple== | ==Simple== | ||
The | The "Simple" set of equations is less accurate and is retained for backward compatibility. However, for Water Density, Water Cp, and Steam Cp, it is recommended to use the new [[#Fast IF97 Sat|Fast IF97 Sat]] equations - these equations are faster, more accurate, and valid for larger temperature ranges. | ||
{| class="wikitable" | {| class="wikitable" | ||
Line 113: | Line 111: | ||
|<math> \cfrac{a_0+a_1T+a_2T^2+a_3T^3+a_4T^4+a_5T^5}{1+b_1T} </math> | |<math> \cfrac{a_0+a_1T+a_2T^2+a_3T^3+a_4T^4+a_5T^5}{1+b_1T} </math> | ||
| | | | ||
* T | * T = Temperature in Celcius | ||
*The valid temperature range is -30°C to 300°C. | *The valid temperature range is -30°C to 300°C. | ||
* If the temperature is outside of the valid range given above, then the minimum or maximum temperatures are used to calculate the density. | * If the temperature is outside of the valid range given above, then the minimum or maximum temperatures are used to calculate the density. | ||
Line 126: | Line 124: | ||
| | | | ||
where | where | ||
<div style="column-count: | <div style="column-count:4 > | ||
:a = 16.749 | :a = 16.749 | ||
:b = 62.120 | :b = 62.120 | ||
Line 132: | Line 130: | ||
:d = 90.391 | :d = 90.391 | ||
</div> | </div> | ||
:T = | :T = Temperature in Kelvin | ||
:This equation is used in the range from 273.15 to 500K. | :This equation is used in the range from 273.15 to 500K. | ||
Line 141: | Line 139: | ||
| | | | ||
where | where | ||
<div style="column-count: | <div style="column-count:4 > | ||
:P = Pressure | :P = Pressure | ||
:R = | :R = Universal gas constant | ||
:T = Temperature in Kelvin | :T = Temperature in Kelvin | ||
:M = | :M = Molecular weight | ||
</div> | </div> | ||
The density is calculated using a rearrangement of the Ideal gas law. | The density is calculated using a rearrangement of the Ideal gas law. | ||
Line 156: | Line 154: | ||
<math> C_p = a + 10^{-3}b T + \cfrac{10^5c}{T^2} + 10^{-6}dT^2\, </math> | <math> C_p = a + 10^{-3}b T + \cfrac{10^5c}{T^2} + 10^{-6}dT^2\, </math> | ||
|where | |where | ||
<div style="column-count: | <div style="column-count:4 > | ||
:a = -38.14955 | :a = -38.14955 | ||
:b = 463.4602 | :b = 463.4602 | ||
Line 162: | Line 160: | ||
:d = -762.3604 | :d = -762.3604 | ||
</div> | </div> | ||
T = | T = Temperature in Kelvin | ||
This equation is used in the range from 298.15 to 500K. | This equation is used in the range from 298.15 to 500K. | ||
|- | |- | ||
Line 170: | Line 168: | ||
| | | | ||
where | where | ||
<div style="column-count: | <div style="column-count:4 > | ||
:A = -3433.74 | :A = -3433.74 | ||
:B = -12.0063 | :B = -12.0063 | ||
Line 184: | Line 182: | ||
Some of the values that are available to users are: | Some of the values that are available to users are: | ||
* Saturated Temperature | * Saturated Temperature | ||
* Saturated Pressure | * Saturated Pressure | ||
* Density | * Density | ||
* Heats of Formation | * Heats of Formation | ||
* Enthalpy | * Enthalpy |
Revision as of 14:17, 23 February 2024
Navigation: User Guide ➔ Stream Properties ➔ Water and Steam Properties
Standard Species Model | Stream Properties | Special Cases | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
Standard Species Model | Stream Properties Calculations | Density | Heat of Formation | Heat of Dilution | Specific Heat (Cp) | Solubility | Boiling Point Elevation | pH | Charge | Steam and Water | Sulfuric Acid |
Related Links: Species Table, Standard Species Model, Plant Model - Species Property Overrides, Fast IF97 Sat Properties
Introduction
As H2O (water / steam) is a common and important component of most models, the property equations for both water and steam are built into SysCAD. These properties are hardcoded and cannot be changed.
However, the user may select the properties of H2O to be calculated using either High Fidelity (IF97), Fast IF97 Sat or Simple equations (discussed in more detail below). High Fidelity equations will return the most accurate results for T and P, but the calculations do require more computing power, slowing the solver speed.
Notes:
- Properties for solid H2O (ice, or water trapped in a solid crystalline structure) are not included. To include these, add a species to the Species Database in the normal way with Phase = 'Solid' but the compound must NOT be 'H2O' and the Definition must NOT be H2O1, rather use O1H2.
- If alternate user-defined properties for water or steam are required, it is possible to "trick" SysCAD using a species with Definition 'O1H2' and a Component that is not "H2O".
- IMPORTANT: Caution should be used when including alternate properties for steam/water as the hardwired steam/water will also exist in the project! (H2O(l) and H2O(g) cannot be removed from the species database.)
- SysCAD does a range check on streams and many units to determine if liquid water, H2O(l), is present above the critical temperature (Tc) of water. If it detects water above Tc it will generate a warning message. Please also see Plant Model - RC for more information.
High Fidelity (IF97)
IF97
- Full details of the algorithms used to calculate the thermodynamic properties of steam/water, such as heat capacity and enthalpy, saturated temperature/pressure, etc., are not presented here, as they are available in the listed references.
- The implementation in SysCAD has been optimised to maximise performance and solver speed.
- Calculations are based on the IAPWS-IF97 standard, described in a paper by W. Wagner et al. "The IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam", published in ASME J. Eng. Gas Turbines and Power, Vol. 122 (2000).
The valid temperature range is:
- Steam: 273.16 to 3000 K
- Water: 273.16 to 647 K
Some good internet references for further information are:
- The International Association for the Properties of Water and Steam (http://www.iapws.org/)
- Thermodynamic and Transport Properties of Water and Steam (http://www.cheresources.com/iapwsif97.shtml)
IF97 Equation Names
In the SysCAD species database, equation names for full IF97 correlations are:
IF97WaterRho(), IF97WaterCp(), IF97SteamRho() and IF97SteamCp()
Note: Due to differences in reference points, numerical methods and implementation techniques there can be very small (generally much less than 0.1%) differences between IAPWS results and SysCAD results.
Water Properties below 0°C
A set of equations for water Cp, H, S and density between 235K and 273.15K are available from Build 139.
To enable these equations, please check the tick boxes for H2O(l).UseNewLowT.Cp and or H2O(l).UseNewLowT.Rho on the View - Plant Model - Species Tab.
- If these options are selected, then water Cp, H, S and Density between 235K and 273.15K are calculated based on references below. For temperatures below 235K (-38.15°C), the properties are extrapolated.
- If UseNewLowT options are not selected, then water Cp and Density are extrapolated below 0°C. (This is original method used in Build 138. and earlier, retained here for backward compatibility purposes.)
- UseNewLowT options are auto selected for any new project or existing project with a minimum temperature >0°C.
- For projects with below 0°C temperatures, user can select this option to improve the accuracy of the enthalpy and density values between 235K and 273.15K.
- If these new UseNewLowT options are not selected in a project, a warning is given in Message Window - Project settings Tab: PlantModel - Using old Water Density extrapolation below 0C. Review PlantModel.H2O(l).UseNewLowT.xxx. (where xxx = Cp or Density)
References:
- V. Holten, C. E. Bertrand, M. A. Anisimov, and J. V. Sengers (2012) "Thermodynamics of Supercooled Water", J. Chem. Phys. 136, 094507
- Smith, Van Ness, Abbott (1996) “Introduction to Chemical Engineering Thermodynamics", 5th Ed
IF97 Equation Version Change Notes
From Build 139.27886, a single equation, "IF97" replaces the old equations ("IF97" and "IF97_2") that were available in previous versions of SysCAD. Projects created using older builds will automatically upgrade to always use the updated "IF97" form of the equation. However, it is important to note that this change may affect results if the project was previously using the old equation with steam at high temperatures.
Prior to Build 139, there were two equations for Steam Density and Steam Cp: "IF97" and "IF97_2".
From Build 139.27886, the "IF97" and "IF97_2" equations have been merged into a single equation "IF97" (equivalent to the old "IF97_2").
Notes for project upgrade:
|
Fast IF97 Sat
A set of lower fidelity equations called "FastIF97Sat" for water and steam Density, Cp, Enthalpy and Entropy is available from Build 139. These equations are based on accurate curve fits to the IF97 data along the Saturation pressure line. The equations are therefore a function of temperature only, pressure is ignored and it is assumed that the pressure is the saturation pressure. These equations are significantly faster and for many operating conditions, especially for water, are suitable for most simulation requirements where the full accuracy of IF97 is not required. Please refer to Fast IF97 Sat Properties for more detail.
Function | Equation | Notes |
---|---|---|
Water Density | Special density function LiqH2ODensity() See Fast IF97 Sat Properties. |
The equation for water density is an equation fitted to IF97 water density along the saturation line. The valid temperature range is 260K to 647.096K. If the temperature is outside of the valid range then the minimum or maximum temperatures are used to calculate the density. In the SysCAD species database, any species can be assigned this density by use of a special density function LiqH2ODensity(). |
Water Cp | Special Cp function LiqH2OCp(). See Fast IF97 Sat Properties. |
The equation for water Cp is an equation fitted to IF97 water heat capacity along the saturation line. The valid temperature range is 260K to 573.15K. In the SysCAD species database, any species can be assigned this Cp by use of a special Cp function LiqH2OCp(). |
Steam Density | Special Density function VapH2ODensity(). See Fast IF97 Sat Properties. |
The density is calculated using the IdealGasDensity equation with the temperature ranged between 260K and 4000K. In the SysCAD species database, any species can be assigned this Density by use of a special Density function VapH2ODensity(). |
Steam Cp | Special Cp function VapH2OCp(). See Fast IF97 Sat Properties. |
The equation for steam Cp is an equation fitted to IF97 steam heat capacity along the saturation line. In the SysCAD species database, any species can be assigned this Cp by use of a special Cp function VapH2OCp(). |
Simple
The "Simple" set of equations is less accurate and is retained for backward compatibility. However, for Water Density, Water Cp, and Steam Cp, it is recommended to use the new Fast IF97 Sat equations - these equations are faster, more accurate, and valid for larger temperature ranges.
Function | Equation | Notes |
---|---|---|
Water Density (less accurate) |
[math]\displaystyle{ \cfrac{a_0+a_1T+a_2T^2+a_3T^3+a_4T^4+a_5T^5}{1+b_1T} }[/math] |
In the SysCAD species database, any species can be assigned this density by use of a special density function LiqH2ORho(). |
Water Cp (less accurate) |
HSC_Cp(a,b,c,d) [math]\displaystyle{ C_p = a + b.10^{-3} T + \cfrac{c.10^5}{T^2} + d.10^{-6} T^2\, }[/math] |
where
In the SysCAD species database, a species can be assigned this heat capacity as HSC_Cp(16.749, 62.12, 32.798, 90.391):Range(K, 273.15, 500). |
Steam Density (less accurate) |
[math]\displaystyle{ \rho = \cfrac{PM}{RT} }[/math] |
where
The density is calculated using a rearrangement of the Ideal gas law. This equation is used for all temperatures and pressures. |
Steam Cp (less accurate) |
CRC_Cp(a,b,c,d) [math]\displaystyle{ C_p = a + 10^{-3}b T + \cfrac{10^5c}{T^2} + 10^{-6}dT^2\, }[/math] |
where
T = Temperature in Kelvin This equation is used in the range from 298.15 to 500K. |
vapour pressure of H2O (less accurate) |
[math]\displaystyle{ \rho = 0.1333224\times10^{\tfrac A T +B\log_{10} \;(T)+CT+D} }[/math]
(0.1333224 converts pressure from mm Hg to kPa.) |
where
|
Water and Steam Properties in PGM
The user may obtain some water and steam values in PGM code using the Species Database Class.
Some of the values that are available to users are:
- Saturated Temperature
- Saturated Pressure
- Density
- Heats of Formation
- Enthalpy