[annan]

Dear foamers,

I have built a scalar transport solver to solve the following equation :

fvScalarMatrix psiEqn
{
fvm::ddt(psi)
+fvm::div(phi, psi)
+fvm::div(phiCf, psi)
-fvm::laplacian(D, psi)
};

phiCf is an additional flux, the solver gives me the results I was expecting. The problem is that I need the total diffusive flux which is equal to (D*div(psi)-div(phiCf*psi)), and I thought of using psiEqn.flux() but I don't know what this function gives as a results.

using psiEqn.flux() :
- do I get only convective flux (div(phi,psi)) ?
- do I get only diffusive flux (D*div(psi)) ?
- do I get the total diffusive flux (D*div(psi)-div(phiCf*psi)) ?
- or the total flux (D*div(psi)-div(phiCf*psi)-div(phi,psi)) ?

Thank you in advance for your help.
Annan

[mAlletto]

See this post

Description of flux() method

[annan]

Originally Posted by mAlletto

See this post

Description of flux() method

Thank you very much mAlletto for your answer
I have in fact seen that post, and I see that pEqn.flux() computes the diffusive flux, but as it was applied to poisson equation, I don't know if in my case it takes account of the convective fluxes also or not ...

[mAlletto]

Flux is a public member Funktion of the fvMatrix class