From d2d3e8710fe2e55f67f1d9e9cd6da0441fd31fe0 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Fri, 14 Apr 2023 10:18:23 +0100 Subject: [PATCH] compressibleVoF,compressibleMultiphaseVoF: Moved Uf correction outside the PISO loop for efficiency --- .../modules/compressibleMultiphaseVoF/pressureCorrector.C | 6 +++--- .../solvers/modules/compressibleVoF/pressureCorrector.C | 6 +++--- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/applications/solvers/modules/compressibleMultiphaseVoF/pressureCorrector.C b/applications/solvers/modules/compressibleMultiphaseVoF/pressureCorrector.C index d928824c18..ee305af28a 100644 --- a/applications/solvers/modules/compressibleMultiphaseVoF/pressureCorrector.C +++ b/applications/solvers/modules/compressibleMultiphaseVoF/pressureCorrector.C @@ -152,9 +152,6 @@ void Foam::solvers::compressibleMultiphaseVoF::pressureCorrector() } } - // Correct Uf if the mesh is moving - fvc::correctUf(Uf, U, fvc::absolute(phi, U), MRF); - // Update densities from change in p_rgh mixture.correctRho(p_rgh - p_rgh_0); mixture.correct(); @@ -164,6 +161,9 @@ void Foam::solvers::compressibleMultiphaseVoF::pressureCorrector() p_rgh.correctBoundaryConditions(); } + // Correct Uf if the mesh is moving + fvc::correctUf(Uf, U, fvc::absolute(phi, U), MRF); + K = 0.5*magSqr(U); clearrAU(); diff --git a/applications/solvers/modules/compressibleVoF/pressureCorrector.C b/applications/solvers/modules/compressibleVoF/pressureCorrector.C index 18289f08ae..945f45f1c2 100644 --- a/applications/solvers/modules/compressibleVoF/pressureCorrector.C +++ b/applications/solvers/modules/compressibleVoF/pressureCorrector.C @@ -211,9 +211,6 @@ void Foam::solvers::compressibleVoF::pressureCorrector() } } - // Correct Uf if the mesh is moving - fvc::correctUf(Uf, U, fvc::absolute(phi, U), MRF); - // Update densities from change in p_rgh mixture_.thermo1().correctRho(psi1*(p_rgh - p_rgh_0)); mixture_.thermo2().correctRho(psi2*(p_rgh - p_rgh_0)); @@ -224,6 +221,9 @@ void Foam::solvers::compressibleVoF::pressureCorrector() p_rgh.correctBoundaryConditions(); } + // Correct Uf if the mesh is moving + fvc::correctUf(Uf, U, fvc::absolute(phi, U), MRF); + K = 0.5*magSqr(U); clearrAU();