00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035 #ifndef Ar_H
00036 #define Ar_H
00037
00038 #include "liquid.H"
00039 #include "NSRDSfunc0.H"
00040 #include "NSRDSfunc1.H"
00041 #include "NSRDSfunc2.H"
00042 #include "NSRDSfunc3.H"
00043 #include "NSRDSfunc4.H"
00044 #include "NSRDSfunc5.H"
00045 #include "NSRDSfunc6.H"
00046 #include "NSRDSfunc7.H"
00047 #include "NSRDSfunc14.H"
00048 #include "APIdiffCoefFunc.H"
00049
00050
00051
00052 namespace Foam
00053 {
00054
00055
00056
00057
00058
00059 class Ar
00060 :
00061 public liquid
00062 {
00063
00064
00065 NSRDSfunc5 rho_;
00066 NSRDSfunc1 pv_;
00067 NSRDSfunc6 hl_;
00068 NSRDSfunc0 cp_;
00069 NSRDSfunc0 h_;
00070 NSRDSfunc0 cpg_;
00071 NSRDSfunc4 B_;
00072 NSRDSfunc1 mu_;
00073 NSRDSfunc2 mug_;
00074 NSRDSfunc0 K_;
00075 NSRDSfunc2 Kg_;
00076 NSRDSfunc6 sigma_;
00077 APIdiffCoefFunc D_;
00078
00079
00080 public:
00081
00082
00083 TypeName("Ar");
00084
00085
00086
00087
00088
00089 Ar()
00090 :
00091 liquid(39.948, 150.86, 4.8981e+6, 0.07459, 0.291, 83.78, 6.88e+4, 87.28, 0.0, 0.0, 1.4138e+4),
00092 rho_(151.922244, 0.286, 150.86, 0.2984),
00093 pv_(39.233, -1051.7, -3.5895, 5.0444e-05, 2),
00094 hl_(150.86, 218509.061780314, 0.352, 0, 0, 0),
00095 cp_(4562.43116050866, -70.7770101131471, 0.367477721037349, 0, 0, 0),
00096
00097
00098
00099 h_(-1460974.49982473, 4562.43116050866, -35.3885050565735, 0.122492573679116, 0, 0),
00100 cpg_(520.326424351657, 0, 0, 0, 0, 0),
00101 B_(0.000952488234705117, -0.379993992189847, -2022.62941824372, 4633523580654.85, -302893761890458.0),
00102 mu_(-8.868, 204.3, -0.3831, -1.3e-22, 10),
00103 mug_(8.386e-07, 0.6175, 75.377, -432.5),
00104 K_(0.1819, -0.0003176, -4.11e-06, 0, 0, 0),
00105 Kg_(0.0001236, 0.8262, -132.8, 16000),
00106 sigma_(150.86, 0.03823, 1.2927, 0, 0, 0),
00107 D_(147.18, 20.1, 39.948, 28)
00108 {}
00109 Ar
00110 (
00111 const liquid& l,
00112 const NSRDSfunc5& density,
00113 const NSRDSfunc1& vapourPressure,
00114 const NSRDSfunc6& heatOfVapourisation,
00115 const NSRDSfunc0& heatCapacity,
00116 const NSRDSfunc0& enthalpy,
00117 const NSRDSfunc0& idealGasHeatCapacity,
00118 const NSRDSfunc4& secondVirialCoeff,
00119 const NSRDSfunc1& dynamicViscosity,
00120 const NSRDSfunc2& vapourDynamicViscosity,
00121 const NSRDSfunc0& thermalConductivity,
00122 const NSRDSfunc2& vapourThermalConductivity,
00123 const NSRDSfunc6& surfaceTension,
00124 const APIdiffCoefFunc& vapourDiffussivity
00125 )
00126 :
00127 liquid(l),
00128 rho_(density),
00129 pv_(vapourPressure),
00130 hl_(heatOfVapourisation),
00131 cp_(heatCapacity),
00132 h_(enthalpy),
00133 cpg_(idealGasHeatCapacity),
00134 B_(secondVirialCoeff),
00135 mu_(dynamicViscosity),
00136 mug_(vapourDynamicViscosity),
00137 K_(thermalConductivity),
00138 Kg_(vapourThermalConductivity),
00139 sigma_(surfaceTension),
00140 D_(vapourDiffussivity)
00141 {}
00142
00143
00144 Ar(Istream& is)
00145 :
00146 liquid(is),
00147 rho_(is),
00148 pv_(is),
00149 hl_(is),
00150 cp_(is),
00151 h_(is),
00152 cpg_(is),
00153 B_(is),
00154 mu_(is),
00155 mug_(is),
00156 K_(is),
00157 Kg_(is),
00158 sigma_(is),
00159 D_(is)
00160 {}
00161
00162
00163
00164
00165
00166 scalar rho(scalar p, scalar T) const
00167 {
00168 return rho_.f(p, T);
00169 }
00170
00171
00172 scalar pv(scalar p, scalar T) const
00173 {
00174 return pv_.f(p, T);
00175 }
00176
00177
00178 scalar hl(scalar p, scalar T) const
00179 {
00180 return hl_.f(p, T);
00181 }
00182
00183
00184 scalar cp(scalar p, scalar T) const
00185 {
00186 return cp_.f(p, T);
00187 }
00188
00189
00190 scalar h(scalar p, scalar T) const
00191 {
00192 return h_.f(p, T);
00193 }
00194
00195
00196 scalar cpg(scalar p, scalar T) const
00197 {
00198 return cpg_.f(p, T);
00199 }
00200
00201
00202 scalar B(scalar p, scalar T) const
00203 {
00204 return B_.f(p, T);
00205 }
00206
00207
00208 scalar mu(scalar p, scalar T) const
00209 {
00210 return mu_.f(p, T);
00211 }
00212
00213
00214 scalar mug(scalar p, scalar T) const
00215 {
00216 return mug_.f(p, T);
00217 }
00218
00219
00220 scalar K(scalar p, scalar T) const
00221 {
00222 return K_.f(p, T);
00223 }
00224
00225
00226 scalar Kg(scalar p, scalar T) const
00227 {
00228 return Kg_.f(p, T);
00229 }
00230
00231
00232 scalar sigma(scalar p, scalar T) const
00233 {
00234 return sigma_.f(p, T);
00235 }
00236
00237
00238 scalar D(scalar p, scalar T) const
00239 {
00240 return D_.f(p, T);
00241 }
00242
00243
00244
00245 void writeData(Ostream& os) const
00246 {
00247 liquid::writeData(os); os << nl;
00248 rho_.writeData(os); os << nl;
00249 pv_.writeData(os); os << nl;
00250 hl_.writeData(os); os << nl;
00251 cp_.writeData(os); os << nl;
00252 cpg_.writeData(os); os << nl;
00253 B_.writeData(os); os << nl;
00254 mu_.writeData(os); os << nl;
00255 mug_.writeData(os); os << nl;
00256 K_.writeData(os); os << nl;
00257 Kg_.writeData(os); os << nl;
00258 sigma_.writeData(os); os << nl;
00259 D_.writeData(os); os << endl;
00260 }
00261
00262
00263
00264
00265 friend Ostream& operator<<(Ostream& os, const Ar& l)
00266 {
00267 l.writeData(os);
00268 return os;
00269 }
00270 };
00271
00272
00273
00274
00275 }
00276
00277
00278
00279 #endif
00280
00281