diff --git a/dumux/decoupled/2p/diffusion/fvmpfa/omethod/fvmpfao2dpressurevelocity2p.hh b/dumux/decoupled/2p/diffusion/fvmpfa/omethod/fvmpfao2dpressurevelocity2p.hh index ef8b6f8b70e577adde0e751e087b3debc166a59b..c9c1db2b25660e5f1815f3150035a1cb59062f4e 100644 --- a/dumux/decoupled/2p/diffusion/fvmpfa/omethod/fvmpfao2dpressurevelocity2p.hh +++ b/dumux/decoupled/2p/diffusion/fvmpfa/omethod/fvmpfao2dpressurevelocity2p.hh @@ -464,9 +464,6 @@ void FvMpfaO2dPressureVelocity2p<TypeTag>::calculateVelocityOnBoundary(const Int pressNwBound = pressBound; } - //get temperature at current position - Scalar temperature = problem_.temperature(*element); - Scalar lambdaWBound = MaterialLaw::krw(problem_.spatialParams().materialLawParams(*element), satW) / viscosity_[wPhaseIdx]; Scalar lambdaNwBound = MaterialLaw::krn(problem_.spatialParams().materialLawParams(*element), satW) @@ -476,9 +473,8 @@ void FvMpfaO2dPressureVelocity2p<TypeTag>::calculateVelocityOnBoundary(const Int Scalar potentialDiffNw = cellData.fluxData().upwindPotential(nPhaseIdx, isIndex); //calculate potential gradient - - potentialDiffW = (cellData.pressure(wPhaseIdx) - pressWBound); - potentialDiffNw = (cellData.pressure(nPhaseIdx) - pressNwBound); + potentialDiffW = (cellData.pressure(wPhaseIdx) - pressWBound); + potentialDiffNw = (cellData.pressure(nPhaseIdx) - pressNwBound); potentialDiffW += density_[wPhaseIdx] * (distVec * problem_.gravity()); potentialDiffNw += density_[nPhaseIdx] * (distVec * problem_.gravity());