A generic compositional multi-phase model using primary-variable switching. More...
#include <pvsmodel.hh>
Public Member Functions | |
PvsModel (Simulator &simulator) | |
std::string | primaryVarName (unsigned pvIdx) const |
Given an primary variable index, return a human readable name. More... | |
std::string | eqName (unsigned eqIdx) const |
Given an equation index, return a human readable name. More... | |
void | updateFailed () |
Called by the update() method if it was unsuccessful. More... | |
void | updateBegin () |
Called by the update() method before it tries to apply the newton method. More... | |
Scalar | primaryVarWeight (unsigned globalDofIdx, unsigned pvIdx) const |
Returns the relative weight of a primary variable for calculating relative errors. More... | |
Scalar | eqWeight (unsigned globalDofIdx, unsigned eqIdx) const |
Returns the relative weight of an equation. More... | |
void | advanceTimeLevel () |
Called by the problem if a time integration was successful, post processing of the solution is done and the result has been written to disk. More... | |
bool | switched () const |
Return true if the primary variables were switched for at least one vertex after the last timestep. | |
template<class DofEntity > | |
void | serializeEntity (std::ostream &outstream, const DofEntity &dofEntity) |
Write the current solution for a degree of freedom to a restart file. More... | |
template<class DofEntity > | |
void | deserializeEntity (std::istream &instream, const DofEntity &dofEntity) |
Reads the current solution variables for a degree of freedom from a restart file. More... | |
void | switchPrimaryVars_ () |
template<class FluidState > | |
void | printSwitchedPhases_ (const ElementContext &elemCtx, unsigned dofIdx, const FluidState &fs, short oldPhasePresence, const PrimaryVariables &newPv) const |
void | registerOutputModules_ () |
![]() | |
MultiPhaseBaseModel (Simulator &simulator) | |
void | finishInit () |
Apply the initial conditions to the model. More... | |
bool | phaseIsConsidered (unsigned phaseIdx OPM_UNUSED) const |
Returns true iff a fluid phase is used by the model. More... | |
void | globalPhaseStorage (EqVector &storage, unsigned phaseIdx) |
Compute the total storage inside one phase of all conservation quantities. More... | |
void | registerOutputModules_ () |
Static Public Member Functions | |
static void | registerParameters () |
Register all run-time parameters for the PVS compositional model. | |
static std::string | name () |
![]() | |
static void | registerParameters () |
Register all run-time parameters for the immiscible model. | |
Public Attributes | |
Scalar | referencePressure_ |
unsigned | numSwitched_ |
int | verbosity_ |
A generic compositional multi-phase model using primary-variable switching.
This model assumes a flow of fluid phases
, each of which is assumed to be a mixture
chemical species
.
By default, the standard multi-phase Darcy approach is used to determine the velocity, i.e.
although the actual approach which is used can be specified via the FluxModule
property. For example, the velocity model can by changed to the Forchheimer approach by
The core of the model is the conservation mass of each component by means of the equation
To close the system mathematically, model equations are also required. This model uses the primary variable switching assumptions, which are given by:
To make this approach applicable, a pseudo primary variable phase presence has to be introduced. Its purpose is to specify for each phase whether it is present or not. It is a pseudo primary variable because it is not directly considered when linearizing the system in the Newton method, but after each Newton iteration, it gets updated like the "real" primary variables. The following rules are used for this update procedure:
In all other cases don't modify the phase presence for phase .
The model always requires primary variables, but their interpretation is dependent on the phase presence:
The first primary variable is always interpreted as the pressure of the phase with the lowest index .
Then, "switching primary variables" follow, which are interpreted depending in the presence of the first
phases: If phase
is present, its saturation
is used as primary variable; if it is not present, the mole fraction
of the component with index
in the phase with the lowest index that is present
is used instead.
Finally, the mole fractions of the components with the largest index in the phase with the lowest index that is present
are used as primary variables.
|
inline |
Called by the problem if a time integration was successful, post processing of the solution is done and the result has been written to disk.
This should prepare the model for the next time integration.
|
inline |
Reads the current solution variables for a degree of freedom from a restart file.
instream | The stream from which the vertex data should be deserialized from |
dof | The Dune entity which's data should be deserialized |
|
inline |
Given an equation index, return a human readable name.
eqIdx | The index of the conservation equation of interest. |
|
inline |
Returns the relative weight of an equation.
globalVertexIdx | The global index of the vertex |
eqIdx | The index of the equation |
|
inlinestatic |
|
inline |
Given an primary variable index, return a human readable name.
pvIdx | The index of the primary variable of interest. |
|
inline |
Returns the relative weight of a primary variable for calculating relative errors.
globalDofIdx | The global index of the degree of freedom |
pvIdx | The index of the primary variable |
|
inline |
Write the current solution for a degree of freedom to a restart file.
outstream | The stream into which the vertex data should be serialized to |
dof | The Dune entity which's data should be serialized |
|
inline |
Called by the update() method before it tries to apply the newton method.
This is primary a hook which the actual model can overload.
|
inline |
Called by the update() method if it was unsuccessful.
This is primary a hook which the actual model can overload.