[han-1]

I was simulation pure species transport using reactingFoam(turn reactions off)
First, the divSchemes were set as follows:
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,Yi_h) Gauss upwind;
div(phi,h) Gauss limitedLinear 1;
div(phi,K) Gauss upwind;
div(phi,Yi) Gauss upwind;
div(phid,p) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,k) Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}
The solution diverges in a few time steps.Then i swapped the position of 'div(phi,Yi_h)' and 'div(phi,Yi)’. It works !!!!!
WHY? Someone can give an explanation to this?

[Alczem]

Hey


Can you post the whole case? It might help to reproduce this issue and seek out an explanation.

[hjasak]

For consistency, you should use the same cenvection scheme on all species and the energy equation - this is why it says Yi (species) _ h (enthalpy).


This is called a multi-variate scheme. Maybe the format has changed ot somebudy messes with your solver.


Hrv