[LuisAlberto]

Hello everyone!

Looking at pisoFoam solver, I am seeing now that, if I am not mistaken, the transport equations for the turbulent kinetic energy (k) and dissipation rate (epsilon) are solved outside piso algorithm loops, i.e., the sentence turbulence->correct() is called upon only when the velocity and pressure fields have already converged.

This implies that the equations are written using the nut field computed in the previous time step, and when using wall functions, the term kp used is also from the previous time step.

Is this correct or there is something I am not understanding correctly? Could the transport equations be solved inside the loop and checked for convergence too, as in Verteeg and Malasakera?