[DuarteMagalhaes]

Hello everyone,

I am using reactingFoam solver and am trying to fully understand the Ueqn file, which is a subject I have never found thoroughly explained:

"
fvm::ddt(rho, U)
+ fvm::div(phi, U)
+ turbulence->divDevRhoReff(U)
==
rho*g
+ fvOptions(rho, U)
"

First entry is for time variation, second for convection, and what about third and last entries, turbulence->divDevRhoReff(U), and fvOptions(rho, U), what is it doing exactly?

And how does it include diffusion when solving in laminar regime? It would be expected to see a laplacian term, I guess.

Any help is much appreciated!
Thanks in advance