The non linear stresses from Baglietto & al 2006 and as written now in Code_Saturne correspond to rho*Rij, the rho is taken account directly in the function divrij in base/predvv.f90 because itypfl is set to 1.
Therefore to obtain the Reynold stresses and not rho*Reynold Stresses, a division by rho is needed
The non linear stresses from Baglietto & al 2006 and as written now in Code_Saturne correspond to rho*Rij, the rho is taken account directly in the function divrij in base/predvv.f90 because itypfl is set to 1.
Therefore to obtain the Reynold stresses and not rho*Reynold Stresses, a division by rho is needed