[Bazinga]

In the pEqn.H of pimpleFoam, face fluxes are corrected after the momentum correction in case of moving meshes in the lines 64ff
https://cpp.openfoam.org/v6/incompre...8H_source.html

// Correct Uf if the mesh is moving
 fvc::correctUf(Uf, U, phi);
 
 // Make the fluxes relative to the mesh motion
 fvc::makeRelative(phi, U);

It is my understanding that mesh motion is calculated before the momentum predictor based on the old flow field and the influence of the moving mesh on face fluxes is accounted for by space conservation law in pimpleFoam.C line 80-96.
https://cpp.openfoam.org/v6/pimpleFoam_8C_source.html

                 mesh.update();
 
                 if (mesh.changing())
                 {
                     MRF.update();
 
                     if (correctPhi)
                     {
                         // Calculate absolute flux
                         // from the mapped surface velocity
                         phi = mesh.Sf() & Uf();
 
                         #include "correctPhi.H"
 
                         // Make the flux relative to the mesh motion
                         fvc::makeRelative(phi, U);
                     }

It is not quite clear to me why the face fluxes are corrected again in pEqn.H. Can someone enlighten me?