From fb332ca5ec1b6d258780c974e197c7d22e7350cd Mon Sep 17 00:00:00 2001 From: andy Date: Thu, 28 Oct 2010 18:12:06 +0100 Subject: [PATCH] ENH: Added stabilisation to parcel analaytical integration scheme --- .../intermediate/IntegrationScheme/Analytical/Analytical.C | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C index 7075f2152a..fa96bdec0a 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C +++ b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C @@ -69,7 +69,8 @@ Foam::Analytical::integrate const scalar expTerm = exp(min(50, -beta*dt)); - retValue.average() = alpha + (phi - alpha)*(1 - expTerm)/(beta*dt); + retValue.average() = + alpha + (phi - alpha)*(1 - expTerm)/(beta*dt + ROOTVSMALL); retValue.value() = alpha + (phi - alpha)*expTerm; return retValue;