zhyang-dev/OpenFOAM-6
1
0
Fork 0
mirror of https://github.com/OpenFOAM/OpenFOAM-6.git synced 2025-12-08 06:57:46 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

reactingMultiphaseEulerFoam: Bug fix to divergence handling

Browse Source
This commit is contained in:
Will Bainbridge
2018-03-14 17:06:35 +00:00
parent ea00c1a4c5
commit 26f2ce9c9e
1 changed files with 3 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C
View File

@ -189,8 +189,6 @@ void Foam::multiphaseSystem::solveAlphas()
); );
volScalarField divU(fvc::div(fvc::absolute(phi_, phases().first().U())));
forAll(phases(), phasei) forAll(phases(), phasei)
{ {
phaseModel& phase = phases()[phasei]; phaseModel& phase = phases()[phasei];
@ -209,20 +207,13 @@ void Foam::multiphaseSystem::solveAlphas()
mesh_ mesh_
), ),
mesh_, mesh_,
dimensionedScalar("Sp", divU.dimensions(), 0.0) dimensionedScalar("Sp", dimless/dimTime, 0)
); );
volScalarField::Internal Su volScalarField::Internal Su
(
IOobject
( (
"Su", "Su",
mesh_.time().timeName(), min(alpha, scalar(1))*fvc::div(fvc::absolute(phi_, phase.U()))
mesh_
),
// Divergence term is handled explicitly to be
// consistent with the explicit transport solution
divU*min(alpha, scalar(1))
); );
if (phase.divU().valid()) if (phase.divU().valid())