BUG: Taking out contErr contribution to TEq for divUp. Fixes #1813.

contErr is based on rhoPhi not on phi. The term contErr/rho*p in TEq makes a large correction which occasionally makes T unbounded
2020-11-10 08:52:31 -08:00 · 2020-11-10 08:52:31 -08:00 · fbcf7111be
commit fbcf7111be
parent 4cf44850a7
1 changed files with 1 additions and 1 deletions
--- a/applications/solvers/multiphase/compressibleInterFoam/TEqn.H
+++ b/applications/solvers/multiphase/compressibleInterFoam/TEqn.H
@ -4,7 +4,7 @@
        fvm::ddt(rho, T) + fvm::div(rhoPhi, T) - fvm::Sp(contErr, T)
      - fvm::laplacian(turbulence.alphaEff(), T)
      + (
-             divUp - contErr/rho*p
+             divUp()// - contErr/rho*p
          +  (fvc::ddt(rho, K) + fvc::div(rhoPhi, K))() - contErr*K
        )
       *(