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
00036 #ifndef C6H6_H
00037 #define C6H6_H
00038
00039 #include "liquid.H"
00040 #include "NSRDSfunc0.H"
00041 #include "NSRDSfunc1.H"
00042 #include "NSRDSfunc2.H"
00043 #include "NSRDSfunc3.H"
00044 #include "NSRDSfunc4.H"
00045 #include "NSRDSfunc5.H"
00046 #include "NSRDSfunc6.H"
00047 #include "NSRDSfunc7.H"
00048 #include "NSRDSfunc14.H"
00049 #include "APIdiffCoefFunc.H"
00050
00051
00052
00053 namespace Foam
00054 {
00055
00056
00057
00058
00059
00060 class C6H6
00061 :
00062 public liquid
00063 {
00064
00065
00066 NSRDSfunc5 rho_;
00067 NSRDSfunc1 pv_;
00068 NSRDSfunc6 hl_;
00069 NSRDSfunc0 cp_;
00070 NSRDSfunc0 h_;
00071 NSRDSfunc7 cpg_;
00072 NSRDSfunc4 B_;
00073 NSRDSfunc1 mu_;
00074 NSRDSfunc2 mug_;
00075 NSRDSfunc0 K_;
00076 NSRDSfunc2 Kg_;
00077 NSRDSfunc6 sigma_;
00078 APIdiffCoefFunc D_;
00079
00080
00081 public:
00082
00083
00084 TypeName("C6H6");
00085
00086
00087
00088
00089
00090 C6H6()
00091 :
00092 liquid(78.114, 562.16, 4.898e+6, 0.25894, 0.271, 278.68, 4.7961e+3, 353.24, 0.0, 0.2108, 1.8706e+4),
00093 rho_(80.5511568, 0.2667, 562.16, 0.2818),
00094 pv_(78.05, -6275.5, -8.4443, 6.26e-06, 2),
00095 hl_(562.16, 649435.440510024, 0.7616, -0.5052, 0.1564, 0),
00096 cp_(1386.69124612745, -0.416058581048212, 0.00542796425736744, 0, 0, 0),
00097
00098
00099
00100 h_(186141.395065592, 1386.69124612745, -0.208029290524106, 0.00180932141912248, 0, 0),
00101 cpg_(568.656066774202, 2970.65826868423, 1494.6, 2203.57426325627, -678.15),
00102 B_(0.00184089919860716, -2.30176408838364, -309176.332027549, -5.12072099751645e+15, -2.90216862534245e+19),
00103 mu_(6.764, 336.4, -2.687, 0, 0),
00104 mug_(3.134e-08, 0.9676, 7.9, 0),
00105 K_(0.2407, -0.0003202, 0, 0, 0, 0),
00106 Kg_(1.652e-05, 1.3117, 491, 0),
00107 sigma_(562.16, 0.07195, 1.2389, 0, 0, 0),
00108 D_(147.18, 20.1, 78.114, 28)
00109 {}
00110 C6H6
00111 (
00112 const liquid& l,
00113 const NSRDSfunc5& density,
00114 const NSRDSfunc1& vapourPressure,
00115 const NSRDSfunc6& heatOfVapourisation,
00116 const NSRDSfunc0& heatCapacity,
00117 const NSRDSfunc0& enthalpy,
00118 const NSRDSfunc7& idealGasHeatCapacity,
00119 const NSRDSfunc4& secondVirialCoeff,
00120 const NSRDSfunc1& dynamicViscosity,
00121 const NSRDSfunc2& vapourDynamicViscosity,
00122 const NSRDSfunc0& thermalConductivity,
00123 const NSRDSfunc2& vapourThermalConductivity,
00124 const NSRDSfunc6& surfaceTension,
00125 const APIdiffCoefFunc& vapourDiffussivity
00126 )
00127 :
00128 liquid(l),
00129 rho_(density),
00130 pv_(vapourPressure),
00131 hl_(heatOfVapourisation),
00132 cp_(heatCapacity),
00133 h_(enthalpy),
00134 cpg_(idealGasHeatCapacity),
00135 B_(secondVirialCoeff),
00136 mu_(dynamicViscosity),
00137 mug_(vapourDynamicViscosity),
00138 K_(thermalConductivity),
00139 Kg_(vapourThermalConductivity),
00140 sigma_(surfaceTension),
00141 D_(vapourDiffussivity)
00142 {}
00143
00144
00145 C6H6(Istream& is)
00146 :
00147 liquid(is),
00148 rho_(is),
00149 pv_(is),
00150 hl_(is),
00151 cp_(is),
00152 h_(is),
00153 cpg_(is),
00154 B_(is),
00155 mu_(is),
00156 mug_(is),
00157 K_(is),
00158 Kg_(is),
00159 sigma_(is),
00160 D_(is)
00161 {}
00162
00163
00164
00165
00166
00167 scalar rho(scalar p, scalar T) const
00168 {
00169 return rho_.f(p, T);
00170 }
00171
00172
00173 scalar pv(scalar p, scalar T) const
00174 {
00175 return pv_.f(p, T);
00176 }
00177
00178
00179 scalar hl(scalar p, scalar T) const
00180 {
00181 return hl_.f(p, T);
00182 }
00183
00184
00185 scalar cp(scalar p, scalar T) const
00186 {
00187 return cp_.f(p, T);
00188 }
00189
00190
00191 scalar h(scalar p, scalar T) const
00192 {
00193 return h_.f(p, T);
00194 }
00195
00196
00197 scalar cpg(scalar p, scalar T) const
00198 {
00199 return cpg_.f(p, T);
00200 }
00201
00202
00203 scalar B(scalar p, scalar T) const
00204 {
00205 return B_.f(p, T);
00206 }
00207
00208
00209 scalar mu(scalar p, scalar T) const
00210 {
00211 return mu_.f(p, T);
00212 }
00213
00214
00215 scalar mug(scalar p, scalar T) const
00216 {
00217 return mug_.f(p, T);
00218 }
00219
00220
00221 scalar K(scalar p, scalar T) const
00222 {
00223 return K_.f(p, T);
00224 }
00225
00226
00227 scalar Kg(scalar p, scalar T) const
00228 {
00229 return Kg_.f(p, T);
00230 }
00231
00232
00233 scalar sigma(scalar p, scalar T) const
00234 {
00235 return sigma_.f(p, T);
00236 }
00237
00238
00239 scalar D(scalar p, scalar T) const
00240 {
00241 return D_.f(p, T);
00242 }
00243
00244
00245
00246 void writeData(Ostream& os) const
00247 {
00248 liquid::writeData(os); os << nl;
00249 rho_.writeData(os); os << nl;
00250 pv_.writeData(os); os << nl;
00251 hl_.writeData(os); os << nl;
00252 cp_.writeData(os); os << nl;
00253 cpg_.writeData(os); os << nl;
00254 B_.writeData(os); os << nl;
00255 mu_.writeData(os); os << nl;
00256 mug_.writeData(os); os << nl;
00257 K_.writeData(os); os << nl;
00258 Kg_.writeData(os); os << nl;
00259 sigma_.writeData(os); os << nl;
00260 D_.writeData(os); os << endl;
00261 }
00262
00263
00264
00265
00266 friend Ostream& operator<<(Ostream& os, const C6H6& l)
00267 {
00268 l.writeData(os);
00269 return os;
00270 }
00271 };
00272
00273
00274
00275
00276 }
00277
00278
00279
00280 #endif
00281
00282