From 143800ab3d27b8d9379d8724267b4e17072e2a33 Mon Sep 17 00:00:00 2001
From: Andrew Heather <a.heather@opencfd.co.uk>
Date: Tue, 26 Jul 2016 15:33:10 +0100
Subject: [PATCH] BUG: Turbulence DFSEM - protect sigFpe for reconstruct - see
#194
---
.../turbulentDFSEMInlet/turbulentDFSEMInletFvPatchVectorField.C | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/finiteVolume/fields/fvPatchFields/derived/turbulentDFSEMInlet/turbulentDFSEMInletFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/turbulentDFSEMInlet/turbulentDFSEMInletFvPatchVectorField.C
index b0df476fe0d..59646f3ad86 100644
--- a/src/finiteVolume/fields/fvPatchFields/derived/turbulentDFSEMInlet/turbulentDFSEMInletFvPatchVectorField.C
+++ b/src/finiteVolume/fields/fvPatchFields/derived/turbulentDFSEMInlet/turbulentDFSEMInletFvPatchVectorField.C
@@ -848,7 +848,7 @@ turbulentDFSEMInletFvPatchVectorField
eddy::debug = debug;
// Set UMean as patch area average value
- UMean_ = gSum(U_*patch().magSf())/gSum(patch().magSf());
+ UMean_ = gSum(U_*patch().magSf())/(gSum(patch().magSf()) + ROOTVSMALL);
}
--
GitLab