zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

reactingMultiphaseEulerFoam: Prevent dereference of non-existent drag fields

Browse Source
This commit is contained in:
Will Bainbridge
2018-03-12 12:18:07 +00:00
parent f46fad4343
commit dce56828e7
2 changed files with 16 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pU/pEqn.H
View File

@ -137,8 +137,14 @@ while (pimple.correct())
forAll(phases, phasei)
{
HbyAs[phasei] -= KdUByAs[phasei];
phiHbyAs[phasei] -= phiKdPhis[phasei];
if (KdUByAs.set(phasei))
{
HbyAs[phasei] -= KdUByAs[phasei];
}
if (phiKdPhis.set(phasei))
{
phiHbyAs[phasei] -= phiKdPhis[phasei];
}
}
}
@ -169,7 +175,10 @@ while (pimple.correct())
forAll(phases, phasei)
{
phiHbyA -= alphafs[phasei]*phiKdPhis[phasei];
if (phiKdPhis.set(phasei))
{
phiHbyA -= alphafs[phasei]*phiKdPhis[phasei];
}
}
}

5
applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pUf/pEqn.H
View File

@ -169,7 +169,10 @@ while (pimple.correct())
forAll(phases, phasei)
{
phiHbyA -= alphafs[phasei]*phiKdPhifs[phasei];
if (phiKdPhifs.set(phasei))
{
phiHbyA -= alphafs[phasei]*phiKdPhifs[phasei];
}
}
}