compressibleMultiphaseInterFoam: Added contErr to avoid unboundedness in U and T

2022-02-08 16:25:44 +00:00
parent 13b78294ff
commit e22870f508
1 changed files with 15 additions and 0 deletions
--- a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/contErr.H
+++ 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);
+}