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: move calcInjCoeff to compile unit

Browse Source
This commit is contained in:
Arne Morten Kvarving
2023-03-06 08:49:25 +01:00
parent 7672e14890
commit 11602bbdda
2 changed files with 42 additions and 33 deletions
Download Patch File Download Diff File Expand all files Collapse all files

41
opm/simulators/wells/RateConverter.cpp
View File

@@ -24,6 +24,8 @@
#include <opm/material/fluidsystems/BlackOilDefaultIndexTraits.hpp> #include <opm/material/fluidsystems/BlackOilDefaultIndexTraits.hpp>
#include <opm/material/fluidsystems/BlackOilFluidSystem.hpp> #include <opm/material/fluidsystems/BlackOilFluidSystem.hpp>
#include <algorithm>
namespace Opm { namespace Opm {
namespace RateConverter { namespace RateConverter {
@@ -55,11 +57,50 @@ sumRates(std::unordered_map<RegionId,Attributes>& attributes_hpv,
} }
} }
template <class FluidSystem, class Region>
template <class Coeff>
void SurfaceToReservoirVoidage<FluidSystem,Region>::
calcInjCoeff(const RegionId r, const int pvtRegionIdx, Coeff& coeff) const
{
const auto& pu = phaseUsage_;
const auto& ra = attr_.attributes(r);
const double p = ra.pressure;
const double T = ra.temperature;
const double saltConcentration = ra.saltConcentration;
const int iw = RegionAttributeHelpers::PhasePos::water(pu);
const int io = RegionAttributeHelpers::PhasePos::oil (pu);
const int ig = RegionAttributeHelpers::PhasePos::gas (pu);
std::fill(& coeff[0], & coeff[0] + phaseUsage_.num_phases, 0.0);
if (RegionAttributeHelpers::PhaseUsed::water(pu)) {
// q[w]_r = q[w]_s / bw
const double bw = FluidSystem::waterPvt().inverseFormationVolumeFactor(pvtRegionIdx, T, p, 0.0, saltConcentration);
coeff[iw] = 1.0 / bw;
}
if (RegionAttributeHelpers::PhaseUsed::oil(pu)) {
const double bo = FluidSystem::oilPvt().inverseFormationVolumeFactor(pvtRegionIdx, T, p, 0.0);
coeff[io] += 1.0 / bo;
}
if (RegionAttributeHelpers::PhaseUsed::gas(pu)) {
const double bg = FluidSystem::gasPvt().inverseFormationVolumeFactor(pvtRegionIdx, T, p, 0.0, 0.0);
coeff[ig] += 1.0 / bg;
}
}
using FS = BlackOilFluidSystem<double,BlackOilDefaultIndexTraits>; using FS = BlackOilFluidSystem<double,BlackOilDefaultIndexTraits>;
template void SurfaceToReservoirVoidage<FS,std::vector<int>>:: template void SurfaceToReservoirVoidage<FS,std::vector<int>>::
sumRates(std::unordered_map<int,Attributes>&, sumRates(std::unordered_map<int,Attributes>&,
std::unordered_map<int,Attributes>&, std::unordered_map<int,Attributes>&,
Parallel::Communication); Parallel::Communication);
template void SurfaceToReservoirVoidage<FS,std::vector<int>>::
calcInjCoeff<std::vector<double>>(const int, const int, std::vector<double>&) const;
} // namespace RateConverter } // namespace RateConverter
} // namespace Opm } // namespace Opm

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

@@ -110,7 +110,6 @@ namespace Opm {
ra.saltConcentration = 0.0; ra.saltConcentration = 0.0;
ra.rsw = 0.0; ra.rsw = 0.0;
ra.rvw = 0.0; ra.rvw = 0.0;
} }
// quantities for pore volume average // quantities for pore volume average
@@ -329,38 +328,7 @@ namespace Opm {
template <class Coeff> template <class Coeff>
void void
calcInjCoeff(const RegionId r, const int pvtRegionIdx, Coeff& coeff) const calcInjCoeff(const RegionId r, const int pvtRegionIdx, Coeff& coeff) const;
{
const auto& pu = phaseUsage_;
const auto& ra = attr_.attributes(r);
const double p = ra.pressure;
const double T = ra.temperature;
const double saltConcentration = ra.saltConcentration;
const int iw = RegionAttributeHelpers::PhasePos::water(pu);
const int io = RegionAttributeHelpers::PhasePos::oil (pu);
const int ig = RegionAttributeHelpers::PhasePos::gas (pu);
std::fill(& coeff[0], & coeff[0] + phaseUsage_.num_phases, 0.0);
if (RegionAttributeHelpers::PhaseUsed::water(pu)) {
// q[w]_r = q[w]_s / bw
const double bw = FluidSystem::waterPvt().inverseFormationVolumeFactor(pvtRegionIdx, T, p, 0.0, saltConcentration);
coeff[iw] = 1.0 / bw;
}
if (RegionAttributeHelpers::PhaseUsed::oil(pu)) {
const double bo = FluidSystem::oilPvt().inverseFormationVolumeFactor(pvtRegionIdx, T, p, 0.0);
coeff[io] += 1.0 / bo;
}
if (RegionAttributeHelpers::PhaseUsed::gas(pu)) {
const double bg = FluidSystem::gasPvt().inverseFormationVolumeFactor(pvtRegionIdx, T, p, 0.0, 0.0);
coeff[ig] += 1.0 / bg;
}
}
/** /**
* Convert surface volume flow rates to reservoir voidage flow * Convert surface volume flow rates to reservoir voidage flow