28#ifndef OPM_BLACKOIL_DIFFUSION_MODULE_HH
29#define OPM_BLACKOIL_DIFFUSION_MODULE_HH
33#include <opm/material/common/Valgrind.hpp>
35#include <dune/common/fvector.hh>
47template <
class TypeTag,
bool enableDiffusion>
50template <
class TypeTag,
bool enableDiffusion>
56template <
class TypeTag>
84 template <
class Context>
95template <
class TypeTag>
104 enum { numPhases = FluidSystem::numPhases };
105 enum { numComponents = FluidSystem::numComponents };
106 enum { conti0EqIdx = Indices::conti0EqIdx };
120 use_mole_fraction_ = eclState.getTableManager().diffMoleFraction();
152 template <
class Context>
157 if(!FluidSystem::enableDiffusion())
162 const auto& diffusivity =
extQuants.diffusivity();
163 const auto& effectiveDiffusionCoefficient =
extQuants.effectiveDiffusionCoefficient();
167 template<
class Flu
idState,
class EvaluationArray>
168 static void addDiffusiveFlux(RateVector&
flux,
171 const Evaluation& diffusivity,
172 const EvaluationArray& effectiveDiffusionCoefficient)
174 unsigned pvtRegionIndex =
fluidStateI.pvtRegionIndex();
176 if (!FluidSystem::phaseIsActive(
phaseIdx)) {
181 if (!FluidSystem::enableDissolvedGasInWater() && FluidSystem::waterPhaseIdx ==
phaseIdx) {
186 if (FluidSystem::gasPhaseIdx ==
phaseIdx && !FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx)) {
199 Evaluation
diffR = 0.0;
200 if (FluidSystem::enableDissolvedGas() && FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx) &&
phaseIdx == FluidSystem::oilPhaseIdx) {
205 if (FluidSystem::enableVaporizedOil() && FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) &&
phaseIdx == FluidSystem::gasPhaseIdx) {
207 convFactor = toFractionGasOil(pvtRegionIndex) / (1.0 +
rvAvg*toFractionGasOil(pvtRegionIndex));
210 if (FluidSystem::enableDissolvedGasInWater() &&
phaseIdx == FluidSystem::waterPhaseIdx) {
215 if (FluidSystem::enableVaporizedWater() &&
phaseIdx == FluidSystem::gasPhaseIdx) {
217 convFactor = toFractionGasWater(pvtRegionIndex)/ (1.0 +
rvAvg*toFractionGasWater(pvtRegionIndex));
244 static Scalar toFractionGasOil (
unsigned regionIdx) {
245 Scalar
mMOil = use_mole_fraction_? FluidSystem::molarMass(FluidSystem::oilCompIdx,
regionIdx) : 1;
246 Scalar
rhoO = FluidSystem::referenceDensity(FluidSystem::oilPhaseIdx,
regionIdx);
247 Scalar
mMGas = use_mole_fraction_? FluidSystem::molarMass(FluidSystem::gasCompIdx,
regionIdx) : 1;
248 Scalar
rhoG = FluidSystem::referenceDensity(FluidSystem::gasPhaseIdx,
regionIdx);
251 static Scalar toFractionGasWater (
unsigned regionIdx) {
252 Scalar
mMWater = use_mole_fraction_? FluidSystem::molarMass(FluidSystem::waterCompIdx,
regionIdx) : 1;
253 Scalar
rhoW = FluidSystem::referenceDensity(FluidSystem::waterPhaseIdx,
regionIdx);
254 Scalar
mMGas = use_mole_fraction_? FluidSystem::molarMass(FluidSystem::gasCompIdx,
regionIdx) : 1;
255 Scalar
rhoG = FluidSystem::referenceDensity(FluidSystem::gasPhaseIdx,
regionIdx);
259 static bool use_mole_fraction_;
262template <
typename TypeTag>
264BlackOilDiffusionModule<TypeTag, true>::use_mole_fraction_;
273template <
class TypeTag,
bool enableDiffusion>
279template <
class TypeTag>
293 throw std::logic_error(
"Method tortuosity() does not make sense "
294 "if diffusion is disabled");
303 throw std::logic_error(
"Method diffusionCoefficient() does not "
304 "make sense if diffusion is disabled");
313 throw std::logic_error(
"Method effectiveDiffusionCoefficient() "
314 "does not make sense if diffusion is disabled");
322 template <
class Flu
idState>
325 const ElementContext&,
334template <
class TypeTag>
342 enum { numPhases = FluidSystem::numPhases };
343 enum { numComponents = FluidSystem::numComponents };
355 if (
this == &rhs)
return *
this;
357 if (FluidSystem::enableDiffusion()) {
358 std::copy(rhs.tortuosity_, rhs.tortuosity_ + numPhases, tortuosity_);
359 for (
size_t i = 0; i < numPhases; ++i) {
360 std::copy(rhs.diffusionCoefficient_[i],
361 rhs.diffusionCoefficient_[i]+numComponents,
362 diffusionCoefficient_[i]);
401 template <
class Flu
idState>
404 const ElementContext& elemCtx,
409 if(!FluidSystem::enableDiffusion())
411 const auto& intQuants = elemCtx.intensiveQuantities(dofIdx,
timeIdx);
415 template<
class Flu
idState>
416 void update_(FluidState& fluidState,
418 const IntensiveQuantities& intQuants) {
422 if (!FluidSystem::phaseIsActive(
phaseIdx)) {
427 if (!FluidSystem::enableDissolvedGasInWater() && FluidSystem::waterPhaseIdx ==
phaseIdx) {
434 constexpr double myeps = 0.0001;
435 const Evaluation&
base =
438 * intQuants.fluidState().saturation(
phaseIdx));
440 1.0 / (intQuants.porosity() * intQuants.porosity())
441 * Toolbox::pow(
base, 10.0/3.0);
445 FluidSystem::diffusionCoefficient(fluidState,
454 Evaluation tortuosity_[numPhases];
455 Evaluation diffusionCoefficient_[numPhases][numComponents];
464template <
class TypeTag,
bool enableDiffusion>
465class BlackOilDiffusionExtensiveQuantities;
470template <
class TypeTag>
487 template <
class Context,
class Flu
idState>
488 void updateBoundary_(
const Context&,
501 throw std::logic_error(
"The method diffusivity() does not "
502 "make sense if diffusion is disabled.");
515 throw std::logic_error(
"The method effectiveDiffusionCoefficient() "
516 "does not make sense if diffusion is disabled.");
524template <
class TypeTag>
535 enum { dimWorld = GridView::dimensionworld };
539 using DimVector = Dune::FieldVector<Scalar, dimWorld>;
540 using DimEvalVector = Dune::FieldVector<Evaluation, dimWorld>;
542 using EvaluationArray = Evaluation[numPhases][numComponents];
551 if(!FluidSystem::enableDiffusion())
554 const auto& stencil = elemCtx.stencil(
timeIdx);
555 const auto&
face = stencil.interiorFace(faceIdx);
560 const Scalar diffusivity = elemCtx.problem().diffusivity(elemCtx,
face.interiorIndex(),
face.exteriorIndex());
561 const Scalar faceArea =
face.area();
562 diffusivity_ = diffusivity / faceArea;
564 Valgrind::CheckDefined(diffusivity_);
568 static void update(EvaluationArray& effectiveDiffusionCoefficient,
573 if (!FluidSystem::phaseIsActive(
phaseIdx)) {
577 if (!FluidSystem::enableDissolvedGasInWater() && FluidSystem::waterPhaseIdx ==
phaseIdx) {
586 Valgrind::CheckDefined(effectiveDiffusionCoefficient[
phaseIdx][
compIdx]);
591 template <
class Context,
class Flu
idState>
592 void updateBoundary_(
const Context&,
597 throw std::runtime_error(
"Not implemented: Diffusion across boundary not implemented for blackoil");
608 {
return diffusivity_; }
620 const auto& effectiveDiffusionCoefficient()
const{
621 return effectiveDiffusionCoefficient_;
626 EvaluationArray effectiveDiffusionCoefficient_;
const Evaluation & effectiveDiffusionCoefficient(unsigned, unsigned) const
The effective diffusion coeffcient of a component in a fluid phase at the face's integration point.
Definition blackoildiffusionmodule.hh:512
const Scalar & diffusivity() const
The diffusivity the face.
Definition blackoildiffusionmodule.hh:499
void update_(const ElementContext &, unsigned, unsigned)
Update the quantities required to calculate the diffusive mass fluxes.
Definition blackoildiffusionmodule.hh:482
Provides the quantities required to calculate diffusive mass fluxes.
Definition blackoildiffusionmodule.hh:526
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 blackoildiffusionmodule.hh:617
void update_(const ElementContext &elemCtx, unsigned faceIdx, unsigned timeIdx)
Update the quantities required to calculate the diffusive mass fluxes.
Definition blackoildiffusionmodule.hh:548
const Scalar & diffusivity() const
The diffusivity of the face.
Definition blackoildiffusionmodule.hh:607
Provides the quantities required to calculate diffusive mass fluxes.
Definition blackoildiffusionmodule.hh:51
Scalar effectiveDiffusionCoefficient(unsigned, unsigned) const
Returns the effective molecular diffusion coefficient of the porous medium for a component in a phase...
Definition blackoildiffusionmodule.hh:311
Scalar diffusionCoefficient(unsigned, unsigned) const
Returns the molecular diffusion coefficient for a component in a phase.
Definition blackoildiffusionmodule.hh:301
void update_(FluidState &, typename FluidSystem::template ParameterCache< typename FluidState::Scalar > &, const ElementContext &, unsigned, unsigned)
Update the quantities required to calculate diffusive mass fluxes.
Definition blackoildiffusionmodule.hh:323
Scalar tortuosity(unsigned) const
Returns the tortuousity of the sub-domain of a fluid phase in the porous medium.
Definition blackoildiffusionmodule.hh:291
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 blackoildiffusionmodule.hh:402
Evaluation diffusionCoefficient(unsigned phaseIdx, unsigned compIdx) const
Returns the molecular diffusion coefficient for a component in a phase.
Definition blackoildiffusionmodule.hh:371
Evaluation tortuosity(unsigned phaseIdx) const
Returns the tortuousity of the sub-domain of a fluid phase in the porous medium.
Definition blackoildiffusionmodule.hh:378
Evaluation effectiveDiffusionCoefficient(unsigned phaseIdx, unsigned compIdx) const
Returns the effective molecular diffusion coefficient of the porous medium for a component in a phase...
Definition blackoildiffusionmodule.hh:385
Provides the volumetric quantities required for the calculation of molecular diffusive fluxes.
Definition blackoildiffusionmodule.hh:274
static void registerParameters()
Register all run-time parameters for the diffusion module.
Definition blackoildiffusionmodule.hh:77
static void addDiffusiveFlux(RateVector &, const Context &, unsigned, unsigned)
Adds the diffusive mass flux flux to the flux vector over a flux integration point.
Definition blackoildiffusionmodule.hh:85
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 integration point.
Definition blackoildiffusionmodule.hh:153
static void registerParameters()
Register all run-time parameters for the diffusion module.
Definition blackoildiffusionmodule.hh:127
Provides the auxiliary methods required for consideration of the diffusion equation.
Definition blackoildiffusionmodule.hh:48
Declare the properties used by the infrastructure code of the finite volume discretizations.
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