27#ifndef OPM_MESITYLENE_HPP
28#define OPM_MESITYLENE_HPP
45template <
class Scalar>
54 static std::string_view
name()
55 {
return "mesitylene"; }
85 {
throw std::runtime_error(
"Not implemented: tripleTemperature for mesitylene"); }
91 {
throw std::runtime_error(
"Not implemented: triplePressure for mesitylene"); }
100 template <
class Evaluation>
103 const Scalar A = 7.07638;
104 const Scalar B = 1571.005;
105 const Scalar C = 209.728;
107 const Evaluation& T = temperature - 273.15;
109 return 100 * 1.334 * pow(10.0, A - (B / (T + C)));
119 template <
class Evaluation>
120 static Evaluation
liquidEnthalpy(
const Evaluation& temperature,
const Evaluation& pressure)
145 template <
class Evaluation>
146 static Evaluation
heatVap(
const Evaluation& temperature,
const Evaluation& )
156 const Scalar DH_v_boil =
158 * (3.978 * Tr1 - 3.958 + 1.555*std::log(p_crit * 1e-5 ) )
163 const Scalar n = 0.375;
164 const Evaluation& DH_vap = DH_v_boil * pow(((1.0 - Tr2)/(1.0 - Tr1)), n);
179 template <
class Evaluation>
180 static Evaluation
gasEnthalpy(
const Evaluation& temperature,
const Evaluation& pressure)
191 template <
class Evaluation>
192 static Evaluation
gasDensity(
const Evaluation& temperature,
const Evaluation& pressure)
201 template <
class Evaluation>
202 static Evaluation
liquidDensity(
const Evaluation& temperature,
const Evaluation& )
230 template <
class Evaluation>
231 static Evaluation
gasViscosity(Evaluation temperature,
const Evaluation& ,
bool =
true)
233 temperature = min(temperature, 500.0);
234 temperature = max(temperature, 250.0);
241 const Evaluation& eta_xi =
242 Fp0*(0.807*pow(Tr,0.618)
243 - 0.357*exp(-0.449*Tr)
244 + 0.34*exp(-4.058*Tr)
247 return eta_xi/xi/1e7;
256 template <
class Evaluation>
259 temperature = min(temperature, 500.0);
260 temperature = max(temperature, 250.0);
262 const Scalar A = -6.749;
263 const Scalar B = 2010.0;
265 return exp(A + B/temperature)*1e-3;
277 template <
class Evaluation>
284 Evaluation H, CH3, C6H5;
285 if(temperature<298.) {
287 H = 13.4 + 1.2*(temperature-273.0)/25.;
288 CH3 = 40.0 + 1.6*(temperature-273.0)/25.;
289 C6H5 = 113.0 + 4.2*(temperature-273.0)/25.;
291 else if((temperature>=298.0)&&(temperature<323.)){
292 H = 14.6 + 0.9*(temperature-298.0)/25.;
293 CH3 = 41.6 + 1.9*(temperature-298.0)/25.;
294 C6H5 = 117.2 + 6.2*(temperature-298.0)/25.;
296 else if((temperature>=323.0)&&(temperature<348.)){
297 H = 15.5 + 1.2*(temperature-323.0)/25.;
298 CH3 = 43.5 + 2.3*(temperature-323.0)/25.;
299 C6H5 = 123.4 + 6.3*(temperature-323.0)/25.;
302 assert(temperature>=348.0);
305 H = 16.7+2.1*(temperature-348.0)/25.;
306 CH3 = 45.8+2.5*(temperature-348.0)/25.;
307 C6H5 = 129.7+6.3*(temperature-348.0)/25.;
322 template <
class Evaluation>
325 temperature = min(temperature, 500.0);
326 temperature = max(temperature, 250.0);
328 const Scalar Z_RA = 0.2556;
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:44
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
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 Mesitylene.
Definition Mesitylene.hpp:47
static bool liquidIsCompressible()
Returns true iff the liquid phase is assumed to be compressible.
Definition Mesitylene.hpp:220
static Evaluation liquidHeatCapacity(const Evaluation &temperature, const Evaluation &)
Specific heat cap of liquid mesitylene .
Definition Mesitylene.hpp:278
static Scalar triplePressure()
Returns the pressure at mesitylene's triple point.
Definition Mesitylene.hpp:90
static Evaluation vaporPressure(const Evaluation &temperature)
The saturation vapor pressure in of pure mesitylene at a given temperature according to Antoine afte...
Definition Mesitylene.hpp:101
static Evaluation liquidEnthalpy(const Evaluation &temperature, const Evaluation &pressure)
Specific enthalpy of liquid mesitylene .
Definition Mesitylene.hpp:120
static Scalar molarMass()
The molar mass in of mesitylene.
Definition Mesitylene.hpp:60
static Scalar criticalTemperature()
Returns the critical temperature of mesitylene.
Definition Mesitylene.hpp:66
static Evaluation gasDensity(const Evaluation &temperature, const Evaluation &pressure)
The density of pure mesitylene vapor at a given pressure and temperature .
Definition Mesitylene.hpp:192
static Evaluation heatVap(const Evaluation &temperature, const Evaluation &)
Latent heat of vaporization for mesitylene .
Definition Mesitylene.hpp:146
static bool gasIsIdeal()
Returns true iff the gas phase is assumed to be ideal.
Definition Mesitylene.hpp:214
static Evaluation gasEnthalpy(const Evaluation &temperature, const Evaluation &pressure)
Specific enthalpy of mesitylene vapor .
Definition Mesitylene.hpp:180
static Scalar tripleTemperature()
Returns the temperature at mesitylene's triple point.
Definition Mesitylene.hpp:84
static bool gasIsCompressible()
Returns true iff the gas phase is assumed to be compressible.
Definition Mesitylene.hpp:208
static Evaluation liquidDensity(const Evaluation &temperature, const Evaluation &)
The density of pure mesitylene at a given pressure and temperature .
Definition Mesitylene.hpp:202
static Evaluation liquidViscosity(Evaluation temperature, const Evaluation &)
The dynamic viscosity of pure mesitylene.
Definition Mesitylene.hpp:257
static Evaluation molarLiquidDensity_(Evaluation temperature)
The molar density of pure mesitylene at a given pressure and temperature .
Definition Mesitylene.hpp:323
static Scalar boilingTemperature()
Returns the temperature at mesitylene's boiling point (1 atm).
Definition Mesitylene.hpp:78
static Scalar criticalPressure()
Returns the critical pressure of mesitylene.
Definition Mesitylene.hpp:72
static Evaluation gasViscosity(Evaluation temperature, const Evaluation &, bool=true)
The dynamic viscosity of mesitylene vapor.
Definition Mesitylene.hpp:231
static std::string_view name()
A human readable name for the mesitylene.
Definition Mesitylene.hpp:54
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30