From e22870f50851d64f721d89ae2c1618b68bcee953 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Tue, 8 Feb 2022 16:25:44 +0000
Subject: [PATCH] compressibleMultiphaseInterFoam: Added contErr to avoid
 unboundedness in U and T

---
 .../compressibleMultiphaseInterFoam/contErr.H     | 15 +++++++++++++++
 1 file changed, 15 insertions(+)
 create mode 100644 applications/solvers/multiphase/compressibleMultiphaseInterFoam/contErr.H

diff --git a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/contErr.H b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/contErr.H
new file mode 100644
index 0000000000..68cd06ed56
--- /dev/null
+++ b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/contErr.H
@@ -0,0 +1,15 @@
+volScalarField::Internal contErr
+(
+    (fvc::ddt(rho) + fvc::div(mixture.rhoPhi()))()
+);
+
+forAllConstIter
+(
+    PtrDictionary<phaseModel>,
+    mixture.phases(),
+    phase
+)
+{
+    const volScalarField& rho = phase().thermo().rho()();
+    contErr -= (fvModels.source(phase(), rho)&rho);
+}