Determines the pressures and saturations of all fluid phases given the total mass of all components. More...
#include <opm/common/ErrorMacros.hpp>#include <opm/common/Exceptions.hpp>#include <opm/material/common/Valgrind.hpp>#include <dune/common/fvector.hh>#include <dune/common/fmatrix.hh>#include <limits>#include <iostream>
Go to the source code of this file.
Classes | |
| class | Opm::ImmiscibleFlash< Scalar, FluidSystem > |
| Determines the pressures and saturations of all fluid phases given the total mass of all components. More... | |
Namespaces | |
| Opm | |
Detailed Description
Determines the pressures and saturations of all fluid phases given the total mass of all components.
In a N-phase, N-component context, we have the following unknowns if assuming immiscibility:
- N pressures
- N saturations
This sums up to 2*N unknowns. On the equations side of things, we have:
- N total component molarities
- 1 The sum of all saturations is 1
- N-1 Relations from capillary pressure
this also sums up to 2*N. We include the capillary pressures and the sum of the saturations explicitly. This means that we only solve for the first pressure and N-1 saturations.
If a fluid phase is incompressible, the pressure cannot determined by this, though. In this case the original pressure is kept, and the saturation of the phase is calculated by dividing the global molarity of the component by the phase density.