[mgao1]

Hi there,

I am recently modifying the denseParticleFoam (DPFoam) and I want to convert it to a conserved solver.

The original code is:
1 fvVectorMatrix UcEqn
2 (
3 fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)
4 - fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc)
5 + continuousPhaseTurbulence->divDevTau(Uc)
6 ==
7 (1.0/rhoc)*cloudSU
8 );

In my opinion, the 4th line will make the momentum equation a non-consered form. Therefore, when I modified the solver, I think I should remove this line, and the code will look like:

alphaRhoc = rhoc*alphac;
alphaRhocf = fvc::interpolate(alphaRhoc);
alphaRhoPhic = alphaRhocf*phic;

9 fvm::ddt(alphaRhoc, Uc) + fvm::div(alphaRhoPhic, Uc)
10 + turbulence->divDevTau(alphaRhoc,Uc)
11 ==
12 cloudSU

However, this will return me different results in all the fields (such as velocity, pressure, etc...) due to the 4th will not be a zero filed.
My question is, should I remove the 4th line of code to make it conserved form? If yes, is there any other code I should add to get a similar result?

Any help will be much appreciated!

Regards,
Mingze