#include <dune/common/fmatrix.hh>#include <dune/common/fvector.hh>#include <dune/grid/yaspgrid.hh>#include <dune/grid/io/file/dgfparser/dgfyasp.hh>#include <opm/material/constraintsolvers/ComputeFromReferencePhase.hpp>#include <opm/material/fluidmatrixinteractions/LinearMaterial.hpp>#include <opm/material/fluidmatrixinteractions/RegularizedBrooksCorey.hpp>#include <opm/material/fluidmatrixinteractions/EffToAbsLaw.hpp>#include <opm/material/fluidmatrixinteractions/MaterialTraits.hpp>#include <opm/material/fluidstates/ImmiscibleFluidState.hpp>#include <opm/material/fluidstates/CompositionalFluidState.hpp>#include <opm/material/fluidsystems/H2OAirFluidSystem.hpp>#include <opm/material/thermal/ConstantSolidHeatCapLaw.hpp>#include <opm/material/thermal/SomertonThermalConductionLaw.hpp>#include <opm/models/common/multiphasebaseparameters.hh>#include <opm/models/discretization/common/fvbasefdlocallinearizer.hh>#include <opm/models/pvs/pvsproperties.hh>#include <opm/simulators/linalg/parallelistlbackend.hh>#include <sstream>#include <string>
Go to the source code of this file.
Classes | |
| struct | Opm::Properties::TTag::WaterAirBaseProblem |
| struct | Opm::Properties::Grid< TypeTag, TTag::WaterAirBaseProblem > |
| struct | Opm::Properties::Problem< TypeTag, TTag::WaterAirBaseProblem > |
| struct | Opm::Properties::MaterialLaw< TypeTag, TTag::WaterAirBaseProblem > |
| struct | Opm::Properties::ThermalConductionLaw< TypeTag, TTag::WaterAirBaseProblem > |
| struct | Opm::Properties::SolidEnergyLaw< TypeTag, TTag::WaterAirBaseProblem > |
| struct | Opm::Properties::FluidSystem< TypeTag, TTag::WaterAirBaseProblem > |
| struct | Opm::Properties::LinearSolverSplice< TypeTag, TTag::WaterAirBaseProblem > |
| struct | Opm::Properties::LinearSolverWrapper< TypeTag, TTag::WaterAirBaseProblem > |
| struct | Opm::Properties::PreconditionerWrapper< TypeTag, TTag::WaterAirBaseProblem > |
| class | Opm::WaterAirProblem< TypeTag > |
| Non-isothermal gas injection problem where a air is injected into a fully water saturated medium. More... | |
Namespaces | |
| namespace | Opm |
| namespace | Opm::Properties |
| namespace | Opm::Properties::TTag |
| The generic type tag for problems using the immiscible multi-phase model. | |
Detailed Description
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.