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

Problem in fvschemes divSchemes cannot use Gauss linearUpwind

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

Like Tree13Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   November 16, 2011, 04:02
Default
  #21
New Member
 
Felipe
Join Date: Apr 2010
Posts: 13
Rep Power: 7
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, 14:34
Default
  #22
uli
New Member
 
Join Date: Jun 2012
Posts: 25
Rep Power: 5
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, 14:40
Default
  #23
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
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 "Advanced Process Simulation of
Solidification and Melting"

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

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   August 2, 2012, 15:12
Default
  #24
uli
New Member
 
Join Date: Jun 2012
Posts: 25
Rep Power: 5
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, 15:18
Default
  #25
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
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 "Advanced Process Simulation of
Solidification and Melting"

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

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   August 21, 2012, 18:34
Default
  #26
uli
New Member
 
Join Date: Jun 2012
Posts: 25
Rep Power: 5
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, 01:54
Default
  #27
Member
 
adarsh tiwari
Join Date: Feb 2014
Location: Bangalore
Posts: 48
Rep Power: 3
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, 05:29
Default
  #28
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 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

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
DivSchemes Gauss limitedLinearV and interfaceCompression nicasch OpenFOAM Running, Solving & CFD 1 July 12, 2010 10: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 08:51.