[Balakrshnan Ramakrishnan]

Dear Foamers,
I am trying to understand the code for simpleFoam and I came across the wiki for the same

http://openfoamwiki.net/index.php/The_SIMPLE_algorithm_in_OpenFOAM

In the above link , in the "Implementation of the SIMPLE algorithm in OpenFOAM" section,

the second equation says that the

tmp<fvVectorMatrix> UEqn
(
fvm::div(phi, U) - laplacian(nu,U)
);

But in the actual file UEqn.H in the /solvers folder , the code is as below

tmp<fvVectorMatrix> UEqn
(
fvm::div(phi, U)
+ turbulence->divDevReff(U)
);

Is the laplacian term replaced with a turbulence term??

Can someone help to figure this out?

Your help is highly appreciated.

Thanks

[santiagomarquezd]

Hi, simpleFoam manages both laminar and turbulent cases via turbulence->divDevReff(U). This term reduces to the laplacian one in laminar cases.

Regards.

[Balakrshnan Ramakrishnan]

Thanks a lot for your reply