[fredo490]

Dear everyone,
I'm looking for a way of converting my continuity equation in an implicit form. This equation correspond to a particle tracking using an Eulerian approach. "rho" is not a gas, it is my particle concentration.

Right now I have :
Continuity equation :

	Info<< "Solve Continuity" << endl;
	tmp<fvScalarMatrix> CEqn
	(
	  fvm::ddt(rho)
	  + fvc::div(phi)
	);
	CEqn().solve();

I want to convert "fvc::div(phi)" in an implicit form.

I'm a little bit confused by the OpenFoam documentation : http://www.foamcfd.org/Nabla/guides/...sGuidese9.html

It says that the convection term can be Imp/Exp with the following form: div(psi,scheme)* where psi is a surfaceScalarField (just like my "phi").

However, if I try the following code, I get a mistake while compiling:

 word scheme("div(phi)");

	Info<< "Solve Continuity" << endl;
	tmp<fvScalarMatrix> CEqn
	(
	  fvm::ddt(rho)
	  + fvm::div(phi, scheme)
	);
	CEqn().solve();

Does anybody has an idea ?
Thanks in advance.