From 5c3c481e85eed71657f33ea01c862f64238d7590 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Thu, 11 Jun 2015 15:56:06 +0100
Subject: [PATCH] Field bounding after topo change

---
 .../compressible/sonicDyMFoam/sonicDyMFoam.C       | 14 +++++++++-----
 1 file changed, 9 insertions(+), 5 deletions(-)

diff --git a/applications/solvers/compressible/sonicDyMFoam/sonicDyMFoam.C b/applications/solvers/compressible/sonicDyMFoam/sonicDyMFoam.C
index 524e4168e..4fb6e2572 100644
--- a/applications/solvers/compressible/sonicDyMFoam/sonicDyMFoam.C
+++ b/applications/solvers/compressible/sonicDyMFoam/sonicDyMFoam.C
@@ -25,12 +25,15 @@ Application
     sonicDyMFoam
 
 Description
-    Transient solver for trans-sonic/supersonic for laminar or flow of a
-    compressible gas with support for mesh motion and topological changes
+    Transient solver for trans-sonic/supersonic for laminar or turbulent
+    flow of a compressible gas with support for mesh motion and
+    topological changes
 
     Uses the flexible PIMPLE (PISO-SIMPLE) solution for time-resolved and
     pseudo-transient simulations.
 
+    Turbulence modelling is generic, i.e. laminar, RAS or LES may be selected.
+
     Updated from sonicFoamAutoMotion by Hrvoje Jasak
 
 Author
@@ -82,9 +85,10 @@ int main(int argc, char *argv[])
         // Mesh motion update
         if (meshChanged)
         {
-            T.correctBoundaryConditions();
-            p.correctBoundaryConditions();
-            e.correctBoundaryConditions();
+            T = max(T, TMin);
+            p = max(p, pMin);
+            e = max(e, thermo.Cv()*TMin);
+
             thermo.correct();
             rho = thermo.rho();