Species Table 9.3 - Density
Navigation: User Guide -> Species Table 9.3 -> Species Table 9.3 - Density
Species Table | Heat of Formation and Entropy |
Density | Specific Heat (Cp) | Phase Change (solubility) |
BPE & Acid/Base Ka/b |
Vapour Properties | Viscosity |
---|
Related Links: Editing User Species Database 9.3, Species Properties ($SDB)
This page is valid for SysCAD 9.3. For SysCAD 9.2, please see Species Table 9.2 and Species Table 9.2 - Solution Data.
Contents
Introduction
The Density field is optional, but it should be filled in to ensure that the density calculations for any stream containing this compound are correct.
This is the density of the species in the defined phase. The unit for density is kg/m^{3}.
Notes:
- SysCAD works with mass flows and uses the density to convert mass flows to volume flows. Thus in order to get accurate volume flows, the user must specify accurate density values for their species.
- The user may enter a constant value for species in all phases, but in the cases of aqueous and gas species it is preferable to use functions, as defined further on.
- If this field is left blank SysCAD will assume the following:
- A constant value of 2000 kg/m^{3} for solids,
- A constant value of 1000 kg/m^{3} for liquids, and
- Ideal gas density for gases. (Note: if a constant is entered, Linear Gas Density will be assumed)
The density provided here will also be used in volume calculations. The equation used is Volume = mass / density.
Refer to Stream Density for an example of how these individual densities are used to determine the density of a stream.
The list of density functions that are available will change if the species is a Solid, a Liquid or a Gas. The image below shows the different functions that are available for the different phases:
Constant Density
The user may enter a constant value for the density for a solid, liquid or gas species, by selecting 'Constant Value'.
In this case the density of the species will not change as a function of temperature or as the mass fraction of the species changes in a stream.
The constant density option is normally used for Solid species and pure Liquid species.
Density as a Function of Temperature
This option is only available for Solid or Liquid species.
If the user has data for changing density as a function of temperature, then they may enter that using:
- Density as a function of T : Polynomial; OR
- Spline under Tension.
Polynomial Function
If the user has a polynomial function that describes density as a function of temperature, then they can select the option:
- Density as a function of T : Polynomial
and enter the function in the following format: Poly_T(C0, C1, C2, C3, C4, C5), where
- [math]Density = C0 + C1T + C2T^2 + C3T^3 + C4T^4 + C5T^5\,[/math]
NOTES:
- T is Temperature in K.
- The user does not need to have all the constants, any that are not required may be set to zero. So, for example if you have a 3rd order equation, you would only need to enter C0, C1, C2 and C3. C4 and C5 would be zero.
Table of Data
If the user has a table with Density versus Temperature, then they may select the option:
- Spline under Tension.
And enter the data in 2 columns. SysCAD will interpolate the data using the TSpline method. The easiest way to enter the data is to copy the data from a spreadsheet.
Note: A fitted equation will be processed more quickly within SysCAD and hence is preferred to entering a table of data.
Density as a Function of Mass Fraction
This option is only available for Solid or Liquid species.
If the user has a polynomial function that describes density as a function of mass fraction, then they can select the option:
- Density as a function of MF : Polynomial
and enter the function in the following format: Poly_MF(C0, C1, C2, C3, C4, C5), where
- [math]Density = C0 + C1*MF + C2*MF^2 + C3*MF^3 + C4*MF^4 + C5*MF^5\,[/math]
NOTES:
- MF is Mass Fraction in %.
- The user does not need to have all the constants, any that are not required may be set to zero. So, for example if you have a 2nd order equation, you would only need to enter C0, C1 and C2. C3, C4 and C5 would be zero.
Liquid Density Correction Functions
The following methods are only available for liquid species, and are normally used for dissolved aqueous species, for example, NaCl(aq).
The density of a solution containing aqueous or ionic species changes according to the mass fraction of dissolved species in solution. The Density Correction functions describe the solution density as a function of solute mass fraction.
Emulate Water Density
If the user has defined a species that will normally exist in the aqueous form, but the user does not have a density correction function, then it is recommended that the density of the species be set to the water density using the special density function LiqH2ORho(), by choosing the 'Water Density Function' option. This will ensure that the species has the same density as water and hence the species will not change the density of the solution. For further info on this equation see Water Density.
Polynomial Function
If the user has a polynomial function that describes density correction factors as a function of mass fraction, then they can select the option:
- Density Correction as a function of MF : Polynomial
and enter the function in the following format: Poly_MF(a, b, c, d, e, f)
Once the data has been entered in the user data dialog box, the data will be summarised in the Species Database as:
Poly(a, b, c, d, e, f), Limit(Limit Fraction, Pure Value, Generate Warnings)
This consists of 2 parts:
- The Polynomial, which represents the density correction function:
- [math] DensCorrFn(MF)_i = a + b.MF + c.MF^2 + d.MF^3 + e.MF^4 + f.MF^5 \,[/math]
- Where MF is the mass fraction of the species in solution.
- Notes:
- SysCAD ignores the a value in the above polynomial, as it represents the pure solution density based on the data for the solvent - usually water. (The value for 'a' is expected because this is the normal form of the equation that users will obtain when carrying out curve fitting). Please see Implementation of the Density correction Function for further information.
- The user only needs to enter the required number of parameters. For example, if the polynomial is 2nd order you need only enter a,b and c.
- As a minimum, the user must enter values for a and b, as this will then be first order with respect to MF.
- If the user only enters a value for a, then the species will have the same density as the solvent, usually water.
- [math] DensCorrFn(MF)_i = a + b.MF + c.MF^2 + d.MF^3 + e.MF^4 + f.MF^5 \,[/math]
- The Limiting Values, which are described below:
- Limit Fraction is the maximum mass fraction of the solute for which the equation is valid.
- If the mass fraction of the solute in the solvent is above the Limit Fraction value, then SysCAD calculates the density of the solution by linearly interpolating between the density at the Limit Fraction and the density of the pure solute (Pure Value).
- Pure Value is only used if MF is greater than Limit Fraction. The purpose of Pure Value is to give a linear extension of the curve, so that reasonable densities are calculated when MF is greater than Limit Fraction. The Pure Value used is typically the density of the pure solute in kg/m^3, but can be any suitable value for linear interpolation between the Limit Fraction and a MF of 100%.
- Generate Warnings is either On or Off. If it is On (recommended), SysCAD will warn the user if the mass fraction of the solute is above the specified Limit Fraction.
Example:
For FeSO4(aq): Poly(0.998,0.951,0.62), Limit(0.2, 2200, On)
Notes:
- If the user does not have a function relating the change in solvent density with solute concentration, then SysCAD will use the density in the species database in a pure mass weighted mean calculation. However, this does not usually produce the correct liquid density or volumetric flow rates.
- Please see Density Correction Calculations for the implementation method for density correction, a description of the 2 methods of calculating the species Mass Fraction, MFi and Density Calculation Examples.
- The user may view the corrected density of a solution consisting of ONLY the solvent and the solute on the Species Properties ($SDB) 9.3 access window.
Table of Data
If the user has a table with density correction factors versus mass fraction in solution, then they may select the option:
- Density Correction as a function of MF : Spline under Tension
And enter the data in 2 columns. SysCAD will interpolate the data using the TSpline method. The easiest way to enter the data is to copy the data from a spreadsheet.
Note: A fitted equation will be processed more quickly within SysCAD and hence is preferred to entering a table of data.
Laliberte´ Function
This method is only valid for aqueous species.
If the user has appropriate constants for the Laliberte´ equation, then they can select the option:
- Density Correction as a function of MF : Laliberte values
and enter the function in the following format: Laliberte_Rho(c_{0}, c_{1}, c_{2}, c_{3}, c_{4})
The solution density is calculated using the water density, p_{w} and the solutes apparent density using following equation:
- [math]\mathbf{\mathit{p_m=\frac{1}{\frac{m_w}{p_w}+\sum{\frac{m_i}{p_{app,i}}}}}}[/math]
The Apparent density of each solute in aqueous solution is calculated from:
- [math]\mathbf{\mathit{p_{app,i}=\frac{(c_0(1-m_w)+c_1)*e^{(0.000001(T+c_4)^2)}}{(1-m_w)+c_2+c_3*T}}}[/math]
Where:
m_{w} | = | mass fraction of water |
m_{i} | = | mass fraction of solute species i |
p_{w} | = | density of water (at stream temperature and pressure), kg/m^{3} |
p_{app,i} | = | solute i apparent density, kg/m^{3} |
p_{m} | = | solution density, kg/m^{3} |
T | = | Temperature in °C |
c_{0} to c_{4} | = | dimensionless empirical constants for each solute species. |
Notes:
- If the user wishes to use the Laliberte´ method, then it is recommended that this method is used for ALL aqueous species in a project, as a mixture of Laliberte´ and other aqueous density methods may result in inaccurate solution density values.
- It is important to note that the constants for many of the aqueous species are valid for temperatures between 0 and approximately 100°C.
- If the unit temperature is outside of the species temperature range, then SysCAD will use the values at the temperature limit.
- Water density is calculated in SysCAD as described here: Water and Steam Properties.
Reference
Laliberte´ M. and Cooper W.E. Model for Calculating the Density of Aqueous Electrolyte Solutions J. Chem. Eng. Data 2004, 49.
Gas Density
For gases the user may use one of the following three input formats (the formulation for Ideal Gas is also shown here):
- Constant Value; or
- Ideal Gas Density - The density value will be calculated based on the Ideal gas law. Equations used are:
- (1) [math] \mathbf {\mathrm{Density_{T,P} = \frac{m}{V}}}[/math] and (2) [math] \mathbf {\mathrm{PV = nRT}}[/math] and (3) [math] \mathbf {\mathrm{m = nM}}[/math]
- Rearranging the above equations will give:
- [math] \mathbf {\mathrm{Density_{T,P} = \frac{PM}{RT}}}[/math]
- Where:
- m = mass of compound
- V = Volume of compound
- P = Partial Pressure of species
- R = Universal Gas Constant = 8.314 472 J/mol.K (Reference: National Institute of Standards and Technology)
- T = Temperature in Kelvin
- n = number of moles of compound
- M = molecular weight of compound
- Linear Gas Density - the density value provided (in brackets) is expected to be at 0°C and Std. Pressure. Density @ T, P will be corrected based on:
- [math] \mathbf {\mathrm{Density_{T,P} = Density_{0,StdP}*\frac{P}{StdP}*\frac{273(K)}{T(K)}}}[/math]; or