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

fix analytical aquifers + energy

Browse Source
This commit is contained in:
Tor Harald Sandve
2022-08-29 10:59:20 +02:00
parent 949bb6ecb5
commit 3ad04d8326
3 changed files with 25 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -158,6 +158,25 @@ public:
rates[BlackoilIndices::conti0EqIdx + compIdx_()]
+= this->Qai_[idx] / model.dofTotalVolume(cellIdx);
if constexpr (enableEnergy) {
auto fs = intQuantsPtr->fluidState();
if (this->Ta0_.has_value() && this->Qai_[idx] > 0)
{
fs.setTemperature(this->Ta0_.value());
typedef typename std::decay<decltype(fs)>::type::Scalar FsScalar;
typename FluidSystem::template ParameterCache<FsScalar> paramCache;
const unsigned pvtRegionIdx = intQuantsPtr->pvtRegionIndex();
paramCache.setRegionIndex(pvtRegionIdx);
paramCache.setMaxOilSat(this->ebos_simulator_.problem().maxOilSaturation(cellIdx));
paramCache.updatePhase(fs, this->phaseIdx_());
const auto& h = FluidSystem::enthalpy(fs, paramCache, this->phaseIdx_());
fs.setEnthalpy(this->phaseIdx_(), h);
}
rates[BlackoilIndices::contiEnergyEqIdx]
+= this->Qai_[idx] *fs.enthalpy(this->phaseIdx_()) * FluidSystem::referenceDensity( this->phaseIdx_(), intQuantsPtr->pvtRegionIndex()) / model.dofTotalVolume(cellIdx);
}
}
std::size_t size() const
@@ -404,6 +423,7 @@ protected:
Scalar Tc_{}; // Time constant
Scalar pa0_{}; // initial aquifer pressure
std::optional<Scalar> Ta0_{}; // initial aquifer temperature
Scalar rhow_{};
Eval W_flux_;

3
opm/simulators/aquifers/AquiferCarterTracy.hpp
View File

@@ -215,6 +215,9 @@ protected:
}
this->pa0_ = this->aquct_data_.initial_pressure.value();
if (this->aquct_data_.initial_temperature.has_value())
this->Ta0_ = this->aquct_data_.initial_temperature.value();
if(this->co2store_()) {
const auto press = this->aquct_data_.initial_pressure.value();

2
opm/simulators/aquifers/AquiferFetkovich.hpp
View File

@@ -163,6 +163,8 @@ protected:
this->rhow_ = this->aqufetp_data_.waterDensity();
this->pa0_ = this->aqufetp_data_.initial_pressure.value();
if (this->aqufetp_data_.initial_temperature.has_value())
this->Ta0_ = this->aqufetp_data_.initial_temperature.value();
this->aquifer_pressure_ = this->pa0_;
}