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

Convergence problem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 11, 2013, 10:53
Default Convergence problem
  #1
Member
 
xxxxx
Join Date: Feb 2013
Posts: 34
Rep Power: 4
xxxx is on a distinguished road
Hello everybody,
I'm trying to run a case with simpleFoam. It's a stationary simulation with k-espilon turbulence model.
Initially it seems the simulation is converging but after some step suddenly the residuals show a strange behaviour and from this on the results have no phisical sense.
I attach the residuals vs iteration plot.
Is there anyone who can help me?
I tryied to make modification to the mesh, to the relaxation factor and to the fvScheme but I wasn't able to solve the problem.
Attached Images
File Type: jpg residuals.jpg (34.6 KB, 26 views)
xxxx is offline   Reply With Quote

Old   February 11, 2013, 12:52
Default
  #2
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
i have a problem like you.whats your case?do you have some region of high gradients?
immortality is offline   Reply With Quote

Old   February 11, 2013, 13:54
Default
  #3
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 8
danvica is on a distinguished road
Can you send your controlDict / fvSolution and fvScheme ?
__________________
Daniele Vicario

blueCFD2.1 - Windows 7
danvica is offline   Reply With Quote

Old   February 11, 2013, 15:47
Default
  #4
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
I use compressible solvers like rhoPimpleFoam but now changing to sonicFoam.
immortality is offline   Reply With Quote

Old   February 11, 2013, 16:00
Default
  #5
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
Code:
application     rhoPimpleFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         .02;//0.020708089

deltaT          5.6e-9;

writeControl    adjustableRunTime;

writeInterval   .000001;

purgeWrite      0;

writeFormat     binary;

writePrecision  10;

writeCompression on;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  yes;

maxCo           0.1;

maxDeltaT       .0001;

libs ( "libOpenFOAM.so" "libgroovyBC.so" ) ;
Code:
ddtSchemes
{
    default         backward;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
    grad(U)         Gauss linear;//cellLimited leastSquares 1.0
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss upwind;
    div(phid,p)     Gauss upwind;// upwind
    div(phi,K)      Gauss upwind;
    div(phi,h)      Gauss upwind;
    div(phi,k)      Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(phi,R)      Gauss upwind;
    div(phi,omega)  Gauss upwind;
    div((rho*R))    Gauss limitedLinear 1;
    div(R)          Gauss limitedLinear 1;
    div(U)          Gauss limitedLinear 1;
    div((muEff*dev2(T(grad(U))))) Gauss linear;
    div(tauMC)  Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited .5;//limited .5
    /*laplacian(muEff,U) Gauss linear corrected;
    laplacian(mut,U) Gauss linear corrected;
    laplacian(DkEff,k) Gauss linear corrected;
    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
    laplacian(DREff,R) Gauss linear corrected;
    laplacian(DomegaEff,omega) Gauss linear corrected;
    laplacian((rho*(1|A(U))),p) Gauss linear corrected;
    laplacian(alphaEff,h) Gauss linear corrected;
    laplacian(k,T)   Gauss linear corrected;  
    laplacian(alpha,e) Gauss linear corrected;
    laplacian(alphaEff,e)  Gauss linear corrected;*/
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default       limited .5;//corrected
}

fluxRequired
{
    default         no;
    p               ;
}
Code:
solvers
{
    p
    {
       solver PCG;
       preconditioner DIC;
       tolerance  1e-13;
       relTol 0;
       maxIter 25000;

     /*solver          GAMG;
        smoother  DICGaussSeidel;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 200;
        nPreSweeps      0;
        nPostSweeps     2;
        nFinestSweeps   2;
        agglomerator  faceAreaPair;
        mergeLevels 2;
        tolerance       1e-012;
        relTol          0;*/
    }

    pFinal
    {
        $p;
        relTol          0;
    }

    "rho.*"
    {
        $p;
        tolerance       1e-11;
        relTol          0;
    }

    "(U|e|h|R|k|epsilon|omega)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-11;
        relTol          0;
        maxIter 25000;
    }

    "(U|h|R|k|epsilon|omega)Final"
    {
        $U;
        relTol          0;
    }
}


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

residualControl
    {
        "(U|k|omega)"
        {
            relTol          0;
            tolerance       0.00000000001;
        }
    }
}
immortality is offline   Reply With Quote

Old   February 11, 2013, 17:22
Default
  #6
Member
 
xxxxx
Join Date: Feb 2013
Posts: 34
Rep Power: 4
xxxx is on a distinguished road
Thanks for your answers.
My case is steasy state incompressible flow in a duct of quite complex geometry. I tried to simplify the geometry by removing chamfers and small details. The mesh is unstructured, with boundary layer, about 500'000 nodes.
I noticed that, since some iterations before the computation failed, the mass flow continuity started not to be satisfied.
Then in some region of the mesh the values of velocity or pressure rise suddenly. When I change the mesh and restart another simulation the point in which the anomalous values appears mooves.
When i do the checkMesh i don't find any error or warning. The skew value is around 0.7 . I tried to change the fvScheme. I set 1st order scheme everywere. This allowed me to make the simulation run for 1 hundred iteration more before to have the same problem.
Does anyone have an idea of the causes?
xxxx is offline   Reply With Quote

Old   February 13, 2013, 16:36
Default
  #7
Member
 
xxxxx
Join Date: Feb 2013
Posts: 34
Rep Power: 4
xxxx is on a distinguished road
Please, can someone answer to some dubts?
1- when computing a simulation, should one let it run without modify relaxation factors or is it possible/advised to modify them? It happens that changing a relaxation factor completely changes the convergence. But is it right to do so?
2- how can i understand when the simulation is converged? Sometimes it happens that although the residuals are low, the results have no physical sense. On the other hand, sometimes although the residuals are high macroscopic quantities ,pressure drop for example, do not change much through time steps and results seem to be quite good. Is there any way to be quite sure that simulation has reached convergence?
xxxx is offline   Reply With Quote

Old   February 15, 2013, 07:16
Default
  #8
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 8
danvica is on a distinguished road
Without posting the necessary info (How to give enough info to get help) it's difficult to help you.

I had a similar problem regarding mass flow rate in a tet mesh... just a blind hint: in your fvsolution use relTol=0 for p.
__________________
Daniele Vicario

blueCFD2.1 - Windows 7
danvica 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
convergence problem when use pisoFoam, LES for wind tunnel case Forrest_Lei OpenFOAM 3 July 19, 2011 06:00
convergence problem commonyue Main CFD Forum 1 December 1, 2009 04:54
Convergence of CFX field in FSI analysis nasdak CFX 2 June 29, 2009 01:17
3D Fluid Flow Convergence problem Emily FLUENT 2 March 21, 2007 23:18
Non Convergence of 3D Heat transfer cfd problem Balraj Main CFD Forum 3 December 9, 2004 01:24


All times are GMT -4. The time now is 20:13.