Commit aa3eb436 authored by Bernd Flemisch's avatar Bernd Flemisch
Browse files

implicit 3p3cni: fix matrixHeatFlux routines such that they also work on boundary faces

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@10725 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 8cd228bd
......@@ -249,8 +249,8 @@ public:
static const Scalar lSn1 = 0.65;
// arithmetic mean of the liquid saturation and the porosity
const int i = fvGeometry.subContVolFace[scvfIdx].i;
const int j = fvGeometry.subContVolFace[scvfIdx].j;
const int i = fluxDat.face().i;
const int j = fluxDat.face().j;
Scalar Sw = std::max(0.0, (elemVolVars[i].saturation(wPhaseIdx) +
elemVolVars[j].saturation(wPhaseIdx)) / 2);
Scalar Sn = std::max(0.0, (elemVolVars[i].saturation(nPhaseIdx) +
......
......@@ -263,8 +263,8 @@ public:
static const Scalar lSn1 = 0.65;
// arithmetic mean of the liquid saturation and the porosity
const int i = fvGeometry.subContVolFace[faceIdx].i;
const int j = fvGeometry.subContVolFace[faceIdx].j;
const int i = fluxDat.face().i;
const int j = fluxDat.face().j;
Scalar Sw = std::max(0.0, (elemVolVars[i].saturation(wPhaseIdx) +
elemVolVars[j].saturation(wPhaseIdx)) / 2);
Scalar Sn = std::max(0.0, (elemVolVars[i].saturation(nPhaseIdx) +
......
......@@ -258,8 +258,8 @@ public:
static const Scalar lSn1 = 0.65;
// arithmetic mean of the liquid saturation and the porosity
const int i = fvGeometry.subContVolFace[faceIdx].i;
const int j = fvGeometry.subContVolFace[faceIdx].j;
const int i = fluxDat.face().i;
const int j = fluxDat.face().j;
Scalar Sw = std::max(0.0, (elemVolVars[i].saturation(wPhaseIdx) +
elemVolVars[j].saturation(wPhaseIdx)) / 2);
Scalar Sn = std::max(0.0, (elemVolVars[i].saturation(nPhaseIdx) +
......
......@@ -216,8 +216,8 @@ public:
// static const Scalar lGranite = 2.8; // [W / (m K ) ]
//
// // arithmetic mean of the liquid saturation and the porosity
// const int i = fvElemGeom.subContVolFace[faceIdx].i;
// const int j = fvElemGeom.subContVolFace[faceIdx].j;
// const int i = fluxDat.face().i;
// const int j = fluxDat.face().j;
// Scalar Sl = std::max(0.0, (vDat[i].saturation(wPhaseIdx) +
// vDat[j].saturation(wPhaseIdx)) / 2);
// Scalar poro = (porosity(element, fvElemGeom, i) +
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment