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

Merge pull request #863 from totto82/fixThermalAqu

Browse Source 
avoid divition by zero for case without conductivity
This commit is contained in:
Kai Bao
2024-01-16 22:59:44 +01:00
committed by GitHub
parent e956bfb029 ef198675f1
commit 75ced1df73
1 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
opm/models/blackoil/blackoilenergymodules.hh
View File

@@ -567,13 +567,13 @@ public:
exLambda = decay<Scalar>(exIq.totalThermalConductivity());
Evaluation H;
if (inLambda > 0.0 && exLambda > 0.0) {
const Evaluation& inH = inLambda*inAlpha;
const Evaluation& exH = exLambda*outAlpha;
if (inH > 0 && exH > 0) {
// compute the "thermal transmissibility". In contrast to the normal
// transmissibility this cannot be done as a preprocessing step because the
// average thermal conductivity is analogous to the permeability but
// depends on the solution.
const Evaluation& inH = inLambda*inAlpha;
const Evaluation& exH = exLambda*outAlpha;
H = 1.0/(1.0/inH + 1.0/exH);
}
else