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

multiphaseEuler::phaseSystem::correctBoundaryFlux: Update the phase fluxes from the phase face-velocity and make relative

Browse Source 
This update is necessary for mesh topology changes and mesh-to-mesh mapping
which update the phase face-velocity fields but not the fluxes.
This commit is contained in:
Henry Weller
2024-05-23 15:41:24 +01:00
parent f9d03d6ad9
commit d444cb168d
1 changed files with 17 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
applications/modules/multiphaseEuler/phaseSystem/phaseSystem/phaseSystem.C
View File

@ -689,9 +689,9 @@ void Foam::phaseSystem::meshUpdate()
void Foam::phaseSystem::correctBoundaryFlux() void Foam::phaseSystem::correctBoundaryFlux()
{ {
forAll(movingPhases(), movingPhasei) forAll(movingPhaseModels_, movingPhasei)
{ {
phaseModel& phase = movingPhases()[movingPhasei]; phaseModel& phase = movingPhaseModels_[movingPhasei];
tmp<volVectorField> tU(phase.U()); tmp<volVectorField> tU(phase.U());
const volVectorField::Boundary& UBf = tU().boundaryField(); const volVectorField::Boundary& UBf = tU().boundaryField();
@ -726,9 +726,18 @@ void Foam::phaseSystem::correctPhi
nonOrthogonalSolutionControl& pimple nonOrthogonalSolutionControl& pimple
) )
{ {
forAll(movingPhases(), movingPhasei) // Update the phase fluxes from the phase face-velocity and make relative
forAll(movingPhaseModels_, movingPhasei)
{ {
phaseModel& phase = movingPhases()[movingPhasei]; phaseModel& phase = movingPhaseModels_[movingPhasei];
phase.phiRef() = mesh_.Sf() & phase.UfRef();
MRF_.makeRelative(phase.phiRef());
fvc::makeRelative(phase.phiRef(), phase.U());
}
forAll(movingPhaseModels_, movingPhasei)
{
phaseModel& phase = movingPhaseModels_[movingPhasei];
volVectorField::Boundary& Ubf = phase.URef().boundaryFieldRef(); volVectorField::Boundary& Ubf = phase.URef().boundaryFieldRef();
surfaceVectorField::Boundary& UfBf = phase.UfRef().boundaryFieldRef(); surfaceVectorField::Boundary& UfBf = phase.UfRef().boundaryFieldRef();
@ -756,9 +765,9 @@ void Foam::phaseSystem::correctPhi
phi_ = Zero; phi_ = Zero;
PtrList<surfaceScalarField> alphafs(phaseModels_.size()); PtrList<surfaceScalarField> alphafs(phaseModels_.size());
forAll(movingPhases(), movingPhasei) forAll(movingPhaseModels_, movingPhasei)
{ {
phaseModel& phase = movingPhases()[movingPhasei]; phaseModel& phase = movingPhaseModels_[movingPhasei];
const label phasei = phase.index(); const label phasei = phase.index();
const volScalarField& alpha = phase; const volScalarField& alpha = phase;
@ -774,7 +783,7 @@ void Foam::phaseSystem::correctPhi
fv::correctPhi fv::correctPhi
( (
phi_, phi_,
movingPhases()[0].U(), movingPhaseModels_[0].U(),
p_rgh, p_rgh,
autoPtr<volScalarField>(), autoPtr<volScalarField>(),
divU, divU,
@ -815,7 +824,7 @@ void Foam::phaseSystem::correctPhi
// Make the flux relative to the mesh motion // Make the flux relative to the mesh motion
MRF_.makeRelative(phi_); MRF_.makeRelative(phi_);
fvc::makeRelative(phi_, movingPhases()[0].U()); fvc::makeRelative(phi_, movingPhaseModels_[0].U());
setMixturePhi(alphafs, phi_); setMixturePhi(alphafs, phi_);
} }