PaSR: ensure the sqrt operates on a positive number

2015-03-06 22:38:12 +00:00
parent 6884300a00
commit 0003b4b008
1 changed files with 1 additions and 1 deletions
--- a/src/combustionModels/PaSR/PaSR.C
+++ b/src/combustionModels/PaSR/PaSR.C
@ -84,7 +84,7 @@ void Foam::combustionModels::PaSR<Type>::correct()
            forAll(epsilon, i)
            {
                scalar tk =
-                    Cmix_*Foam::sqrt(muEff[i]/rho[i]/(epsilon[i] + SMALL));
+                    Cmix_*sqrt(max(muEff[i]/rho[i]/(epsilon[i] + SMALL), 0));

                if (tk > SMALL)
                {