From 8599a55d40716ddb9ae411110ea47b5df97473d2 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 7554e0811e6..9e14ce1d27a 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)
{
--
GitLab