[qi.yang@polimi.it]

Hi guys, I need to add a term into the mass conservation equation for twophaseeulerfoam, as shown in following codes which are in the p-equation. I hope I added it in a right position. I think I should use plus sign for this term however the simulation crashed when I use "plus". While I use "minus". It goes well without problem. Do you have any suggestion?
pEqnComp1 =
(
contErr1
- fvc::Sp(fvc::ddt(alpha1) + fvc::div(alphaPhi1), rho1)
)/rho1
+ (alpha1*psi1/rho1)*correction(fvm::ddt(p_rgh)) + fvc::laplacian(phase2.turbulence().nut()/mysigma,alpha1);


pEqnComp2 =
(
contErr2
- fvc::Sp(fvc::ddt(alpha2) + fvc::div(alphaPhi2), rho2)
)/rho2
+ (alpha2*psi2/rho2)*correction(fvm::ddt(p_rgh)) + fvc::laplacian(phase2.turbulence().nut()/mysigma,alpha2);