BUG: (kklOmega.C) limiting exp to avoid under flow

2025-11-28 03:28:01 +00:00 · 2012-05-24 14:30:52 +01:00
parent 78b5630da7
commit 2395d810d4
1 changed files with 17 additions and 8 deletions
--- a/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C
+++ b/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C
@ -155,16 +155,25 @@ tmp<volScalarField> kkLOmega::fOmega
 tmp<volScalarField> kkLOmega::gammaBP(const volScalarField& omega) const
 {
    return
+    (
+        min
        (
            max
            (
                kt_/nu()
             / (
                    omega
-              + dimensionedScalar("ROTVSMALL", omega.dimensions(), ROOTVSMALL)
+                  + dimensionedScalar
+                    (
+                        "ROTVSMALL",
+                        omega.dimensions(),
+                        ROOTVSMALL
+                    )
                )
              - CbpCrit_,
                scalar(0)
+            ),
+            scalar(50.0)
        )
    );
 }