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

- changed pointer to refernace for IQ for possibility of using rvalues

Browse Source
This commit is contained in:
hnil 2023-03-28 09:37:00 +02:00
parent 88e5b235b8
commit f0b5e6a670
1 changed files with 5 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
opm/simulators/aquifers/AquiferAnalytical.hpp
View File

@ -150,26 +150,23 @@ public:
if (idx < 0) if (idx < 0)
return; return;
const auto* intQuantsPtr = &model.intensiveQuantities(cellIdx, timeIdx); const auto& intQuants = model.intensiveQuantities(cellIdx, timeIdx);
if (intQuantsPtr == nullptr) {
throw std::logic_error("Invalid intensive quantities cache detected in AquiferAnalytical::addToSource()");
}
// This is the pressure at td + dt // This is the pressure at td + dt
this->updateCellPressure(this->pressure_current_, idx, *intQuantsPtr); this->updateCellPressure(this->pressure_current_, idx, intQuants);
this->calculateInflowRate(idx, this->ebos_simulator_); this->calculateInflowRate(idx, this->ebos_simulator_);
rates[BlackoilIndices::conti0EqIdx + compIdx_()] rates[BlackoilIndices::conti0EqIdx + compIdx_()]
+= this->Qai_[idx] / model.dofTotalVolume(cellIdx); += this->Qai_[idx] / model.dofTotalVolume(cellIdx);
if constexpr (enableEnergy) { if constexpr (enableEnergy) {
auto fs = intQuantsPtr->fluidState(); auto fs = intQuants.fluidState();
if (this->Ta0_.has_value() && this->Qai_[idx] > 0) if (this->Ta0_.has_value() && this->Qai_[idx] > 0)
{ {
fs.setTemperature(this->Ta0_.value()); fs.setTemperature(this->Ta0_.value());
typedef typename std::decay<decltype(fs)>::type::Scalar FsScalar; typedef typename std::decay<decltype(fs)>::type::Scalar FsScalar;
typename FluidSystem::template ParameterCache<FsScalar> paramCache; typename FluidSystem::template ParameterCache<FsScalar> paramCache;
const unsigned pvtRegionIdx = intQuantsPtr->pvtRegionIndex(); const unsigned pvtRegionIdx = intQuants.pvtRegionIndex();
paramCache.setRegionIndex(pvtRegionIdx); paramCache.setRegionIndex(pvtRegionIdx);
paramCache.setMaxOilSat(this->ebos_simulator_.problem().maxOilSaturation(cellIdx)); paramCache.setMaxOilSat(this->ebos_simulator_.problem().maxOilSaturation(cellIdx));
paramCache.updatePhase(fs, this->phaseIdx_()); paramCache.updatePhase(fs, this->phaseIdx_());
@ -177,7 +174,7 @@ public:
fs.setEnthalpy(this->phaseIdx_(), h); fs.setEnthalpy(this->phaseIdx_(), h);
} }
rates[BlackoilIndices::contiEnergyEqIdx] rates[BlackoilIndices::contiEnergyEqIdx]
+= this->Qai_[idx] *fs.enthalpy(this->phaseIdx_()) * FluidSystem::referenceDensity( this->phaseIdx_(), intQuantsPtr->pvtRegionIndex()) / model.dofTotalVolume(cellIdx); += this->Qai_[idx] *fs.enthalpy(this->phaseIdx_()) * FluidSystem::referenceDensity( this->phaseIdx_(), intQuants.pvtRegionIndex()) / model.dofTotalVolume(cellIdx);
} }
} }