27#ifndef OPM_H2O_N2_FLUID_SYSTEM_HPP
28#define OPM_H2O_N2_FLUID_SYSTEM_HPP
50template <
class Scalar>
58 typedef ::Opm::IdealGas<Scalar>
IdealGas;
65 template <
class Evaluation>
83 static const char* name[] = {
89 return name[phaseIdx];
154 static const char* name[] = {
160 return name[compIdx];
167 return (compIdx ==
H2OIdx)
181 return (compIdx ==
H2OIdx)
195 return (compIdx ==
H2OIdx)
209 return (compIdx ==
H2OIdx)
250 if (H2O::isTabulated) {
252 pressMin, pressMax, nPress);
259 template <
class Flu
idState,
class LhsEval =
typename Flu
idState::Scalar,
class ParamCacheEval = LhsEval>
260 static LhsEval
density(
const FluidState& fluidState,
266 const auto& T = decay<LhsEval>(fluidState.temperature(phaseIdx));
267 const auto& p = decay<LhsEval>(fluidState.pressure(phaseIdx));
269 LhsEval sumMoleFrac = 0;
270 for (
unsigned compIdx = 0; compIdx <
numComponents; ++compIdx)
271 sumMoleFrac += decay<LhsEval>(fluidState.moleFraction(phaseIdx, compIdx));
289 const auto& xgH2O = decay<LhsEval>(fluidState.moleFraction(
gasPhaseIdx,
H2OIdx));
290 const auto& xgN2 = decay<LhsEval>(fluidState.moleFraction(
gasPhaseIdx,
N2Idx));
293 return (rho_gH2O + rho_gN2)/max(1e-5, sumMoleFrac);
297 template <
class Flu
idState,
class LhsEval =
typename Flu
idState::Scalar,
class ParamCacheEval = LhsEval>
304 const auto& T = decay<LhsEval>(fluidState.temperature(phaseIdx));
305 const auto& p = decay<LhsEval>(fluidState.pressure(phaseIdx));
321 LhsEval muResult = 0;
328 for (
unsigned compIdx = 0; compIdx <
numComponents; ++compIdx)
329 sumx += decay<LhsEval>(fluidState.moleFraction(phaseIdx, compIdx));
330 sumx = max(1e-10, sumx);
339 decay<LhsEval>(fluidState.moleFraction(phaseIdx, j))
343 decay<LhsEval>(fluidState.moleFraction(phaseIdx, i))
350 template <
class Flu
idState,
class LhsEval =
typename Flu
idState::Scalar,
class ParamCacheEval = LhsEval>
359 const auto& T = decay<LhsEval>(fluidState.temperature(phaseIdx));
360 const auto& p = decay<LhsEval>(fluidState.pressure(phaseIdx));
377 template <
class Flu
idState,
class LhsEval =
typename Flu
idState::Scalar,
class ParamCacheEval = LhsEval>
384 const auto& T = decay<LhsEval>(fluidState.temperature(phaseIdx));
385 const auto& p = decay<LhsEval>(fluidState.pressure(phaseIdx));
397 template <
class Flu
idState,
class LhsEval =
typename Flu
idState::Scalar,
class ParamCacheEval = LhsEval>
398 static LhsEval
enthalpy(
const FluidState& fluidState,
402 const auto& T = decay<LhsEval>(fluidState.temperature(phaseIdx));
403 const auto& p = decay<LhsEval>(fluidState.pressure(phaseIdx));
404 Valgrind::CheckDefined(T);
405 Valgrind::CheckDefined(p);
418 const auto& XgH2O = decay<LhsEval>(fluidState.massFraction(
gasPhaseIdx,
H2OIdx));
419 const auto& XgN2 = decay<LhsEval>(fluidState.massFraction(
gasPhaseIdx,
N2Idx));
427 template <
class Flu
idState,
class LhsEval =
typename Flu
idState::Scalar,
class ParamCacheEval = LhsEval>
434 const auto& T = decay<LhsEval>(fluidState.temperature(phaseIdx));
435 const auto& p = decay<LhsEval>(fluidState.pressure(phaseIdx));
443 const auto& xH2O = decay<LhsEval>(fluidState.moleFraction(phaseIdx,
H2OIdx));
444 const auto& xN2 = decay<LhsEval>(fluidState.moleFraction(phaseIdx,
N2Idx));
451 return lambdaN2 + lambdaH2O;
455 template <
class Flu
idState,
class LhsEval =
typename Flu
idState::Scalar,
class ParamCacheEval = LhsEval>
460 const auto& T = decay<LhsEval>(fluidState.temperature(phaseIdx));
461 const auto& p = decay<LhsEval>(fluidState.pressure(phaseIdx));
462 const auto& xAlphaH2O = decay<LhsEval>(fluidState.moleFraction(phaseIdx,
H2OIdx));
463 const auto& xAlphaN2 = decay<LhsEval>(fluidState.moleFraction(phaseIdx,
N2Idx));
464 const auto& XAlphaH2O = decay<LhsEval>(fluidState.massFraction(phaseIdx,
H2OIdx));
465 const auto& XAlphaN2 = decay<LhsEval>(fluidState.massFraction(phaseIdx,
N2Idx));
481 return XAlphaH2O*c_pH2O + XAlphaN2*c_pN2;
The base class for all fluid systems.
Material properties of pure water .
Binary coefficients for water and nitrogen.
Relations valid for an ideal gas.
Properties of pure molecular nitrogen .
A parameter cache which does nothing.
A simple version of pure water.
A generic class which tabulates all thermodynamic properties of a given component.
Some templates to wrap the valgrind client request macros.
The base class for all fluid systems.
Definition BaseFluidSystem.hpp:46
Scalar Scalar
The type used for scalar quantities.
Definition BaseFluidSystem.hpp:51
static Evaluation liquidDiffCoeff(const Evaluation &temperature, const Evaluation &)
Diffusion coefficent for molecular nitrogen in liquid water.
Definition H2O_N2.hpp:102
static Evaluation gasDiffCoeff(const Evaluation &temperature, const Evaluation &pressure)
Binary diffusion coefficent for molecular water and nitrogen.
Definition H2O_N2.hpp:70
static Evaluation henry(const Evaluation &temperature)
Henry coefficent for molecular nitrogen in liquid water.
Definition H2O_N2.hpp:52
A two-phase fluid system with water and nitrogen as components.
Definition H2ON2FluidSystem.hpp:53
static const int H2OIdx
The component index of water.
Definition H2ON2FluidSystem.hpp:139
static Scalar acentricFactor(unsigned compIdx)
The acentric factor of a component [].
Definition H2ON2FluidSystem.hpp:207
static LhsEval thermalConductivity(const FluidState &fluidState, const ParameterCache< ParamCacheEval > &, unsigned phaseIdx)
Thermal conductivity of a fluid phase [W/(m K)].
Definition H2ON2FluidSystem.hpp:428
static Scalar criticalTemperature(unsigned compIdx)
Critical temperature of a component [K].
Definition H2ON2FluidSystem.hpp:179
TabulatedH2O H2O
The component for pure water.
Definition H2ON2FluidSystem.hpp:144
static void init()
Initialize the fluid system's static parameters.
Definition H2ON2FluidSystem.hpp:226
static const char * phaseName(unsigned phaseIdx)
Return the human readable name of a fluid phase.
Definition H2ON2FluidSystem.hpp:81
static void init(Scalar tempMin, Scalar tempMax, unsigned nTemp, Scalar pressMin, Scalar pressMax, unsigned nPress)
Initialize the fluid system's static parameters using problem specific temperature and pressure range...
Definition H2ON2FluidSystem.hpp:247
static LhsEval density(const FluidState &fluidState, const ParameterCache< ParamCacheEval > &, unsigned phaseIdx)
Calculate the density [kg/m^3] of a fluid phase.
Definition H2ON2FluidSystem.hpp:260
static bool isIdealMixture(unsigned)
Returns true if and only if a fluid phase is assumed to be an ideal mixture.
Definition H2ON2FluidSystem.hpp:122
static bool isLiquid(unsigned phaseIdx)
Return whether a phase is liquid.
Definition H2ON2FluidSystem.hpp:93
static const int liquidPhaseIdx
Index of the liquid phase.
Definition H2ON2FluidSystem.hpp:76
static const int numPhases
Number of fluid phases in the fluid system.
Definition H2ON2FluidSystem.hpp:73
static const int gasPhaseIdx
Index of the gas phase.
Definition H2ON2FluidSystem.hpp:78
static LhsEval heatCapacity(const FluidState &fluidState, const ParameterCache< ParamCacheEval > &, unsigned phaseIdx)
Specific isobaric heat capacity of a fluid phase [J/kg].
Definition H2ON2FluidSystem.hpp:456
static LhsEval diffusionCoefficient(const FluidState &fluidState, const ParameterCache< ParamCacheEval > &, unsigned phaseIdx, unsigned)
Calculate the binary molecular diffusion coefficient for a component in a fluid phase [mol^2 * s / (k...
Definition H2ON2FluidSystem.hpp:378
static Scalar criticalPressure(unsigned compIdx)
Critical pressure of a component [Pa].
Definition H2ON2FluidSystem.hpp:193
static LhsEval enthalpy(const FluidState &fluidState, const ParameterCache< ParamCacheEval > &, unsigned phaseIdx)
Given a phase's composition, temperature, pressure and density, calculate its specific enthalpy [J/kg...
Definition H2ON2FluidSystem.hpp:398
static const char * componentName(unsigned compIdx)
Return the human readable name of a component.
Definition H2ON2FluidSystem.hpp:152
static const int numComponents
Number of chemical species in the fluid system.
Definition H2ON2FluidSystem.hpp:136
static Scalar molarMass(unsigned compIdx)
Return the molar mass of a component in [kg/mol].
Definition H2ON2FluidSystem.hpp:164
static bool isIdealGas(unsigned phaseIdx)
Returns true if and only if a fluid phase is assumed to be an ideal gas.
Definition H2ON2FluidSystem.hpp:111
static LhsEval viscosity(const FluidState &fluidState, const ParameterCache< ParamCacheEval > &, unsigned phaseIdx)
Calculate the dynamic viscosity of a fluid phase [Pa*s].
Definition H2ON2FluidSystem.hpp:298
SimpleN2 N2
The component for pure nitrogen.
Definition H2ON2FluidSystem.hpp:149
static const int N2Idx
The component index of molecular nitrogen.
Definition H2ON2FluidSystem.hpp:141
static bool isCompressible(unsigned phaseIdx)
Returns true if and only if a fluid phase is assumed to be compressible.
Definition H2ON2FluidSystem.hpp:100
static LhsEval fugacityCoefficient(const FluidState &fluidState, const ParameterCache< ParamCacheEval > &, unsigned phaseIdx, unsigned compIdx)
Calculate the fugacity coefficient [Pa] of an individual component in a fluid phase.
Definition H2ON2FluidSystem.hpp:351
Material properties of pure water .
Definition H2O.hpp:65
Relations valid for an ideal gas.
Definition IdealGas.hpp:38
Properties of pure molecular nitrogen .
Definition N2.hpp:49
static Evaluation gasThermalConductivity(const Evaluation &, const Evaluation &)
Specific heat conductivity of steam .
Definition N2.hpp:303
static Evaluation gasDensity(const Evaluation &temperature, const Evaluation &pressure)
The density of gas at a given pressure and temperature.
Definition N2.hpp:145
static bool gasIsIdeal()
Returns true iff the gas phase is assumed to be ideal.
Definition N2.hpp:160
static Scalar criticalPressure()
Returns the critical pressure of molecular nitrogen.
Definition N2.hpp:74
static Scalar criticalTemperature()
Returns the critical temperature of molecular nitrogen.
Definition N2.hpp:68
static Evaluation gasViscosity(const Evaluation &temperature, const Evaluation &)
The dynamic viscosity of at a given pressure and temperature.
Definition N2.hpp:267
static const char * name()
A human readable name for nitrogen.
Definition N2.hpp:56
static Scalar acentricFactor()
Acentric factor of .
Definition N2.hpp:85
static Scalar molarMass()
The molar mass in of molecular nitrogen.
Definition N2.hpp:62
static Evaluation gasHeatCapacity(const Evaluation &temperature, const Evaluation &)
Specific isobaric heat capacity of pure nitrogen gas.
Definition N2.hpp:236
static Evaluation gasEnthalpy(const Evaluation &temperature, const Evaluation &)
Specific enthalpy of pure nitrogen gas.
Definition N2.hpp:186
A parameter cache which does nothing.
Definition NullParameterCache.hpp:40
A generic class which tabulates all thermodynamic properties of a given component.
Definition TabulatedComponent.hpp:56
static Evaluation gasThermalConductivity(const Evaluation &temperature, const Evaluation &pressure)
The thermal conductivity of gaseous water .
Definition TabulatedComponent.hpp:495
static Evaluation gasEnthalpy(const Evaluation &temperature, const Evaluation &pressure)
Specific enthalpy of the gas .
Definition TabulatedComponent.hpp:282
static Evaluation liquidHeatCapacity(const Evaluation &temperature, const Evaluation &pressure)
Specific isobaric heat capacity of the liquid .
Definition TabulatedComponent.hpp:333
static void init(Scalar tempMin, Scalar tempMax, unsigned nTemp, Scalar pressMin, Scalar pressMax, unsigned nPress)
Initialize the tables.
Definition TabulatedComponent.hpp:72
static Scalar criticalTemperature()
Returns the critical temperature in of the component.
Definition TabulatedComponent.hpp:227
static Scalar criticalPressure()
Returns the critical pressure in of the component.
Definition TabulatedComponent.hpp:233
static Scalar molarMass()
The molar mass in of the component.
Definition TabulatedComponent.hpp:221
static bool liquidIsCompressible()
Returns true iff the liquid phase is assumed to be compressible.
Definition TabulatedComponent.hpp:408
static Evaluation liquidViscosity(const Evaluation &temperature, const Evaluation &pressure)
The dynamic viscosity of liquid.
Definition TabulatedComponent.hpp:478
static Evaluation gasDensity(const Evaluation &temperature, const Evaluation &pressure)
The density of gas at a given pressure and temperature .
Definition TabulatedComponent.hpp:426
static Evaluation gasHeatCapacity(const Evaluation &temperature, const Evaluation &pressure)
Specific isobaric heat capacity of the gas .
Definition TabulatedComponent.hpp:316
static bool gasIsIdeal()
Returns true iff the gas phase is assumed to be ideal.
Definition TabulatedComponent.hpp:414
static Scalar acentricFactor()
Returns the acentric factor of the component.
Definition TabulatedComponent.hpp:239
static Evaluation liquidDensity(const Evaluation &temperature, const Evaluation &pressure)
The density of liquid at a given pressure and temperature .
Definition TabulatedComponent.hpp:444
static Evaluation vaporPressure(const Evaluation &temperature)
The vapor pressure in of the component at a given temperature.
Definition TabulatedComponent.hpp:267
static Evaluation liquidThermalConductivity(const Evaluation &temperature, const Evaluation &pressure)
The thermal conductivity of liquid water .
Definition TabulatedComponent.hpp:512
static Evaluation gasViscosity(const Evaluation &temperature, const Evaluation &pressure)
The dynamic viscosity of gas.
Definition TabulatedComponent.hpp:461
static Evaluation liquidEnthalpy(const Evaluation &temperature, const Evaluation &pressure)
Specific enthalpy of the liquid .
Definition TabulatedComponent.hpp:299
static const char * name()
A human readable name for the component.
Definition TabulatedComponent.hpp:215
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30