mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
127 lines
4.1 KiB
C++
127 lines
4.1 KiB
C++
/**/
|
|
|
|
#ifndef OPM_FULLYIMPLICITTWOPHASESOLVER_HEADER_INCLUDED
|
|
#define OPM_FULLYIMPLICITTWOPHASESOLVER_HEADER_INCLUDED
|
|
|
|
#include <opm/polymer/fullyimplicit/AutoDiffBlock.hpp>
|
|
#include <opm/polymer//fullyimplicit/AutoDiffHelpers.hpp>
|
|
#include <opm/polymer/fullyimplicit/IncompPropsAdInterface.hpp>
|
|
#include <opm/core/pressure/tpfa/trans_tpfa.h>
|
|
#include <opm/core/wells/WellsManager.hpp>
|
|
#include <opm/core/simulator/WellState.hpp>
|
|
|
|
struct UnstructuredGrid;
|
|
namespace Opm {
|
|
// struct HelperOps;
|
|
class LinearSolverInterface;
|
|
class TwophaseState;
|
|
|
|
|
|
class FullyImplicitTwoPhaseSolver
|
|
{
|
|
public:
|
|
FullyImplicitTwoPhaseSolver(const UnstructuredGrid& grid,
|
|
const IncompPropsAdInterface& fluid,
|
|
const Wells& wells,
|
|
const LinearSolverInterface& linsolver,
|
|
const double* gravity);
|
|
|
|
void step(const double dt,
|
|
TwophaseState& state,
|
|
const std::vector<double>& src,
|
|
WellState& wstate);
|
|
private:
|
|
typedef AutoDiffBlock<double> ADB;
|
|
typedef ADB::V V;
|
|
typedef ADB::M M;
|
|
typedef Eigen::Array<double,
|
|
Eigen::Dynamic,
|
|
Eigen::Dynamic,
|
|
Eigen::RowMajor> DataBlock;
|
|
struct SolutionState {
|
|
SolutionState(const int np);
|
|
ADB pressure;
|
|
std::vector<ADB> saturation;
|
|
ADB bhp;
|
|
};
|
|
/*
|
|
struct Source {
|
|
Wells& wells;
|
|
std::vector<double> src;
|
|
} source;
|
|
*/
|
|
struct WellOps {
|
|
WellOps(const Wells& wells);
|
|
M w2p; // well->perf
|
|
M p2w; // perf->well
|
|
};
|
|
|
|
const UnstructuredGrid& grid_;
|
|
const IncompPropsAdInterface& fluid_;
|
|
const Wells& wells_;
|
|
const LinearSolverInterface& linsolver_;
|
|
const double* grav_;
|
|
const std::vector<int> cells_;
|
|
HelperOps ops_;
|
|
const WellOps wops_;
|
|
|
|
std::vector<ADB> mob_;
|
|
|
|
struct {
|
|
std::vector<ADB> mass_balance;
|
|
ADB well_flux_eq;
|
|
ADB well_eq;
|
|
} residual_;
|
|
|
|
|
|
SolutionState
|
|
constantState(const TwophaseState& x,
|
|
const WellState& xw);
|
|
SolutionState
|
|
variableState(const TwophaseState& x,
|
|
const WellState& xw);
|
|
void
|
|
assemble(const V& pvdt,
|
|
const SolutionState& old_state,
|
|
const TwophaseState& x,
|
|
const WellState& xw,
|
|
const std::vector<double>& src);
|
|
V solveJacobianSystem() const;
|
|
void updateState(const V& dx,
|
|
TwophaseState& x,
|
|
WellState& xw)const;
|
|
std::vector<ADB>
|
|
computeRelPerm(const SolutionState& state) const;
|
|
V
|
|
transmissibility() const;
|
|
ADB
|
|
computeFracFlow(const int phase) const;
|
|
ADB
|
|
accumSource(const int phase,
|
|
const std::vector<ADB>& kr,
|
|
const std::vector<double>& src) const;
|
|
ADB
|
|
computeMassFlux(const int phase,
|
|
const V& trans,
|
|
const std::vector<ADB>& kr,
|
|
const SolutionState& state);
|
|
double
|
|
residualNorm() const;
|
|
|
|
ADB
|
|
rockPorosity(const ADB& p) const;
|
|
ADB
|
|
rockPermeability(const ADB& p) const;
|
|
const double
|
|
fluidDensity(const int phase) const;
|
|
ADB
|
|
transMult(const ADB& p) const;
|
|
|
|
ADB
|
|
fluidDensity(const int phase,
|
|
const ADB& p) const;
|
|
const double* gravity() const { return grav_; }
|
|
};
|
|
} // namespace Opm
|
|
#endif// OPM_FULLYIMPLICITTWOPHASESOLVER_HEADER_INCLUDED
|