From fbcf7111be013e22382e97b88dea75c7a1cb433c Mon Sep 17 00:00:00 2001
From: sergio <s.ferraris@opencfd.co.uk>
Date: Tue, 10 Nov 2020 08:52:31 -0800
Subject: [PATCH] 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
---
applications/solvers/multiphase/compressibleInterFoam/TEqn.H | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/applications/solvers/multiphase/compressibleInterFoam/TEqn.H b/applications/solvers/multiphase/compressibleInterFoam/TEqn.H
index 29ba87088f0..ed933ccd236 100644
--- 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
)
*(
--
GitLab