Continuity equation in implicit form
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 : Code:
Info<< "Solve Continuity" << endl; 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: Code:
word scheme("div(phi)"); Thanks in advance. 
Did you already try fvm::div(phi,rho)? I think that should do the trick for you.

Hi Frédéric,
Defining two different phi's would indeed be the way to go. Otherwise, you would end up with two unknowns in the continuity equation ( and ). You would not be able to solve for both. I would be careful with the definition of phiRhoU. It might be better to do an interpolation of U * Rho, but this is something you could test of course. How do you solve the momentum equation? Is given from solving the continuity equation? I have one small remark. You seem to be using both rho and alpha for the volume fraction of the particles. Correct me if I am wrong. I think it would be better to be consistent to avoid confusion. Good luck. 
I have applied your method and it seems to work perfectly. I made a short run and the results look similar to my explicit solver.
To gives some details, here is roughly what happen in my code Code:
//// Thanks again ;) 
What happens if you try
Code:
phiR = fvc::interpolate(rho*U) & mesh.Sf() 
It also works and it simplify my code :D
I have another question related to the field management. Right now I have a RelativeReynolds number stored in a volScalarField and a drag value in another volScalarField. Until now, I use the following code : Code:
// Compute the Relative Reynolds number field Quote:

All times are GMT 4. The time now is 06:08. 