The spatial parameters for the 2p2c chemical nonequilibrium problem. More...
#include <test/porousmediumflow/2p2c/implicit/chemicalnonequilibrium/spatialparams.hh>
The spatial parameters for the 2p2c chemical nonequilibrium problem.
Public Types | |
using | PermeabilityType = Scalar |
Export permeability type. More... | |
using | MaterialLaw = EffToAbsLaw< EffectiveLaw > |
Export the type used for the material law. More... | |
using | MaterialLawParams = typename MaterialLaw::Params |
using | EffectiveIALawAwn = AwnSurfacePcMaxFct< Scalar > |
using | AwnSurface = EffToAbsLawIA< EffectiveIALawAwn, MaterialLawParams > |
using | AwnSurfaceParams = typename AwnSurface::Params |
using | AwsSurfaceParams = Scalar |
using | AnsSurfaceParams = Scalar |
Public Member Functions | |
TwoPTwoCChemicalNonequilibriumSpatialParams (std::shared_ptr< const GridGeometry > gridGeometry) | |
template<class ElementSolution > | |
PermeabilityType | permeability (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Scalar | porosityAtPos (const GlobalPosition &globalPos) const |
Defines the porosity \([-]\) of the soil. More... | |
const MaterialLawParams & | materialLawParamsAtPos (const GlobalPosition &globalPos) const |
Function for defining the parameters needed by constitutive relationships (kr-sw, pc-sw, etc.). More... | |
template<class ElementSolution > | |
const AwnSurfaceParams & | aWettingNonWettingSurfaceParams (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Returns a reference to the container object for the parametrization of the surface between wetting and non-Wetting phase. More... | |
template<class ElementSolution > | |
const Scalar | pcMax (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Returns the maximum capillary pressure for the given pc-Sw curve. More... | |
const Scalar | characteristicLengthAtPos (const GlobalPosition &globalPos) const |
Returns the characteristic length for the mass transfer. More... | |
const Scalar | factorMassTransferAtPos (const GlobalPosition &globalPos) const |
Return the pre factor the the energy transfer. More... | |
template<class FluidSystem > | |
int | wettingPhaseAtPos (const GlobalPosition &globalPos) const |
Function for defining which phase is to be considered as the wetting phase. More... | |
const Scalar | characteristicLength (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Return the characteristic length for the mass transfer. More... | |
const Scalar | factorEnergyTransfer (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Return the pre-factor the the energy transfer. More... | |
const Scalar | factorEnergyTransferAtPos (const GlobalPosition &globalPos) const |
Return the pre factor the the energy transfer. More... | |
const Scalar | factorMassTransfer (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Return the pre-factor the the mass transfer. More... | |
decltype(auto) | materialLawParams (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Function for defining the parameters needed by constitutive relationships (kr-sw, pc-sw, etc.). More... | |
int | wettingPhase (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Function for defining which phase is to be considered as the wetting phase. More... | |
const GlobalPosition & | gravity (const GlobalPosition &pos) const |
Returns the acceleration due to gravity \(\mathrm{[m/s^2]}\). More... | |
Scalar | harmonicMean (const Scalar T1, const Scalar T2, const GlobalPosition &normal) const |
Harmonic average of a discontinuous scalar field at discontinuity interface (for compatibility reasons with the function below) More... | |
DimWorldMatrix | harmonicMean (const DimWorldMatrix &T1, const DimWorldMatrix &T2, const GlobalPosition &normal) const |
Harmonic average of a discontinuous tensorial field at discontinuity interface. More... | |
decltype(auto) | permeability (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Function for defining the (intrinsic) permeability \([m^2]\). More... | |
Scalar | porosity (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Function for defining the porosity. That is possibly solution dependent. More... | |
Scalar | inertVolumeFraction (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol, int compIdx) const |
Function for defining the solid volume fraction. That is possibly solution dependent. More... | |
Scalar | inertVolumeFraction (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol, int compIdx) const |
Scalar | inertVolumeFraction (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol, int compIdx) const |
Scalar | beaversJosephCoeffAtPos (const GlobalPosition &globalPos) const |
Function for defining the Beavers-Joseph coefficient for multidomain problems \(\mathrm{[-]}\). More... | |
Scalar | forchCoeff (const SubControlVolumeFace &scvf) const |
Apply the Forchheimer coefficient for inertial forces calculation. More... | |
const GridGeometry & | fvGridGeometry () const |
The finite volume grid geometry. More... | |
const GridGeometry & | gridGeometry () const |
The finite volume grid geometry. More... | |
Static Public Member Functions | |
static constexpr bool | evaluatePermeabilityAtScvfIP () |
If the permeability should be evaluated directly at the scvf integration point (for convergence tests with analytical and continuous perm functions) or is evaluated at the scvs (for permeability fields with discontinuities) -> default. More... | |
Protected Member Functions | |
TwoPTwoCChemicalNonequilibriumSpatialParams< GridGeometry, Scalar > & | asImp_ () |
const TwoPTwoCChemicalNonequilibriumSpatialParams< GridGeometry, Scalar > & | asImp_ () const |
|
inherited |
using Dumux::TwoPTwoCChemicalNonequilibriumSpatialParams< GridGeometry, Scalar >::AwnSurface = EffToAbsLawIA<EffectiveIALawAwn, MaterialLawParams> |
using Dumux::TwoPTwoCChemicalNonequilibriumSpatialParams< GridGeometry, Scalar >::AwnSurfaceParams = typename AwnSurface::Params |
|
inherited |
using Dumux::TwoPTwoCChemicalNonequilibriumSpatialParams< GridGeometry, Scalar >::EffectiveIALawAwn = AwnSurfacePcMaxFct<Scalar> |
using Dumux::TwoPTwoCChemicalNonequilibriumSpatialParams< GridGeometry, Scalar >::MaterialLaw = EffToAbsLaw<EffectiveLaw> |
Export the type used for the material law.
using Dumux::TwoPTwoCChemicalNonequilibriumSpatialParams< GridGeometry, Scalar >::MaterialLawParams = typename MaterialLaw::Params |
using Dumux::TwoPTwoCChemicalNonequilibriumSpatialParams< GridGeometry, Scalar >::PermeabilityType = Scalar |
Export permeability type.
|
inline |
|
inlineprotectedinherited |
|
inlineprotectedinherited |
|
inline |
Returns a reference to the container object for the parametrization of the surface between wetting and non-Wetting phase.
The position is determined based on the coordinate of the vertex belonging to the considered sub-control volume.
element | The finite element |
scv | The sub-control volume |
elemSol | The element solution |
|
inlineinherited |
Function for defining the Beavers-Joseph coefficient for multidomain problems \(\mathrm{[-]}\).
globalPos | The global position |
|
inlineinherited |
Return the characteristic length for the mass transfer.
The position is determined based on the coordinate of the vertex belonging to the considered sub control volume.
|
inline |
Returns the characteristic length for the mass transfer.
globalPos | The position in global coordinates. |
|
inlinestaticconstexprinherited |
If the permeability should be evaluated directly at the scvf integration point (for convergence tests with analytical and continuous perm functions) or is evaluated at the scvs (for permeability fields with discontinuities) -> default.
|
inlineinherited |
Return the pre-factor the the energy transfer.
The position is determined based on the coordinate of the vertex belonging to the considered sub control volume.
|
inlineinherited |
Return the pre factor the the energy transfer.
globalPos | The position in global coordinates. |
|
inlineinherited |
Return the pre-factor the the mass transfer.
The position is determined based on the coordinate of the vertex belonging to the considered sub control volume.
|
inline |
Return the pre factor the the energy transfer.
globalPos | The position in global coordinates. |
|
inlineinherited |
Apply the Forchheimer coefficient for inertial forces calculation.
scvf | The sub-control volume face where the intrinsic velocity ought to be calculated. |
|
inlineinherited |
The finite volume grid geometry.
|
inlineinherited |
Returns the acceleration due to gravity \(\mathrm{[m/s^2]}\).
The default behaviour is a constant gravity vector; if the Problem.EnableGravity
parameter is true, \(\boldsymbol{g} = ( 0,\dots,\ -9.81)^T \), else \(\boldsymbol{g} = ( 0,\dots, 0)^T \).
pos | the spatial position at which to evaulate the gravity vector |
|
inlineinherited |
The finite volume grid geometry.
|
inlineinherited |
Harmonic average of a discontinuous tensorial field at discontinuity interface.
T1 | first tensor |
T2 | second tensor |
normal | The unit normal vector of the interface |
|
inlineinherited |
Harmonic average of a discontinuous scalar field at discontinuity interface (for compatibility reasons with the function below)
T1 | first scalar parameter |
T2 | second scalar parameter |
normal | The unit normal vector of the interface |
|
inlineinherited |
Function for defining the solid volume fraction. That is possibly solution dependent.
element | The current element |
scv | The sub-control volume inside the element. |
elemSol | The solution at the dofs connected to the element. |
compIdx | The solid component index |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
Function for defining the parameters needed by constitutive relationships (kr-sw, pc-sw, etc.).
element | The current element |
scv | The sub-control volume inside the element. |
elemSol | The solution at the dofs connected to the element. |
|
inline |
Function for defining the parameters needed by constitutive relationships (kr-sw, pc-sw, etc.).
globalPos | The global position of the sub-control volume. |
|
inline |
Returns the maximum capillary pressure for the given pc-Sw curve.
Of course physically there is no such thing as a maximum capillary pressure. The parametrization (VG/BC) goes to infinity and physically there is only one pressure for single phase conditions. Here, this is used for fitting the interfacial area surface: the capillary pressure, where the interfacial area is zero. Technically this value is obtained as the capillary pressure of saturation zero. This value of course only exists for the case of a regularized pc-Sw relation.
element | The finite element |
scv | The sub-control volume |
elemSol | The element solution |
|
inlineinherited |
Function for defining the (intrinsic) permeability \([m^2]\).
element | The current element |
scv | The sub-control volume inside the element. |
elemSol | The solution at the dofs connected to the element. |
|
inline |
|
inlineinherited |
Function for defining the porosity. That is possibly solution dependent.
element | The current element |
scv | The sub-control volume inside the element. |
elemSol | The solution at the dofs connected to the element. |
|
inline |
Defines the porosity \([-]\) of the soil.
globalPos | The global position of the sub-control volume. |
|
inlineinherited |
Function for defining which phase is to be considered as the wetting phase.
element | The current element |
scv | The sub-control volume inside the element. |
elemSol | The solution at the dofs connected to the element. |
|
inline |
Function for defining which phase is to be considered as the wetting phase.
globalPos | The global position of the sub-control volume. |