/usr/include/freefoam/liquids/IDEA.H is in libfreefoam-dev 0.1.0+dfsg-1build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 | /*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
Class
Foam::IDEA
Description
The IDEA fuel is constructed by adding 30% alphaMethylNaphthalene
with 70% n-decane.
The new properties have been calculated by adding the values in these
proportions and making a least square fit, using the same NSRDS-eq.
so that Y = 0.3*Y_naphthalene + 0.7*Y_decane
The valid Temperature range for n-decane is normally 243.51 - 617.70 K
and for the naphthalene it is 242.67 - 772.04 K
The least square fit was done in the interval 244 - 617 K
The critical temperature was taken to be 618.074 K, since this
is the 'c'-value in the rho-equation, which corresponds to Tcrit,
This value was then used in the fit for the NSRDS6-eq, which uses Tcrit.
(important for the latent heat and surface tension)
The molecular weights are 142.20 and 142.285 and for the IDEA fuel
it is thus 142.26 ( approximately 0.3*142.2 + 0.7*142.285 )
Critical pressure was set to the lowest one (n-Decane)
Critical volume... also the lowest one (naphthalene) 0.523 m^3/kmol
Second Virial Coefficient is n-Decane
SourceFiles
IDEA.C
\*---------------------------------------------------------------------------*/
#ifndef IDEA_H
#define IDEA_H
#include <liquids/liquid.H>
#include <thermophysicalFunctions/NSRDSfunc0.H>
#include <thermophysicalFunctions/NSRDSfunc1.H>
#include <thermophysicalFunctions/NSRDSfunc2.H>
#include <thermophysicalFunctions/NSRDSfunc3.H>
#include <thermophysicalFunctions/NSRDSfunc4.H>
#include <thermophysicalFunctions/NSRDSfunc5.H>
#include <thermophysicalFunctions/NSRDSfunc6.H>
#include <thermophysicalFunctions/NSRDSfunc7.H>
#include <thermophysicalFunctions/APIdiffCoefFunc.H>
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
namespace Foam
{
/*---------------------------------------------------------------------------*\
Class IDEA Declaration
\*---------------------------------------------------------------------------*/
class IDEA
:
public liquid
{
// Private data
NSRDSfunc5 rho_;
NSRDSfunc1 pv_;
NSRDSfunc6 hl_;
NSRDSfunc0 cp_;
NSRDSfunc0 h_;
NSRDSfunc7 cpg_;
NSRDSfunc4 B_;
NSRDSfunc1 mu_;
NSRDSfunc2 mug_;
NSRDSfunc0 K_;
NSRDSfunc2 Kg_;
NSRDSfunc6 sigma_;
APIdiffCoefFunc D_;
public:
//- Runtime type information
TypeName("IDEA");
// Constructors
//- Construct null
IDEA();
// Construct from components
IDEA
(
const liquid& l,
const NSRDSfunc5& density,
const NSRDSfunc1& vapourPressure,
const NSRDSfunc6& heatOfVapourisation,
const NSRDSfunc0& heatCapacity,
const NSRDSfunc0& enthalpy,
const NSRDSfunc7& idealGasHeatCapacity,
const NSRDSfunc4& secondVirialCoeff,
const NSRDSfunc1& dynamicViscosity,
const NSRDSfunc2& vapourDynamicViscosity,
const NSRDSfunc0& thermalConductivity,
const NSRDSfunc2& vapourThermalConductivity,
const NSRDSfunc6& surfaceTension,
const APIdiffCoefFunc& vapourDiffussivity
);
//- Construct from Istream
IDEA(Istream& is);
// Member Functions
//- Liquid density [kg/m^3]
inline scalar rho(scalar p, scalar T) const;
//- Vapour pressure [Pa]
inline scalar pv(scalar p, scalar T) const;
//- Heat of vapourisation [J/kg]
inline scalar hl(scalar p, scalar T) const;
//- Liquid heat capacity [J/(kg K)]
inline scalar cp(scalar p, scalar T) const;
//- Liquid Enthalpy [J/(kg)]
inline scalar h(scalar p, scalar T) const;
//- Ideal gas heat capacity [J/(kg K)]
inline scalar cpg(scalar p, scalar T) const;
//- Second Virial Coefficient [m^3/kg]
inline scalar B(scalar p, scalar T) const;
//- Liquid viscosity [Pa s]
inline scalar mu(scalar p, scalar T) const;
//- Vapour viscosity [Pa s]
inline scalar mug(scalar p, scalar T) const;
//- Liquid thermal conductivity [W/(m K)]
inline scalar K(scalar p, scalar T) const;
//- Vapour thermal conductivity [W/(m K)]
inline scalar Kg(scalar p, scalar T) const;
//- Surface tension [N/m]
inline scalar sigma(scalar p, scalar T) const;
//- Vapour diffussivity [m2/s]
inline scalar D(scalar p, scalar T) const;
//- Vapour diffussivity [m2/s] with specified binary pair
inline scalar D(scalar p, scalar T, scalar Wb) const;
// I-O
//- Write the function coefficients
void writeData(Ostream& os) const
{
liquid::writeData(os); os << nl;
rho_.writeData(os); os << nl;
pv_.writeData(os); os << nl;
hl_.writeData(os); os << nl;
cp_.writeData(os); os << nl;
cpg_.writeData(os); os << nl;
mu_.writeData(os); os << nl;
mug_.writeData(os); os << nl;
K_.writeData(os); os << nl;
Kg_.writeData(os); os << nl;
sigma_.writeData(os); os << nl;
D_.writeData(os); os << endl;
}
//- Ostream Operator
friend Ostream& operator<<(Ostream& os, const IDEA& l)
{
l.writeData(os);
return os;
}
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace Foam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#include "IDEAI.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#endif
// ************************ vim: set sw=4 sts=4 et: ************************ //
|