From 2885c7efcfac4566db68ab228c054a009a1b83ff Mon Sep 17 00:00:00 2001
From: DennisGlaeser <dennis.glaeser@iws.uni-stuttgart.de>
Date: Tue, 23 Feb 2016 14:12:36 +0100
Subject: [PATCH] [2pdfm] mobility update in the volume variables
After updating the matrix saturations in the interface condition the mobilities in the matrix have to be updated. This is essential in the dfm approach and was not implemented.
---
.../2pdfm/implicit/volumevariables.hh | 14 ++++++++++++++
1 file changed, 14 insertions(+)
diff --git a/dumux/porousmediumflow/2pdfm/implicit/volumevariables.hh b/dumux/porousmediumflow/2pdfm/implicit/volumevariables.hh
index 4cc216aea5..1102f67961 100644
--- a/dumux/porousmediumflow/2pdfm/implicit/volumevariables.hh
+++ b/dumux/porousmediumflow/2pdfm/implicit/volumevariables.hh
@@ -211,6 +211,20 @@ public:
porosityFracture_ = problem.spatialParams().porosityFracture(element,
fvGeometry,
scvIdx);
+
+ // After modifying the Matrix saturations we have to update the fluid state and the matrix mobilities
+ PrimaryVariables updatedMatrixPV;
+ updatedMatrixPV[pressureIdx] = priVars[pressureIdx];
+ updatedMatrixPV[saturationIdx] = satWMatrix_;
+ this->completeFluidState(priVars, problem, element, fvGeometry, scvIdx, fluidState_);
+
+ mobilityMatrix_[wPhaseIdx] =
+ MaterialLaw::krw(materialParamsMatrix, satWMatrix_)
+ / fluidState_.viscosity(wPhaseIdx);
+
+ mobilityMatrix_[nPhaseIdx] =
+ MaterialLaw::krn(materialParamsMatrix, satWMatrix_)
+ / fluidState_.viscosity(nPhaseIdx);
}
/*!
--
GitLab