27 #ifndef OPM_XYLENE_HPP
28 #define OPM_XYLENE_HPP
45 template <
class Scalar>
49 typedef ::Opm::IdealGas<Scalar>
IdealGas;
80 {
throw std::runtime_error(
"Not implemented: tripleTemperature for xylene"); }
86 {
throw std::runtime_error(
"Not implemented: triplePressure for xylene"); }
94 template <
class Evaluation>
97 const Scalar A = 7.00909;;
98 const Scalar B = 1462.266;;
99 const Scalar C = 215.110;;
101 return 100*1.334*pow(10.0, (A - (B/(temperature - 273.15 + C))));
112 template <
class Evaluation>
115 Evaluation CH3,C6H5,H;
120 if(temperature < 298.0){
121 H = 13.4 + 1.2*(temperature - 273.0)/25.0;
122 CH3 = 40.0 + 1.6*(temperature - 273.0)/25.0;
123 C6H5 = 113.0 + 4.2*(temperature - 273.0)/25.0;
125 else if(temperature < 323.0){
126 H = 14.6 + 0.9*(temperature - 298.0)/25.0;
127 CH3 = 41.6 + 1.9*(temperature - 298.0)/25.0;
128 C6H5 = 117.2 + 6.2*(temperature - 298.0)/25.0;
130 else if(temperature < 348.0){
131 H = 15.5 + 1.2*(temperature - 323.0)/25.0;
132 CH3 = 43.5 + 2.3*(temperature - 323.0)/25.0;
133 C6H5 = 123.4 + 6.3*(temperature - 323.0)/25.0;
136 H = 16.7 + 2.1*(temperature - 348.0)/25.0;
137 CH3 = 45.8 + 2.5*(temperature - 348.0)/25.0;
138 C6H5 = 129.7 + 6.3*(temperature - 348.0)/25.0;
148 template <
class Evaluation>
149 static Evaluation
liquidEnthalpy(
const Evaluation& temperature,
const Evaluation& pressure)
180 template <
class Evaluation>
181 static Evaluation
heatVap(Evaluation temperature,
185 temperature = max(temperature, 0.0);
192 const Scalar DH_v_boil =
Consts::R * T_crit * Tr1
193 * (3.978 * Tr1 - 3.958 + 1.555*std::log(p_crit * 1e-5 ) )
198 const Scalar n = 0.375;
199 const Evaluation& DH_vap = DH_v_boil * pow(((1.0 - Tr2)/(1.0 - Tr1)), n);
210 template <
class Evaluation>
211 static Evaluation
gasEnthalpy(
const Evaluation& temperature,
const Evaluation& pressure)
219 template <
class Evaluation>
220 static Evaluation
gasDensity(
const Evaluation& temperature,
const Evaluation& pressure)
231 template <
class Evaluation>
232 static Evaluation
molarGasDensity(
const Evaluation& temperature,
const Evaluation& pressure)
246 template <
class Evaluation>
252 temperature = min(temperature, 500.0);
253 temperature = max(temperature, 250.0);
255 const Scalar A1 = 0.25919;
256 const Scalar A2 = 0.0014569;
258 return 1.0/(A2*pow(A1, expo));
264 template <
class Evaluation>
265 static Evaluation
liquidDensity(
const Evaluation& temperature,
const Evaluation& pressure)
289 template <
class Evaluation>
290 static Evaluation
gasViscosity(Evaluation temperature,
const Evaluation& )
292 temperature = min(temperature, 500.0);
293 temperature = max(temperature, 250.0);
296 const Scalar Fp0 = 1.0;
297 const Scalar xi = 0.004623;
298 const Evaluation& eta_xi = Fp0*(0.807*pow(Tr, 0.618)
299 - 0.357*exp(-0.449*Tr)
300 + 0.34*exp(-4.058*Tr)
302 return eta_xi/xi / 1e7;
308 template <
class Evaluation>
311 temperature = min(temperature, 500.0);
312 temperature = max(temperature, 250.0);
314 const Scalar A = -3.82;
315 const Scalar B = 1027.0;
316 const Scalar C = -6.38e-4;
317 const Scalar D = 4.52e-7;
322 + D*temperature*temperature);
Abstract base class of a pure chemical species.
A central place for various physical constants occuring in some equations.
Relations valid for an ideal gas.
Abstract base class of a pure chemical species.
Definition: Component.hpp:42
A central place for various physical constants occuring in some equations.
Definition: Constants.hpp:41
static const Scalar R
The ideal gas constant [J/(mol K)].
Definition: Constants.hpp:45
Relations valid for an ideal gas.
Definition: IdealGas.hpp:38
static Evaluation density(const Evaluation &avgMolarMass, const Evaluation &temperature, const Evaluation &pressure)
The density of the gas in , depending on pressure, temperature and average molar mass of the gas.
Definition: IdealGas.hpp:48
Component for Xylene.
Definition: Xylene.hpp:47
static Evaluation liquidViscosity(Evaluation temperature, const Evaluation &)
The dynamic liquid viscosity of the pure component.
Definition: Xylene.hpp:309
static Evaluation gasDensity(const Evaluation &temperature, const Evaluation &pressure)
The density in of the component at a given pressure in and temperature in .
Definition: Xylene.hpp:220
static Evaluation vaporPressure(const Evaluation &temperature)
The saturation vapor pressure in of pure xylene at a given temperature according to Antoine after Be...
Definition: Xylene.hpp:95
static Evaluation gasEnthalpy(const Evaluation &temperature, const Evaluation &pressure)
Specific enthalpy of the pure component in gas.
Definition: Xylene.hpp:211
static const char * name()
A human readable name for the xylene.
Definition: Xylene.hpp:55
static bool liquidIsCompressible()
Returns true iff the liquid phase is assumed to be compressible.
Definition: Xylene.hpp:283
static Evaluation molarLiquidDensity(Evaluation temperature, const Evaluation &)
The molar density of pure xylene at a given pressure and temperature .
Definition: Xylene.hpp:247
static Evaluation spHeatCapLiquidPhase(const Evaluation &temperature, const Evaluation &)
Specific heat cap of liquid xylene .
Definition: Xylene.hpp:113
static Scalar criticalTemperature()
Returns the critical temperature of xylene.
Definition: Xylene.hpp:67
static Evaluation liquidEnthalpy(const Evaluation &temperature, const Evaluation &pressure)
Specific enthalpy of the pure component in liquid.
Definition: Xylene.hpp:149
static Scalar tripleTemperature()
Returns the temperature at xylene's triple point.
Definition: Xylene.hpp:79
static Scalar molarMass()
The molar mass in of xylene.
Definition: Xylene.hpp:61
static Evaluation heatVap(Evaluation temperature, const Evaluation &)
Latent heat of vaporization for xylene .
Definition: Xylene.hpp:181
static Evaluation gasViscosity(Evaluation temperature, const Evaluation &)
The dynamic viscosity of the pure component at a given pressure in and temperature in .
Definition: Xylene.hpp:290
static Evaluation molarGasDensity(const Evaluation &temperature, const Evaluation &pressure)
The density of xylene gas at a given pressure and temperature.
Definition: Xylene.hpp:232
static Scalar triplePressure()
Returns the pressure at xylene's triple point.
Definition: Xylene.hpp:85
static Scalar boilingTemperature()
Returns the temperature at xylene's boiling point (1 atm).
Definition: Xylene.hpp:169
static Evaluation liquidDensity(const Evaluation &temperature, const Evaluation &pressure)
The density of the liquid component at a given pressure in and temperature in .
Definition: Xylene.hpp:265
static bool gasIsIdeal()
Returns true iff the gas phase is assumed to be ideal.
Definition: Xylene.hpp:277
static bool gasIsCompressible()
Returns true iff the gas phase is assumed to be compressible.
Definition: Xylene.hpp:271
static Scalar criticalPressure()
Returns the critical pressure of xylene.
Definition: Xylene.hpp:73