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

Merge pull request #4937 from totto82/addOutputResidual

Browse Source 
Add output of the residual at the end of a time step
This commit is contained in:
Bård Skaflestad
2023-10-24 15:26:15 +02:00
committed by GitHub
parent 922fa5b9bf 81549a0477
commit 7d4ec59415
3 changed files with 45 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

32
ebos/eclgenericoutputblackoilmodule.cc
View File

@@ -551,6 +551,28 @@ assignToSolution(data::Solution& sol)
data::TargetType::RESTART_OPM_EXTENDED); data::TargetType::RESTART_OPM_EXTENDED);
} }
if (FluidSystem::phaseIsActive(waterPhaseIdx) &&
! this->residual_[waterPhaseIdx].empty())
{
sol.insert("RES_WAT", UnitSystem::measure::liquid_surface_volume,
std::move(this->residual_[waterPhaseIdx]),
data::TargetType::RESTART_OPM_EXTENDED);
}
if (FluidSystem::phaseIsActive(gasPhaseIdx) &&
! this->residual_[gasPhaseIdx].empty())
{
sol.insert("RES_GAS", UnitSystem::measure::gas_surface_volume,
std::move(this->residual_[gasPhaseIdx]),
data::TargetType::RESTART_OPM_EXTENDED);
}
if (FluidSystem::phaseIsActive(oilPhaseIdx) &&
! this->residual_[oilPhaseIdx].empty())
{
sol.insert("RES_OIL", UnitSystem::measure::liquid_surface_volume,
std::move(this->residual_[oilPhaseIdx]),
data::TargetType::RESTART_OPM_EXTENDED);
}
// Fluid in place // Fluid in place
if (this->outputFipRestart_) { if (this->outputFipRestart_) {
for (const auto& phase : Inplace::phases()) { for (const auto& phase : Inplace::phases()) {
@@ -1135,6 +1157,16 @@ doAllocBuffers(const unsigned bufferSize,
} }
} }
if (rstKeywords["RESIDUAL"] > 0) {
rstKeywords["RESIDUAL"] = 0;
for (int phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
{
if (FluidSystem::phaseIsActive(phaseIdx)) {
this->residual_[phaseIdx].resize(bufferSize, 0.0);
}
}
}
// ROCKC // ROCKC
if (rstKeywords["ROCKC"] > 0) { if (rstKeywords["ROCKC"] > 0) {
rstKeywords["ROCKC"] = 0; rstKeywords["ROCKC"] = 0;

3
ebos/eclgenericoutputblackoilmodule.hh
View File

@@ -470,6 +470,9 @@ protected:
std::vector<ScalarBuffer> tracerConcentrations_; std::vector<ScalarBuffer> tracerConcentrations_;
std::array<ScalarBuffer, numPhases> residual_;
std::array<ScalarBuffer, numPhases> flowsi_; std::array<ScalarBuffer, numPhases> flowsi_;
std::array<ScalarBuffer, numPhases> flowsj_; std::array<ScalarBuffer, numPhases> flowsj_;
std::array<ScalarBuffer, numPhases> flowsk_; std::array<ScalarBuffer, numPhases> flowsk_;

10
ebos/ecloutputblackoilmodule.hh
View File

@@ -278,6 +278,7 @@ public:
return; return;
const auto& problem = elemCtx.simulator().problem(); const auto& problem = elemCtx.simulator().problem();
const auto& ebosResid = elemCtx.simulator().model().linearizer().residual();
for (unsigned dofIdx = 0; dofIdx < elemCtx.numPrimaryDof(/*timeIdx=*/0); ++dofIdx) { for (unsigned dofIdx = 0; dofIdx < elemCtx.numPrimaryDof(/*timeIdx=*/0); ++dofIdx) {
const auto& intQuants = elemCtx.intensiveQuantities(dofIdx, /*timeIdx=*/0); const auto& intQuants = elemCtx.intensiveQuantities(dofIdx, /*timeIdx=*/0);
const auto& fs = intQuants.fluidState(); const auto& fs = intQuants.fluidState();
@@ -626,6 +627,15 @@ public:
tracerModel.tracerConcentration(tracerIdx, globalDofIdx); tracerModel.tracerConcentration(tracerIdx, globalDofIdx);
} }
} }
// output residual
for ( int phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx )
{
if (!this->residual_[phaseIdx].empty()) {
const unsigned activeCompIdx = Indices::canonicalToActiveComponentIndex(FluidSystem::solventComponentIndex(phaseIdx));
this->residual_[phaseIdx][globalDofIdx] = ebosResid[globalDofIdx][activeCompIdx];
}
}
} }
} }