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

ECL peaceman well: use the reservoir cell's phase density for injected fluids

Browse Source 
(instead of calculating it from the hole pressure. this seems "less
correct" but the way Eclipse does things...)
This commit is contained in:
Andreas Lauser 2015-01-12 11:53:10 +01:00
parent a9e40ce36a
commit 1233e12e76
1 changed files with 2 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
applications/ebos/eclpeacemanwell.hh
View File

@ -1139,7 +1139,6 @@ protected:
// gravity constant
Scalar g = simulator_.problem().gravity()[dimWorld - 1];
typename FluidSystem::ParameterCache paramCache;
for (int phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
// well model due to Peaceman; see Chen et al., p. 449
@ -1147,11 +1146,10 @@ protected:
Scalar p = dofVars.pressure[phaseIdx];
// density and mobility of fluid phase
Scalar rho;
Scalar rho = dofVars.density[phaseIdx];
Scalar lambda;
if (wellType_ == Producer) {
//assert(p < pbh);
rho = dofVars.density[phaseIdx];
lambda = dofVars.mobility[phaseIdx];
}
else if (wellType_ == Injector) {
@ -1159,21 +1157,14 @@ protected:
if (phaseIdx != injectedPhaseIdx_)
continue;
injectionFluidState_.setPressure(phaseIdx, p);
typename FluidSystem::ParameterCache paramCache;
paramCache.updateAll(injectionFluidState_);
rho = FluidSystem::density(injectionFluidState_, paramCache, phaseIdx);
// use the total mobility, i.e. the sum of all phase mobilities at the
// injector cell. this seems a bit weird: at the wall of the borehole,
// there should only be injected phase present, so its mobility should be
// 1/viscosity...
lambda = 0.0;
for (int phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
for (int phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
lambda += dofVars.mobility[phaseIdx];
}
}
else
OPM_THROW(std::logic_error,
"Type of well \"" << name() << "\" is undefined");