1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 2004-2011 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
11 OpenFOAM is free software: you can redistribute it and/or modify it
12 under the terms of the GNU General Public License as published by
13 the Free Software Foundation, either version 3 of the License, or
14 (at your option) any later version.
16 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 You should have received a copy of the GNU General Public License
22 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
33 \*---------------------------------------------------------------------------*/
38 #include "liquidProperties.H"
39 #include "NSRDSfunc0.H"
40 #include "NSRDSfunc1.H"
41 #include "NSRDSfunc2.H"
42 #include "NSRDSfunc3.H"
43 #include "NSRDSfunc4.H"
44 #include "NSRDSfunc5.H"
45 #include "NSRDSfunc6.H"
46 #include "NSRDSfunc7.H"
47 #include "NSRDSfunc14.H"
48 #include "APIdiffCoefFunc.H"
50 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
55 /*---------------------------------------------------------------------------*\
56 Class C2H5OH Declaration
57 \*---------------------------------------------------------------------------*/
61 public liquidProperties
82 //- Runtime type information
91 //- Construct from components
94 const liquidProperties& l,
95 const NSRDSfunc5& density,
96 const NSRDSfunc1& vapourPressure,
97 const NSRDSfunc6& heatOfVapourisation,
98 const NSRDSfunc0& heatCapacity,
99 const NSRDSfunc0& enthalpy,
100 const NSRDSfunc7& idealGasHeatCapacity,
101 const NSRDSfunc4& secondVirialCoeff,
102 const NSRDSfunc1& dynamicViscosity,
103 const NSRDSfunc2& vapourDynamicViscosity,
104 const NSRDSfunc0& thermalConductivity,
105 const NSRDSfunc2& vapourThermalConductivity,
106 const NSRDSfunc6& surfaceTension,
107 const APIdiffCoefFunc& vapourDiffussivity
110 //- Construct from Istream
113 //- Construct from dictionary
114 C2H5OH(const dictionary& dict);
117 C2H5OH(const C2H5OH& liq);
119 //- Construct and return clone
120 virtual autoPtr<liquidProperties> clone() const
122 return autoPtr<liquidProperties>(new C2H5OH(*this));
128 //- Liquid density [kg/m^3]
129 inline scalar rho(scalar p, scalar T) const;
131 //- Vapour pressure [Pa]
132 inline scalar pv(scalar p, scalar T) const;
134 //- Heat of vapourisation [J/kg]
135 inline scalar hl(scalar p, scalar T) const;
137 //- Liquid heat capacity [J/(kg K)]
138 inline scalar Cp(scalar p, scalar T) const;
140 //- Liquid Enthalpy [J/(kg)]
141 inline scalar h(scalar p, scalar T) const;
143 //- Ideal gas heat capacity [J/(kg K)]
144 inline scalar Cpg(scalar p, scalar T) const;
146 //- Second Virial Coefficient [m^3/kg]
147 inline scalar B(scalar p, scalar T) const;
149 //- Liquid viscosity [Pa s]
150 inline scalar mu(scalar p, scalar T) const;
152 //- Vapour viscosity [Pa s]
153 inline scalar mug(scalar p, scalar T) const;
155 //- Liquid thermal conductivity [W/(m K)]
156 inline scalar K(scalar p, scalar T) const;
158 //- Vapour thermal conductivity [W/(m K)]
159 inline scalar Kg(scalar p, scalar T) const;
161 //- Surface tension [N/m]
162 inline scalar sigma(scalar p, scalar T) const;
164 //- Vapour diffussivity [m2/s]
165 inline scalar D(scalar p, scalar T) const;
167 //- Vapour diffussivity [m2/s] with specified binary pair
168 inline scalar D(scalar p, scalar T, scalar Wb) const;
173 //- Write the function coefficients
174 void writeData(Ostream& os) const
176 liquidProperties::writeData(os); os << nl;
177 rho_.writeData(os); os << nl;
178 pv_.writeData(os); os << nl;
179 hl_.writeData(os); os << nl;
180 Cp_.writeData(os); os << nl;
181 Cpg_.writeData(os); os << nl;
182 B_.writeData(os); os << nl;
183 mu_.writeData(os); os << nl;
184 mug_.writeData(os); os << nl;
185 K_.writeData(os); os << nl;
186 Kg_.writeData(os); os << nl;
187 sigma_.writeData(os); os << nl;
188 D_.writeData(os); os << endl;
192 friend Ostream& operator<<(Ostream& os, const C2H5OH& l)
200 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
202 } // End namespace Foam
204 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
208 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
212 // ************************************************************************* //