CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   bounded Gauss upwind Scheme (https://www.cfd-online.com/Forums/openfoam-solving/148973-bounded-gauss-upwind-scheme.html)

deepinheart February 23, 2015 05:35

bounded Gauss upwind Scheme
 
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'?


Code:

/*--------------------------------*- 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 February 23, 2015 05:57

Quote:

Originally Posted by deepinheart (Post 532993)
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'?


Code:

/*--------------------------------*- 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:

Code:

Gauss upwind Gauss linear
This will work

-
Best Luck!


All times are GMT -4. The time now is 04:54.