fix issue with voidfraction slightly > 1 in SyamlalThermCond

This commit is contained in:
danielque
2017-08-30 13:18:31 +02:00
parent 91b84ca607
commit d1c8707c3e

View File

@ -104,9 +104,9 @@ tmp<volScalarField> 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_)
{