CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

bounded Gauss upwind Scheme

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   February 23, 2015, 06:35
Default bounded Gauss upwind Scheme
  #1
New Member
 
Join Date: Apr 2009
Posts: 11
Rep Power: 8
deepinheart is on a distinguished road
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               ;
}


// ************************************************************************* //
deepinheart is offline   Reply With Quote

Old   February 23, 2015, 06:57
Default
  #2
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 318
Blog Entries: 1
Rep Power: 8
Tushar@cfd is on a distinguished road
Quote:
Originally Posted by deepinheart View Post
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:

fvSchemes in FOAM extend 3.1

As suggested by H. Jasak use:

Code:
Gauss upwind Gauss linear
This will work

-
Best Luck!
__________________
Mistakes Are Inevitable, You Can Only Grow Stronger If You Lean From Them
Tushar@cfd is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Free Surface Ship Flow timfranke OpenFOAM Running, Solving & CFD 302 April 21, 2015 08:34
Implementation of 2nd order upwind scheme jaason OpenFOAM Running, Solving & CFD 4 February 6, 2015 18:40
Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPO adarsh tiwari OpenFOAM Running, Solving & CFD 8 February 7, 2014 05:52
bounded Gauss and normal Gauss question Tobi OpenFOAM 2 June 18, 2013 00:35
solution diverges when linear upwind interpolation scheme is used subash OpenFOAM 0 May 29, 2010 01:23


All times are GMT -4. The time now is 10:42.