OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-10 08:23:04 -06:00
Code Issues Packages Projects Releases Wiki Activity

ECL peaceman well: remove some artifact code from the development process

Browse Source 
during debugging this was used to make sure that the rate limits are
correct. the issue is that this is slow if many wells are involved
because it uses a full iteration over the grid for each well to
calculate these rates...
This commit is contained in:
Andreas Lauser 2014-08-06 23:58:24 +02:00
parent 6e478f6f18
commit 4e59104a80
1 changed files with 1 additions and 46 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

47
ewoms/wells/eclpeacemanwell.hh
View File

@ -566,52 +566,7 @@ public:
* accumulation callback.
*/
void beginIterationPostProcess()
{
for (int phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx)
currentSurfaceRates_[phaseIdx] = 0.0;
std::array<Scalar, numPhases> dofReservoirRates;
std::array<Scalar, numPhases> dofSurfaceRates;
// call the accumulation routines
ElementContext elemCtx(simulator_);
auto elemIt = simulator_.gridManager().gridView().template begin</*codim=*/0>();
const auto &elemEndIt = simulator_.gridManager().gridView().template end</*codim=*/0>();
for (; elemIt != elemEndIt; ++elemIt) {
elemCtx.updateStencil(*elemIt);
elemCtx.updateIntensiveQuantities(/*timeIdx=*/0);
for (int dofIdx = 0; dofIdx < elemCtx.numPrimaryDof(/*timeIdx=*/0); ++dofIdx) {
int globalDofIdx = elemCtx.globalSpaceIndex(dofIdx, /*timeIdx=*/0);
if (!applies(globalDofIdx))
continue;
const DofVariables &dofVars = dofVariables_.at(globalDofIdx);
computeUnconstraintVolumetricDofRates_(dofReservoirRates,
dofVars,
elemCtx,
dofIdx,
/*timeIdx=*/0);
limitVolumetricReservoirRates_(dofReservoirRates,
dofVars,
elemCtx,
dofIdx,
/*timeIdx=*/0);
computeSurfaceRates_(dofSurfaceRates,
dofReservoirRates,
elemCtx.intensiveQuantities(dofIdx, /*timeIdx=*/0).fluidState());
for (int phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx)
currentSurfaceRates_[phaseIdx] += dofSurfaceRates[phaseIdx];
}
}
Valgrind::CheckDefined(currentSurfaceRates_[oilPhaseIdx]);
Valgrind::CheckDefined(currentSurfaceRates_[gasPhaseIdx]);
Valgrind::CheckDefined(currentSurfaceRates_[waterPhaseIdx]);
}
{ }
/*!
* \brief Called by the simulator after each Newton-Raphson iteration.