Provides the auxiliary methods required for consideration of the diffusion equation. More...

#include <blackoildiffusionmodule.hh>

Public Types
using	ExtensiveQuantities = BlackOilDiffusionExtensiveQuantities<TypeTag,true>

Static Public Member Functions
static void	initFromState (const EclipseState &eclState)
	Initialize all internal data structures needed by the diffusion module.
static void	registerParameters ()
	Register all run-time parameters for the diffusion module.
template<class Context>
static OPM_HOST_DEVICE void	addDiffusiveFlux (RateVector &flux, const Context &context, unsigned spaceIdx, unsigned timeIdx)
	Adds the mass flux due to molecular diffusion to the flux vector over the integration point.
template<class IntensiveQuantities, class EvaluationArray, class RateVectorT>
static OPM_HOST_DEVICE void	addDiffusiveFlux (RateVectorT &flux, const IntensiveQuantities &inIq, const IntensiveQuantities &exIq, const Evaluation &diffusivity, const EvaluationArray &effectiveDiffusionCoefficient)
template<class IntensiveQuantities, class EvaluationArray>
static OPM_HOST_DEVICE void	addBioDiffFlux (RateVector &flux, const IntensiveQuantities &inIq, const IntensiveQuantities &exIq, const Evaluation &diffusivity, const EvaluationArray &effectiveBioDiffCoefficient)

Detailed Description

template<class TypeTag>
class Opm::BlackOilDiffusionModule< TypeTag, true >

Provides the auxiliary methods required for consideration of the diffusion equation.

Member Function Documentation

◆ addDiffusiveFlux()

template<class TypeTag>

template<class Context>

OPM_HOST_DEVICE void Opm::BlackOilDiffusionModule< TypeTag, true >::addDiffusiveFlux	(	RateVector &	flux,
		const Context &	context,
		unsigned	spaceIdx,
		unsigned	timeIdx )

inlinestatic

Adds the mass flux due to molecular diffusion to the flux vector over the integration point.

Following the notation in blackoilmodel.hh, the diffusive flux for component $\kappa$ in phase $\alpha$ is given by: $-\phi b_\alpha S_\alpha D \mathbf{grad}X_\alpha^\kappa$ , where $b_\alpha$ is the shrinkage/expansion factor [-], $S_\alpha$ is the saturation [-] D is the diffusion coefficient [L/T^2] and $X_\alpha^\kappa$ the component mass fraction [-] or molar fraction [-], depending on the input use_mole_fraction_ (default true) Each component mass/molar fraction can be computed using $R_s,\;R_v,\;R_{sw},\;R_{vw}$ . For example the mass fraction are given by, $X_w^G=\frac{R_{sw}}{R_{sw}+\rho_w/\rho_g}$ , where $\rho_w$ and $\rho_g$ are the reference densities. Considering the water phase and gas component as an example, for cells i and j, the discrete version of the diffusive flux at the face's integration point is given by $-b_{w,ij}S_{w,ij}D_{w,ij}(\frac{1}{R_{sw,ij}+\rho_w/\rho_g})DT_{ij}(R_{sw,i}-R_{sw,j})$ where $b_{w,ij}$ , $S_{w,ij}$ , $D_{w,ij}$ , and $R_{sw,ij}$ are computed using the arithmetic mean, and the ratio $\frac{1}{R_{sw,ij}+\rho_w/\rho_g}$ is denoted as conversion factor. The diffusivity $DT_{ij}$ is computed in ecltransmissibility_impl.hh, using the cells porosity, face area and distance between cell center and the integration point. For mol fraction the convertion factor is given by $\frac{1}{R_{sw,ij}+(Mm_g\rho_w)/(Mm_w\rho_g)}$