From bcb7c3e5d979562ad3da4c285ded245d9f244242 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Fri, 25 Feb 2011 12:12:56 +0000
Subject: [PATCH] Added evaluation of turbulence boundary conditions after topo
 change: attach/detach bug

---
 src/turbulenceModels/compressible/RAS/LRR/LRR.C        | 10 +++++++++-
 .../RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C          | 10 +++++++++-
 .../compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C |  8 ++++++++
 .../compressible/RAS/RNGkEpsilon/RNGkEpsilon.C         | 10 +++++++++-
 .../compressible/RAS/SpalartAllmaras/SpalartAllmaras.C |  7 +++++++
 .../compressible/RAS/kEpsilon/kEpsilon.C               | 10 +++++++++-
 .../compressible/RAS/kOmegaSST/kOmegaSST.C             |  8 ++++++++
 .../compressible/RAS/kOmegaSST_LowRe/kOmegaSST_LowRe.C |  8 ++++++++
 .../compressible/RAS/realizableKE/realizableKE.C       | 10 +++++++++-
 src/turbulenceModels/incompressible/RAS/LRR/LRR.C      | 10 +++++++++-
 .../incompressible/RAS/LamBremhorstKE/LamBremhorstKE.C |  8 ++++++++
 .../RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C          | 10 +++++++++-
 .../RAS/LaunderSharmaKE/LaunderSharmaKE.C              |  8 ++++++++
 .../incompressible/RAS/LienCubicKE/LienCubicKE.C       | 10 +++++++++-
 .../RAS/LienCubicKELowRe/LienCubicKELowRe.C            |  8 ++++++++
 .../RAS/LienLeschzinerLowRe/LienLeschzinerLowRe.C      |  8 ++++++++
 .../RAS/NonlinearKEShih/NonlinearKEShih.C              |  8 ++++++++
 .../incompressible/RAS/RNGkEpsilon/RNGkEpsilon.C       | 10 +++++++++-
 .../RAS/SpalartAllmaras/SpalartAllmaras.C              |  5 +++++
 .../incompressible/RAS/kEpsilon/kEpsilon.C             | 10 +++++++++-
 .../incompressible/RAS/kOmega/kOmega.C                 |  8 ++++++++
 .../incompressible/RAS/kOmegaSST/kOmegaSST.C           |  8 ++++++++
 .../RAS/kOmegaSST_LowRe/kOmegaSST_LowRe.C              |  8 ++++++++
 src/turbulenceModels/incompressible/RAS/qZeta/qZeta.C  |  8 ++++++++
 .../incompressible/RAS/realizableKE/realizableKE.C     |  8 ++++++++
 25 files changed, 206 insertions(+), 10 deletions(-)

diff --git a/src/turbulenceModels/compressible/RAS/LRR/LRR.C b/src/turbulenceModels/compressible/RAS/LRR/LRR.C
index e524d9214..62dec7f61 100644
--- a/src/turbulenceModels/compressible/RAS/LRR/LRR.C
+++ b/src/turbulenceModels/compressible/RAS/LRR/LRR.C
@@ -321,6 +321,14 @@ bool LRR::read()
 
 void LRR::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     if (!turbulence_)
     {
         // Re-calculate viscosity
@@ -339,7 +347,7 @@ void LRR::correct()
     volSymmTensorField P = -twoSymm(R_ & fvc::grad(U_));
     volScalarField G("RASModel::G", 0.5*mag(tr(P)));
 
-    // Update espsilon and G at the wall
+    // Update epsilon and G at the wall
     epsilon_.boundaryField().updateCoeffs();
 
     // Dissipation equation
diff --git a/src/turbulenceModels/compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C b/src/turbulenceModels/compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C
index 61cc2d626..d88343101 100644
--- a/src/turbulenceModels/compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C
+++ b/src/turbulenceModels/compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C
@@ -354,6 +354,14 @@ bool LaunderGibsonRSTM::read()
 
 void LaunderGibsonRSTM::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     if (!turbulence_)
     {
         // Re-calculate viscosity
@@ -377,7 +385,7 @@ void LaunderGibsonRSTM::correct()
     volSymmTensorField P = -twoSymm(R_ & fvc::grad(U_));
     volScalarField G("RASModel::G", 0.5*mag(tr(P)));
 
-    // Update espsilon and G at the wall
+    // Update epsilon and G at the wall
     epsilon_.boundaryField().updateCoeffs();
 
     // Dissipation equation
diff --git a/src/turbulenceModels/compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C b/src/turbulenceModels/compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C
index 08761f521..d93775510 100644
--- a/src/turbulenceModels/compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C
+++ b/src/turbulenceModels/compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C
@@ -272,6 +272,14 @@ bool LaunderSharmaKE::read()
 
 void LaunderSharmaKE::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     if (!turbulence_)
     {
         // Re-calculate viscosity
diff --git a/src/turbulenceModels/compressible/RAS/RNGkEpsilon/RNGkEpsilon.C b/src/turbulenceModels/compressible/RAS/RNGkEpsilon/RNGkEpsilon.C
index 7d9f09a65..24c1f8def 100644
--- a/src/turbulenceModels/compressible/RAS/RNGkEpsilon/RNGkEpsilon.C
+++ b/src/turbulenceModels/compressible/RAS/RNGkEpsilon/RNGkEpsilon.C
@@ -273,6 +273,14 @@ bool RNGkEpsilon::read()
 
 void RNGkEpsilon::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     if (!turbulence_)
     {
         // Re-calculate viscosity
@@ -307,7 +315,7 @@ void RNGkEpsilon::correct()
     volScalarField R =
         ((eta*(-eta/eta0_ + scalar(1)))/(beta_*eta3 + scalar(1)));
 
-    // Update espsilon and G at the wall
+    // Update epsilon and G at the wall
     epsilon_.boundaryField().updateCoeffs();
 
     // Dissipation equation
diff --git a/src/turbulenceModels/compressible/RAS/SpalartAllmaras/SpalartAllmaras.C b/src/turbulenceModels/compressible/RAS/SpalartAllmaras/SpalartAllmaras.C
index 95a869049..12d34c0bf 100644
--- a/src/turbulenceModels/compressible/RAS/SpalartAllmaras/SpalartAllmaras.C
+++ b/src/turbulenceModels/compressible/RAS/SpalartAllmaras/SpalartAllmaras.C
@@ -328,6 +328,13 @@ bool SpalartAllmaras::read()
 
 void SpalartAllmaras::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(nuTilda_, dimensionedScalar("0", nuTilda_.dimensions(), 0.0));
+    }
+
     if (!turbulence_)
     {
         // Re-calculate viscosity
diff --git a/src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C b/src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C
index 58953d267..ded440a8d 100644
--- a/src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C
+++ b/src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C
@@ -254,6 +254,14 @@ bool kEpsilon::read()
 
 void kEpsilon::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     if (!turbulence_)
     {
         // Re-calculate viscosity
@@ -280,7 +288,7 @@ void kEpsilon::correct()
     volScalarField G("RASModel::G", mut_*(tgradU() && dev(twoSymm(tgradU()))));
     tgradU.clear();
 
-    // Update espsilon and G at the wall
+    // Update epsilon and G at the wall
     epsilon_.boundaryField().updateCoeffs();
 
     // Dissipation equation
diff --git a/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C b/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C
index 4eb2d0db5..7ec82bced 100644
--- a/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C
+++ b/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C
@@ -358,6 +358,14 @@ bool kOmegaSST::read()
 
 void kOmegaSST::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(omega_, omega0_);
+    }
+
     if (!turbulence_)
     {
         // Re-calculate viscosity
diff --git a/src/turbulenceModels/compressible/RAS/kOmegaSST_LowRe/kOmegaSST_LowRe.C b/src/turbulenceModels/compressible/RAS/kOmegaSST_LowRe/kOmegaSST_LowRe.C
index cbe1077f0..41ebfeb22 100644
--- a/src/turbulenceModels/compressible/RAS/kOmegaSST_LowRe/kOmegaSST_LowRe.C
+++ b/src/turbulenceModels/compressible/RAS/kOmegaSST_LowRe/kOmegaSST_LowRe.C
@@ -378,6 +378,14 @@ bool kOmegaSST_LowRe::read()
 
 void kOmegaSST_LowRe::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(omega_, omega0_);
+    }
+
     if (!turbulence_)
     {
         // Re-calculate viscosity
diff --git a/src/turbulenceModels/compressible/RAS/realizableKE/realizableKE.C b/src/turbulenceModels/compressible/RAS/realizableKE/realizableKE.C
index f38b5a3a7..15be7d502 100644
--- a/src/turbulenceModels/compressible/RAS/realizableKE/realizableKE.C
+++ b/src/turbulenceModels/compressible/RAS/realizableKE/realizableKE.C
@@ -287,6 +287,14 @@ bool realizableKE::read()
 
 void realizableKE::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     if (!turbulence_)
     {
         // Re-calculate viscosity
@@ -318,7 +326,7 @@ void realizableKE::correct()
 
     volScalarField G("RASModel::G", mut_*(gradU && dev(twoSymm(gradU))));
 
-    // Update espsilon and G at the wall
+    // Update epsilon and G at the wall
     epsilon_.boundaryField().updateCoeffs();
 
     // Dissipation equation
diff --git a/src/turbulenceModels/incompressible/RAS/LRR/LRR.C b/src/turbulenceModels/incompressible/RAS/LRR/LRR.C
index c89704675..8d5d0e45b 100644
--- a/src/turbulenceModels/incompressible/RAS/LRR/LRR.C
+++ b/src/turbulenceModels/incompressible/RAS/LRR/LRR.C
@@ -289,6 +289,14 @@ bool LRR::read()
 
 void LRR::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        R_.correctBoundaryConditions();
+        bound(epsilon_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
@@ -299,7 +307,7 @@ void LRR::correct()
     volSymmTensorField P = -twoSymm(R_ & fvc::grad(U_));
     volScalarField G("RASModel::G", 0.5*mag(tr(P)));
 
-    // Update espsilon and G at the wall
+    // Update epsilon and G at the wall
     epsilon_.boundaryField().updateCoeffs();
 
     // Dissipation equation
diff --git a/src/turbulenceModels/incompressible/RAS/LamBremhorstKE/LamBremhorstKE.C b/src/turbulenceModels/incompressible/RAS/LamBremhorstKE/LamBremhorstKE.C
index 5597df705..af206f733 100644
--- a/src/turbulenceModels/incompressible/RAS/LamBremhorstKE/LamBremhorstKE.C
+++ b/src/turbulenceModels/incompressible/RAS/LamBremhorstKE/LamBremhorstKE.C
@@ -220,6 +220,14 @@ bool LamBremhorstKE::read()
 
 void LamBremhorstKE::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
diff --git a/src/turbulenceModels/incompressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C b/src/turbulenceModels/incompressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C
index 94a037f50..ef38dcde7 100644
--- a/src/turbulenceModels/incompressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C
+++ b/src/turbulenceModels/incompressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C
@@ -325,6 +325,14 @@ bool LaunderGibsonRSTM::read()
 
 void LaunderGibsonRSTM::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        R_.correctBoundaryConditions();
+        bound(epsilon_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
@@ -340,7 +348,7 @@ void LaunderGibsonRSTM::correct()
     volSymmTensorField P = -twoSymm(R_ & fvc::grad(U_));
     volScalarField G("RASModel::G", 0.5*mag(tr(P)));
 
-    // Update espsilon and G at the wall
+    // Update epsilon and G at the wall
     epsilon_.boundaryField().updateCoeffs();
 
     // Dissipation equation
diff --git a/src/turbulenceModels/incompressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C b/src/turbulenceModels/incompressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C
index 0d1713fd8..259ae725a 100644
--- a/src/turbulenceModels/incompressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C
+++ b/src/turbulenceModels/incompressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C
@@ -226,6 +226,14 @@ bool LaunderSharmaKE::read()
 
 void LaunderSharmaKE::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilonTilda_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
diff --git a/src/turbulenceModels/incompressible/RAS/LienCubicKE/LienCubicKE.C b/src/turbulenceModels/incompressible/RAS/LienCubicKE/LienCubicKE.C
index 4b298520e..ef889dc1e 100644
--- a/src/turbulenceModels/incompressible/RAS/LienCubicKE/LienCubicKE.C
+++ b/src/turbulenceModels/incompressible/RAS/LienCubicKE/LienCubicKE.C
@@ -313,6 +313,14 @@ bool LienCubicKE::read()
 
 void LienCubicKE::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
@@ -331,7 +339,7 @@ void LienCubicKE::correct()
         Cmu_*sqr(k_)/epsilon_*S2 - (nonlinearStress_ && gradU_)
     );
 
-    // Update espsilon and G at the wall
+    // Update epsilon and G at the wall
     epsilon_.boundaryField().updateCoeffs();
 
     // Dissipation equation
diff --git a/src/turbulenceModels/incompressible/RAS/LienCubicKELowRe/LienCubicKELowRe.C b/src/turbulenceModels/incompressible/RAS/LienCubicKELowRe/LienCubicKELowRe.C
index 820a8569c..23c36385a 100644
--- a/src/turbulenceModels/incompressible/RAS/LienCubicKELowRe/LienCubicKELowRe.C
+++ b/src/turbulenceModels/incompressible/RAS/LienCubicKELowRe/LienCubicKELowRe.C
@@ -384,6 +384,14 @@ bool LienCubicKELowRe::read()
 
 void LienCubicKELowRe::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
diff --git a/src/turbulenceModels/incompressible/RAS/LienLeschzinerLowRe/LienLeschzinerLowRe.C b/src/turbulenceModels/incompressible/RAS/LienLeschzinerLowRe/LienLeschzinerLowRe.C
index 91ea566b7..bdb01c5bc 100644
--- a/src/turbulenceModels/incompressible/RAS/LienLeschzinerLowRe/LienLeschzinerLowRe.C
+++ b/src/turbulenceModels/incompressible/RAS/LienLeschzinerLowRe/LienLeschzinerLowRe.C
@@ -270,6 +270,14 @@ bool LienLeschzinerLowRe::read()
 
 void LienLeschzinerLowRe::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
diff --git a/src/turbulenceModels/incompressible/RAS/NonlinearKEShih/NonlinearKEShih.C b/src/turbulenceModels/incompressible/RAS/NonlinearKEShih/NonlinearKEShih.C
index 5ac9f08ef..ae8f48107 100644
--- a/src/turbulenceModels/incompressible/RAS/NonlinearKEShih/NonlinearKEShih.C
+++ b/src/turbulenceModels/incompressible/RAS/NonlinearKEShih/NonlinearKEShih.C
@@ -306,6 +306,14 @@ bool NonlinearKEShih::read()
 
 void NonlinearKEShih::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
diff --git a/src/turbulenceModels/incompressible/RAS/RNGkEpsilon/RNGkEpsilon.C b/src/turbulenceModels/incompressible/RAS/RNGkEpsilon/RNGkEpsilon.C
index 234b1f99e..bb2b0aabe 100644
--- a/src/turbulenceModels/incompressible/RAS/RNGkEpsilon/RNGkEpsilon.C
+++ b/src/turbulenceModels/incompressible/RAS/RNGkEpsilon/RNGkEpsilon.C
@@ -239,6 +239,14 @@ bool RNGkEpsilon::read()
 
 void RNGkEpsilon::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
@@ -254,7 +262,7 @@ void RNGkEpsilon::correct()
     volScalarField R =
         ((eta*(scalar(1) - eta/eta0_))/(scalar(1) + beta_*eta*sqr(eta)));
 
-    // Update espsilon and G at the wall
+    // Update epsilon and G at the wall
     epsilon_.boundaryField().updateCoeffs();
 
     // Dissipation equation
diff --git a/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C b/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C
index 9616c9ebd..a32117c23 100644
--- a/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C
+++ b/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C
@@ -348,6 +348,11 @@ bool SpalartAllmaras::read()
 
 void SpalartAllmaras::correct()
 {
+    if (mesh_.changing())
+    {
+        bound(nuTilda_, dimensionedScalar("0", nuTilda_.dimensions(), 0.0));
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
diff --git a/src/turbulenceModels/incompressible/RAS/kEpsilon/kEpsilon.C b/src/turbulenceModels/incompressible/RAS/kEpsilon/kEpsilon.C
index f14d00284..9706982f2 100644
--- a/src/turbulenceModels/incompressible/RAS/kEpsilon/kEpsilon.C
+++ b/src/turbulenceModels/incompressible/RAS/kEpsilon/kEpsilon.C
@@ -208,6 +208,14 @@ bool kEpsilon::read()
 
 void kEpsilon::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
@@ -217,7 +225,7 @@ void kEpsilon::correct()
 
     volScalarField G("RASModel::G", nut_*2*magSqr(symm(fvc::grad(U_))));
 
-    // Update espsilon and G at the wall
+    // Update epsilon and G at the wall
     epsilon_.boundaryField().updateCoeffs();
 
     // Dissipation equation
diff --git a/src/turbulenceModels/incompressible/RAS/kOmega/kOmega.C b/src/turbulenceModels/incompressible/RAS/kOmega/kOmega.C
index 787cf2771..217f73137 100644
--- a/src/turbulenceModels/incompressible/RAS/kOmega/kOmega.C
+++ b/src/turbulenceModels/incompressible/RAS/kOmega/kOmega.C
@@ -217,6 +217,14 @@ bool kOmega::read()
 
 void kOmega::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(omega_, omega0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
diff --git a/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C b/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C
index 45989de6d..088281c0f 100644
--- a/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C
+++ b/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C
@@ -325,6 +325,14 @@ bool kOmegaSST::read()
 
 void kOmegaSST::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(omega_, omega0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
diff --git a/src/turbulenceModels/incompressible/RAS/kOmegaSST_LowRe/kOmegaSST_LowRe.C b/src/turbulenceModels/incompressible/RAS/kOmegaSST_LowRe/kOmegaSST_LowRe.C
index 37b8562de..6cb1130b3 100644
--- a/src/turbulenceModels/incompressible/RAS/kOmegaSST_LowRe/kOmegaSST_LowRe.C
+++ b/src/turbulenceModels/incompressible/RAS/kOmegaSST_LowRe/kOmegaSST_LowRe.C
@@ -346,6 +346,14 @@ bool kOmegaSST_LowRe::read()
 
 void kOmegaSST_LowRe::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(omega_, omega0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
diff --git a/src/turbulenceModels/incompressible/RAS/qZeta/qZeta.C b/src/turbulenceModels/incompressible/RAS/qZeta/qZeta.C
index 2a3e569c0..c953c608a 100644
--- a/src/turbulenceModels/incompressible/RAS/qZeta/qZeta.C
+++ b/src/turbulenceModels/incompressible/RAS/qZeta/qZeta.C
@@ -274,6 +274,14 @@ bool qZeta::read()
 
 void qZeta::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(q_, sqrt(k0_));
+        bound(zeta_, epsilon0_/(2*sqrt(k0_)));
+    }
+
     RASModel::correct();
 
     if (!turbulence_)
diff --git a/src/turbulenceModels/incompressible/RAS/realizableKE/realizableKE.C b/src/turbulenceModels/incompressible/RAS/realizableKE/realizableKE.C
index 238af0651..d473b3b6f 100644
--- a/src/turbulenceModels/incompressible/RAS/realizableKE/realizableKE.C
+++ b/src/turbulenceModels/incompressible/RAS/realizableKE/realizableKE.C
@@ -262,6 +262,14 @@ bool realizableKE::read()
 
 void realizableKE::correct()
 {
+    // Bound in case of topological change
+    // HJ, 22/Aug/2007
+    if (mesh_.changing())
+    {
+        bound(k_, k0_);
+        bound(epsilon_, epsilon0_);
+    }
+
     RASModel::correct();
 
     if (!turbulence_)