[PicklER]

Hi
In shockTube tutorial test case, the divSchemes looks like follows:

divSchemes
{
    default         none;
    div(tauMC)      Gauss linear;
}

And in solver rhoCentralFoam, solving for density:

solve(fvm::ddt(rho) + fvc::div(phi));

My question: what scheme is used to solve for phi? As I understand, if a scheme is not specified, OpenFOAM gives an error and suggests the possible schemes.

Thanks

[elones]

Hi,
look into rhoCentralFoam.C, the solver is based on Kurganov and Tadmor central schemes, ie, Riemann-solver-free approach.

[PicklER]

Hi
Thank you for the quick reply. I think I understand it now. phi is calculated by the interpolated variables of the cell faces and therefore solved. I noticed that for n number of cells in a straight line, like in the shocktube test case, phi has only n-1 values, which confirms that phi only exists on the faces. And the obvious definition of phi

surfaceScalarField phi("phi", mesh.Sf() & fvc::interpolate(rhoU));

The same is done for the momentum and energy equations. I found a good article on rhoCenralFoam and it explains how the convective terms in the conservation equations are handled:
here
I am busy writing up my understanding of it in a way I can better explain it.
Thanks again