[sergio]

Hi there,
I am following the equation for Psi for turbulent deflagration. It looks as:

solve
(
fvm::ddt(rho, Xi)
+ fvm::div(phi + phiXi, Xi, "div(phiXi,Xi)")
- fvm::Sp(fvc::div(phiXi), Xi)
==
rho*R
- fvm::Sp(rho*(R - G), Xi)
- fvm::Sp
(
rho*max
(
sigmat - sigmas,
dimensionedScalar("0", sigmat.dimensions(), 0)
),
Xi
)
);

I quiet understand all the terms but:

-fvm::Sp(fvc::div(phiXi), Xi)

This term is div(phiXi,Xi) but it is already included in fvm::div(phi + phiXi, Xi, "div(phiXi,Xi)")

and is after substracted?. There is something not clear to me about how this term is included in the equation.

thanks

Sergio