1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2008 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 the
13 Free Software Foundation; either version 2 of the License, or (at your
14 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, write to the Free Software Foundation,
23 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
34 \*---------------------------------------------------------------------------*/
40 #include "NSRDSfunc0.H"
41 #include "NSRDSfunc1.H"
42 #include "NSRDSfunc2.H"
43 #include "NSRDSfunc3.H"
44 #include "NSRDSfunc4.H"
45 #include "NSRDSfunc5.H"
46 #include "NSRDSfunc6.H"
47 #include "NSRDSfunc7.H"
48 #include "NSRDSfunc14.H"
49 #include "APIdiffCoefFunc.H"
51 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
56 /*---------------------------------------------------------------------------*\
57 Class bC10H7CH3 Declaration
58 \*---------------------------------------------------------------------------*/
83 //- Runtime type information
92 liquid(142.2, 761.0, 3.25e+6, 0.507, 0.260, 307.73, 1.7374e+1, 514.20, 1.4010e-30, 0.3459, 1.987e+4),
93 rho_(67.36014, 0.23843, 761, 0.2559),
94 pv_(134.31, -12103, -16.195, 6.9659e-06, 2),
95 hl_(761.0, 513150.492264416, 0.4044, 0, 0, 0),
96 cp_(811.322081575246, 2.30225035161744, 0.000862869198312236, 0, 0, 0),
97 // NN: enthalpy, h_, is not used in the sprayModel.
98 // For consistency, the enthalpy is derived from hlat and hl.
99 // It is, however, convenient to have it available.
100 h_(45001.2311880177, 811.322081575246, 1.15112517580872, 0.000287623066104079, 0, 0),
101 cpg_(760.126582278481, 2699.08579465542, 1564.1, 1994.51476793249, 727.49),
102 B_(0.00229430379746835, -3.53720112517581, -1067158.93108298, 2.29746835443038e+18, -2.68438818565401e+21),
103 mu_(-63.276, 4219, 7.5549, 0, 0),
104 mug_(2.1791e-06, 0.3717, 712.53, 0),
105 K_(0.1962, -0.00018414, 0, 0, 0, 0),
106 Kg_(0.4477, -0.1282, -345.89, 2340100),
107 sigma_(761.0, 0.066442, 1.2634, 0, 0, 0),
108 D_(147.18, 20.1, 142.2, 28) // NN: Same as nHeptane
113 const NSRDSfunc5& density,
114 const NSRDSfunc1& vapourPressure,
115 const NSRDSfunc6& heatOfVapourisation,
116 const NSRDSfunc0& heatCapacity,
117 const NSRDSfunc0& enthalpy,
118 const NSRDSfunc7& idealGasHeatCapacity,
119 const NSRDSfunc4& secondVirialCoeff,
120 const NSRDSfunc1& dynamicViscosity,
121 const NSRDSfunc2& vapourDynamicViscosity,
122 const NSRDSfunc0& thermalConductivity,
123 const NSRDSfunc2& vapourThermalConductivity,
124 const NSRDSfunc6& surfaceTension,
125 const APIdiffCoefFunc& vapourDiffussivity
131 hl_(heatOfVapourisation),
134 cpg_(idealGasHeatCapacity),
135 B_(secondVirialCoeff),
136 mu_(dynamicViscosity),
137 mug_(vapourDynamicViscosity),
138 K_(thermalConductivity),
139 Kg_(vapourThermalConductivity),
140 sigma_(surfaceTension),
141 D_(vapourDiffussivity)
144 //- Construct from Istream
145 bC10H7CH3(Istream& is)
166 //- Liquid density [kg/m^3]
167 scalar rho(scalar p, scalar T) const
172 //- Vapour pressure [Pa]
173 scalar pv(scalar p, scalar T) const
178 //- Heat of vapourisation [J/kg]
179 scalar hl(scalar p, scalar T) const
184 //- Liquid heat capacity [J/(kg K)]
185 scalar cp(scalar p, scalar T) const
190 //- Liquid Enthalpy [J/(kg)]
191 scalar h(scalar p, scalar T) const
196 //- Ideal gas heat capacity [J/(kg K)]
197 scalar cpg(scalar p, scalar T) const
202 //- Second Virial Coefficient [m^3/kg]
203 scalar B(scalar p, scalar T) const
208 //- Liquid viscosity [Pa s]
209 scalar mu(scalar p, scalar T) const
214 //- Vapour viscosity [Pa s]
215 scalar mug(scalar p, scalar T) const
220 //- Liquid thermal conductivity [W/(m K)]
221 scalar K(scalar p, scalar T) const
226 //- Vapour thermal conductivity [W/(m K)]
227 scalar Kg(scalar p, scalar T) const
232 //- Surface tension [N/m]
233 scalar sigma(scalar p, scalar T) const
235 return sigma_.f(p, T);
238 //- Vapour diffussivity [m2/s]
239 scalar D(scalar p, scalar T) const
245 //- Write the function coefficients
246 void writeData(Ostream& os) const
248 liquid::writeData(os); os << nl;
249 rho_.writeData(os); os << nl;
250 pv_.writeData(os); os << nl;
251 hl_.writeData(os); os << nl;
252 cp_.writeData(os); os << nl;
253 cpg_.writeData(os); os << nl;
254 B_.writeData(os); os << nl;
255 mu_.writeData(os); os << nl;
256 mug_.writeData(os); os << nl;
257 K_.writeData(os); os << nl;
258 Kg_.writeData(os); os << nl;
259 sigma_.writeData(os); os << nl;
260 D_.writeData(os); os << endl;
266 friend Ostream& operator<<(Ostream& os, const bC10H7CH3& l)
274 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
276 } // End namespace Foam
278 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
282 // ************************************************************************* //