Description

template<class TypeTag>
class Dumux::ScheideggersDispersionTensor< TypeTag >

Scheidegger's dispersion tensor

This class calculates the dispersion tensor for compositional and thermal models using Scheidegger's model. The dispersion tensor is given by:

\[D = \frac{\mathbf{v} \mathbf{v}^T}{\|\mathbf{v}\|} \cdot (\alpha_L - \alpha_T) + \|\mathbf{v}\| \cdot \alpha_T \cdot \mathbf{I} \]

where \(\mathbf{v}\) is the velocity vector, \(\alpha_L\) and \(\alpha_T\) are the longitudinal and transverse dispersivities, respectively, and \(\mathbf{I}\) is the identity matrix. The velocity is either taken from the spatial parameters or, for the box method, evaluated with the potential gradient of the advection type (AdvectionType::potentialGradient), i.e. with the same density in the gravity term as the advective flux. The dispersivities are taken from the spatial parameters.

Static Public Member Functions
template<class ElementFluxVariablesCache>
static DimWorldMatrix	compositionalDispersionTensor (const Problem &problem, const SubControlVolumeFace &scvf, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const int phaseIdx, const int compIdx)
template<class ElementFluxVariablesCache>
static DimWorldMatrix	thermalDispersionTensor (const Problem &problem, const SubControlVolumeFace &scvf, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const int phaseIdx)

Member Function Documentation

◆ compositionalDispersionTensor()

template<class TypeTag>

template<class ElementFluxVariablesCache>

DimWorldMatrix Dumux::ScheideggersDispersionTensor< TypeTag >::compositionalDispersionTensor	(	const Problem &	problem,
		const SubControlVolumeFace &	scvf,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const int	phaseIdx,
		const int	compIdx )

inlinestatic

◆ thermalDispersionTensor()