[sourav8016]

I want to code this equation in interFoam solver.


where c is the sediment volumetric concentration, is sediment settling velocity defined by


is gravitional accleration, is kinematic viscosity and is particle diameter and is relative density of sediment.



is sediment diffusivity and
is eddy viscosity, is turbulent prandtl schimdt number.

Any idea how I can code this interFoam?

I was trying to include a cEqn.H file in interFoam solver as

volScalarField epsilonEff
(
"epsilonEff",
4*turbulence->nut()/Prt
);

fvScalarMatrix cEqn
(
fvm::ddt(c)
+ fvm::div(phi, c)
- fvm::laplacian(epsilonEff, c)
);

cEqn.relax();

cEqn.solve();
}

But I don't know how to modify phi and include here as a closuse.

please help me. Thanks in advance.