[virginie_e]

Hello,

I have a question concerning the SIMPLE loop in interTrackFoam, especially about the velocity equation.

In interTrackFoam, UEqn is implemented this way:

tmp<fvvectormatrix> UEqn
(
fvm::ddt(rho, U)
+ fvm::div(phiNet, U)
- fvm::laplacian(mu, U)
);

UEqn().relax();

solve(UEqn() == - fvc::grad(p));


whereas in other solvers using the SIMPLE algorithm, it is implemented this way:

tmp<fvvectormatrix> UEqn
(
fvm::ddt(U)
+fvm::div(phi,U)
-fvm::laplacian(nu,U)
);
UEqn.relax();
solve (UEqn == -fvc::grad(p));


I understand that the difference is that the interTrackFoam equation is the second one multiplied by rho, but in that case, should the second term -fvc::grad(p) not be multiplied by rho as well? so that the equation should become something like:

tmp<fvvectormatrix> UEqn
(
fvm::ddt(rho, U)
+ fvm::div(phiNet, U)
- fvm::laplacian(mu, U)
);

UEqn().relax();

solve(UEqn() == - rho()*fvc::grad(p));

or something like that?

If what I said is wrong, is there something that I misunderstood in the equation solving in OpenFOAM?

Otherwise, how should I write the equivalent of the line

solve(UEqn() == - rho()*fvc::grad(p));

?

Thank you in advance