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

reactingMultiphaseEulerFoam: Simplified the velocity and flux correction after partial elimination

Browse Source
This commit is contained in:
Henry Weller
2020-07-10 12:04:08 +01:00
parent 2e62bfe5bc
commit 5c1f24054d
2 changed files with 21 additions and 106 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
applications/solvers/multiphase/reactingEulerFoam/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
View File

@ -1210,7 +1210,7 @@ void Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::partialElimination
movingPhasei++
)
{
Um.ref() +=
Um +=
this->movingPhases()[movingPhasei]
*this->movingPhases()[movingPhasei].U();
}

125
applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseSystem/phaseSystem.C
View File

@ -145,55 +145,22 @@ Foam::tmp<Foam::volScalarField> Foam::phaseSystem::sumAlphaMoving() const
}
void Foam::phaseSystem::setMixtureU(const volVectorField& Um)
void Foam::phaseSystem::setMixtureU(const volVectorField& Um0)
{
phaseSystem::phaseModelList& phases = phaseModels_;
// Calculate the mean velocity difference with respect to Um0
// from the current velocity of the moving phases
volVectorField dUm(Um0);
PtrList<volVectorField> U0s(phases.size());
forAll(phases, i)
forAll(movingPhaseModels_, movingPhasei)
{
if (!phases[i].stationary())
{
U0s.set
(
i,
volVectorField::New
(
IOobject::groupName("U0", phases[i].name()),
phases[i].URef()
)
);
}
dUm -=
movingPhaseModels_[movingPhasei]
*movingPhaseModels_[movingPhasei].U();
}
forAll(phases, i)
forAll(movingPhaseModels_, movingPhasei)
{
if (!phases[i].stationary())
{
phases[i].URef() = Um;
forAll(phases, j)
{
if (i != j && !phases[j].stationary())
{
phases[i].URef() += phases[j]*(U0s[i] - U0s[j]);
}
}
}
}
if (!stationaryPhaseModels_.empty())
{
const volScalarField sumAlphaMoving(this->sumAlphaMoving());
forAll(phases, i)
{
if (!phases[i].stationary())
{
phases[i].URef() /= sumAlphaMoving;
}
}
movingPhaseModels_[movingPhasei].URef() += dUm;
}
}
@ -201,75 +168,23 @@ void Foam::phaseSystem::setMixtureU(const volVectorField& Um)
void Foam::phaseSystem::setMixturePhi
(
const PtrList<surfaceScalarField>& alphafs,
const surfaceScalarField& phim
const surfaceScalarField& phim0
)
{
phaseSystem::phaseModelList& phases = phaseModels_;
// Calculate the mean flux difference with respect to phim0
// from the current flux of the moving phases
surfaceScalarField dphim(phim0);
PtrList<surfaceScalarField> phi0s(phases.size());
forAll(phases, i)
forAll(movingPhaseModels_, movingPhasei)
{
if (!phases[i].stationary())
{
phi0s.set
(
i,
surfaceScalarField::New
(
IOobject::groupName("phi0", phases[i].name()),
phases[i].phiRef()
)
);
}
dphim -=
alphafs[movingPhaseModels_[movingPhasei].index()]
*movingPhaseModels_[movingPhasei].phi();
}
forAll(phases, i)
forAll(movingPhaseModels_, movingPhasei)
{
if (!phases[i].stationary())
{
phases[i].phiRef() = phim;
forAll(phases, j)
{
if (i != j && !phases[j].stationary())
{
phases[i].phiRef() += alphafs[j]*(phi0s[i] - phi0s[j]);
}
}
}
}
if (!stationaryPhaseModels_.empty())
{
surfaceScalarField sumAlphafMoving
(
surfaceScalarField::New
(
"sumAlphafMoving",
alphafs[movingPhaseModels_[0].index()],
calculatedFvsPatchScalarField::typeName
)
);
for
(
label movingPhasei=1;
movingPhasei<movingPhaseModels_.size();
movingPhasei++
)
{
sumAlphafMoving +=
alphafs[movingPhaseModels_[movingPhasei].index()];
}
forAll(phases, i)
{
if (!phases[i].stationary())
{
phases[i].phiRef() /= sumAlphafMoving;
}
}
movingPhaseModels_[movingPhasei].phiRef() += dphim;
}
}