Non-isothermal gas injection problem where a air is injected into a fully water saturated medium.
More...
|
| WaterAirProblem (Simulator &simulator) |
|
void | finishInit () |
| Called by the Ewoms::Simulator in order to initialize the problem. More...
|
|
|
std::string | name () const |
| The problem name. More...
|
|
void | endTimeStep () |
| Called by the simulator after each time integration. More...
|
|
template<class Context > |
const DimMatrix & | intrinsicPermeability (const Context &context, unsigned spaceIdx, unsigned timeIdx) const |
|
template<class Context > |
Scalar | porosity (const Context &context, unsigned spaceIdx, unsigned timeIdx) const |
|
template<class Context > |
const MaterialLawParams & | materialLawParams (const Context &context, unsigned spaceIdx, unsigned timeIdx) const |
|
template<class Context > |
Scalar | heatCapacitySolid (const Context &context OPM_UNUSED, unsigned spaceIdx OPM_UNUSED, unsigned timeIdx OPM_UNUSED) const |
|
template<class Context > |
const HeatConductionLawParams & | heatConductionParams (const Context &context, unsigned spaceIdx, unsigned timeIdx) const |
|
|
template<class Context > |
void | boundary (BoundaryRateVector &values, const Context &context, unsigned spaceIdx, unsigned timeIdx) const |
| Evaluate the boundary conditions for a boundary segment. More...
|
|
|
template<class Context > |
void | initial (PrimaryVariables &values, const Context &context, unsigned spaceIdx, unsigned timeIdx) const |
| Evaluate the initial value for a control volume. More...
|
|
template<class Context > |
void | source (RateVector &rate, const Context &context OPM_UNUSED, unsigned spaceIdx OPM_UNUSED, unsigned timeIdx OPM_UNUSED) const |
| Evaluate the source term for all phases within a given sub-control-volume. More...
|
|
template<class TypeTag>
class Ewoms::WaterAirProblem< TypeTag >
Non-isothermal gas injection problem where a air is injected into a fully water saturated medium.
During buoyancy driven upward migration, the gas passes a rectangular high temperature area. This decreases the temperature of the high-temperature area and accelerates gas infiltration due to the lower viscosity of the gas. (Be aware that the pressure of the gas is approximately constant within the lens, so the density of the gas is reduced. This more than off-sets the viscosity increase of the gas at constant density.)
The domain is sized 40 m times 40 m. The rectangular area with increased temperature (380 K) starts at (20 m, 5 m) and ends at (30 m, 35 m).
For the mass conservation equation, no-flow boundary conditions are used on the top and on the bottom of the domain, while free-flow conditions apply on the left and the right boundary. Gas is injected at bottom from 15 m to 25 m at a rate of 0.001 kg/(s m^2) by means if a forced inflow boundary condition.
At the free-flow boundaries, the initial condition for the bulk part of the domain is assumed, i. e. hydrostatic pressure, a gas saturation of zero and a geothermal temperature gradient of 0.03 K/m.