28#ifndef OPM_DIFFUSION_MODULE_HH
29#define OPM_DIFFUSION_MODULE_HH
31#include <dune/common/fvector.hh>
33#include <opm/material/common/Valgrind.hpp>
47template <
class TypeTag,
bool enableDiffusion>
53template <
class TypeTag>
71 template <
class Context>
82template <
class TypeTag>
91 enum { numPhases = FluidSystem::numPhases };
92 enum { numComponents = FluidSystem::numComponents };
93 enum { conti0EqIdx = Indices::conti0EqIdx };
108 template <
class Context>
140template <
class TypeTag,
bool enableDiffusion>
146template <
class TypeTag>
160 throw std::logic_error(
"Method tortuosity() does not make sense "
161 "if diffusion is disabled");
170 throw std::logic_error(
"Method diffusionCoefficient() does not "
171 "make sense if diffusion is disabled");
180 throw std::logic_error(
"Method effectiveDiffusionCoefficient() "
181 "does not make sense if diffusion is disabled");
189 template <
class Flu
idState>
192 const ElementContext&,
201template <
class TypeTag>
209 enum { numPhases = FluidSystem::numPhases };
210 enum { numComponents = FluidSystem::numComponents };
239 template <
class Flu
idState>
257 const Evaluation&
base =
263 * Toolbox::pow(
base, 10.0/3.0);
267 FluidSystem::diffusionCoefficient(fluidState,
276 Evaluation tortuosity_[numPhases];
277 Evaluation diffusionCoefficient_[numPhases][numComponents];
286template <
class TypeTag,
bool enableDiffusion>
292template <
class TypeTag>
309 template <
class Context,
class Flu
idState>
310 void updateBoundary_(
const Context&,
326 throw std::logic_error(
"The method moleFractionGradient() does not "
327 "make sense if diffusion is disabled.");
340 throw std::logic_error(
"The method effectiveDiffusionCoefficient() "
341 "does not make sense if diffusion is disabled.");
348template <
class TypeTag>
356 enum { dimWorld = GridView::dimensionworld };
360 using DimVector = Dune::FieldVector<Scalar, dimWorld>;
361 using DimEvalVector = Dune::FieldVector<Evaluation, dimWorld>;
374 const auto& normal =
face.normal();
395 for (
unsigned i = 0; i < normal.size(); ++i)
398 Opm::Valgrind::CheckDefined(moleFractionGradientNormal_[
phaseIdx][
compIdx]);
408 Opm::Valgrind::CheckDefined(effectiveDiffusionCoefficient_[
phaseIdx][
compIdx]);
413 template <
class Context,
class Flu
idState>
414 void updateBoundary_(
const Context& context,
417 const FluidState& fluidState)
419 const auto& stencil = context.stencil(
timeIdx);
420 const auto&
face = stencil.boundaryFace(
bfIdx);
422 const auto&
elemCtx = context.elementContext();
431 distVec -= context.element().geometry().global(
insideScv.localGeometry().center());
451 Opm::Valgrind::CheckDefined(moleFractionGradientNormal_[
phaseIdx][
compIdx]);
458 Opm::Valgrind::CheckDefined(effectiveDiffusionCoefficient_[
phaseIdx][
compIdx]);
484 Evaluation moleFractionGradientNormal_[numPhases][numComponents];
485 Evaluation effectiveDiffusionCoefficient_[numPhases][numComponents];
void update_(const ElementContext &, unsigned, unsigned)
Update the quantities required to calculate the diffusive mass fluxes.
Definition diffusionmodule.hh:304
const Evaluation & effectiveDiffusionCoefficient(unsigned, unsigned) const
The effective diffusion coeffcient of a component in a fluid phase at the face's integration point.
Definition diffusionmodule.hh:337
const Evaluation & moleFractionGradientNormal(unsigned, unsigned) const
The the gradient of the mole fraction times the face normal.
Definition diffusionmodule.hh:323
const Evaluation & moleFractionGradientNormal(unsigned phaseIdx, unsigned compIdx) const
The the gradient of the mole fraction times the face normal.
Definition diffusionmodule.hh:470
const Evaluation & effectiveDiffusionCoefficient(unsigned phaseIdx, unsigned compIdx) const
The effective diffusion coeffcient of a component in a fluid phase at the face's integration point.
Definition diffusionmodule.hh:480
void update_(const ElementContext &elemCtx, unsigned faceIdx, unsigned timeIdx)
Update the quantities required to calculate the diffusive mass fluxes.
Definition diffusionmodule.hh:368
Provides the quantities required to calculate diffusive mass fluxes.
Definition diffusionmodule.hh:287
Scalar diffusionCoefficient(unsigned, unsigned) const
Returns the molecular diffusion coefficient for a component in a phase.
Definition diffusionmodule.hh:168
Scalar effectiveDiffusionCoefficient(unsigned, unsigned) const
Returns the effective molecular diffusion coefficient of the porous medium for a component in a phase...
Definition diffusionmodule.hh:178
void update_(FluidState &, typename FluidSystem::template ParameterCache< typename FluidState::Scalar > &, const ElementContext &, unsigned, unsigned)
Update the quantities required to calculate diffusive mass fluxes.
Definition diffusionmodule.hh:190
Scalar tortuosity(unsigned) const
Returns the tortuousity of the sub-domain of a fluid phase in the porous medium.
Definition diffusionmodule.hh:158
Evaluation tortuosity(unsigned phaseIdx) const
Returns the tortuousity of the sub-domain of a fluid phase in the porous medium.
Definition diffusionmodule.hh:224
void update_(FluidState &fluidState, typename FluidSystem::template ParameterCache< typename FluidState::Scalar > ¶mCache, const ElementContext &elemCtx, unsigned dofIdx, unsigned timeIdx)
Update the quantities required to calculate diffusive mass fluxes.
Definition diffusionmodule.hh:240
Evaluation effectiveDiffusionCoefficient(unsigned phaseIdx, unsigned compIdx) const
Returns the effective molecular diffusion coefficient of the porous medium for a component in a phase...
Definition diffusionmodule.hh:231
Evaluation diffusionCoefficient(unsigned phaseIdx, unsigned compIdx) const
Returns the molecular diffusion coefficient for a component in a phase.
Definition diffusionmodule.hh:217
Provides the volumetric quantities required for the calculation of molecular diffusive fluxes.
Definition diffusionmodule.hh:141
static void addDiffusiveFlux(RateVector &, const Context &, unsigned, unsigned)
Adds the diffusive mass flux flux to the flux vector over a flux integration point.
Definition diffusionmodule.hh:72
static void registerParameters()
Register all run-time parameters for the diffusion module.
Definition diffusionmodule.hh:64
static 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 flux integration point.
Definition diffusionmodule.hh:109
static void registerParameters()
Register all run-time parameters for the diffusion module.
Definition diffusionmodule.hh:101
Provides the auxiliary methods required for consideration of the diffusion equation.
Definition diffusionmodule.hh:48
Callback class for a mole fraction of a component in a phase.
Definition quantitycallbacks.hh:425
Declare the properties used by the infrastructure code of the finite volume discretizations.
Defines the common properties required by the porous medium multi-phase models.
This file contains a set of helper functions used by VFPProd / VFPInj.
Definition blackoilboundaryratevector.hh:37
constexpr auto getPropValue()
get the value data member of a property
Definition propertysystem.hh:242
typename Properties::Detail::GetPropImpl< TypeTag, Property >::type::type GetPropType
get the type alias defined in the property (equivalent to old macro GET_PROP_TYPE(....
Definition propertysystem.hh:235
This method contains all callback classes for quantities that are required by some extensive quantiti...