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

adding regression tests for filtration cake model

Browse Source
This commit is contained in:
Kai Bao 2023-07-06 16:11:15 +02:00
parent 1981a3bf4d
commit 22269c92c3
3 changed files with 24 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
opm/simulators/wells/WellInterfaceGeneric.cpp
View File

@ -740,7 +740,7 @@ updateFiltrationParticleVolume(const double dt, const size_t water_index,
auto& ws = well_state.well(this->index_of_well_); auto& ws = well_state.well(this->index_of_well_);
const auto& connection_rates = ws.perf_data.phase_rates; const auto& connection_rates = ws.perf_data.phase_rates;
const size_t np = well_state.numPhases(); const std::size_t np = well_state.numPhases();
for (int perf = 0; perf < this->number_of_perforations_; ++perf) { for (int perf = 0; perf < this->number_of_perforations_; ++perf) {
// not considering the production water // not considering the production water
const double water_rates = std::max(0., connection_rates[perf * np + water_index]); const double water_rates = std::max(0., connection_rates[perf * np + water_index]);
@ -768,9 +768,8 @@ updateInjFCMult(const std::vector<double>& filtration_particle_volume) {
const auto cskinfactor = connection.skinFactor(); const auto cskinfactor = connection.skinFactor();
const double K = connection.Kh() / connection.connectionLength(); const double K = connection.Kh() / connection.connectionLength();
const double factor = filter_cake.sf_multiplier; const double factor = filter_cake.sf_multiplier;
// compute a multiplier for the transmissibility // compute a multiplier for the well connection transmissibility
if (filter_cake.geometry == FilterCake::FilterCakeGeometry::LINEAR) { if (filter_cake.geometry == FilterCake::FilterCakeGeometry::LINEAR) {
// but we are using this form just for first prototype
const double skin_factor = thickness / rw * K / perm * factor; const double skin_factor = thickness / rw * K / perm * factor;
const auto denom = std::log(cr0 / std::min(crw, cr0)) + cskinfactor; const auto denom = std::log(cr0 / std::min(crw, cr0)) + cskinfactor;
const auto denom2 = denom + skin_factor; const auto denom2 = denom + skin_factor;

8
parallelTests.cmake
View File

@ -181,6 +181,14 @@ add_test_compare_parallel_simulation(CASENAME winjmult_msw
DIR winjmult DIR winjmult
TEST_ARGS --enable-tuning=true) TEST_ARGS --enable-tuning=true)
add_test_compare_parallel_simulation(CASENAME winjdam_msw
FILENAME WINJDAM_MSW
SIMULATOR flow
ABS_TOL ${abs_tol}
REL_TOL 0.1
DIR winjdam
TEST_ARGS --enable-tuning=true)
add_test_compare_parallel_simulation(CASENAME 3_a_mpi_multflt_mod2 add_test_compare_parallel_simulation(CASENAME 3_a_mpi_multflt_mod2
FILENAME 3_A_MPI_MULTFLT_SCHED_MODEL2 FILENAME 3_A_MPI_MULTFLT_SCHED_MODEL2
SIMULATOR flow SIMULATOR flow

14
regressionTests.cmake
View File

@ -1213,6 +1213,20 @@ add_test_compareECLFiles(CASENAME winjmult_msw
REL_TOL ${rel_tol} REL_TOL ${rel_tol}
DIR winjmult DIR winjmult
TEST_ARGS --enable-tuning=true) TEST_ARGS --enable-tuning=true)
add_test_compareECLFiles(CASENAME winjdam_stdw
FILENAME WINJDAM_STDW
SIMULATOR flow
ABS_TOL ${abs_tol}
REL_TOL ${rel_tol}
DIR winjdam
TEST_ARGS --enable-tuning=true)
add_test_compareECLFiles(CASENAME winjdam_msw
FILENAME WINJDAM_MSW
SIMULATOR flow
ABS_TOL ${abs_tol}
REL_TOL ${rel_tol}
DIR winjdam
TEST_ARGS --enable-tuning=true)
add_test_compareECLFiles(CASENAME 01_vappars add_test_compareECLFiles(CASENAME 01_vappars
FILENAME VAPPARS-01 FILENAME VAPPARS-01
SIMULATOR flow SIMULATOR flow