27 #ifndef OPM_MESITYLENE_HPP
28 #define OPM_MESITYLENE_HPP
43 template <
class Scalar>
53 {
return "mesitylene"; }
83 {
throw std::runtime_error(
"Not implemented: tripleTemperature for mesitylene"); }
89 {
throw std::runtime_error(
"Not implemented: triplePressure for mesitylene"); }
98 template <
class Evaluation>
101 const Scalar A = 7.07638;
102 const Scalar B = 1571.005;
103 const Scalar C = 209.728;
105 const Evaluation& T = temperature - 273.15;
107 return 100 * 1.334 * pow(10.0, A - (B / (T + C)));
117 template <
class Evaluation>
118 static Evaluation
liquidEnthalpy(
const Evaluation& temperature,
const Evaluation& pressure)
143 template <
class Evaluation>
144 static Evaluation
heatVap(
const Evaluation& temperature,
const Evaluation& )
154 const Scalar DH_v_boil =
156 * (3.978 * Tr1 - 3.958 + 1.555*std::log(p_crit * 1e-5 ) )
161 const Scalar n = 0.375;
162 const Evaluation& DH_vap = DH_v_boil * pow(((1.0 - Tr2)/(1.0 - Tr1)), n);
177 template <
class Evaluation>
178 static Evaluation
gasEnthalpy(
const Evaluation& temperature,
const Evaluation& pressure)
189 template <
class Evaluation>
190 static Evaluation
gasDensity(
const Evaluation& temperature,
const Evaluation& pressure)
199 template <
class Evaluation>
200 static Evaluation
liquidDensity(
const Evaluation& temperature,
const Evaluation& )
228 template <
class Evaluation>
229 static Evaluation
gasViscosity(Evaluation temperature,
const Evaluation& ,
bool =
true)
231 temperature = min(temperature, 500.0);
232 temperature = max(temperature, 250.0);
239 const Evaluation& eta_xi =
240 Fp0*(0.807*pow(Tr,0.618)
241 - 0.357*exp(-0.449*Tr)
242 + 0.34*exp(-4.058*Tr)
245 return eta_xi/xi/1e7;
254 template <
class Evaluation>
257 temperature = min(temperature, 500.0);
258 temperature = max(temperature, 250.0);
260 const Scalar A = -6.749;
261 const Scalar B = 2010.0;
263 return exp(A + B/temperature)*1e-3;
275 template <
class Evaluation>
282 Evaluation H, CH3, C6H5;
283 if(temperature<298.) {
285 H = 13.4 + 1.2*(temperature-273.0)/25.;
286 CH3 = 40.0 + 1.6*(temperature-273.0)/25.;
287 C6H5 = 113.0 + 4.2*(temperature-273.0)/25.;
289 else if((temperature>=298.0)&&(temperature<323.)){
290 H = 14.6 + 0.9*(temperature-298.0)/25.;
291 CH3 = 41.6 + 1.9*(temperature-298.0)/25.;
292 C6H5 = 117.2 + 6.2*(temperature-298.0)/25.;
294 else if((temperature>=323.0)&&(temperature<348.)){
295 H = 15.5 + 1.2*(temperature-323.0)/25.;
296 CH3 = 43.5 + 2.3*(temperature-323.0)/25.;
297 C6H5 = 123.4 + 6.3*(temperature-323.0)/25.;
300 assert(temperature>=348.0);
303 H = 16.7+2.1*(temperature-348.0)/25.;
304 CH3 = 45.8+2.5*(temperature-348.0)/25.;
305 C6H5 = 129.7+6.3*(temperature-348.0)/25.;
320 template <
class Evaluation>
323 temperature = min(temperature, 500.0);
324 temperature = max(temperature, 250.0);
326 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: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
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:45
static bool liquidIsCompressible()
Returns true iff the liquid phase is assumed to be compressible.
Definition: Mesitylene.hpp:218
static Evaluation liquidHeatCapacity(const Evaluation &temperature, const Evaluation &)
Specific heat cap of liquid mesitylene .
Definition: Mesitylene.hpp:276
static Scalar triplePressure()
Returns the pressure at mesitylene's triple point.
Definition: Mesitylene.hpp:88
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:99
static Evaluation liquidEnthalpy(const Evaluation &temperature, const Evaluation &pressure)
Specific enthalpy of liquid mesitylene .
Definition: Mesitylene.hpp:118
static Scalar molarMass()
The molar mass in of mesitylene.
Definition: Mesitylene.hpp:58
static Scalar criticalTemperature()
Returns the critical temperature of mesitylene.
Definition: Mesitylene.hpp:64
static Evaluation gasDensity(const Evaluation &temperature, const Evaluation &pressure)
The density of pure mesitylene vapor at a given pressure and temperature .
Definition: Mesitylene.hpp:190
static Evaluation heatVap(const Evaluation &temperature, const Evaluation &)
Latent heat of vaporization for mesitylene .
Definition: Mesitylene.hpp:144
static bool gasIsIdeal()
Returns true iff the gas phase is assumed to be ideal.
Definition: Mesitylene.hpp:212
static Evaluation gasEnthalpy(const Evaluation &temperature, const Evaluation &pressure)
Specific enthalpy of mesitylene vapor .
Definition: Mesitylene.hpp:178
static Scalar tripleTemperature()
Returns the temperature at mesitylene's triple point.
Definition: Mesitylene.hpp:82
static bool gasIsCompressible()
Returns true iff the gas phase is assumed to be compressible.
Definition: Mesitylene.hpp:206
static Evaluation liquidDensity(const Evaluation &temperature, const Evaluation &)
The density of pure mesitylene at a given pressure and temperature .
Definition: Mesitylene.hpp:200
static Evaluation liquidViscosity(Evaluation temperature, const Evaluation &)
The dynamic viscosity of pure mesitylene.
Definition: Mesitylene.hpp:255
static Evaluation molarLiquidDensity_(Evaluation temperature)
The molar density of pure mesitylene at a given pressure and temperature .
Definition: Mesitylene.hpp:321
static Scalar boilingTemperature()
Returns the temperature at mesitylene's boiling point (1 atm).
Definition: Mesitylene.hpp:76
static const char * name()
A human readable name for the mesitylene.
Definition: Mesitylene.hpp:52
static Scalar criticalPressure()
Returns the critical pressure of mesitylene.
Definition: Mesitylene.hpp:70
static Evaluation gasViscosity(Evaluation temperature, const Evaluation &, bool=true)
The dynamic viscosity of mesitylene vapor.
Definition: Mesitylene.hpp:229