26#ifndef DUMUX_INJECTION_SPATIAL_PARAMETERS_HH
27#define DUMUX_INJECTION_SPATIAL_PARAMETERS_HH
43template<
class Gr
idGeometry,
class Scalar>
46 DissolutionSpatialParams<GridGeometry, Scalar>>
48 using GridView =
typename GridGeometry::GridView;
49 using FVElementGeometry =
typename GridGeometry::LocalView;
50 using SubControlVolume =
typename FVElementGeometry::SubControlVolume;
51 using Element =
typename GridView::template Codim<0>::Entity;
58 using GlobalPosition =
typename SubControlVolume::GlobalPosition;
70 solubilityLimit_ = getParam<Scalar>(
"SpatialParams.SolubilityLimit", 0.26);
71 referencePorosity_ = getParam<Scalar>(
"SpatialParams.referencePorosity", 0.11);
72 referencePermeability_ = getParam<Scalar>(
"SpatialParams.referencePermeability", 2.23e-14);
73 irreducibleLiqSat_ = getParam<Scalar>(
"SpatialParams.IrreducibleLiqSat", 0.2);
74 irreducibleGasSat_ = getParam<Scalar>(
"SpatialParams.IrreducibleGasSat", 1e-3);
75 pEntry1_ = getParam<Scalar>(
"SpatialParams.Pentry1", 500);
76 bcLambda1_ = getParam<Scalar>(
"SpatialParams.BCLambda1", 2);
79 materialParams_.setSwr(irreducibleLiqSat_);
80 materialParams_.setSnr(irreducibleGasSat_);
83 materialParams_.setPe(pEntry1_);
84 materialParams_.setLambda(bcLambda1_);
102 template<
class Sol
idSystem>
104 {
return 1.0-referencePorosity_; }
116 {
return referencePorosity_; }
127 template<
class ElementSolution>
129 const SubControlVolume& scv,
130 const ElementSolution& elemSol)
const
132 auto priVars =
evalSolution(element, element.geometry(), elemSol, scv.center(),
true);
134 Scalar sumPrecipitates = priVars[3];
137 const auto poro = max(1e-5, referencePorosity_ - sumPrecipitates);
138 return permLaw_.evaluatePermeability(referencePermeability_, referencePorosity_, poro);
142 {
return solubilityLimit_; }
144 Scalar
theta(
const SubControlVolume &scv)
const
149 {
return materialParams_; }
152 template<
class Flu
idSystem>
154 {
return FluidSystem::H2OIdx; }
162 Scalar solubilityLimit_;
163 Scalar referencePorosity_;
165 Scalar irreducibleLiqSat_;
166 Scalar irreducibleGasSat_;
Linear capillary pressure and relative permeability <-> saturation relations.
Implementation of a regularized version of the Brooks-Corey capillary pressure / relative permeabilit...
The Kozeny-Carman relationship for the calculation of a porosity-dependent permeability.
Class for the evaluation of the porosity subject to precipitation.
The base class for spatial parameters of multi-phase problems using a fully implicit discretization m...
PrimaryVariables evalSolution(const Element &element, const typename Element::Geometry &geometry, const typename FVElementGeometry::GridGeometry &gridGeometry, const BoxElementSolution< FVElementGeometry, PrimaryVariables > &elemSol, const typename Element::Geometry::GlobalCoordinate &globalPos, bool ignoreState=false)
Interpolates a given box element solution at a given global position. Uses the finite element cache o...
Definition: evalsolution.hh:95
make the local view function available whenever we use the grid geometry
Definition: adapt.hh:29
This material law takes a material law defined for effective saturations and converts it to a materia...
Definition: 2p/efftoabslaw.hh:60
AbsParamsT Params
Definition: 2p/efftoabslaw.hh:64
Implementation of the regularized Brooks-Corey capillary pressure / relative permeability <-> saturat...
Definition: regularizedbrookscorey.hh:62
The Kozeny-Carman relationship for the calculation of a porosity-dependent permeability....
Definition: permeabilitykozenycarman.hh:42
The base class for spatial parameters of multi-phase problems using a fully implicit discretization m...
Definition: fv.hh:57
const GridGeometry & gridGeometry() const
The finite volume grid geometry.
Definition: fv1p.hh:334
Spatial parameters for the dissolution problem where water is injected in a for flushing precipitated...
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:47
Scalar solubilityLimit() const
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:141
int wettingPhaseAtPos(const GlobalPosition &globalPos) const
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:153
PermeabilityType permeability(const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:128
Scalar inertVolumeFractionAtPos(const GlobalPosition &globalPos, int compIdx) const
Defines the volume fraction of the inert component.
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:103
Scalar referencePorosity(const Element &element, const SubControlVolume &scv) const
Defines the reference porosity distribution.
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:115
typename MaterialLaw::Params MaterialLawParams
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:65
Scalar PermeabilityType
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:62
const MaterialLawParams & materialLawParamsAtPos(const GlobalPosition &globalPos) const
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:148
Scalar theta(const SubControlVolume &scv) const
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:144
DissolutionSpatialParams(std::shared_ptr< const GridGeometry > gridGeometry)
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:67
Scalar minimalPorosity(const Element &element, const SubControlVolume &scv) const
Defines the minimum porosity distribution.
Definition: porousmediumflow/2pncmin/implicit/isothermal/spatialparams.hh:93
This material law takes a material law defined for effective saturations and converts it to a materia...