The fluid system for the oil, gas and water phases of the SPE5 problem. More...

#include <Spe5FluidSystem.hpp>

Inheritance diagram for Opm::Spe5FluidSystem< Scalar >:

Classes
struct	ParameterCache
	The type of the fluid system's parameter cache. More...

Public Types
typedef ::Opm::H2O< Scalar >	H2O
	The component for pure water to be used. More...

typedef Scalar	Scalar
	The type used for scalar quantities. More...

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

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

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

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

static bool	isIdealMixture (unsigned phaseIdx)
	Returns true if and only if a fluid phase is assumed to be an ideal mixture. 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...

static Scalar	criticalTemperature (unsigned compIdx)
	Critical temperature of a component [K]. More...

static Scalar	criticalPressure (unsigned compIdx)
	Critical pressure of a component [Pa]. More...

static Scalar	criticalMolarVolume (unsigned compIdx)
	Molar volume of a component at the critical point [m^3/mol]. More...

static Scalar	acentricFactor (unsigned compIdx)
	The acentric factor of a component []. More...

static Scalar	interactionCoefficient (unsigned comp1Idx, unsigned comp2Idx)
	Returns the interaction coefficient for two components. More...

static void	init (Scalar minT=273.15, Scalar maxT=373.15, Scalar minP=1e4, Scalar maxP=100e6)
	Initialize the fluid system's static parameters. More...

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

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

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

static void	init ()
	Initialize the fluid system's static parameters. More...

static LhsEval	density (const FluidState &, const ParamCache &, unsigned)
	Calculate the density [kg/m^3] of a fluid phase. More...

static LhsEval	fugacityCoefficient (const FluidState &, ParamCache &, unsigned, unsigned)
	Calculate the fugacity coefficient [Pa] of an individual component in a fluid phase. More...

static LhsEval	viscosity (const FluidState &, ParamCache &, unsigned)
	Calculate the dynamic viscosity of a fluid phase [Pa*s]. More...

static LhsEval	diffusionCoefficient (const FluidState &, ParamCache &, unsigned, unsigned)
	Calculate the binary molecular diffusion coefficient for a component in a fluid phase [mol^2 * s / (kg*m^3)]. More...

static LhsEval	enthalpy (const FluidState &, ParamCache &, unsigned)
	Given a phase's composition, temperature, pressure and density, calculate its specific enthalpy [J/kg]. More...

static LhsEval	thermalConductivity (const FluidState &, ParamCache &, unsigned)
	Thermal conductivity of a fluid phase [W/(m K)]. More...

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

static unsigned	phaseIsActive (unsigned)
	Returns whether a fluid phase is active. More...

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

static const int	gasPhaseIdx = 0
	Index of the gas phase. More...

static const int	waterPhaseIdx = 1
	Index of the water phase. More...

static const int	oilPhaseIdx = 2
	Index of the oil phase. More...

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

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

static const int	C1Idx = 1
	Index of the C1 component. More...

static const int	C3Idx = 2
	Index of the C3 component. More...

static const int	C6Idx = 3
	Index of the C6 component. More...

static const int	C10Idx = 4
	Index of the C10 component. More...

static const int	C15Idx = 5
	Index of the C15 component. More...

static const int	C20Idx = 6
	Index of the C20 component. More...

Static Protected Member Functions
template<class LhsEval >
static LhsEval	henryCoeffWater_ (unsigned compIdx, const LhsEval &temperature)

Detailed Description

template<class Scalar>
class Opm::Spe5FluidSystem< Scalar >

The fluid system for the oil, gas and water phases of the SPE5 problem.

This problem comprises $H_2O$ , $C_1$ , $C_3$ , $C_6$ , $C_10$ , $C_15$ and $C_20$ as components.

See:

J.E. Killough, et al.: Fifth Comparative Solution Project: Evaluation of Miscible Flood Simulators, Ninth SPE Symposium on Reservoir Simulation, 1987

Member Typedef Documentation

◆ H2O

template<class Scalar >

typedef ::Opm::H2O<Scalar> Opm::Spe5FluidSystem< Scalar >::H2O

The component for pure water to be used.

◆ Scalar

typedef Scalar Opm::BaseFluidSystem< Scalar , Spe5FluidSystem< Scalar > >::Scalar

inherited

The type used for scalar quantities.

Member Function Documentation

◆ acentricFactor()

template<class Scalar >

static Scalar Opm::Spe5FluidSystem< Scalar >::acentricFactor ( unsigned compIdx )

inlinestatic

The acentric factor of a component [].

References Opm::H2O< Scalar >::acentricFactor(), Opm::Spe5FluidSystem< Scalar >::C10Idx, Opm::Spe5FluidSystem< Scalar >::C15Idx, Opm::Spe5FluidSystem< Scalar >::C1Idx, Opm::Spe5FluidSystem< Scalar >::C20Idx, Opm::Spe5FluidSystem< Scalar >::C3Idx, Opm::Spe5FluidSystem< Scalar >::C6Idx, and Opm::Spe5FluidSystem< Scalar >::H2OIdx.

◆ componentName()

template<class Scalar >

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

inlinestatic

Return the human readable name of a component.

References Opm::H2O< Scalar >::name(), and Opm::Spe5FluidSystem< Scalar >::numComponents.

◆ criticalMolarVolume()

template<class Scalar >

static Scalar Opm::Spe5FluidSystem< Scalar >::criticalMolarVolume ( unsigned compIdx )

inlinestatic

Molar volume of a component at the critical point [m^3/mol].

References Opm::Spe5FluidSystem< Scalar >::C10Idx, Opm::Spe5FluidSystem< Scalar >::C15Idx, Opm::Spe5FluidSystem< Scalar >::C1Idx, Opm::Spe5FluidSystem< Scalar >::C20Idx, Opm::Spe5FluidSystem< Scalar >::C3Idx, Opm::Spe5FluidSystem< Scalar >::C6Idx, Opm::H2O< Scalar >::criticalMolarVolume(), Opm::Spe5FluidSystem< Scalar >::criticalPressure(), Opm::Spe5FluidSystem< Scalar >::criticalTemperature(), and Opm::Spe5FluidSystem< Scalar >::H2OIdx.

◆ criticalPressure()

template<class Scalar >

static Scalar Opm::Spe5FluidSystem< Scalar >::criticalPressure ( unsigned compIdx )

inlinestatic

Critical pressure of a component [Pa].

References Opm::Spe5FluidSystem< Scalar >::C10Idx, Opm::Spe5FluidSystem< Scalar >::C15Idx, Opm::Spe5FluidSystem< Scalar >::C1Idx, Opm::Spe5FluidSystem< Scalar >::C20Idx, Opm::Spe5FluidSystem< Scalar >::C3Idx, Opm::Spe5FluidSystem< Scalar >::C6Idx, Opm::H2O< Scalar >::criticalPressure(), and Opm::Spe5FluidSystem< Scalar >::H2OIdx.

Referenced by Opm::Spe5FluidSystem< Scalar >::criticalMolarVolume().

◆ criticalTemperature()

template<class Scalar >

static Scalar Opm::Spe5FluidSystem< Scalar >::criticalTemperature ( unsigned compIdx )

inlinestatic

Critical temperature of a component [K].

References Opm::Spe5FluidSystem< Scalar >::C10Idx, Opm::Spe5FluidSystem< Scalar >::C15Idx, Opm::Spe5FluidSystem< Scalar >::C1Idx, Opm::Spe5FluidSystem< Scalar >::C20Idx, Opm::Spe5FluidSystem< Scalar >::C3Idx, Opm::Spe5FluidSystem< Scalar >::C6Idx, Opm::H2O< Scalar >::criticalTemperature(), and Opm::Spe5FluidSystem< Scalar >::H2OIdx.

Referenced by Opm::Spe5FluidSystem< Scalar >::criticalMolarVolume().

◆ density() [1/2]

static LhsEval Opm::BaseFluidSystem< Scalar , Spe5FluidSystem< Scalar > >::density	(	const FluidState &	,
		const ParamCache &	,
		unsigned
	)

inlinestaticinherited

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

◆ density() [2/2]

template<class Scalar >

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

static LhsEval Opm::Spe5FluidSystem< Scalar >::density	(	const FluidState &	fluidState,
		const ParameterCache< ParamCacheEval > &	paramCache,
		unsigned	phaseIdx
	)

inlinestatic

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

References Opm::Spe5ParameterCache< Scalar, FluidSystem >::molarVolume(), and Opm::Spe5FluidSystem< Scalar >::numPhases.

◆ diffusionCoefficient()

static LhsEval Opm::BaseFluidSystem< Scalar , Spe5FluidSystem< Scalar > >::diffusionCoefficient	(	const FluidState &	,
		ParamCache &	,
		unsigned	,
		unsigned
	)

inlinestaticinherited

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.

◆ enthalpy()

static LhsEval Opm::BaseFluidSystem< Scalar , Spe5FluidSystem< Scalar > >::enthalpy	(	const FluidState &	,
		ParamCache &	,
		unsigned
	)

inlinestaticinherited

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

◆ fugacityCoefficient() [1/2]

static LhsEval Opm::BaseFluidSystem< Scalar , Spe5FluidSystem< Scalar > >::fugacityCoefficient	(	const FluidState &	,
		ParamCache &	,
		unsigned	,
		unsigned
	)

inlinestaticinherited

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}$

◆ fugacityCoefficient() [2/2]

template<class Scalar >

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

static LhsEval Opm::Spe5FluidSystem< Scalar >::fugacityCoefficient	(	const FluidState &	fluidState,
		const ParameterCache< ParamCacheEval > &	paramCache,
		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 Opm::PengRobinsonMixture< Scalar, StaticParameters >::computeFugacityCoefficient(), Opm::Spe5FluidSystem< Scalar >::gasPhaseIdx, Opm::Spe5FluidSystem< Scalar >::henryCoeffWater_(), Opm::Spe5FluidSystem< Scalar >::numComponents, Opm::Spe5FluidSystem< Scalar >::numPhases, Opm::Spe5FluidSystem< Scalar >::oilPhaseIdx, and Opm::Spe5FluidSystem< Scalar >::waterPhaseIdx.

◆ heatCapacity()

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

inlinestaticinherited

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

◆ henryCoeffWater_()

template<class Scalar >

template<class LhsEval >

static LhsEval Opm::Spe5FluidSystem< Scalar >::henryCoeffWater_	(	unsigned	compIdx,
		const LhsEval &	temperature
	)

inlinestaticprotected

References Opm::Spe5FluidSystem< Scalar >::C10Idx, Opm::Spe5FluidSystem< Scalar >::C15Idx, Opm::Spe5FluidSystem< Scalar >::C1Idx, Opm::Spe5FluidSystem< Scalar >::C20Idx, Opm::Spe5FluidSystem< Scalar >::C3Idx, Opm::Spe5FluidSystem< Scalar >::C6Idx, Opm::Spe5FluidSystem< Scalar >::H2OIdx, and Opm::H2O< Scalar >::vaporPressure().

Referenced by Opm::Spe5FluidSystem< Scalar >::fugacityCoefficient().

◆ init() [1/2]

static void Opm::BaseFluidSystem< Scalar , Spe5FluidSystem< Scalar > >::init ( )

inlinestaticinherited

Initialize the fluid system's static parameters.

◆ init() [2/2]

template<class Scalar >

static void Opm::Spe5FluidSystem< Scalar >::init	(	Scalar	minT = `273.15`,
		Scalar	maxT = `373.15`,
		Scalar	minP = `1e4`,
		Scalar	maxP = `100e6`
	)

inlinestatic

Initialize the fluid system's static parameters.

Parameters

minT	The minimum temperature possibly encountered during the simulation
maxT	The maximum temperature possibly encountered during the simulation
minP	The minimum pressure possibly encountered during the simulation
maxP	The maximum pressure possibly encountered during the simulation

References Opm::PengRobinsonParams< Scalar >::a(), Opm::PengRobinsonParams< Scalar >::b(), Opm::Spe5FluidSystem< Scalar >::gasPhaseIdx, Opm::PengRobinson< Scalar >::init(), Opm::max(), Opm::min(), Opm::Spe5FluidSystem< Scalar >::numComponents, Opm::PengRobinsonParamsMixture< Scalar, FluidSystem, phaseIdx, useSpe5Relations >::pureParams(), and Opm::PengRobinsonParamsMixture< Scalar, FluidSystem, phaseIdx, useSpe5Relations >::updatePure().

◆ interactionCoefficient()

template<class Scalar >

static Scalar Opm::Spe5FluidSystem< Scalar >::interactionCoefficient	(	unsigned	comp1Idx,
		unsigned	comp2Idx
	)

inlinestatic

Returns the interaction coefficient for two components.

The values are given by the SPE5 paper.

References Opm::Spe5FluidSystem< Scalar >::C15Idx, Opm::Spe5FluidSystem< Scalar >::C1Idx, Opm::Spe5FluidSystem< Scalar >::C20Idx, Opm::Spe5FluidSystem< Scalar >::C3Idx, Opm::max(), and Opm::min().

◆ isCompressible()

template<class Scalar >

static bool Opm::Spe5FluidSystem< Scalar >::isCompressible ( unsigned )

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.

In the SPE-5 problems all fluids are compressible...

◆ isIdealGas()

template<class Scalar >

static bool Opm::Spe5FluidSystem< Scalar >::isIdealGas ( unsigned )

inlinestatic

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

◆ isIdealMixture()

template<class Scalar >

static bool Opm::Spe5FluidSystem< Scalar >::isIdealMixture ( unsigned phaseIdx )

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.

References Opm::Spe5FluidSystem< Scalar >::waterPhaseIdx.