Added evaluation of turbulence boundary conditions after topo change: attach/detach bug

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

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

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

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

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

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

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

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

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

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

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_)

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

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_)

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

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_)

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_)

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_)

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

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_)

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

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_)

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_)

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_)

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_)

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_)