20 #ifndef OPM_COMPRESSIBLETPFA_HEADER_INCLUDED 21 #define OPM_COMPRESSIBLETPFA_HEADER_INCLUDED 26 struct UnstructuredGrid;
35 class BlackoilPropertiesInterface;
36 class RockCompressibility;
37 class LinearSolverInterface;
68 const double residual_tol,
69 const double change_tol,
71 const double* gravity,
80 void solve(
const double dt,
93 virtual void computePerSolveDynamicData(
const double dt,
96 void computePerIterationDynamicData(
const double dt,
99 virtual void computeCellDynamicData(
const double dt,
102 void computeFaceDynamicData(
const double dt,
105 void computeWellDynamicData(
const double dt,
108 void assemble(
const double dt,
111 void solveIncrement();
112 double residualNorm()
const;
113 double incrementNorm()
const;
120 const UnstructuredGrid& grid_;
124 const double residual_tol_;
125 const double change_tol_;
127 const double* gravity_;
129 std::vector<double> htrans_;
130 std::vector<double> trans_ ;
131 std::vector<int> allcells_;
137 std::vector<double> wellperf_wdp_;
138 std::vector<double> initial_porevol_;
141 std::vector<double> cell_A_;
142 std::vector<double> cell_dA_;
143 std::vector<double> cell_viscosity_;
144 std::vector<double> cell_phasemob_;
145 std::vector<double> cell_voldisc_;
146 std::vector<double> face_A_;
147 std::vector<double> face_phasemob_;
148 std::vector<double> face_gravcap_;
149 std::vector<double> wellperf_A_;
150 std::vector<double> wellperf_phasemob_;
151 std::vector<double> porevol_;
152 std::vector<double> rock_comp_;
154 std::vector<double> pressure_increment_;
165 #endif // OPM_COMPRESSIBLETPFA_HEADER_INCLUDED Result structure that presents the fully assembled system of linear equations, linearised around the ...
Definition: cfs_tpfa_residual.h:98
bool singularPressure() const
After solve(), was the resulting pressure singular.
Definition: CompressibleTpfa.cpp:201
void solve(const double dt, BlackoilState &state, WellState &well_state)
Solve the pressure equation by Newton-Raphson scheme.
Definition: CompressibleTpfa.cpp:122
CompressibleTpfa(const UnstructuredGrid &grid, const BlackoilPropertiesInterface &props, const RockCompressibility *rock_comp_props, const LinearSolverInterface &linsolver, const double residual_tol, const double change_tol, const int maxiter, const double *gravity, const Wells *wells)
Construct solver.
Definition: CompressibleTpfa.cpp:62
Abstract interface for linear solvers.
Definition: LinearSolverInterface.hpp:32
virtual ~CompressibleTpfa()
Destructor.
Definition: CompressibleTpfa.cpp:113
Definition: AnisotropicEikonal.cpp:446
Data structure aggregating static information about all wells in a scenario.
Definition: wells.h:50
Definition: RockCompressibility.hpp:32
Simulator state for a blackoil simulator.
Definition: BlackoilState.hpp:40
void computeWellPotentials(const BlackoilState &state)
Compute well potentials.
Definition: CompressibleTpfa.cpp:211
Abstract base class for blackoil fluid and reservoir properties.
Definition: BlackoilPropertiesInterface.hpp:37
Encapsulating a tpfa pressure solver for the compressible-fluid case.
Definition: CompressibleTpfa.hpp:44
The state of a set of wells.
Definition: WellState.hpp:39