CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Bugs

Problem in fvschemes divSchemes cannot use Gauss linearUpwind

Register Blogs Community New Posts Updated Threads Search

Like Tree21Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 16, 2011, 04:02
Default
  #21
New Member
 
Felipe
Join Date: Apr 2010
Posts: 13
Rep Power: 16
fippo_dk is on a distinguished road
how do your results compare to measurements? are you under or over predicting?

In my case I over predict from 10 to 20 degrees
fippo_dk is offline   Reply With Quote

Old   August 2, 2012, 15:34
Default
  #22
uli
New Member
 
Join Date: Jun 2012
Posts: 25
Rep Power: 14
uli is on a distinguished road
hi all

I noticed that both in both tutorial cases and posted fvschemes there is always

Code:
div(R) Gauss linear;
Is that because there is no flux (phi = rho U) in that term and should it be left like this?

I tried

Code:
div(R)          Gauss limitedLinear 1;
div(R)          Gauss QUICK;
div(R)          Gauss MUSCL;
and it gave me similar results without any problems/errors.

thanks
Uli
uli is offline   Reply With Quote

Old   August 2, 2012, 15:40
Default
  #23
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 252
Blog Entries: 1
Rep Power: 19
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Hi! It is for other turbulence model where you use Reynolds stress
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Metallurgical Applications of Magnetohydrodynamics"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   August 2, 2012, 16:12
Default
  #24
uli
New Member
 
Join Date: Jun 2012
Posts: 25
Rep Power: 14
uli is on a distinguished road
Thanks for your reply.

So using a k-Epsilon model does not involve R? But how is it possible that I get different results after switching the scheme for div(R)?
uli is offline   Reply With Quote

Old   August 2, 2012, 16:18
Default
  #25
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 252
Blog Entries: 1
Rep Power: 19
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
You wrote you get the SAME results. You can easily check which schemes you solver uses, just set "default none;", and you will get errors until you define all schemes.
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Metallurgical Applications of Magnetohydrodynamics"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   August 21, 2012, 19:34
Default
  #26
uli
New Member
 
Join Date: Jun 2012
Posts: 25
Rep Power: 14
uli is on a distinguished road
hi, sorry for the late reply.

I wrote "similar" results since the plots are slightly different. I removed the schemes and it still worked, thanks for your help.
uli is offline   Reply With Quote

Old   April 9, 2014, 02:54
Default
  #27
Member
 
adarsh tiwari's Avatar
 
adarsh tiwari
Join Date: Feb 2014
Location: Bangalore
Posts: 42
Blog Entries: 5
Rep Power: 12
adarsh tiwari is on a distinguished road
Hi all,

I have to simulate the case of time varying BC 'uniformValue' for p and T

I tried various ways listed here and also from other open-sources, for schemes and solvers, to solve the 'p' with rhoSimpleFoam.

whatever may be the conditions applied, the solver is unable to run after third time-step. here I am posting my schemes and solvers. I am thinking that the problem is with these files only because I have cross-checked the p, T and U files and am also able to run 'rhoPimpleFoam' successfully.

Please let me know what went wrong.

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.2                                 |
|   \\  /    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         cellLimited leastSquares 1;//Gauss linear;
//    grad(U)         Gauss linear;
    grad(P)         Gauss linear;
}

divSchemes
{
    default         none;//Gauss linear corrected; //
    div(phi,U)      bounded Gauss upwind;
    div((muEff*dev2(T(grad(U))))) Gauss linear;
    div(phi,e)      bounded Gauss upwind;
    div(phi,h)      bounded Gauss upwind;
    div(phi,epsilon) bounded Gauss upwind;
    div(phi,omega)   bounded Gauss upwind;
    div(phi,k)      bounded Gauss upwind;
    div(phi,K)      bounded Gauss upwind;
    div(phi,Ekp)    bounded Gauss upwind;
    div(phid,p)     Gauss upwind;//bounded Gauss upwind;
    div(U,p)      bounded Gauss upwind;
}

laplacianSchemes
{
    default         Gauss linear orthogonal;//Gauss linear corrected;
//    laplacian(muEff,U) Gauss linear corrected;
//    laplacian(alphaEff,e) Gauss linear corrected;
//    laplacian(alphaEff,h) Gauss linear corrected;
//    laplacian((rho*(1|A(U))),p) Gauss linear corrected;
//    laplacian((rho|A(U)),p) Gauss linear corrected;
//    laplacian((rho*rAU),p) Gauss linear corrected;
//    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
//   laplacian(DomegaEff,omega) Gauss linear corrected;
//    laplacian(DkEff,k) Gauss linear corrected;
//   laplacian(1,p) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
    div(U,p)        upwind phi;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}


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

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

solvers
{
    p
    {
        solver         PCG;//GAMG;//
        preconditioner  none;//DIC;//FDIC;//GAMG;//diagonal;//
        agglomerator    faceAreaPair;
       tolerance       1e-07;
        relTol          0;//0.05;
        smoother        GaussSeidel;
        cacheAgglomeration true;//off;
        nCellsInCoarsestLevel 10;//20;
        mergeLevels     1;
    } //temperoraily not using 


    "(U|e|h|R|k|epsilon|omega)" //"(U|p|e|h|R|k|epsilon|omega)"
    {
        solver          smoothSolver;//solver          diagonal;//
        smoother        GaussSeidel;
        nSweeps         2;
        tolerance       1e-06;
//        relTol          0; //0.1;
//        solver          PBiCG;
//        preconditioner  DILU;
//        tolerance       1e-06;
//        relTol          0; //0.1;
    }

//    "(k|epsilon|omega)"
//    {
//        $U;
//        tolerance       1e-05;
//        relTol          0; //0.1;
//    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    rhoMin          rhoMin [ 1 -3 0 0 0 ] 0.4;
    rhoMax          rhoMax [ 1 -3 0 0 0 ] 1.5;

   residualControl
    {
        p               1e-2;
        U               1e-4;
        e               1e-3;

        // possibly check turbulence fields
        "(k|epsilon|omega)" 1e-3;
    }                       // residual control temporarily not used
}

PIMPLE
{
    momentumPredictor yes;
    nOuterCorrectors  1;
    nCorrectors       1;
    nNonOrthogonalCorrectors 0;
    rhoMin            rhoMin [ 1 -3 0 0 0 ] 0.5;
    rhoMax            rhoMax [ 1 -3 0 0 0 ] 2.0;

    maxCo             0.2;
    rDeltaTSmoothingCoeff 0.1;
    rDeltaTDampingCoeff 1;
    maxDeltaT         1;
}

relaxationFactors
{
    fields
    {
        p               0.03;
        rho             0.05;
    }
    equations
    {
        U               0.07;
        "(k|epsilon|omega)"   0.04;
        e               0.05;
        h               0.01;
    }
}


// ************************************************************************* //
adarsh tiwari is offline   Reply With Quote

Old   June 3, 2014, 06:29
Default
  #28
Senior Member
 
T. Chourushi
Join Date: Jul 2009
Posts: 321
Blog Entries: 1
Rep Power: 18
Tushar@cfd is on a distinguished road
Quote:
Originally Posted by adarsh tiwari View Post
Hi all,
relaxationFactors
{
fields
{
p 0.03;
rho 0.05;
}
equations
{
U 0.07;
"(k|epsilon|omega)" 0.04;
e 0.05;
h 0.01;
}
}

// ************************************************** *********************** //[/CODE]
Is it a serious issue (bug)?

A user can edit those files.

Use the following conditions for relaxation factors:

relaxationFactors
{
fields
{
p 0.3;
rho 1;
}
equations
{
U 0.7;
"(k|epsilon|omega)" 0.7;
e 0.7;
h 0.7;
}
}

User has an option to edit the "fvSchemes" as well depending on problem statement.

Best Luck!
Tushar@cfd is offline   Reply With Quote

Old   August 9, 2018, 07:52
Default
  #29
New Member
 
Sunil
Join Date: Aug 2016
Posts: 7
Rep Power: 10
sunilkonatham is on a distinguished road
Quote:
Originally Posted by makaveli_lcf View Post
Hi!

Code:
gradSchemes
{
    default         cellLimited leastSquares 1;
}

divSchemes
{
    default        none;
    div(phi,epsilon) Gauss linearUpwind grad(epsilon);  // OF20x syntax
    div(phi,k)       Gauss linearUpwind grad(k);  // OF20x syntax
    div(phi,U)       Gauss linearUpwind grad(U);  // OF20x syntax
    div((nuEff*dev(T(grad(U)))))     Gauss linear;  // OF20x syntax
}
My example which I use with pimpleFoam.
It has been a while since this thread is posted, can anyone please answer my question?

In this case, does it mean that divergence scheme Gauss linearUpwind uses the gradient scheme of respective variables?


Thanks,
Sunil.
sunilkonatham is offline   Reply With Quote

Old   August 9, 2018, 08:23
Default
  #30
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 252
Blog Entries: 1
Rep Power: 19
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Quote:
Originally Posted by sunilkonatham View Post
It has been a while since this thread is posted, can anyone please answer my question?

In this case, does it mean that divergence scheme Gauss linearUpwind uses the gradient scheme of respective variables?


Thanks,
Sunil.

Yes, you are right. It can either use specific scheme listed in the gradient section or just use default scheme. You can also specify your name for the scheme to use it with several variables, for example like this:


Code:
gradSchemes
{
   default      none;
   grad(DDD)    cellLimited Gauss linear 1;
   grad(p)      leastSquares;
}

divSchemes
{
    default        none;
    div(phi,epsilon) Gauss linearUpwind grad(DDD);
    div(phi,k)       Gauss linearUpwind grad(DDD);
    div(phi,U)       Gauss linearUpwind grad(DDD);
}
lukasf likes this.
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Metallurgical Applications of Magnetohydrodynamics"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   August 9, 2018, 08:39
Default
  #31
New Member
 
Sunil
Join Date: Aug 2016
Posts: 7
Rep Power: 10
sunilkonatham is on a distinguished road
Great. Thanks a lot Dr. Vakhrushev.
sunilkonatham is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
DivSchemes Gauss limitedLinearV and interfaceCompression nicasch OpenFOAM Running, Solving & CFD 1 July 12, 2010 11:26
FvSchemes sega OpenFOAM Running, Solving & CFD 2 February 15, 2010 12:07
DivSchemes limitedLinearV and interfaceCompression nicasch OpenFOAM 0 February 28, 2008 11:33
Divschemes limitedLinearV 1 and interfaceCompression nicasch OpenFOAM 0 February 27, 2008 13:24
FvSchemes from 12 to 11 fedegavo OpenFOAM Running, Solving & CFD 1 January 20, 2006 14:49


All times are GMT -4. The time now is 06:55.