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

Refined test to have multi-perforation wells.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen
2013-05-22 08:34:50 +02:00
parent 89bd40f81e
commit 476220c84b
1 changed files with 7 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
examples/test_impestpfa_ad.cpp
View File

@@ -108,7 +108,7 @@ int
main(int argc, char* argv[]) main(int argc, char* argv[])
{ {
const Opm::parameter::ParameterGroup param(argc, argv, false); const Opm::parameter::ParameterGroup param(argc, argv, false);
const Opm::GridManager gm(20, 1); const Opm::GridManager gm(5, 5);
const UnstructuredGrid* g = gm.c_grid(); const UnstructuredGrid* g = gm.c_grid();
const int nc = g->number_of_cells; const int nc = g->number_of_cells;
@@ -121,14 +121,14 @@ main(int argc, char* argv[])
typedef Opm::BlackoilPropsAd BOFluid; typedef Opm::BlackoilPropsAd BOFluid;
typedef Opm::ImpesTPFAAD<BOFluid, GeoProps> PSolver; typedef Opm::ImpesTPFAAD<BOFluid, GeoProps> PSolver;
Wells* wells = create_wells(2, 2, 2); Wells* wells = create_wells(2, 2, 5);
const double inj_frac[] = { 1.0, 0.0 }; const double inj_frac[] = { 1.0, 0.0 };
const double prod_frac[] = { 0.0, 0.0 }; const double prod_frac[] = { 0.0, 0.0 };
const int inj_cell = 0; const int inj_cells[] = { 0, 1, 2 };
const int prod_cell = g->number_of_cells - 1; const int prod_cells[] = { 20, 21 };
const double WI = 1e-14; const double WI[3] = { 1e-14 };
bool ok = add_well(INJECTOR, 0.0, 1, inj_frac, &inj_cell, &WI, "Inj", wells); bool ok = add_well(INJECTOR, 0.0, 1, inj_frac, inj_cells, WI, "Inj", wells);
ok = ok && add_well(PRODUCER, 0.0, 1, prod_frac, &prod_cell, &WI, "Prod", wells); ok = ok && add_well(PRODUCER, 0.0, 1, prod_frac, prod_cells, WI, "Prod", wells);
ok = ok && append_well_controls(BHP, 500.0*Opm::unit::barsa, 0, 0, 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); // ok = ok && append_well_controls(BHP, 200.0*Opm::unit::barsa, 0, 1, wells);
double oildistr[2] = { 0.0, 1.0 }; double oildistr[2] = { 0.0, 1.0 };