From d1c8707c3e1e6261064320f9780135cae73a057f Mon Sep 17 00:00:00 2001 From: danielque Date: Wed, 30 Aug 2017 13:18:31 +0200 Subject: [PATCH] fix issue with voidfraction slightly > 1 in SyamlalThermCond --- .../thermCondModel/SyamlalThermCond/SyamlalThermCond.C | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/lagrangian/cfdemParticle/subModels/thermCondModel/SyamlalThermCond/SyamlalThermCond.C b/src/lagrangian/cfdemParticle/subModels/thermCondModel/SyamlalThermCond/SyamlalThermCond.C index 7a08b85b..0e2b7e9d 100644 --- a/src/lagrangian/cfdemParticle/subModels/thermCondModel/SyamlalThermCond/SyamlalThermCond.C +++ b/src/lagrangian/cfdemParticle/subModels/thermCondModel/SyamlalThermCond/SyamlalThermCond.C @@ -104,9 +104,9 @@ tmp SyamlalThermCond::thermCond() const ) ); volScalarField& svf = tvf.ref(); - - svf = (1-sqrt(1-voidfraction_)) / (voidfraction_) * kf0_; - + + svf = (1-sqrt(1-voidfraction_+SMALL)) / (voidfraction_) * kf0_; + // if a wallQFactor field is present, use it to scale heat transport through a patch if (hasWallQFactor_) {