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

Merge pull request #4370 from totto82/computeCVIR

Browse Source 
compute CVIR and CVPR
This commit is contained in:
Bård Skaflestad 2023-02-22 17:38:33 +01:00 committed by GitHub
commit 4e8381b6f9
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
3 changed files with 29 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

34
opm/simulators/wells/WellInterfaceFluidSystem.cpp
View File

@ -64,7 +64,7 @@ WellInterfaceFluidSystem(const Well& well,
{ {
} }
template<typename FluidSystem> template <typename FluidSystem>
void void
WellInterfaceFluidSystem<FluidSystem>:: WellInterfaceFluidSystem<FluidSystem>::
calculateReservoirRates(SingleWellState& ws) const calculateReservoirRates(SingleWellState& ws) const
@ -72,14 +72,32 @@ calculateReservoirRates(SingleWellState& ws) const
const int fipreg = 0; // not considering the region for now const int fipreg = 0; // not considering the region for now
const int np = number_of_phases_; const int np = number_of_phases_;
std::vector<double> surface_rates(np, 0.0); this->rateConverter_
for (int p = 0; p < np; ++p) { .calcReservoirVoidageRates(fipreg,
surface_rates[p] = ws.surface_rates[p]; this->pvtRegionIdx_,
} ws.surface_rates,
ws.reservoir_rates);
std::vector<double> voidage_rates(np, 0.0); // Compute total connection reservoir rate CVPR/CVIR
rateConverter_.calcReservoirVoidageRates(fipreg, pvtRegionIdx_, surface_rates, voidage_rates); auto& perf_data = ws.perf_data;
ws.reservoir_rates = voidage_rates; const auto num_perf_well = perf_data.size();
const auto& surf_perf_rates = perf_data.phase_rates;
for (auto i = 0*num_perf_well; i < num_perf_well; ++i) {
const auto surface_rates_perf = std::vector<double>
{ surf_perf_rates.begin() + (i + 0)*np ,
surf_perf_rates.begin() + (i + 1)*np };
std::vector<double> voidage_rates_perf(np, 0.0);
this->rateConverter_
.calcReservoirVoidageRates(fipreg,
this->pvtRegionIdx_,
surface_rates_perf,
voidage_rates_perf);
perf_data.rates[i] =
std::accumulate(voidage_rates_perf.begin(),
voidage_rates_perf.end(), 0.0);
}
} }
template <typename FluidSystem> template <typename FluidSystem>

4
parallelTests.cmake
View File

@ -48,7 +48,7 @@ add_test_compare_parallel_simulation(CASENAME spe3
SIMULATOR flow SIMULATOR flow
ABS_TOL ${abs_tol_parallel} ABS_TOL ${abs_tol_parallel}
REL_TOL ${coarse_rel_tol_parallel} REL_TOL ${coarse_rel_tol_parallel}
TEST_ARGS --linear-solver-reduction=1e-7 --tolerance-cnv=5e-6 --tolerance-mb=1e-8) TEST_ARGS --linear-solver-reduction=1e-7 --tolerance-cnv=5e-6 --tolerance-mb=1e-8 --tolerance-wells=1e-7)
add_test_compare_parallel_simulation(CASENAME spe1_solvent add_test_compare_parallel_simulation(CASENAME spe1_solvent
FILENAME SPE1CASE2_SOLVENT FILENAME SPE1CASE2_SOLVENT
@ -92,7 +92,7 @@ add_test_compare_parallel_simulation(CASENAME spe1_water
ABS_TOL ${abs_tol} ABS_TOL ${abs_tol}
REL_TOL ${rel_tol} REL_TOL ${rel_tol}
DIR spe1 DIR spe1
TEST_ARGS --linear-solver-reduction=1e-7 --tolerance-cnv=5e-6 --tolerance-mb=1e-8) TEST_ARGS --linear-solver-reduction=1e-7 --tolerance-cnv=5e-6 --tolerance-mb=1e-8 --tolerance-wells=1e-7)
add_test_compare_parallel_simulation(CASENAME spe1_brine add_test_compare_parallel_simulation(CASENAME spe1_brine
FILENAME SPE1CASE1_BRINE FILENAME SPE1CASE1_BRINE

2
restartTests.cmake
View File

@ -69,7 +69,7 @@ add_test_compare_restarted_simulation(CASENAME msw_3d_hfa
ABS_TOL ${abs_tol_restart_msw} ABS_TOL ${abs_tol_restart_msw}
REL_TOL ${rel_tol_restart_msw} REL_TOL ${rel_tol_restart_msw}
RESTART_STEP 10 RESTART_STEP 10
TEST_ARGS --enable-adaptive-time-stepping=false --sched-restart=true) TEST_ARGS --enable-adaptive-time-stepping=false --sched-restart=true --tolerance-wells=1e-7)
# Basic restart tests which only compare the summary output, this test driver should # Basic restart tests which only compare the summary output, this test driver should