A generic compositional multi-phase model using primary-variable switching. More...
#include <opm/material/densead/Math.hpp>
#include "pvsproperties.hh"
#include "pvslocalresidual.hh"
#include "pvsnewtonmethod.hh"
#include "pvsprimaryvariables.hh"
#include "pvsratevector.hh"
#include "pvsboundaryratevector.hh"
#include "pvsintensivequantities.hh"
#include "pvsextensivequantities.hh"
#include "pvsindices.hh"
#include <ewoms/models/common/multiphasebasemodel.hh>
#include <ewoms/models/common/diffusionmodule.hh>
#include <ewoms/models/common/energymodule.hh>
#include <ewoms/io/vtkcompositionmodule.hh>
#include <ewoms/io/vtkenergymodule.hh>
#include <ewoms/io/vtkdiffusionmodule.hh>
#include <opm/material/fluidmatrixinteractions/NullMaterial.hpp>
#include <opm/material/fluidmatrixinteractions/MaterialTraits.hpp>
#include <opm/common/ErrorMacros.hpp>
#include <opm/common/Exceptions.hpp>
#include <iostream>
#include <sstream>
#include <string>
#include <vector>
Go to the source code of this file.
Classes | |
class | Ewoms::PvsModel< TypeTag > |
A generic compositional multi-phase model using primary-variable switching. More... | |
class | Ewoms::PvsModel< TypeTag > |
A generic compositional multi-phase model using primary-variable switching. More... | |
Functions | |
Ewoms::Properties::NEW_TYPE_TAG (PvsModel, INHERITS_FROM(MultiPhaseBaseModel, VtkPhasePresence, VtkComposition, VtkEnergy, VtkDiffusion)) | |
The type tag for the isothermal single phase problems. | |
Ewoms::Properties::SET_TYPE_PROP (PvsModel, LocalResidual, Ewoms::PvsLocalResidual< TypeTag >) | |
Use the PVS local jacobian operator for the PVS model. | |
Ewoms::Properties::SET_TYPE_PROP (PvsModel, NewtonMethod, Ewoms::PvsNewtonMethod< TypeTag >) | |
Use the PVS specific newton method for the PVS model. | |
Ewoms::Properties::SET_TYPE_PROP (PvsModel, Model, Ewoms::PvsModel< TypeTag >) | |
the Model property | |
Ewoms::Properties::SET_TYPE_PROP (PvsModel, PrimaryVariables, Ewoms::PvsPrimaryVariables< TypeTag >) | |
the PrimaryVariables property | |
Ewoms::Properties::SET_TYPE_PROP (PvsModel, RateVector, Ewoms::PvsRateVector< TypeTag >) | |
the RateVector property | |
Ewoms::Properties::SET_TYPE_PROP (PvsModel, BoundaryRateVector, Ewoms::PvsBoundaryRateVector< TypeTag >) | |
the BoundaryRateVector property | |
Ewoms::Properties::SET_TYPE_PROP (PvsModel, IntensiveQuantities, Ewoms::PvsIntensiveQuantities< TypeTag >) | |
the IntensiveQuantities property | |
Ewoms::Properties::SET_TYPE_PROP (PvsModel, ExtensiveQuantities, Ewoms::PvsExtensiveQuantities< TypeTag >) | |
the ExtensiveQuantities property | |
Ewoms::Properties::SET_TYPE_PROP (PvsModel, Indices, Ewoms::PvsIndices< TypeTag, 0 >) | |
The indices required by the isothermal PVS model. | |
Ewoms::Properties::SET_INT_PROP (PvsModel, PvsVerbosity, 1) | |
Ewoms::Properties::SET_BOOL_PROP (PvsModel, EnableEnergy, false) | |
Disable the energy equation by default. | |
Ewoms::Properties::SET_BOOL_PROP (PvsModel, EnableDiffusion, false) | |
Ewoms::Properties::SET_SCALAR_PROP (PvsModel, PvsPressureBaseWeight, 1.0) | |
The basis value for the weight of the pressure primary variable. | |
Ewoms::Properties::SET_SCALAR_PROP (PvsModel, PvsSaturationsBaseWeight, 1.0) | |
The basis value for the weight of the saturation primary variables. | |
Ewoms::Properties::SET_SCALAR_PROP (PvsModel, PvsMoleFractionsBaseWeight, 1.0) | |
The basis value for the weight of the mole fraction primary variables. | |
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.