A fluid system with water, gas and NAPL as phases and water, air and mesitylene (DNAPL) as components. More...

#include <H2OAirMesityleneFluidSystem.hpp>

Inheritance diagram for Opm::FluidSystems::H2OAirMesitylene< Scalar >:

Public Types
typedef NullParameterCache	ParameterCache
	The type of the fluid system's parameter cache. More...

typedef Opm::Mesitylene< Scalar >	NAPL
	The type of the mesithylene/napl component. More...

typedef Opm::Air< Scalar >	Air
	The type of the air component. More...

typedef TabulatedH2O	H2O
	The type of the water component. More...

Static Public Member Functions
static void	init ()
	Initialize the fluid system's static parameters. More...

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 ranges. More...

static bool	isLiquid (unsigned phaseIdx)
	Return whether a phase is liquid. More...

static bool	isIdealGas (unsigned phaseIdx)
	Returns true if and only if a fluid phase is assumed to be an ideal gas. More...

static bool	isCompressible (unsigned phaseIdx)
	Returns true if and only if a fluid phase is assumed to be compressible. More...

static bool	isIdealMixture (unsigned)
	Returns true if and only if a fluid phase is assumed to be an ideal mixture. More...

static const char *	phaseName (unsigned phaseIdx)
	Return the human readable name of a fluid phase. More...

static const char *	componentName (unsigned compIdx)
	Return the human readable name of a component. More...

static Scalar	molarMass (unsigned compIdx)
	Return the molar mass of a component in [kg/mol]. More...

template<class FluidState , class LhsEval = typename FluidState::Scalar>
static LhsEval	density (const FluidState &fluidState, const ParameterCache &, unsigned phaseIdx)
	Calculate the density [kg/m^3] of a fluid phase. More...

template<class FluidState , class LhsEval = typename FluidState::Scalar>
static LhsEval	viscosity (const FluidState &fluidState, const ParameterCache &, unsigned phaseIdx)
	Calculate the dynamic viscosity of a fluid phase [Pa*s]. More...

template<class FluidState , class LhsEval = typename FluidState::Scalar>
static LhsEval	diffusionCoefficient (const FluidState &, const ParameterCache &, unsigned, unsigned)
	Calculate the binary molecular diffusion coefficient for a component in a fluid phase [mol^2 * s / (kg*m^3)]. More...

template<class FluidState , class LhsEval = typename FluidState::Scalar>
static LhsEval	fugacityCoefficient (const FluidState &fluidState, const ParameterCache &, unsigned phaseIdx, unsigned compIdx)
	Calculate the fugacity coefficient [Pa] of an individual component in a fluid phase. More...

template<class FluidState , class LhsEval = typename FluidState::Scalar>
static LhsEval	enthalpy (const FluidState &fluidState, const ParameterCache &, unsigned phaseIdx)
	Given a phase's composition, temperature, pressure and density, calculate its specific enthalpy [J/kg]. More...

template<class FluidState , class LhsEval = typename FluidState::Scalar>
static LhsEval	thermalConductivity (const FluidState &fluidState, const ParameterCache &, unsigned phaseIdx)
	Thermal conductivity of a fluid phase [W/(m K)]. More...

static LhsEval	heatCapacity (const FluidState &, const ParameterCache &, unsigned)
	Specific isobaric heat capacity of a fluid phase [J/kg]. More...

Static Public Attributes
static const int	numPhases = 3
	Number of fluid phases in the fluid system. More...

static const int	numComponents = 3
	Number of chemical species in the fluid system. More...

static const int	waterPhaseIdx = 0
	The index of the water phase. More...

static const int	naplPhaseIdx = 1
	The index of the NAPL phase. More...

static const int	gasPhaseIdx = 2
	The index of the gas phase. More...

static const int	H2OIdx = 0
	The index of the water component. More...

static const int	NAPLIdx = 1
	The index of the NAPL component. More...

static const int	airIdx = 2
	The index of the air pseudo-component. More...

Detailed Description

template<class Scalar>
class Opm::FluidSystems::H2OAirMesitylene< Scalar >

A fluid system with water, gas and NAPL as phases and water, air and mesitylene (DNAPL) as components.

Member Typedef Documentation

template<class Scalar >

typedef Opm::Air<Scalar> Opm::FluidSystems::H2OAirMesitylene< Scalar >::Air

The type of the air component.

template<class Scalar >

typedef TabulatedH2O Opm::FluidSystems::H2OAirMesitylene< Scalar >::H2O

The type of the water component.

template<class Scalar >

typedef Opm::Mesitylene<Scalar> Opm::FluidSystems::H2OAirMesitylene< Scalar >::NAPL

The type of the mesithylene/napl component.

template<class Scalar >

typedef NullParameterCache Opm::FluidSystems::H2OAirMesitylene< Scalar >::ParameterCache

The type of the fluid system's parameter cache.

The parameter cache can be used to avoid re-calculating expensive parameters for multiple quantities. Be aware that what the parameter cache actually does is specific for each fluid system and that it is opaque outside the fluid system.

Member Function Documentation

template<class Scalar >

static const char* Opm::FluidSystems::H2OAirMesitylene< Scalar >::componentName ( unsigned compIdx )

inlinestatic

Return the human readable name of a component.

References Opm::Mesitylene< Scalar >::name(), Opm::Air< Scalar >::name(), and Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::name().

template<class Scalar >

template<class FluidState , class LhsEval = typename FluidState::Scalar>

static LhsEval Opm::FluidSystems::H2OAirMesitylene< Scalar >::density	(	const FluidState &	fluidState,
		const ParameterCache &	,
		unsigned	phaseIdx
	)

inlinestatic

Calculate the density [kg/m^3] of a fluid phase.

References Opm::Air< Scalar >::gasDensity(), Opm::Mesitylene< Scalar >::gasDensity(), Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::gasDensity(), Opm::Mesitylene< Scalar >::liquidDensity(), Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::liquidDensity(), Opm::Mesitylene< Scalar >::liquidIsCompressible(), Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::liquidIsCompressible(), Opm::Mesitylene< Scalar >::molarMass(), Opm::Air< Scalar >::molarMass(), and Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::molarMass().

template<class Scalar >

template<class FluidState , class LhsEval = typename FluidState::Scalar>

static LhsEval Opm::FluidSystems::H2OAirMesitylene< Scalar >::diffusionCoefficient	(	const FluidState &	,
		const ParameterCache &	,
		unsigned	,
		unsigned
	)

inlinestatic

Calculate the binary molecular diffusion coefficient for a component in a fluid phase [mol^2 * s / (kg*m^3)].

Molecular diffusion of a compoent $\kappa$ is caused by a gradient of the mole fraction and follows the law

$J = - D \mathbf{grad} x^\kappa_\alpha$

where $x_\alpha^\kappa$ is the component's mole fraction in phase $\alpha$ , $D$ is the diffusion coefficient and $J$ is the diffusive flux.

References Opm::BinaryCoeff::Air_Mesitylene::gasDiffCoeff(), Opm::BinaryCoeff::H2O_Mesitylene::gasDiffCoeff(), and Opm::BinaryCoeff::H2O_Air::gasDiffCoeff().

template<class Scalar >

template<class FluidState , class LhsEval = typename FluidState::Scalar>

static LhsEval Opm::FluidSystems::H2OAirMesitylene< Scalar >::enthalpy	(	const FluidState &	fluidState,
		const ParameterCache &	,
		unsigned	phaseIdx
	)

inlinestatic

Given a phase's composition, temperature, pressure and density, calculate its specific enthalpy [J/kg].

References Opm::Mesitylene< Scalar >::gasEnthalpy(), Opm::Air< Scalar >::gasEnthalpy(), Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::gasEnthalpy(), Opm::Mesitylene< Scalar >::liquidEnthalpy(), and Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::liquidEnthalpy().

template<class Scalar >

template<class FluidState , class LhsEval = typename FluidState::Scalar>

static LhsEval Opm::FluidSystems::H2OAirMesitylene< Scalar >::fugacityCoefficient	(	const FluidState &	fluidState,
		const ParameterCache &	,
		unsigned	phaseIdx,
		unsigned	compIdx
	)

inlinestatic

Calculate the fugacity coefficient [Pa] of an individual component in a fluid phase.

The fugacity coefficient $\phi_\kappa$ is connected to the fugacity $f_\kappa$ and the component's molarity $x_\kappa$ by means of the relation

$f_\kappa = \phi_\kappa\,x_{\kappa}$

References Valgrind::CheckDefined(), Opm::BinaryCoeff::H2O_Mesitylene::henry(), Opm::BinaryCoeff::H2O_N2::henry(), Opm::Mesitylene< Scalar >::vaporPressure(), and Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::vaporPressure().

static LhsEval Opm::BaseFluidSystem< Scalar, H2OAirMesitylene< Scalar > >::heatCapacity	(	const FluidState &	,
		const ParameterCache &	,
		unsigned
	)

inlinestaticinherited

Specific isobaric heat capacity of a fluid phase [J/kg].

template<class Scalar >

static void Opm::FluidSystems::H2OAirMesitylene< Scalar >::init ( )

inlinestatic

Initialize the fluid system's static parameters.

template<class Scalar >

static void Opm::FluidSystems::H2OAirMesitylene< Scalar >::init	(	Scalar	tempMin,
		Scalar	tempMax,
		unsigned	nTemp,
		Scalar	pressMin,
		Scalar	pressMax,
		unsigned	nPress
	)

inlinestatic

Initialize the fluid system's static parameters using problem specific temperature and pressure ranges.

Parameters

tempMin	The minimum temperature used for tabulation of water [K]
tempMax	The maximum temperature used for tabulation of water [K]
nTemp	The number of ticks on the temperature axis of the table of water
pressMin	The minimum pressure used for tabulation of water [Pa]
pressMax	The maximum pressure used for tabulation of water [Pa]
nPress	The number of ticks on the pressure axis of the table of water

References Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::init(), and Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::isTabulated.

template<class Scalar >

static bool Opm::FluidSystems::H2OAirMesitylene< Scalar >::isCompressible ( unsigned phaseIdx )

inlinestatic

Returns true if and only if a fluid phase is assumed to be compressible.

Compressible means that the partial derivative of the density to the fluid pressure is always larger than zero.

References Opm::Mesitylene< Scalar >::liquidIsCompressible(), Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::liquidIsCompressible(), and Opm::FluidSystems::H2OAirMesitylene< Scalar >::waterPhaseIdx.

template<class Scalar >

static bool Opm::FluidSystems::H2OAirMesitylene< Scalar >::isIdealGas ( unsigned phaseIdx )

inlinestatic

Returns true if and only if a fluid phase is assumed to be an ideal gas.

References Opm::Air< Scalar >::gasIsIdeal(), Opm::Mesitylene< Scalar >::gasIsIdeal(), and Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::gasIsIdeal().

template<class Scalar >

static bool Opm::FluidSystems::H2OAirMesitylene< Scalar >::isIdealMixture ( unsigned )

inlinestatic

Returns true if and only if a fluid phase is assumed to be an ideal mixture.

We define an ideal mixture as a fluid phase where the fugacity coefficients of all components times the pressure of the phase are independent on the fluid composition. This assumption is true if Henry's law and Rault's law apply. If you are unsure what this function should return, it is safe to return false. The only damage done will be (slightly) increased computation times in some cases.

template<class Scalar >

static bool Opm::FluidSystems::H2OAirMesitylene< Scalar >::isLiquid ( unsigned phaseIdx )

inlinestatic

Return whether a phase is liquid.

References Opm::FluidSystems::H2OAirMesitylene< Scalar >::gasPhaseIdx.

template<class Scalar >

static Scalar Opm::FluidSystems::H2OAirMesitylene< Scalar >::molarMass ( unsigned compIdx )

inlinestatic

Return the molar mass of a component in [kg/mol].

References Opm::FluidSystems::H2OAirMesitylene< Scalar >::airIdx, Opm::Mesitylene< Scalar >::molarMass(), Opm::Air< Scalar >::molarMass(), Opm::TabulatedComponent< ScalarT, RawComponent, useVaporPressure >::molarMass(), and Opm::FluidSystems::H2OAirMesitylene< Scalar >::NAPLIdx.

template<class Scalar >

static const char* Opm::FluidSystems::H2OAirMesitylene< Scalar >::phaseName ( unsigned phaseIdx )

inlinestatic

Return the human readable name of a fluid phase.

template<class Scalar >

template<class FluidState , class LhsEval = typename FluidState::Scalar>

static LhsEval Opm::FluidSystems::H2OAirMesitylene< Scalar >::thermalConductivity	(	const FluidState &	fluidState,
		const ParameterCache &	,
		unsigned	phaseIdx
	)