OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

add well rate contribution to polymer mass balance equation.

Browse Source
This commit is contained in:
Liu Ming
2014-09-30 15:22:50 +08:00
parent 9db24a22d2
commit 9beaf4e03b
2 changed files with 24 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
opm/polymer/fullyimplicit/FullyImplicitBlackoilPolymerSolver.hpp
View File

@@ -102,7 +102,8 @@ namespace Opm {
void
step(const double dt ,
PolymerBlackoilState& state ,
WellStateFullyImplicitBlackoil& wstate);
WellStateFullyImplicitBlackoil& wstate,
const std::vector<double>& polymer_inflow);
private:
// Types and enums
@@ -211,7 +212,8 @@ namespace Opm {
void
addWellEq(const SolutionState& state,
WellStateFullyImplicitBlackoil& xw,
V& aliveWells);
V& aliveWells,
const std::vector<double>& polymer_inflow);
void updateWellControls(ADB& bhp,
ADB& well_phase_flow_rate,
@@ -220,7 +222,8 @@ namespace Opm {
void
assemble(const V& dtpv,
const PolymerBlackoilState& x,
WellStateFullyImplicitBlackoil& xw);
WellStateFullyImplicitBlackoil& xw,
const std::vector<double>& polymer_inflow);
V solveJacobianSystem() const;

24
opm/polymer/fullyimplicit/FullyImplicitBlackoilPolymerSolver_impl.hpp
View File

@@ -649,7 +649,7 @@ namespace {
template<class T>
void FullyImplicitBlackoilPolymerSolver<T>::computeWellConnectionPressures(const SolutionState& state,
const WellStateFullyImplicitBlackoil& xw)
const WellStateFullyImplicitBlackoil& xw)
{
using namespace Opm::AutoDiffGrid;
// 1. Compute properties required by computeConnectionPressureDelta().
@@ -725,7 +725,8 @@ namespace {
FullyImplicitBlackoilPolymerSolver<T>::
assemble(const V& pvdt,
const BlackoilState& x ,
WellStateFullyImplicitBlackoil& xw )
WellStateFullyImplicitBlackoil& xw,
const std::vector<double>& polymer_inflow)
{
using namespace Opm::AutoDiffGrid;
// Create the primary variables.
@@ -796,7 +797,7 @@ namespace {
const Opm::PhaseUsage& pu = fluid_.phaseUsage();
// Add polymer equation.
if (has_polymer) {
if (has_polymer_) {
residula_.material_balance_eq[poly_pos_] = pvdt*(rq_[poly_pos_].accum[1] - rq_[poly_pos_].accum[0])
+ ops_.div*rq_[poly_pos_].mflux;
}
@@ -805,7 +806,7 @@ namespace {
// a well control is switched.
updateWellControls(state.bhp, state.qs, xw);
V aliveWells;
addWellEq(state, xw, aliveWells);
addWellEq(state, xw, aliveWells, polymer_inflow);
addWellControlEq(state, xw, aliveWells);
}
@@ -815,8 +816,9 @@ namespace {
template <class T>
void FullyImplicitBlackoilPolymerSolver<T>::addWellEq(const SolutionState& state,
WellStateFullyImplicitBlackoil& xw,
V& aliveWells)
WellStateFullyImplicitBlackoil& xw,
V& aliveWells,
const std::vector<double>& polymer_inflow)
{
const int nc = Opm::AutoDiffGrid::numCells(grid_);
const int np = wells_.number_of_phases;
@@ -998,6 +1000,16 @@ namespace {
residual_.material_balance_eq[phase] -= superset(cq_s[phase],well_cells,nc);
}
// Add well contributions to polymer mass_balance equation
if (has_polmer_) {
const ADB mc = computeMc(state);
const V polyin = Eigen::Map<const V>(&polymer_inflow[0], nc);
const V poly_in_perf = subset(polyin, well_cells);
const V poly_mc_perf = subset(mc, well_cells).value();
residual_.material_balance_eq[poly_pos_] += superset(cq_ps[pu.phase_pos[Water]]*poly_mc_perf
+ (wops_.w2p * mix_s[pu.phase_pos[Water]])*cqt_is*poly_in_perf,
well_cells, nc);
}
// Add WELL EQUATIONS
ADB qs = state.qs;