From 0003b4b0081a53b461a7436a160ea60cdee5c2c4 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Fri, 6 Mar 2015 22:38:12 +0000
Subject: [PATCH] PaSR: ensure the sqrt operates on a positive number

---
 src/combustionModels/PaSR/PaSR.C | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/combustionModels/PaSR/PaSR.C b/src/combustionModels/PaSR/PaSR.C
index 7554e0811..9e14ce1d2 100644
--- 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)
                 {