From e8325bbcecdec4e54e5476d1bc9755478cb86c80 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Sun, 23 Oct 2022 10:13:04 +0100 Subject: [PATCH] combustionModels::EDC: Corrected handling of limits in the v2016 variant Resolves bug-report https://bugs.openfoam.org/view.php?id=3905 --- src/combustionModels/EDC/EDC.C | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/combustionModels/EDC/EDC.C b/src/combustionModels/EDC/EDC.C index f41f11bacb..c06300ede8 100644 --- a/src/combustionModels/EDC/EDC.C +++ b/src/combustionModels/EDC/EDC.C @@ -135,10 +135,10 @@ void Foam::combustionModels::EDC::correct() max(min(sqrt(nu[i]/(epsilon[i] + small))/tc[i], 10), 1e-10); const scalar ReT = sqr(k[i])/(nu[i]*epsilon[i] + small); - const scalar CtauI = min(C1_/(Da*sqrt(ReT + 1)), 2.1377); + const scalar CtauI = min(C1_/(Da*sqrt(ReT + 1)), Ctau_); const scalar CgammaI = - max(min(C2_*sqrt(Da*(ReT + 1)), 5), 0.4082); + max(min(C2_*sqrt(Da*(ReT + 1)), 5), Cgamma_); const scalar gammaL = CgammaI*pow025(nu[i]*epsilon[i]/(sqr(k[i]) + small));