BUG: parcel phase change - bypass calc if phase fraction < SMALL - mantis #1298

This commit is contained in:
andy
2014-05-19 17:19:36 +01:00
committed by Andrew Heather
parent 876df23704
commit 98b9eaac98

View File

@ -60,14 +60,14 @@ void Foam::ReactingParcel<ParcelType>::calcPhaseChange
typedef typename TrackData::cloudType::reactingCloudType reactingCloudType;
PhaseChangeModel<reactingCloudType>& phaseChange = td.cloud().phaseChange();
if (!phaseChange.active())
if (!phaseChange.active() || (YPhase < SMALL))
{
return;
}
scalar Tvap = phaseChange.Tvap(YComponents);
if (T < Tvap || YPhase < SMALL)
if (T < Tvap)
{
return;
}