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

Merge branch 'master' of github.com-OpenFOAM:OpenFOAM/OpenFOAM-dev

Browse Source
This commit is contained in:
Henry Weller
2018-03-20 22:29:57 +00:00
parent a5d6778281 26f2ce9c9e
commit 84dd9a5047
1 changed files with 3 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
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",
( min(alpha, scalar(1))*fvc::div(fvc::absolute(phi_, phase.U()))
"Su",
mesh_.time().timeName(),
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())