22#ifndef OPM_RATECONVERTER_HPP_HEADER_INCLUDED
23#define OPM_RATECONVERTER_HPP_HEADER_INCLUDED
25#include <opm/grid/utility/RegionMapping.hpp>
32#include <dune/grid/common/gridenums.hh>
33#include <dune/grid/common/rangegenerators.hh>
37#include <unordered_map>
52 namespace RateConverter {
69 template <
class Flu
idSystem,
class Region>
72 using Scalar =
typename FluidSystem::Scalar;
81 , attr_ (rmap_, Attributes())
92 template <
typename ElementContext,
class Simulator>
97 for (
const auto& reg : rmap_.activeRegions()) {
98 auto& ra = attr_.attributes(reg);
100 ra.temperature = 0.0;
104 ra.saltConcentration = 0.0;
110 std::unordered_map<RegionId, Attributes> attributes_pv;
113 std::unordered_map<RegionId, Attributes> attributes_hpv;
115 for (
const auto& reg : rmap_.activeRegions()) {
116 attributes_pv.insert({reg, Attributes()});
117 attributes_hpv.insert({reg, Attributes()});
120 ElementContext elemCtx( simulator );
121 const auto& gridView = simulator.
gridView();
122 const auto& comm = gridView.comm();
125 for (
const auto& elem : elements(gridView, Dune::Partitions::interior)) {
126 elemCtx.updatePrimaryStencil(elem);
127 elemCtx.updatePrimaryIntensiveQuantities(0);
128 const unsigned cellIdx = elemCtx.globalSpaceIndex(0, 0);
129 const auto& intQuants = elemCtx.intensiveQuantities(0, 0);
130 const auto& fs = intQuants.fluidState();
133 simulator.
model().dofTotalVolume(simulator.
vanguard().gridEquilIdxToGridIdx(cellIdx))
134 * intQuants.porosity().value();
138 if (FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx)) {
139 hydrocarbon -= fs.saturation(FluidSystem::waterPhaseIdx).value();
142 const int reg = rmap_.region(cellIdx);
146 const Scalar hydrocarbonPV = pv_cell*hydrocarbon;
147 if (hydrocarbonPV > 0.) {
148 auto& attr = attributes_hpv[reg];
149 attr.pv += hydrocarbonPV;
150 if (FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) && FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) {
151 attr.rs += fs.Rs().value() * hydrocarbonPV;
152 attr.rv += fs.Rv().value() * hydrocarbonPV;
154 if (FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx)) {
155 attr.pressure += fs.pressure(FluidSystem::oilPhaseIdx).value() * hydrocarbonPV;
156 attr.temperature += fs.temperature(FluidSystem::oilPhaseIdx).value() * hydrocarbonPV;
158 assert(FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx));
159 attr.pressure += fs.pressure(FluidSystem::gasPhaseIdx).value() * hydrocarbonPV;
160 attr.temperature += fs.temperature(FluidSystem::gasPhaseIdx).value() * hydrocarbonPV;
162 attr.saltConcentration += fs.saltConcentration().value() * hydrocarbonPV;
163 if (FluidSystem::enableDissolvedGasInWater()) {
164 attr.rsw += fs.Rsw().value() * hydrocarbonPV;
166 if (FluidSystem::enableVaporizedWater()) {
167 attr.rvw += fs.Rvw().value() * hydrocarbonPV;
172 auto& attr = attributes_pv[reg];
174 if (FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) && FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) {
175 attr.rs += fs.Rs().value() * pv_cell;
176 attr.rv += fs.Rv().value() * pv_cell;
178 if (FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx)) {
179 attr.pressure += fs.pressure(FluidSystem::oilPhaseIdx).value() * pv_cell;
180 attr.temperature += fs.temperature(FluidSystem::oilPhaseIdx).value() * pv_cell;
181 }
else if (FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) {
182 attr.pressure += fs.pressure(FluidSystem::gasPhaseIdx).value() * pv_cell;
183 attr.temperature += fs.temperature(FluidSystem::gasPhaseIdx).value() * pv_cell;
185 assert(FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx));
186 attr.pressure += fs.pressure(FluidSystem::waterPhaseIdx).value() * pv_cell;
187 attr.temperature += fs.temperature(FluidSystem::waterPhaseIdx).value() * pv_cell;
189 attr.saltConcentration += fs.saltConcentration().value() * pv_cell;
190 if (FluidSystem::enableDissolvedGasInWater()) {
191 attr.rsw += fs.Rsw().value() * pv_cell;
193 if (FluidSystem::enableVaporizedWater()) {
194 attr.rvw += fs.Rvw().value() * pv_cell;
201 this->sumRates(attributes_hpv,
240 template <
class Coeff>
244 template <
class Coeff ,
class Rates>
246 calcCoeff(
const RegionId r,
const int pvtRegionIdx,
const Rates& surface_rates, Coeff& coeff)
const;
248 template <
class Coeff>
257 const Scalar saltConcentration,
260 template <
class Coeff>
286 template <
class Rates>
288 const int pvtRegionIdx,
289 const Rates& surface_rates,
290 Rates& voidage_rates)
const;
326 template <
typename SurfaceRates,
typename Vo
idageRates>
334 const Scalar saltConcentration,
335 const SurfaceRates& surface_rates,
336 VoidageRates& voidage_rates)
const;
338 template <
class Rates>
339 std::pair<Scalar, Scalar>
342 const Rates& surface_rates)
const;
356 template <
class SolventModule>
360 const auto& ra = attr_.attributes(r);
361 const Scalar p = ra.pressure;
362 const Scalar T = ra.temperature;
363 const auto& solventPvt = SolventModule::solventPvt();
364 const Scalar bs = solventPvt.inverseFormationVolumeFactor(pvtRegionIdx, T, p);
372 const RegionMapping<Region> rmap_;
381 , temperature(data[1])
387 , saltConcentration(data[7])
390 Attributes(
const Attributes& rhs)
393 this->data = rhs.data;
396 Attributes& operator=(
const Attributes& rhs)
398 this->data = rhs.data;
402 std::array<Scalar,8> data;
410 Scalar& saltConcentration;
413 void sumRates(std::unordered_map<RegionId,Attributes>& attributes_hpv,
414 std::unordered_map<RegionId,Attributes>& attributes_pv,
417 RegionAttributeHelpers::RegionAttributes<RegionId, Attributes> attr_;
#define OPM_END_PARALLEL_TRY_CATCH(prefix, comm)
Catch exception and throw in a parallel try-catch clause.
Definition: DeferredLoggingErrorHelpers.hpp:192
#define OPM_BEGIN_PARALLEL_TRY_CATCH()
Macro to setup the try of a parallel try-catch.
Definition: DeferredLoggingErrorHelpers.hpp:158
Definition: RateConverter.hpp:70
std::pair< Scalar, Scalar > inferDissolvedVaporisedRatio(const Scalar rsMax, const Scalar rvMax, const Rates &surface_rates) const
SurfaceToReservoirVoidage(const Region ®ion)
Definition: RateConverter.hpp:79
void calcReservoirVoidageRates(const RegionId r, const int pvtRegionIdx, const Rates &surface_rates, Rates &voidage_rates) const
void calcCoeff(const RegionId r, const int pvtRegionIdx, const Rates &surface_rates, Coeff &coeff) const
void calcInjCoeff(const RegionId r, const int pvtRegionIdx, Coeff &coeff) const
void calcCoeff(const RegionId r, const int pvtRegionIdx, Coeff &coeff) const
void calcCoeff(const int pvtRegionIdx, const Scalar p, const Scalar rs, const Scalar rv, const Scalar rsw, const Scalar rvw, const Scalar T, const Scalar saltConcentration, Coeff &coeff) const
void defineState(const Simulator &simulator)
Definition: RateConverter.hpp:93
typename RegionMapping< Region >::RegionId RegionId
Definition: RateConverter.hpp:211
void calcReservoirVoidageRates(const int pvtRegionIdx, const Scalar p, const Scalar rs, const Scalar rv, const Scalar rsw, const Scalar rvw, const Scalar T, const Scalar saltConcentration, const SurfaceRates &surface_rates, VoidageRates &voidage_rates) const
void calcCoeffSolvent(const RegionId r, const int pvtRegionIdx, Scalar &coeff) const
Definition: RateConverter.hpp:358
typename FluidSystem::Scalar Scalar
Definition: RateConverter.hpp:72
Manages the initializing and running of time dependent problems.
Definition: simulator.hh:84
Vanguard & vanguard()
Return a reference to the grid manager of simulation.
Definition: simulator.hh:234
const GridView & gridView() const
Return the grid view for which the simulation is done.
Definition: simulator.hh:246
Model & model()
Return the physical model used in the simulation.
Definition: simulator.hh:252
Dune::Communication< MPIComm > Communication
Definition: ParallelCommunication.hpp:30
Definition: blackoilboundaryratevector.hh:39
int RegionId
Definition: WellConstraints.hpp:37