mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
Add support for PBDV and PDDV in initstateequil.hh
This commit is contained in:
Tor Harald Sandve
parent
b8e714d60a
commit
46fc4b742b
152
tests/data/equil_pbvd_and_pdvd.DATA
Normal file
152
tests/data/equil_pbvd_and_pdvd.DATA
Normal file
@ -0,0 +1,152 @@
|
||||
NOECHO
|
||||
|
||||
RUNSPEC ======
|
||||
|
||||
WATER
|
||||
OIL
|
||||
GAS
|
||||
DISGAS
|
||||
VAPOIL
|
||||
|
||||
TABDIMS
|
||||
1 1 40 20 1 20 /
|
||||
|
||||
DIMENS
|
||||
1 1 20
|
||||
/
|
||||
|
||||
WELLDIMS
|
||||
30 10 2 30 /
|
||||
|
||||
START
|
||||
1 'JAN' 1990 /
|
||||
|
||||
NSTACK
|
||||
25 /
|
||||
|
||||
EQLDIMS
|
||||
-- NTEQUL
|
||||
1 /
|
||||
|
||||
|
||||
FMTOUT
|
||||
FMTIN
|
||||
|
||||
GRID ======
|
||||
|
||||
DX
|
||||
20*1.0
|
||||
/
|
||||
|
||||
DY
|
||||
20*1.0
|
||||
/
|
||||
|
||||
DZ
|
||||
20*5.0
|
||||
/
|
||||
|
||||
|
||||
TOPS
|
||||
0.0
|
||||
/
|
||||
|
||||
PORO
|
||||
20*0.2 /
|
||||
|
||||
|
||||
PERMZ
|
||||
20*1.0
|
||||
/
|
||||
|
||||
PERMY
|
||||
20*100.0
|
||||
/
|
||||
|
||||
PERMX
|
||||
20*100.0
|
||||
/
|
||||
|
||||
BOX
|
||||
1 1 1 1 1 1 /
|
||||
|
||||
|
||||
|
||||
PROPS ======
|
||||
|
||||
PVTO
|
||||
-- Rs Pbub Bo Vo
|
||||
0 1. 1.0000 1.20 /
|
||||
20 40. 1.0120 1.17 /
|
||||
40 80. 1.0255 1.14 /
|
||||
60 120. 1.0380 1.11 /
|
||||
80 160. 1.0510 1.08 /
|
||||
100 200. 1.0630 1.06 /
|
||||
120 240. 1.0750 1.03 /
|
||||
140 280. 1.0870 1.00 /
|
||||
160 320. 1.0985 .98 /
|
||||
180 360. 1.1100 .95 /
|
||||
200 400. 1.1200 .94
|
||||
500. 1.1189 .94 /
|
||||
/
|
||||
|
||||
PVTG
|
||||
-- Pg Rv Bg Vg
|
||||
100 0.0001 0.010 0.1
|
||||
0.0 0.0104 0.1 /
|
||||
200 0.0004 0.005 0.2
|
||||
0.0 0.0054 0.2 /
|
||||
/
|
||||
|
||||
SWOF
|
||||
0.2 0 1 0.9
|
||||
1 1 0 0.1
|
||||
/
|
||||
|
||||
SGOF
|
||||
0 0 1 0.2
|
||||
0.8 1 0 0.5
|
||||
/
|
||||
|
||||
PVTW
|
||||
--RefPres Bw Comp Vw Cv
|
||||
1. 1.0 4.0E-5 0.96 0.0 /
|
||||
|
||||
|
||||
ROCK
|
||||
--RefPres Comp
|
||||
1. 5.0E-5 /
|
||||
|
||||
DENSITY
|
||||
700 1000 1
|
||||
/
|
||||
|
||||
SOLUTION ======
|
||||
|
||||
EQUIL
|
||||
45 150 50 0.25 45 0.35 1 1 0
|
||||
/
|
||||
|
||||
PBVD
|
||||
0 1.0
|
||||
50 150. /
|
||||
|
||||
PDVD
|
||||
50. 150.
|
||||
100. 100 /
|
||||
|
||||
RPTSOL
|
||||
'PRES' 'PGAS' 'PWAT' 'SOIL' 'SWAT' 'SGAS' 'RS' 'RESTART=2' /
|
||||
|
||||
SUMMARY ======
|
||||
RUNSUM
|
||||
|
||||
SEPARATE
|
||||
|
||||
SCHEDULE ======
|
||||
|
||||
RPTSCHED
|
||||
'PRES' 'PGAS' 'PWAT' 'SOIL' 'SWAT' 'SGAS' 'RS' 'RESTART=3' 'NEWTON=2' /
|
||||
|
||||
|
||||
END
|
||||
@ -849,6 +849,97 @@ void test_DeckWithRSVDAndRVVD()
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
void test_DeckWithPBVDAndPDVD()
|
||||
{
|
||||
typedef typename TTAG(TestEquilTypeTag) TypeTag;
|
||||
auto simulator = initSimulator<TypeTag>("data/equil_pbvd_and_pdvd.DATA");
|
||||
const auto& eclipseState = simulator->gridManager().eclState();
|
||||
Opm::GridManager gm(eclipseState.getInputGrid());
|
||||
const UnstructuredGrid& grid = *(gm.c_grid());
|
||||
|
||||
Ewoms::EQUIL::DeckDependent::InitialStateComputer<TypeTag> comp(*simulator->problem().materialLawManager(), eclipseState, simulator->gridManager().grid(), 9.80665);
|
||||
const auto& pressures = comp.press();
|
||||
REQUIRE(pressures.size() == 3);
|
||||
REQUIRE(int(pressures[0].size()) == grid.number_of_cells);
|
||||
|
||||
const int first = 0, last = grid.number_of_cells - 1;
|
||||
// The relative tolerance is too loose to be very useful,
|
||||
// but the answer we are checking is the result of an ODE
|
||||
// solver, and it is unclear if we should check it against
|
||||
// the true answer or something else.
|
||||
const double reltol = 1.0e-6;
|
||||
CHECK_CLOSE(pressures[0][first], 14821552, reltol);
|
||||
CHECK_CLOSE(pressures[0][last], 15479828, reltol);
|
||||
CHECK_CLOSE(pressures[1][first], 14911552, reltol);
|
||||
CHECK_CLOSE(pressures[1][last], 15489828, reltol);
|
||||
|
||||
const auto& sats = comp.saturation();
|
||||
// std::cout << "Saturations:\n";
|
||||
// for (const auto& sat : sats) {
|
||||
// for (const double s : sat) {
|
||||
// std::cout << s << ' ';
|
||||
// }
|
||||
// std::cout << std::endl;
|
||||
// }
|
||||
|
||||
const std::vector<double> s_opm[3]{ // opm
|
||||
{ 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2426402656423233, 0.5383705390740118, 0.7844998821510003, 0.9152832369551807, 1, 1, 1, 1, 1, 1, 1, 1, 1 },
|
||||
{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.1817779931230221, 0.08471676304481934, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
|
||||
{ 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.7573597343576767, 0.4616294609259882, 0.03372212472597758, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 }
|
||||
};
|
||||
|
||||
for (int phase = 0; phase < 3; ++phase) {
|
||||
REQUIRE(sats[phase].size() == s_opm[phase].size());
|
||||
for (size_t i = 0; i < s_opm[phase].size(); ++i) {
|
||||
//std::cout << std::setprecision(10) << sats[phase][i] << '\n';
|
||||
CHECK_CLOSE(sats[phase][i], s_opm[phase][i], reltol);
|
||||
}
|
||||
}
|
||||
|
||||
const auto& rs = comp.rs();
|
||||
const std::vector<double> rs_opm { // opm
|
||||
74.55776480956456,
|
||||
74.6008507125663,
|
||||
74.6439680789467,
|
||||
74.68711693934459,
|
||||
74.73029732443825,
|
||||
74.77350926494491,
|
||||
74.81675279162118,
|
||||
74.86802321984302,
|
||||
74.96677993174352,
|
||||
75.09034523640406,
|
||||
75, 75, 75,75,75, 75, 75, 75, 75, 75 };
|
||||
|
||||
const auto& rv = comp.rv();
|
||||
const std::vector<double> rv_opm {
|
||||
0.0002488465888573874,
|
||||
0.0002491051042753978,
|
||||
0.0002493638084736803,
|
||||
0.0002496227016360676,
|
||||
0.0002498817839466295,
|
||||
0.00025,
|
||||
0.00025,
|
||||
0.00025,
|
||||
0.00025,
|
||||
0.000251180039180951,
|
||||
0.0002522295187440788,
|
||||
0.0002275000000000001,
|
||||
0.0002125,
|
||||
0.0001975,
|
||||
0.0001825,
|
||||
0.0001675,
|
||||
0.0001525,
|
||||
0.0001375,
|
||||
0.0001225,
|
||||
0.0001075};
|
||||
|
||||
for (size_t i = 0; i < rv_opm.size(); ++i) {
|
||||
CHECK_CLOSE(rs[i], rs_opm[i], reltol);
|
||||
CHECK_CLOSE(rv[i], rv_opm[i], reltol);
|
||||
}
|
||||
}
|
||||
|
||||
void test_DeckWithSwatinit()
|
||||
{
|
||||
#if 0
|
||||
@ -1016,6 +1107,7 @@ int main(int argc, char** argv)
|
||||
test_DeckWithLiveOil();
|
||||
test_DeckWithLiveGas();
|
||||
test_DeckWithRSVDAndRVVD();
|
||||
test_DeckWithPBVDAndPDVD();
|
||||
//test_DeckWithSwatinit();
|
||||
|
||||
return 0;
|
||||
|
||||
Loading…
Reference in New Issue
Block a user