24#ifndef DUMUX_BINARY_COEFF_AIR_MESITYLENE_HH
25#define DUMUX_BINARY_COEFF_AIR_MESITYLENE_HH
31namespace BinaryCoeff {
44 template <
class Scalar>
46 { DUNE_THROW(Dune::NotImplemented,
47 "Henry coefficient of air in mesitylene");
58 template <
class Scalar>
75 const Scalar M_m = 1e3*Mesitylene::molarMass();
76 const Scalar M_a = 1e3*Air::molarMass();
77 const Scalar Tb_m = 437.9;
78 const Scalar sigma_a = 3.711;
79 const Scalar T_scal_a = 78.6;
80 const Scalar V_B_m = 162.6;
83 const Scalar sigma_m = 1.18*cbrt(V_B_m);
84 const Scalar sigma_am = 0.5*(sigma_a + sigma_m);
85 const Scalar T_scal_m = 1.15*Tb_m;
86 const Scalar T_scal_am = sqrt(T_scal_a*T_scal_m);
89 T_star = max(T_star, 1e-5);
91 const Scalar Omega = 1.06036/pow(T_star, 0.1561) + 0.193/exp(T_star*0.47635)
92 + 1.03587/exp(T_star*1.52996) + 1.76474/exp(T_star*3.89411);
93 const Scalar B_ = 0.00217 - 0.0005*sqrt(1.0/M_a + 1.0/M_m);
94 const Scalar Mr = (M_a + M_m)/(M_a*M_m);
96 /(1e-5*
pressure*power(sigma_am, 2) * Omega);
108 template <
class Scalar>
Properties of mesitylene.
A simple class for the air fluid properties.
std::string temperature() noexcept
I/O name of temperature for equilibrium models.
Definition: name.hh:51
std::string pressure(int phaseIdx) noexcept
I/O name of pressure for multiphase systems.
Definition: name.hh:34
Binary coefficients for water and mesitylene.
Definition: air_mesitylene.hh:38
static Scalar liquidDiffCoeff(Scalar temperature, Scalar pressure)
Diffusion coefficient for air and mesitylene in liquid water.
Definition: air_mesitylene.hh:109
static Scalar gasDiffCoeff(Scalar temperature, Scalar pressure)
Binary diffusion coefficient for air and mesitylene. I used the method according to Wilke and Lee se...
Definition: air_mesitylene.hh:59
static Scalar henry(Scalar temperature)
Henry coefficient for mesitylene in air.
Definition: air_mesitylene.hh:45
A class for the air fluid properties.
Definition: air.hh:47
mesitylene
Definition: mesitylene.hh:48