28#ifndef EWOMS_FLASH_MODEL_HH
29#define EWOMS_FLASH_MODEL_HH
31#include <opm/material/constraintsolvers/NcpFlash.hpp>
33#include <opm/material/densead/Math.hpp>
35#include <opm/material/fluidmatrixinteractions/MaterialTraits.hpp>
36#include <opm/material/fluidmatrixinteractions/NullMaterial.hpp>
59template <
class TypeTag>
63namespace Opm::Properties {
67struct FlashModel {
using InheritsFrom = std::tuple<MultiPhaseBaseModel>; };
71template<
class TypeTag>
75template<
class TypeTag>
81template<
class TypeTag>
85template<
class TypeTag>
89template<
class TypeTag>
93template<
class TypeTag>
97template<
class TypeTag>
101template<
class TypeTag>
105template<
class TypeTag>
109template<
class TypeTag>
113template<
class TypeTag>
179template <
class TypeTag>
199 explicit FlashModel(Simulator& simulator)
200 : ParentType(simulator)
219 Parameters::Register<Parameters::FlashTolerance<Scalar>>
220 (
"The maximum tolerance for the flash solver to "
221 "consider the solution converged");
225 Parameters::SetDefault<Parameters::EnableIntensiveQuantityCache>(
true);
230 Parameters::SetDefault<Parameters::EnableThermodynamicHints>(
true);
244 const std::string& tmp = EnergyModule::primaryVarName(
pvIdx);
248 std::ostringstream
oss;
249 if (Indices::cTot0Idx <=
pvIdx &&
pvIdx < Indices::cTot0Idx
251 oss <<
"c_tot," << FluidSystem::componentName(
pvIdx
252 - Indices::cTot0Idx);
264 const std::string& tmp = EnergyModule::eqName(
eqIdx);
268 std::ostringstream
oss;
269 if (Indices::conti0EqIdx <=
eqIdx &&
eqIdx < Indices::conti0EqIdx
272 oss <<
"continuity^" << FluidSystem::componentName(
compIdx);
285 Scalar tmp = EnergyModule::primaryVarWeight(*
this, globalDofIdx,
pvIdx);
296 return FluidSystem::molarMass(
compIdx) / 100.0;
304 Scalar tmp = EnergyModule::eqWeight(*
this, globalDofIdx,
eqIdx);
311 return FluidSystem::molarMass(
compIdx);
314 void registerOutputModules_()
316 ParentType::registerOutputModules_();
Provides the auxiliary methods required for consideration of the energy equation.
Definition energymodule.hh:50
Implements a boundary vector for the fully implicit compositional multi-phase model which is based on...
Definition flashboundaryratevector.hh:45
This template class contains the data which is required to calculate all fluxes of components over a ...
Definition flashextensivequantities.hh:54
Defines the primary variable and equation indices for the compositional multi-phase model based on fl...
Definition flashindices.hh:47
Contains the intensive quantities of the flash-based compositional multi-phase model.
Definition flashintensivequantities.hh:58
Calculates the local residual of the compositional multi-phase model based on flash calculations.
Definition flashlocalresidual.hh:46
A compositional multi-phase model based on flash-calculations.
Definition flashmodel.hh:191
static void registerParameters()
Register all run-time parameters for the immiscible model.
Definition flashmodel.hh:206
Scalar primaryVarWeight(unsigned globalDofIdx, unsigned pvIdx) const
Returns the relative weight of a primary variable for calculating relative errors.
Definition flashmodel.hh:283
static std::string name()
Definition flashmodel.hh:236
Scalar eqWeight(unsigned globalDofIdx, unsigned eqIdx) const
Returns the relative weight of an equation.
Definition flashmodel.hh:302
std::string primaryVarName(unsigned pvIdx) const
Given an primary variable index, return a human readable name.
Definition flashmodel.hh:242
std::string eqName(unsigned eqIdx) const
Given an equation index, return a human readable name.
Definition flashmodel.hh:262
Represents the primary variables used by the compositional flow model based on flash calculations.
Definition flashprimaryvariables.hh:51
Implements a vector representing rates of conserved quantities.
Definition flashratevector.hh:48
A base class for fully-implicit multi-phase porous-media flow models which assume multiple fluid phas...
Definition multiphasebasemodel.hh:153
static void registerParameters()
Register all run-time parameters for the immiscible model.
Definition multiphasebasemodel.hh:179
Manages the initializing and running of time dependent problems.
Definition simulator.hh:97
VTK output module for the fluid composition.
Definition vtkcompositionmodule.hpp:57
static void registerParameters()
Register all run-time parameters for the Vtk output module.
Definition vtkcompositionmodule.hpp:86
VTK output module for quantities which make sense for models which incorperate molecular diffusion.
Definition vtkdiffusionmodule.hpp:58
static void registerParameters()
Register all run-time parameters for the Vtk output module.
Definition vtkdiffusionmodule.hpp:87
VTK output module for quantities which make sense for models which assume thermal equilibrium.
Definition vtkenergymodule.hpp:58
static void registerParameters()
Register all run-time parameters for the Vtk output module.
Definition vtkenergymodule.hpp:86
Contains the classes required to consider energy as a conservation quantity in a multi-phase module.
Defines the primary variable and equation indices for the compositional multi-phase model based on fl...
Contains the intensive quantities of the flash-based compositional multi-phase model.
Calculates the local residual of the compositional multi-phase model based on flash calculations.
Declares the parameters for the compositional multi-phase model based on flash calculations.
Represents the primary variables used by the compositional flow model based on flash calculations.
Implements a boundary vector for the fully implicit compositional multi-phase model which is based on...
This template class contains the data which is required to calculate all fluxes of components over a ...
Declares the properties required by the compositional multi-phase model based on flash calculations.
Implements a vector representing rates of conserved quantities.
A base class for fully-implicit multi-phase porous-media flow models which assume multiple fluid phas...
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
Type of object for specifying boundary conditions.
Definition fvbaseproperties.hh:119
Enable diffusive fluxes?
Definition multiphasebaseproperties.hh:79
Specify whether energy should be considered as a conservation quantity or not.
Definition multiphasebaseproperties.hh:76
Data required to calculate a flux over a face.
Definition fvbaseproperties.hh:149
The type of the flash constraint solver.
Definition flashproperties.hh:39
Enumerations used by the model.
Definition multiphasebaseproperties.hh:48
The secondary variables within a sub-control volume.
Definition fvbaseproperties.hh:133
The type of the local residual function.
Definition fvbaseproperties.hh:94
The type of the model.
Definition basicproperties.hh:88
A vector of primary variables within a sub-control volume.
Definition fvbaseproperties.hh:130
Vector containing volumetric or areal rates of quantities.
Definition fvbaseproperties.hh:116
The type tag for the isothermal single phase problems.
Definition flashmodel.hh:67
VTK output module for the fluid composition.
VTK output module for quantities which make sense for models which incorperate molecular diffusion.
VTK output module for quantities which make sense for models which assume thermal equilibrium.