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

compressibleVoF,compressibleMultiphaseVoF: Moved Uf correction outside the PISO loop

Browse Source 
for efficiency
This commit is contained in:
Henry Weller
2023-04-14 10:18:23 +01:00
parent 5e8d1ff99e
commit d2d3e8710f
2 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
applications/solvers/modules/compressibleMultiphaseVoF/pressureCorrector.C
View File

@ -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 // Update densities from change in p_rgh
mixture.correctRho(p_rgh - p_rgh_0); mixture.correctRho(p_rgh - p_rgh_0);
mixture.correct(); mixture.correct();
@ -164,6 +161,9 @@ void Foam::solvers::compressibleMultiphaseVoF::pressureCorrector()
p_rgh.correctBoundaryConditions(); p_rgh.correctBoundaryConditions();
} }
// Correct Uf if the mesh is moving
fvc::correctUf(Uf, U, fvc::absolute(phi, U), MRF);
K = 0.5*magSqr(U); K = 0.5*magSqr(U);
clearrAU(); clearrAU();

6
applications/solvers/modules/compressibleVoF/pressureCorrector.C
View File

@ -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 // Update densities from change in p_rgh
mixture_.thermo1().correctRho(psi1*(p_rgh - p_rgh_0)); mixture_.thermo1().correctRho(psi1*(p_rgh - p_rgh_0));
mixture_.thermo2().correctRho(psi2*(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(); p_rgh.correctBoundaryConditions();
} }
// Correct Uf if the mesh is moving
fvc::correctUf(Uf, U, fvc::absolute(phi, U), MRF);
K = 0.5*magSqr(U); K = 0.5*magSqr(U);
clearrAU(); clearrAU();