research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

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

Browse source
This commit is contained in:
andy 2014-05-19 17:19:36 +01:00
parent 0eb39ebe0f
commit 71ff5ebd12
1 changed files with 2 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
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;
}