Merge pull request #645 from totto82/fixCond

fix conductivity harmonic mean computation
This commit is contained in:
Bård Skaflestad 2021-03-27 13:58:44 +01:00 committed by GitHub
commit 754c37106d

View File

@ -528,9 +528,10 @@ public:
// transmissibility this cannot be done as a preprocessing step because the
// average thermal thermal conductivity is analogous to the permeability but
// depends on the solution.
Scalar alpha = elemCtx.problem().thermalHalfTransmissibility(elemCtx, scvfIdx, timeIdx);
const Evaluation& inH = inLambda*alpha;
const Evaluation& exH = exLambda*alpha;
Scalar inAlpha = elemCtx.problem().thermalHalfTransmissibilityIn(elemCtx, scvfIdx, timeIdx);
Scalar outAlpha = elemCtx.problem().thermalHalfTransmissibilityOut(elemCtx, scvfIdx, timeIdx);
const Evaluation& inH = inLambda*inAlpha;
const Evaluation& exH = exLambda*outAlpha;
H = 1.0/(1.0/inH + 1.0/exH);
}
else