[deepinheart]

I run a case using simpleFoam in OF 230 and everything goes well. When I turn to OF ext 3.1, I can't find 'bouned Gauss upwind' scheme as before. Instead, I use 'Gauss upwind' but it failed.

Anybody has similar experience? Would you please explain which one is the similar to 'bounded Gauss upwind'?

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

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
    // added for Of ext3.1
    grad(p)         Gauss linear;
    grad(U)         Gauss linear;
}

divSchemes
{
    default         none;
    //div(phi,U)      bounded Gauss upwind;
    //div(phi,k)      bounded Gauss upwind;
    //dv(phi,epsilon) bounded Gauss upwind;
    //div(phi,R)      bounded Gauss upwind;
    //div(R)          Gauss linear;
    //div(phi,nuTilda) bounded Gauss upwind;
    //div((nuEff*dev(T(grad(U))))) Gauss linear;

// added for Of ext3.1
    div(phi,U)      Gauss upwind;
    div(phi,k)      Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(phi,R)      Gauss upwind;
    div(R)          Gauss linear;
    div(phi,nuTilda) Gauss upwind;

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

laplacianSchemes
{
// added for Of ext3.1
    //default         Gauss linear corrected;
    default         none;
    laplacian(nuEff,U) Gauss linear corrected;
    laplacian((1|A(U)),p) Gauss linear corrected;
    laplacian(DkEff,k) Gauss linear corrected;
    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
    laplacian(DREff,R) Gauss linear corrected;
    laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
// added for Of ext3.1
    interpolate(U)  linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}


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

[Tushar@cfd]

Originally Posted by deepinheart

I run a case using simpleFoam in OF 230 and everything goes well. When I turn to OF ext 3.1, I can't find 'bouned Gauss upwind' scheme as before. Instead, I use 'Gauss upwind' but it failed.

Anybody has similar experience? Would you please explain which one is the similar to 'bounded Gauss upwind'?

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

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
    // added for Of ext3.1
    grad(p)         Gauss linear;
    grad(U)         Gauss linear;
}

divSchemes
{
    default         none;
    //div(phi,U)      bounded Gauss upwind;
    //div(phi,k)      bounded Gauss upwind;
    //dv(phi,epsilon) bounded Gauss upwind;
    //div(phi,R)      bounded Gauss upwind;
    //div(R)          Gauss linear;
    //div(phi,nuTilda) bounded Gauss upwind;
    //div((nuEff*dev(T(grad(U))))) Gauss linear;

// added for Of ext3.1
    div(phi,U)      Gauss upwind;
    div(phi,k)      Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(phi,R)      Gauss upwind;
    div(R)          Gauss linear;
    div(phi,nuTilda) Gauss upwind;

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

laplacianSchemes
{
// added for Of ext3.1
    //default         Gauss linear corrected;
    default         none;
    laplacian(nuEff,U) Gauss linear corrected;
    laplacian((1|A(U)),p) Gauss linear corrected;
    laplacian(DkEff,k) Gauss linear corrected;
    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
    laplacian(DREff,R) Gauss linear corrected;
    laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
// added for Of ext3.1
    interpolate(U)  linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}


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

Dear Deepinheart,

I would you suggest you to have a look at the following post:

http://www.cfd-online.com/Forums/ope...end-3-1-a.html

As suggested by H. Jasak use:

Gauss upwind Gauss linear

This will work

-
Best Luck!