zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

src/thermophysicalModels/specie/equationOfState: Added documentation with example usage

Browse Source
This commit is contained in:
Henry Weller
2021-08-26 16:44:11 +01:00
parent 987d5906ef
commit c86c64db72
11 changed files with 182 additions and 34 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
src/thermophysicalModels/specie/equationOfState/Boussinesq/Boussinesq.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2015-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2015-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -35,6 +35,24 @@ Description
Coefficient mixing is very inaccurate and not supported,
so this equation of state is not applicable to mixtures.
Usage
\table
Property | Description
rho0 | Reference density
T0 | Reference temperature
beta | Coefficient of thermal expansion
\endtable
Example specification of the Boussinesq equation of state:
\verbatim
equationOfState
{
rho0 1;
T0 300;
beta 3e-03;
}
\endverbatim
SourceFiles
BoussinesqI.H
Boussinesq.C

29
src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGas.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2014-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2014-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -25,7 +25,32 @@ Class
Foam::PengRobinsonGas
Description
PengRobinsonGas gas equation of state.
PengRobinsonGas cubic equation of state for gases
Coefficient mixing is very inaccurate and not supported, so this equation of
state can only be applied to mixtures by using a thermophysical properties
model which mixes specie property value.
Usage
\table
Property | Description
Tc | Critical temperature [K]
Vc | Critical volume [m^3/kmol]
Pc | Critical pressure [Pa]
omega | Acentric factor [-]
\endtable
Example specification of the PengRobinsonGas equation of state for methane:
\verbatim
equationOfState
{
Tc 190.55;
Vc 0.0285;
Pc 4.595e6;
omega 0.0115;
}
\endverbatim
SourceFiles
PengRobinsonGasI.H

30
src/thermophysicalModels/specie/equationOfState/adiabaticPerfectFluid/adiabaticPerfectFluid.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2013-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2013-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -25,11 +25,37 @@ Class
Foam::adiabaticPerfectFluid
Description
Adiabatic perfect fluid equation of state.
Adiabatic perfect fluid equation of state for liquids:
\verbatim
rho = rho0*pow((p + B)/(p0 + B), 1/gamma)
\endverbatim
Equivalent to the Murnaghan-Tait equation of state with different
coefficient specification.
Coefficient mixing is very inaccurate and not supported,
so this equation of state is not applicable to mixtures.
Usage
\table
Property | Description
rho0 | Reference density
R | Fluid constant
\endtable
Example specification of the adiabaticPerfectFluid equation of state
for LNG:
\verbatim
equationOfState
{
rho0 450;
p0 1e5;
gamma 15.4;
B 50000000
}
\endverbatim
SourceFiles
adiabaticPerfectFluidI.H
adiabaticPerfectFluid.C

20
src/thermophysicalModels/specie/equationOfState/icoPolynomial/icoPolynomial.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -28,8 +28,14 @@ Description
Incompressible, polynomial form of equation of state, using a polynomial
function for density.
Coefficient mixing is very inaccurate and not supported,
so this equation of state is not applicable to mixtures.
Coefficient mixing is very inaccurate and not supported, so this equation of
state is not applicable to mixtures.
The polynomial expression is evaluated from:
\verbatim
rho = 1000 - 0.05*T + 0.003*sqr(T)
\endverbatim
Usage
\table
@ -37,7 +43,7 @@ Usage
rhoCoeffs<8> | Density polynomial coefficients
\endtable
Example of the specification of the equation of state:
Example specification of the icoPolynomial equation of state:
\verbatim
equationOfState
{
@ -45,12 +51,6 @@ Usage
}
\endverbatim
The polynomial expression is evaluated as so:
\f[
\rho = 1000 - 0.05 T + 0.003 T^2
\f]
SourceFiles
icoPolynomialI.H
icoPolynomial.C

4
src/thermophysicalModels/specie/equationOfState/icoTabulated/icoTabulated.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2020-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -34,7 +34,7 @@ Usage
rho | Density vs temperature table
\endtable
Example of the specification of the equation of state:
Example specification of the icoTabulated equation of state:
\verbatim
equationOfState
{

22
src/thermophysicalModels/specie/equationOfState/incompressiblePerfectGas/incompressiblePerfectGas.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2012-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2012-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -27,7 +27,25 @@ Class
Description
Incompressible gas equation of state using a constant reference pressure in
the perfect gas equation of state rather than the local pressure so that the
density only varies with temperature and composition.
density only varies with temperature:
\verbatim
rho = pRef/(R*T)
\endverbatim
Usage
\table
Property | Description
pRef | Reference pressure
\endtable
Example specification of the incompressiblePerfectGas equation of state:
\verbatim
equationOfState
{
pRef 1e5;
}
\endverbatim
SourceFiles
incompressiblePerfectGasI.H

20
src/thermophysicalModels/specie/equationOfState/linear/linear.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2013-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2013-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -25,7 +25,7 @@ Class
Foam::linear
Description
Linear equation of state with constant compressibility
Linear equation of state with constant compressibility:
\verbatim
rho = rho0 + psi*p
@ -34,6 +34,22 @@ Description
Coefficient mixing is very inaccurate and not supported,
so this equation of state is not applicable to mixtures.
Usage
\table
Property | Description
rho0 | Reference density
psi | Constant compressibility
\endtable
Example specification of the linear equation of state:
\verbatim
equationOfState
{
rho0 1000;
psi 1e-5;
}
\endverbatim
SourceFiles
linearI.H
linear.C

28
src/thermophysicalModels/specie/equationOfState/perfectFluid/perfectFluid.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2012-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2012-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -26,14 +26,34 @@ Class
Description
Simple extension of the perfect gas equation of state to liquids by
the addition of a constant density off-set.
the addition of a constant density off-set:
Coefficient mixing is very inaccurate and not supported,
so this equation of state is not applicable to mixtures.
\verbatim
rho = rho0 + p/(R*T)
\endverbatim
Coefficient mixing is very inaccurate and not supported, so this equation of
state is not applicable to mixtures.
This equation of state is rather inaccurate and has been superseded by
rPolynomial which is much more accurate and supports mixtures.
Usage
\table
Property | Description
rho0 | Reference density
R | Fluid constant
\endtable
Example specification of the perfectFluid equation of state:
\verbatim
equationOfState
{
rho0 1000;
R 100;
}
\endverbatim
SourceFiles
perfectFluidI.H
perfectFluid.C

14
src/thermophysicalModels/specie/equationOfState/perfectGas/perfectGas.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -25,7 +25,17 @@ Class
Foam::perfectGas
Description
Perfect gas equation of state.
Perfect gas equation of state:
\verbatim
rho = p/(R*T)
\endverbatim
Usage
The gas constant R used by the perfect gas equation of state is obtained
directly from the specie or mixture molecular weight so there is no need to
provide an equation of \c equationOfState enry in the \c mixture
specification.
SourceFiles
perfectGasI.H

10
src/thermophysicalModels/specie/equationOfState/rPolynomial/rPolynomial.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2019-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2019-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -27,9 +27,9 @@ Class
Description
Reciprocal polynomial equation of state for liquids and solids
\f[
1/\rho = C_0 + C_1 T + C_2 T^2 - C_3 p - C_4 p T
\f]
\verbatim
1/rho = C0 + C1*T + C2*sqr(T) - C3*p - C4*p*T
\endverbatim
This polynomial for the reciprocal of the density provides a much better fit
than the equivalent polynomial for the density and has the advantage that it
@ -42,7 +42,7 @@ Usage
C | Density polynomial coefficients
\endtable
Example of the specification of the equation of state for pure water:
Example specification of the rPolynomial equation of state for pure water:
\verbatim
equationOfState
{

19
src/thermophysicalModels/specie/equationOfState/rhoConst/rhoConst.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2012-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2012-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -25,7 +25,22 @@ Class
Foam::rhoConst
Description
RhoConst (rho = const) of state.
Constant density equations of state
Usage
\table
Property | Description
rho | Constant density
\endtable
Example specification of the rhoConst equation of state:
\verbatim
equationOfState
{
rho 1000;
}
\endverbatim
SourceFiles
rhoConstI.H