[shouvik ghorui]

Hi everyone,
I want to implement this(below) scalar transport equation. where Rz is a source term, expressed as below where u_n and u both are particle velocity only.
I have implemented
fvm::ddt(C)
+fvm::div(phiup,C)
-fvm::laplacian(gammaPsi,C)
==
(1/rp)*gammaPsi*Foam::mag(fvc::grad(C))
-fvm::SuSp(Foam::mag(up)/rp,C)
my confusion is related to the source term, Is there any other way to work with this source term. (taking mag of velocity is okay?)
I am not getting the correct results using this. Any suggestion would be helpful. Thanks.




note:
I am getting the result(wrong) which may be due to the large sink term(-fvm::SuSp(Foam::mag(up)/rp,C).