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

Added wells to solver interface, not used for anything yet.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen
2013-05-12 21:43:54 +02:00
parent 2f6f32bf6e
commit b0cf134479
2 changed files with 35 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
ImpesTPFAAD.hpp
View File

@@ -26,6 +26,7 @@
#include <opm/core/simulator/BlackoilState.hpp>
#include <opm/core/simulator/WellState.hpp>
#include <opm/core/wells.h>
#include <algorithm>
#include <cassert>
@@ -96,7 +97,7 @@ namespace Opm {
A_ .data(), dA_ .data());
fluid_.viscosity(nc_, & p[0], & z[0], & cells_[0],
mu_.data(), dmu_.data());
mu_.data(), /*dmu_.data()*/ 0);
}
ADB
@@ -165,9 +166,11 @@ namespace Opm {
public:
ImpesTPFAAD(const UnstructuredGrid& grid ,
const BOFluid& fluid,
const GeoProps& geo )
const GeoProps& geo ,
const Wells& wells)
: grid_ (grid)
, geo_ (geo)
, wells_ (wells)
, pdepfdata_(grid.number_of_cells, fluid)
, ops_ (grid)
{
@@ -175,11 +178,12 @@ namespace Opm {
void
solve(const double dt,
BlackoilState& state)
BlackoilState& state,
WellState& well_state)
{
pdepfdata_.computeSatQuant(state);
assemble(dt, state);
assemble(dt, state, well_state);
}
private:
@@ -192,12 +196,14 @@ namespace Opm {
const UnstructuredGrid& grid_;
const GeoProps& geo_ ;
const Wells& wells_;
PDepFData pdepfdata_;
HelperOps ops_;
void
assemble(const double dt,
const BlackoilState& state)
const BlackoilState& state,
const WellState& well_state)
{
typedef typename ADB::V V;
@@ -209,6 +215,7 @@ namespace Opm {
const V& pv = geo_.poreVolume();
const int nc = grid_.number_of_cells;
const int np = state.numPhases();
const int nw = wells_.number_of_wells;
pdepfdata_.computePressQuant(state);
@@ -218,6 +225,7 @@ namespace Opm {
const V transi = subset(geo_.transmissibility(),
ops_.internal_faces);
const V p0 = Eigen::Map<const V>(&state.pressure()[0], nc, 1);
const V bhp = Eigen::Map<const V>(&well_state.bhp()[0], nw, 1);
const std::vector<int> bpat(1, nc);
ADB p = ADB::variable(0, p0, bpat);

24
test_impestpfa_ad.cpp
View File

@@ -28,9 +28,13 @@
#include <opm/core/pressure/tpfa/trans_tpfa.h>
#include <opm/core/utility/parameters/ParameterGroup.hpp>
#include <opm/core/utility/Units.hpp>
#include <opm/core/simulator/initState.hpp>
#include <opm/core/wells.h>
// #include <opm/core/WellsManager.hpp>
#include <algorithm>
namespace {
@@ -82,14 +86,30 @@ main(int argc, char* argv[])
typedef Opm::BlackoilPropertiesInterface BOFluid;
typedef Opm::ImpesTPFAAD<BOFluid, GeoProps> PSolver;
Wells* wells = create_wells(2, 2, 2);
const double inj_frac[] = { 1.0, 0.0 };
const double prod_frac[] = { 0.0, 0.0 };
const int inj_cell = 0;
const int prod_cell = g->number_of_cells - 1;
const double WI = 1e-8;
bool ok = add_well(INJECTOR, 0.0, 1, inj_frac, &inj_cell, &WI, "Inj", wells);
ok = ok && add_well(PRODUCER, 0.0, 1, prod_frac, &prod_cell, &WI, "Prod", wells);
ok = ok && append_well_controls(BHP, 500.0*Opm::unit::barsa, 0, 0, wells);
ok = ok && append_well_controls(BHP, 200.0*Opm::unit::barsa, 0, 1, wells);
if (!ok) {
THROW("Something went wrong with well init.");
}
GeoProps geo(*g, props);
PSolver ps (*g, props, geo);
PSolver ps (*g, props, geo, *wells);
Opm::BlackoilState state;
initStateBasic(*g, props, param, 0.0, state);
initBlackoilSurfvol(*g, props, state);
Opm::WellState well_state;
well_state.init(wells, state);
ps.solve(1.0, state);
ps.solve(1.0, state, well_state);
return 0;
}