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

BlackoilWellModel: use Scalar type

Browse Source
This commit is contained in:
Arne Morten Kvarving
2024-02-20 22:22:41 +01:00
parent f0e7f8842b
commit f5d6b69703
2 changed files with 47 additions and 47 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
opm/simulators/wells/BlackoilWellModel.hpp
View File

@@ -264,7 +264,7 @@ namespace Opm {
this->assignWellTracerRates(wsrpt, tracerRates); this->assignWellTracerRates(wsrpt, tracerRates);
BlackoilWellModelGuideRates<double>(*this).assignWellGuideRates(wsrpt, this->reportStepIndex()); BlackoilWellModelGuideRates(*this).assignWellGuideRates(wsrpt, this->reportStepIndex());
this->assignShutConnections(wsrpt, this->reportStepIndex()); this->assignShutConnections(wsrpt, this->reportStepIndex());
return wsrpt; return wsrpt;
@@ -333,7 +333,7 @@ namespace Opm {
WellInterfacePtr getWell(const std::string& well_name) const; WellInterfacePtr getWell(const std::string& well_name) const;
bool hasWell(const std::string& well_name) const; bool hasWell(const std::string& well_name) const;
using PressureMatrix = Dune::BCRSMatrix<Opm::MatrixBlock<double, 1, 1>>; using PressureMatrix = Dune::BCRSMatrix<Opm::MatrixBlock<Scalar, 1, 1>>;
void addWellPressureEquations(PressureMatrix& jacobian, const BVector& weights,const bool use_well_weights) const; void addWellPressureEquations(PressureMatrix& jacobian, const BVector& weights,const bool use_well_weights) const;
@@ -355,8 +355,8 @@ namespace Opm {
void updateWellControlsDomain(DeferredLogger& deferred_logger, const Domain& domain); void updateWellControlsDomain(DeferredLogger& deferred_logger, const Domain& domain);
void logPrimaryVars() const; void logPrimaryVars() const;
std::vector<double> getPrimaryVarsDomain(const Domain& domain) const; std::vector<Scalar> getPrimaryVarsDomain(const Domain& domain) const;
void setPrimaryVarsDomain(const Domain& domain, const std::vector<double>& vars); void setPrimaryVarsDomain(const Domain& domain, const std::vector<Scalar>& vars);
void setupDomains(const std::vector<Domain>& domains); void setupDomains(const std::vector<Domain>& domains);
@@ -389,8 +389,8 @@ namespace Opm {
std::size_t global_num_cells_{}; std::size_t global_num_cells_{};
// the number of the cells in the local grid // the number of the cells in the local grid
std::size_t local_num_cells_{}; std::size_t local_num_cells_{};
double gravity_{}; Scalar gravity_{};
std::vector<double> depth_{}; std::vector<Scalar> depth_{};
bool alternative_well_rate_init_{}; bool alternative_well_rate_init_{};
std::unique_ptr<RateConverterType> rateConverter_{}; std::unique_ptr<RateConverterType> rateConverter_{};
@@ -495,7 +495,7 @@ namespace Opm {
ExceptionType::ExcEnum& exc_type, ExceptionType::ExcEnum& exc_type,
DeferredLogger& deferred_logger) override; DeferredLogger& deferred_logger) override;
const std::vector<double>& wellPerfEfficiencyFactors() const; const std::vector<Scalar>& wellPerfEfficiencyFactors() const;
void calculateProductivityIndexValuesShutWells(const int reportStepIdx, DeferredLogger& deferred_logger) override; void calculateProductivityIndexValuesShutWells(const int reportStepIdx, DeferredLogger& deferred_logger) override;
void calculateProductivityIndexValues(DeferredLogger& deferred_logger) override; void calculateProductivityIndexValues(DeferredLogger& deferred_logger) override;
@@ -539,12 +539,12 @@ namespace Opm {
void calcRates(const int fipnum, void calcRates(const int fipnum,
const int pvtreg, const int pvtreg,
const std::vector<double>& production_rates, const std::vector<Scalar>& production_rates,
std::vector<double>& resv_coeff) override; std::vector<Scalar>& resv_coeff) override;
void calcInjRates(const int fipnum, void calcInjRates(const int fipnum,
const int pvtreg, const int pvtreg,
std::vector<double>& resv_coeff) override; std::vector<Scalar>& resv_coeff) override;
void computeWellTemperature(); void computeWellTemperature();

72
opm/simulators/wells/BlackoilWellModel_impl.hpp
View File

@@ -90,9 +90,9 @@ namespace Opm {
.localCellIndex([this](const std::size_t globalIndex) .localCellIndex([this](const std::size_t globalIndex)
{ return this->compressedIndexForInterior(globalIndex); }) { return this->compressedIndexForInterior(globalIndex); })
.evalCellSource([this](const int localCell, .evalCellSource([this](const int localCell,
PAvgDynamicSourceData::SourceDataSpan<double> sourceTerms) PAvgDynamicSourceData::SourceDataSpan<Scalar> sourceTerms)
{ {
using Item = PAvgDynamicSourceData::SourceDataSpan<double>::Item; using Item = typename PAvgDynamicSourceData::SourceDataSpan<Scalar>::Item;
const auto* intQuants = this->simulator_.model() const auto* intQuants = this->simulator_.model()
.cachedIntensiveQuantities(localCell, /*timeIndex = */0); .cachedIntensiveQuantities(localCell, /*timeIndex = */0);
@@ -511,9 +511,9 @@ namespace Opm {
//update guide rates //update guide rates
const auto& comm = simulator_.vanguard().grid().comm(); const auto& comm = simulator_.vanguard().grid().comm();
std::vector<double> pot(this->numPhases(), 0.0); std::vector<Scalar> pot(this->numPhases(), 0.0);
const Group& fieldGroup = this->schedule().getGroup("FIELD", reportStepIdx); const Group& fieldGroup = this->schedule().getGroup("FIELD", reportStepIdx);
WellGroupHelpers<double>::updateGuideRates(fieldGroup, WellGroupHelpers<Scalar>::updateGuideRates(fieldGroup,
this->schedule(), this->schedule(),
summaryState, summaryState,
this->phase_usage_, this->phase_usage_,
@@ -533,7 +533,7 @@ namespace Opm {
calculator.second->template defineState<ElementContext>(simulator_); calculator.second->template defineState<ElementContext>(simulator_);
} }
const double dt = simulator_.timeStepSize(); const double dt = simulator_.timeStepSize();
WellGroupHelpers<double>::updateGpMaintTargetForGroups(fieldGroup, WellGroupHelpers<Scalar>::updateGpMaintTargetForGroups(fieldGroup,
this->schedule_, this->schedule_,
regionalAveragePressureCalculator_, regionalAveragePressureCalculator_,
reportStepIdx, reportStepIdx,
@@ -594,8 +594,8 @@ namespace Opm {
// some preparation before the well can be used // some preparation before the well can be used
well->init(&this->phase_usage_, depth_, gravity_, local_num_cells_, B_avg_, true); well->init(&this->phase_usage_, depth_, gravity_, local_num_cells_, B_avg_, true);
double well_efficiency_factor = wellEcl.getEfficiencyFactor(); Scalar well_efficiency_factor = wellEcl.getEfficiencyFactor();
WellGroupHelpers<double>::accumulateGroupEfficiencyFactor(this->schedule().getGroup(wellEcl.groupName(), WellGroupHelpers<Scalar>::accumulateGroupEfficiencyFactor(this->schedule().getGroup(wellEcl.groupName(),
timeStepIdx), timeStepIdx),
this->schedule(), this->schedule(),
timeStepIdx, timeStepIdx,
@@ -799,7 +799,7 @@ namespace Opm {
BlackoilWellModel<TypeTag>:: BlackoilWellModel<TypeTag>::
initializeWellState(const int timeStepIdx) initializeWellState(const int timeStepIdx)
{ {
std::vector<double> cellPressures(this->local_num_cells_, 0.0); std::vector<Scalar> cellPressures(this->local_num_cells_, 0.0);
ElementContext elemCtx(simulator_); ElementContext elemCtx(simulator_);
const auto& gridView = simulator_.vanguard().gridView(); const auto& gridView = simulator_.vanguard().gridView();
@@ -811,7 +811,7 @@ namespace Opm {
const auto& fs = elemCtx.intensiveQuantities(/*spaceIdx=*/0, /*timeIdx=*/0).fluidState(); const auto& fs = elemCtx.intensiveQuantities(/*spaceIdx=*/0, /*timeIdx=*/0).fluidState();
// copy of get perfpressure in Standard well except for value // copy of get perfpressure in Standard well except for value
double& perf_pressure = cellPressures[elemCtx.globalSpaceIndex(/*spaceIdx=*/0, /*timeIdx=*/0)]; Scalar& perf_pressure = cellPressures[elemCtx.globalSpaceIndex(/*spaceIdx=*/0, /*timeIdx=*/0)];
if (Indices::oilEnabled) { if (Indices::oilEnabled) {
perf_pressure = fs.pressure(FluidSystem::oilPhaseIdx).value(); perf_pressure = fs.pressure(FluidSystem::oilPhaseIdx).value();
} else if (Indices::waterEnabled) { } else if (Indices::waterEnabled) {
@@ -959,7 +959,7 @@ namespace Opm {
if (it != this->node_pressures_.end()) { if (it != this->node_pressures_.end()) {
// The well belongs to a group which has a network nodal pressure, // The well belongs to a group which has a network nodal pressure,
// set the dynamic THP constraint based on the network nodal pressure // set the dynamic THP constraint based on the network nodal pressure
const double nodal_pressure = it->second; const Scalar nodal_pressure = it->second;
well->setDynamicThpLimit(nodal_pressure); well->setDynamicThpLimit(nodal_pressure);
} }
} }
@@ -1217,9 +1217,9 @@ namespace Opm {
const double simulationTime = simulator_.time(); const double simulationTime = simulator_.time();
const auto& comm = simulator_.vanguard().grid().comm(); const auto& comm = simulator_.vanguard().grid().comm();
const auto& summaryState = simulator_.vanguard().summaryState(); const auto& summaryState = simulator_.vanguard().summaryState();
std::vector<double> pot(this->numPhases(), 0.0); std::vector<Scalar> pot(this->numPhases(), 0.0);
const Group& fieldGroup = this->schedule().getGroup("FIELD", reportStepIdx); const Group& fieldGroup = this->schedule().getGroup("FIELD", reportStepIdx);
WellGroupHelpers<double>::updateGuideRates(fieldGroup, WellGroupHelpers<Scalar>::updateGuideRates(fieldGroup,
this->schedule(), this->schedule(),
summaryState, summaryState,
this->phase_usage_, this->phase_usage_,
@@ -1286,7 +1286,7 @@ namespace Opm {
bool do_glift_optimization = false; bool do_glift_optimization = false;
int num_wells_changed = 0; int num_wells_changed = 0;
const double simulation_time = simulator_.time(); const double simulation_time = simulator_.time();
const double min_wait = simulator_.vanguard().schedule().glo(simulator_.episodeIndex()).min_wait(); const Scalar min_wait = simulator_.vanguard().schedule().glo(simulator_.episodeIndex()).min_wait();
// We only optimize if a min_wait time has past. // We only optimize if a min_wait time has past.
// If all_newton is true we still want to optimize several times pr timestep // If all_newton is true we still want to optimize several times pr timestep
// i.e. we also optimize if check simulation_time == last_glift_opt_time_ // i.e. we also optimize if check simulation_time == last_glift_opt_time_
@@ -1390,13 +1390,13 @@ namespace Opm {
if (num_rates_to_sync > 0) { if (num_rates_to_sync > 0) {
std::vector<int> group_indexes; std::vector<int> group_indexes;
group_indexes.reserve(num_rates_to_sync); group_indexes.reserve(num_rates_to_sync);
std::vector<double> group_alq_rates; std::vector<Scalar> group_alq_rates;
group_alq_rates.reserve(num_rates_to_sync); group_alq_rates.reserve(num_rates_to_sync);
std::vector<double> group_oil_rates; std::vector<Scalar> group_oil_rates;
group_oil_rates.reserve(num_rates_to_sync); group_oil_rates.reserve(num_rates_to_sync);
std::vector<double> group_gas_rates; std::vector<Scalar> group_gas_rates;
group_gas_rates.reserve(num_rates_to_sync); group_gas_rates.reserve(num_rates_to_sync);
std::vector<double> group_water_rates; std::vector<Scalar> group_water_rates;
group_water_rates.reserve(num_rates_to_sync); group_water_rates.reserve(num_rates_to_sync);
if (comm.rank() == i) { if (comm.rank() == i) {
for (auto idx : groups_to_sync) { for (auto idx : groups_to_sync) {
@@ -1897,10 +1897,10 @@ namespace Opm {
bool more_network_update = false; bool more_network_update = false;
if (this->shouldBalanceNetwork(episodeIdx, iterationIdx) || mandatory_network_balance) { if (this->shouldBalanceNetwork(episodeIdx, iterationIdx) || mandatory_network_balance) {
const auto local_network_imbalance = this->updateNetworkPressures(episodeIdx); const auto local_network_imbalance = this->updateNetworkPressures(episodeIdx);
const double network_imbalance = comm.max(local_network_imbalance); const Scalar network_imbalance = comm.max(local_network_imbalance);
const auto& balance = this->schedule()[episodeIdx].network_balance(); const auto& balance = this->schedule()[episodeIdx].network_balance();
constexpr double relaxtion_factor = 10.0; constexpr Scalar relaxation_factor = 10.0;
const double tolerance = relax_network_tolerance ? relaxtion_factor * balance.pressure_tolerance() : balance.pressure_tolerance(); const Scalar tolerance = relax_network_tolerance ? relaxation_factor * balance.pressure_tolerance() : balance.pressure_tolerance();
more_network_update = this->networkActive() && network_imbalance > tolerance; more_network_update = this->networkActive() && network_imbalance > tolerance;
} }
@@ -2068,7 +2068,7 @@ namespace Opm {
BlackoilWellModel<TypeTag>:: BlackoilWellModel<TypeTag>::
makeWellSourceEvaluatorFactory(const std::vector<Well>::size_type wellIdx) const makeWellSourceEvaluatorFactory(const std::vector<Well>::size_type wellIdx) const
{ {
using Span = PAvgDynamicSourceData::SourceDataSpan<double>; using Span = PAvgDynamicSourceData::SourceDataSpan<Scalar>;
using Item = typename Span::Item; using Item = typename Span::Item;
return [wellIdx, this]() -> ParallelWBPCalculation::Evaluator return [wellIdx, this]() -> ParallelWBPCalculation::Evaluator
@@ -2244,7 +2244,7 @@ namespace Opm {
DeferredLogger& deferred_logger) DeferredLogger& deferred_logger)
{ {
const int np = this->numPhases(); const int np = this->numPhases();
std::vector<double> potentials; std::vector<Scalar> potentials;
const auto& well = well_container_[widx]; const auto& well = well_container_[widx];
std::string cur_exc_msg; std::string cur_exc_msg;
auto cur_exc_type = ExceptionType::NONE; auto cur_exc_type = ExceptionType::NONE;
@@ -2528,8 +2528,8 @@ namespace Opm {
BlackoilWellModel<TypeTag>:: BlackoilWellModel<TypeTag>::
calcRates(const int fipnum, calcRates(const int fipnum,
const int pvtreg, const int pvtreg,
const std::vector<double>& production_rates, const std::vector<Scalar>& production_rates,
std::vector<double>& resv_coeff) std::vector<Scalar>& resv_coeff)
{ {
rateConverter_->calcCoeff(fipnum, pvtreg, production_rates, resv_coeff); rateConverter_->calcCoeff(fipnum, pvtreg, production_rates, resv_coeff);
} }
@@ -2539,7 +2539,7 @@ namespace Opm {
BlackoilWellModel<TypeTag>:: BlackoilWellModel<TypeTag>::
calcInjRates(const int fipnum, calcInjRates(const int fipnum,
const int pvtreg, const int pvtreg,
std::vector<double>& resv_coeff) std::vector<Scalar>& resv_coeff)
{ {
rateConverter_->calcInjCoeff(fipnum, pvtreg, resv_coeff); rateConverter_->calcInjCoeff(fipnum, pvtreg, resv_coeff);
} }
@@ -2554,17 +2554,17 @@ namespace Opm {
return; return;
int np = this->numPhases(); int np = this->numPhases();
double cellInternalEnergy; Scalar cellInternalEnergy;
double cellBinv; Scalar cellBinv;
double cellDensity; Scalar cellDensity;
double perfPhaseRate; Scalar perfPhaseRate;
const int nw = this->numLocalWells(); const int nw = this->numLocalWells();
for (auto wellID = 0*nw; wellID < nw; ++wellID) { for (auto wellID = 0*nw; wellID < nw; ++wellID) {
const Well& well = this->wells_ecl_[wellID]; const Well& well = this->wells_ecl_[wellID];
if (well.isInjector()) if (well.isInjector())
continue; continue;
std::array<double,2> weighted{0.0,0.0}; std::array<Scalar,2> weighted{0.0,0.0};
auto& [weighted_temperature, total_weight] = weighted; auto& [weighted_temperature, total_weight] = weighted;
auto& well_info = this->local_parallel_well_info_[wellID].get(); auto& well_info = this->local_parallel_well_info_[wellID].get();
@@ -2579,9 +2579,9 @@ namespace Opm {
const auto& fs = intQuants.fluidState(); const auto& fs = intQuants.fluidState();
// we on only have one temperature pr cell any phaseIdx will do // we on only have one temperature pr cell any phaseIdx will do
double cellTemperatures = fs.temperature(/*phaseIdx*/0).value(); Scalar cellTemperatures = fs.temperature(/*phaseIdx*/0).value();
double weight_factor = 0.0; Scalar weight_factor = 0.0;
for (unsigned phaseIdx = 0; phaseIdx < FluidSystem::numPhases; ++phaseIdx) for (unsigned phaseIdx = 0; phaseIdx < FluidSystem::numPhases; ++phaseIdx)
{ {
if (!FluidSystem::phaseIsActive(phaseIdx)) { if (!FluidSystem::phaseIsActive(phaseIdx)) {
@@ -2611,7 +2611,7 @@ namespace Opm {
for (const auto& w : well_container_) { for (const auto& w : well_container_) {
os << w->name() << ":"; os << w->name() << ":";
auto pv = w->getPrimaryVars(); auto pv = w->getPrimaryVars();
for (const double v : pv) { for (const Scalar v : pv) {
os << ' ' << v; os << ' ' << v;
} }
os << '\n'; os << '\n';
@@ -2622,11 +2622,11 @@ namespace Opm {
template <typename TypeTag> template <typename TypeTag>
std::vector<double> std::vector<typename BlackoilWellModel<TypeTag>::Scalar>
BlackoilWellModel<TypeTag>:: BlackoilWellModel<TypeTag>::
getPrimaryVarsDomain(const Domain& domain) const getPrimaryVarsDomain(const Domain& domain) const
{ {
std::vector<double> ret; std::vector<Scalar> ret;
for (const auto& well : well_container_) { for (const auto& well : well_container_) {
if (well_domain_.at(well->name()) == domain.index) { if (well_domain_.at(well->name()) == domain.index) {
const auto& pv = well->getPrimaryVars(); const auto& pv = well->getPrimaryVars();
@@ -2641,7 +2641,7 @@ namespace Opm {
template <typename TypeTag> template <typename TypeTag>
void void
BlackoilWellModel<TypeTag>:: BlackoilWellModel<TypeTag>::
setPrimaryVarsDomain(const Domain& domain, const std::vector<double>& vars) setPrimaryVarsDomain(const Domain& domain, const std::vector<Scalar>& vars)
{ {
std::size_t offset = 0; std::size_t offset = 0;
for (auto& well : well_container_) { for (auto& well : well_container_) {