20 #ifndef OPM_TRANSPORTSOLVERCOMPRESSIBLETWOPHASEREORDER_HEADER_INCLUDED
21 #define OPM_TRANSPORTSOLVERCOMPRESSIBLETWOPHASEREORDER_HEADER_INCLUDED
31 class BlackoilPropertiesInterface;
59 void solve(
const double* darcyflux,
60 const double* pressure,
62 const double* porevolume0,
63 const double* porevolume,
66 std::vector<double>& saturation,
67 std::vector<double>& surfacevol);
82 void solveGravity(
const std::vector<std::vector<int> >& columns,
84 std::vector<double>& saturation,
85 std::vector<double>& surfacevol);
88 virtual void solveSingleCell(
const int cell);
89 virtual void solveMultiCell(
const int num_cells,
const int* cells);
90 void solveSingleCellGravity(
const std::vector<int>& cells,
92 const double* gravflux);
93 int solveGravityColumn(
const std::vector<int>& cells);
94 void initGravityDynamic();
99 std::vector<int> allcells_;
100 std::vector<double> visc_;
101 std::vector<double> A_;
102 std::vector<double> smin_;
103 std::vector<double> smax_;
107 const double* darcyflux_;
108 const double* surfacevol0_;
109 const double* porevolume0_;
110 const double* porevolume_;
111 const double* source_;
113 std::vector<double> saturation_;
114 std::vector<double> fractionalflow_;
116 const double* gravity_;
117 std::vector<double> trans_;
118 std::vector<double> density_;
119 std::vector<double> gravflux_;
120 std::vector<double> mob_;
121 std::vector<double> s0_;
124 std::vector<int> ia_upw_;
125 std::vector<int> ja_upw_;
126 std::vector<int> ia_downw_;
127 std::vector<int> ja_downw_;
130 double fracFlow(
double s,
int cell)
const;
132 struct GravityResidual;
133 void mobility(
double s,
int cell,
double* mob)
const;
138 #endif // OPM_TRANSPORTMODELCOMPRESSIBLETWOPHASE_HEADER_INCLUDED
Definition: AnisotropicEikonal.hpp:43
void solveGravity(const std::vector< std::vector< int > > &columns, const double dt, std::vector< double > &saturation, std::vector< double > &surfacevol)
TransportSolverCompressibleTwophaseReorder(const UnstructuredGrid &grid, const Opm::BlackoilPropertiesInterface &props, const double tol, const int maxit)
Definition: ReorderSolverInterface.hpp:38
Definition: BlackoilPropertiesInterface.hpp:37
Definition: TransportSolverCompressibleTwophaseReorder.hpp:35
void solve(const double *darcyflux, const double *pressure, const double *temperature, const double *porevolume0, const double *porevolume, const double *source, const double dt, std::vector< double > &saturation, std::vector< double > &surfacevol)
const UnstructuredGrid & grid
Definition: ColumnExtract.hpp:31
void initGravity(const double *grav)
std::vector< double > temperature(const Grid &, const Region &, const CellRange &cells)
Definition: initStateEquil_impl.hpp:667