/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         cellLimited Gauss linear 1;
    grad(p)         cellLimited Gauss linear 1;
    grad(U)         cellLimited Gauss linear 1;
    grad(T)         cellLimited Gauss linear 1;
    grad(e)         cellLimited Gauss linear 1;
    grad(h)         cellLimited Gauss linear 1;
    grad(k)         cellLimited Gauss linear 1;
    grad(omega)     cellLimited Gauss linear 1;
}

divSchemes
{
    default                                 none;

    // div(phi,U)                              bounded Gauss limitedLinearV 0.5;
    // div(phi,e)                              bounded Gauss limitedLinear 0.5;
    // div(phi,h)                              bounded Gauss limitedLinear 0.5;
    // div(phi,he)                             bounded Gauss limitedLinear 0.5;
    // div(phi,K)                              bounded Gauss limitedLinear 0.5;
    // div(phi,Ekp)                            bounded Gauss limitedLinear 0.5;

    // div(phi,k)                              bounded Gauss limitedLinear 0.5;
    // div(phi,omega)                          bounded Gauss limitedLinear 0.5;

    // Dissipative for the first convergence
    div(phi,U)      bounded Gauss upwind;
    div(phi,e)      bounded Gauss upwind;
    div(phi,h)      bounded Gauss upwind;
    div(phi,he)     bounded Gauss upwind;
    div(phi,K)      bounded Gauss upwind;
    div(phi,Ekp)    bounded Gauss upwind;

    div(phi,k)      bounded Gauss upwind;
    div(phi,omega)  bounded Gauss upwind;

    div(phid,p)                             Gauss upwind;
    div(phiv,p)                             Gauss upwind;

    // div(((rho*nuEff)*dev2(T(grad(U)))))     Gauss linear;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; // Dissipative for the first convergence
    div((muEff*dev2(T(grad(U)))))           Gauss linear;
    div((nuEff*dev2(T(grad(U)))))           Gauss linear;

    div((nuEff*dev(T(grad(U)))))            Gauss linear;
    div((muEff*dev(T(grad(U)))))            Gauss linear;
    div(((rho*nuEff)*dev(T(grad(U)))))      Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited 0.5;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.5;
}

wallDist
{
    method          meshWave;
}

// ************************************************************************* //
