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

RateConverter: coalesce global reductions

Browse Source
This commit is contained in:
Arne Morten Kvarving
2023-03-06 14:53:28 +01:00
parent fc06fc77c3
commit 917af4e28e
1 changed files with 45 additions and 54 deletions
Download Patch File Download Diff File Expand all files Collapse all files

99
opm/simulators/wells/RateConverter.hpp
View File

@@ -31,6 +31,8 @@
#include <dune/grid/common/rangegenerators.hh> #include <dune/grid/common/rangegenerators.hh>
#include <algorithm> #include <algorithm>
#include <array>
#include <cassert>
#include <cmath> #include <cmath>
#include <memory> #include <memory>
#include <stdexcept> #include <stdexcept>
@@ -203,49 +205,24 @@ namespace Opm {
OPM_END_PARALLEL_TRY_CATCH("SurfaceToReservoirVoidage::defineState() failed: ", simulator.vanguard().grid().comm()); OPM_END_PARALLEL_TRY_CATCH("SurfaceToReservoirVoidage::defineState() failed: ", simulator.vanguard().grid().comm());
for (const auto& reg : rmap_.activeRegions()) { for (const auto& reg : rmap_.activeRegions()) {
// Calculating first using the hydrocarbon pore volumes
auto& ra = attr_.attributes(reg); auto& ra = attr_.attributes(reg);
const double hpv_sum = comm.sum(attributes_hpv[reg].pv); const auto& attri_hpv = attributes_hpv[reg];
ra.data = attri_hpv.data;
comm.sum(ra.data.data(), ra.data.size());
// TODO: should we have some epsilon here instead of zero? // TODO: should we have some epsilon here instead of zero?
if (hpv_sum > 0.) { // No hydrocarbon pore volume, redo but this time using full pore volumes.
const auto& attri_hpv = attributes_hpv[reg]; if (ra.pv <= 0.) {
const double p_hpv_sum = comm.sum(attri_hpv.pressure);
const double T_hpv_sum = comm.sum(attri_hpv.temperature);
const double rs_hpv_sum = comm.sum(attri_hpv.rs);
const double rv_hpv_sum = comm.sum(attri_hpv.rv);
const double sc_hpv_sum = comm.sum(attri_hpv.saltConcentration);
const double rsw_hpv_sum = comm.sum(attri_hpv.rsw);
const double rvw_hpv_sum = comm.sum(attri_hpv.rvw);
ra.pressure = p_hpv_sum / hpv_sum;
ra.temperature = T_hpv_sum / hpv_sum;
ra.rs = rs_hpv_sum / hpv_sum;
ra.rv = rv_hpv_sum / hpv_sum;
ra.rsw = rsw_hpv_sum / hpv_sum;
ra.rvw = rvw_hpv_sum / hpv_sum;
ra.pv = hpv_sum;
ra.saltConcentration = sc_hpv_sum / hpv_sum;
} else {
// using the pore volume to do the averaging // using the pore volume to do the averaging
const auto& attri_pv = attributes_pv[reg]; const auto& attri_pv = attributes_pv[reg];
const double pv_sum = comm.sum(attri_pv.pv); ra.data = attri_pv.data;
assert(pv_sum > 0.); comm.sum(ra.data.data(), ra.data.size());
const double p_pv_sum = comm.sum(attri_pv.pressure); assert(ra.pv > 0.);
const double T_pv_sum = comm.sum(attri_pv.temperature);
const double rs_pv_sum = comm.sum(attri_pv.rs);
const double rv_pv_sum = comm.sum(attri_pv.rv);
const double rsw_pv_sum = comm.sum(attri_pv.rsw);
const double rvw_pv_sum = comm.sum(attri_pv.rvw);
const double sc_pv_sum = comm.sum(attri_pv.saltConcentration);
ra.pressure = p_pv_sum / pv_sum;
ra.temperature = T_pv_sum / pv_sum;
ra.rs = rs_pv_sum / pv_sum;
ra.rv = rv_pv_sum / pv_sum;
ra.rsw = rsw_pv_sum / pv_sum;
ra.rvw = rvw_pv_sum / pv_sum;
ra.pv = pv_sum;
ra.saltConcentration = sc_pv_sum / pv_sum;
} }
const double pv_sum = ra.pv;
for (double& d : ra.data)
d /= pv_sum;
ra.pv = pv_sum;
} }
} }
@@ -614,24 +591,38 @@ namespace Opm {
*/ */
struct Attributes { struct Attributes {
Attributes() Attributes()
: pressure (0.0) : data{0.0}
, temperature(0.0) , pressure(data[0])
, rs(0.0) , temperature(data[1])
, rv(0.0) , rs(data[2])
, rsw(0.0) , rv(data[3])
, rvw(0.0) , rsw(data[4])
, pv(0.0) , rvw(data[5])
, saltConcentration(0.0) , pv(data[6])
, saltConcentration(data[7])
{} {}
double pressure; Attributes(const Attributes& rhs)
double temperature; : Attributes()
double rs; {
double rv; this->data = rhs.data;
double rsw; }
double rvw;
double pv; Attributes& operator=(const Attributes& rhs)
double saltConcentration; {
this->data = rhs.data;
return *this;
}
std::array<double,8> data;
double& pressure;
double& temperature;
double& rs;
double& rv;
double& rsw;
double& rvw;
double& pv;
double& saltConcentration;
}; };
RegionAttributeHelpers::RegionAttributes<RegionId, Attributes> attr_; RegionAttributeHelpers::RegionAttributes<RegionId, Attributes> attr_;