From 2395d810d415ff109e7b576e7151e136650f2fd6 Mon Sep 17 00:00:00 2001 From: sergio Date: Thu, 24 May 2012 14:30:52 +0100 Subject: [PATCH] BUG: (kklOmega.C) limiting exp to avoid under flow --- .../incompressible/RAS/kkLOmega/kkLOmega.C | 25 +++++++++++++------ 1 file changed, 17 insertions(+), 8 deletions(-) diff --git a/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C b/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C index 7dd7adb80e..612e8e25b4 100644 --- a/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C +++ b/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C @@ -156,15 +156,24 @@ tmp kkLOmega::gammaBP(const volScalarField& omega) const { return ( - max + min ( - kt_/nu() - / ( - omega - + dimensionedScalar("ROTVSMALL", omega.dimensions(), ROOTVSMALL) - ) - - CbpCrit_, - scalar(0) + max + ( + kt_/nu() + / ( + omega + + dimensionedScalar + ( + "ROTVSMALL", + omega.dimensions(), + ROOTVSMALL + ) + ) + - CbpCrit_, + scalar(0) + ), + scalar(50.0) ) ); }