28#ifndef EWOMS_NCP_INTENSIVE_QUANTITIES_HH
29#define EWOMS_NCP_INTENSIVE_QUANTITIES_HH
36#include <opm/material/constraintsolvers/NcpFlash.hpp>
37#include <opm/material/fluidstates/CompositionalFluidState.hpp>
38#include <opm/material/constraintsolvers/CompositionFromFugacities.hpp>
39#include <opm/material/common/Valgrind.hpp>
41#include <dune/common/fvector.hh>
42#include <dune/common/fmatrix.hh>
52template <
class TypeTag>
54 :
public GetPropType<TypeTag, Properties::DiscIntensiveQuantities>
57 ,
public GetPropType<TypeTag, Properties::FluxModule>::FluxIntensiveQuantities
76 enum { fugacity0Idx = Indices::fugacity0Idx };
77 enum { saturation0Idx = Indices::saturation0Idx };
78 enum { pressure0Idx = Indices::pressure0Idx };
79 enum { dimWorld = GridView::dimensionworld };
83 using ComponentVector = Dune::FieldVector<Evaluation, numComponents>;
84 using DimMatrix = Dune::FieldMatrix<Scalar, dimWorld, dimWorld>;
87 using FluxIntensiveQuantities =
typename FluxModule::FluxIntensiveQuantities;
105 ParentType::checkDefined();
117 fluidState_.setSaturation(numPhases - 1, 1.0 -
sumSat);
118 Opm::Valgrind::CheckDefined(
sumSat);
124 const auto& problem =
elemCtx.problem();
131 const Evaluation&
pressure0 = priVars.makeEvaluation(pressure0Idx,
timeIdx);
152 CompositionFromFugacitiesSolver::guessInitial(fluidState_,
phaseIdx,
fug);
161 Opm::Valgrind::CheckDefined(porosity_);
164 MaterialLaw::relativePermeabilities(relativePermeability_,
materialParams, fluidState_);
194 {
return fluidState_; }
200 {
return intrinsicPerm_; }
206 {
return relativePermeability_[
phaseIdx]; }
218 {
return porosity_; }
225#if !defined NDEBUG && HAVE_VALGRIND
226 ParentType::checkDefined();
228 Opm::Valgrind::CheckDefined(porosity_);
229 Opm::Valgrind::CheckDefined(relativePermeability_);
231 fluidState_.checkDefined();
236 DimMatrix intrinsicPerm_;
237 FluidState fluidState_;
238 Evaluation porosity_;
239 Evaluation relativePermeability_[numPhases];
240 Evaluation mobility_[numPhases];
Provides the volumetric quantities required for the calculation of molecular diffusive fluxes.
Definition diffusionmodule.hh:141
Provides the volumetric quantities required for the energy equation.
Definition energymodule.hh:532
Contains the quantities which are are constant within a finite volume in the compositional multi-phas...
Definition ncpintensivequantities.hh:58
void checkDefined() const
IntensiveQuantities::checkDefined.
Definition ncpintensivequantities.hh:223
const FluidState & fluidState() const
ImmiscibleIntensiveQuantities::fluidState.
Definition ncpintensivequantities.hh:193
const Evaluation & mobility(unsigned phaseIdx) const
ImmiscibleIntensiveQuantities::mobility.
Definition ncpintensivequantities.hh:211
const Evaluation & porosity() const
ImmiscibleIntensiveQuantities::porosity.
Definition ncpintensivequantities.hh:217
const Evaluation & relativePermeability(unsigned phaseIdx) const
ImmiscibleIntensiveQuantities::relativePermeability.
Definition ncpintensivequantities.hh:205
void update(const ElementContext &elemCtx, unsigned dofIdx, unsigned timeIdx)
IntensiveQuantities::update.
Definition ncpintensivequantities.hh:100
const DimMatrix & intrinsicPermeability() const
ImmiscibleIntensiveQuantities::intrinsicPermeability.
Definition ncpintensivequantities.hh:199
Classes required for molecular diffusion.
Contains the classes required to consider energy as a conservation quantity in a multi-phase module.
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
Declares the properties required for the NCP compositional multi-phase model.