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

delta T decreases suddenly in pimplefoam

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

Like Tree2Likes
  • 1 Post By Tushar@cfd
  • 1 Post By Tushar@cfd

Reply
 
LinkBack Thread Tools Display Modes
Old   June 26, 2014, 11:33
Default delta T decreases suddenly in pimplefoam
  #1
Member
 
Ali reza
Join Date: Mar 2014
Posts: 95
Rep Power: 3
1988 is on a distinguished road
hi
I have simulated tiny pipe with newtoniam incompressible flow which is laminar by pimplefoam but the time steps are getting smaller and smaller during the solution process(it start from 0.003 and after 10 time steps gets to 10e-23 and more)and finally it is divrged .courant was set as 1 and nothing has changed by increasing the amounts of nOuterCorrectors ,nCorrectors and nNonOrthogonalCorrector.what should I do ?any suggestion?
Code:
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           9078
    faces:            90732
    internal faces:   83908
    cells:            43660
    faces per cell:   4
    boundary patches: 3
    point zones:      0
    face zones:       1
    cell zones:       1

Overall number of cells of each type:
    hexahedra:     0
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    43660
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology                  
    body                6496     3278     ok (non-closed singly connected)  
    outlet              164      98       ok (non-closed singly connected)  
    inlet               164      98       ok (non-closed singly connected)  

Checking geometry...
    Overall domain bounding box (-0.0215399 -0.0115756 -0.002365) (0.0198627 0.00851412 0.002365)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (3.2476e-18 2.81028e-17 -1.76588e-17) OK.
    Max cell openness = 1.42371e-16 OK.
    Max aspect ratio = 6.62732 OK.
    Minimum face area = 4.49238e-08. Maximum face area = 2.90028e-07.  Face area magnitudes OK.
    Min volume = 3.78117e-12. Max volume = 4.92374e-11.  Total volume = 8.03586e-07.  Cell volumes OK.
    Mesh non-orthogonality Max: 56.6451 average: 19.2493
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 0.702453 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End
thanks
1988 is offline   Reply With Quote

Old   June 27, 2014, 06:56
Default
  #2
Member
 
Ali reza
Join Date: Mar 2014
Posts: 95
Rep Power: 3
1988 is on a distinguished road
there is no idea?
1988 is offline   Reply With Quote

Old   June 27, 2014, 07:06
Default
  #3
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
By uploading more information about your case, you can get quick replies.

-
Best Regards!
Tushar@cfd is offline   Reply With Quote

Old   June 27, 2014, 07:10
Default
  #4
Member
 
Ali reza
Join Date: Mar 2014
Posts: 95
Rep Power: 3
1988 is on a distinguished road
thanks for your answer
these are fv solution and scheme
Code:
solvers
{
     p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          1e-05;
    }

    pFinal
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-08;
        relTol          0;
    }

    "(U|k|epsilon)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-05;
        relTol          0.1;
    }

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

PIMPLE
{
    nOuterCorrectors 4;
    nCorrectors     2;
    nNonOrthogonalCorrectors 3;
    pRefCell        0;
    pRefValue       0;
residualControl
    {
        U
        {
            tolerance 1e-2;
            relTol 0;
        }

        p
        {
            tolerance 1e-2;
            relTol 0;
        }
    }
}

relaxationFactors
{
    fields
    {
    }
    equations
    {
        "p.*"           0.3;
        "U.*"           0.7;
        "k.*"           0.4;
        "epsilon.*"     0.5;
    }
}
Code:
ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
    grad(U)         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss limitedLinearV 1;
    div(phi,k)      Gauss limitedLinear 1;
    div(phi,epsilon) Gauss limitedLinear 1;
    div(phi,R)      Gauss limitedLinear 1;
    div(R)          Gauss linear;
    div(phi,nuTilda) Gauss limitedLinear 1;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    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;
    interpolate(U)  linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}
1988 is offline   Reply With Quote

Old   June 27, 2014, 07:23
Default
  #5
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
You're welcome!

In fvscheme file try these instead of mentioned ones. Best is to try with 1st order and check your results.

divSchemes
{
default none;
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(T(grad(U))))) Gauss linear;
}

Let us know about your results. post your reply here

-
Best Regards!
Tushar@cfd is offline   Reply With Quote

Old   June 27, 2014, 07:44
Default
  #6
Member
 
Ali reza
Join Date: Mar 2014
Posts: 95
Rep Power: 3
1988 is on a distinguished road
that works and everyting is ok
thank you
but may I ask you to explain how did you find that I mean what the big differnces are there between these two div schemes?
thanks
1988 is offline   Reply With Quote

Old   June 27, 2014, 07:58
Default
  #7
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
Good!

Actually, there are various ways to avoid such situations. One such problem could be due to the above situation.

You can find the reference for these in any of the CFD books.

These "div schemes" differ in order of accuracy and many other facts.

Quick Hint: Upwind scheme is unconditionally stable of all.
1988 likes this.
Tushar@cfd is offline   Reply With Quote

Old   June 27, 2014, 15:47
Default
  #8
Member
 
Ali reza
Join Date: Mar 2014
Posts: 95
Rep Power: 3
1988 is on a distinguished road
thanks
As you know yhis is a first order scheme, thus too diffusive and I will not get a good solution,am I right?and what is the solution
1988 is offline   Reply With Quote

Old   June 28, 2014, 00:47
Default
  #9
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
Please, excuse for late reply. I was away from the desk.

Yes you're correct. With increased mesh resolution you can obtain very close solutions.

I don't know the details of your problem case. It's always wise to compare the CFD results with the available resources (experimental or analytical).

Anyways, with experience and literature survey you can gain grip over many of these CFD techniques.

Best of Luck for your future endeavour.

-
Best Regards!
Tushar@cfd is offline   Reply With Quote

Old   June 28, 2014, 01:14
Default
  #10
Member
 
Ali reza
Join Date: Mar 2014
Posts: 95
Rep Power: 3
1988 is on a distinguished road
thanks for protect and excuse me for lots of question (I am so curious ) but I promise that this is the last one
I think one way to decrease inaccuracy in using upwind linear in starting the algorithm with upwind and then changing it to more accurate schemes but I dont know when is the right time to do that,I mean is there any clue in the terminal that can aware us?
thanks alot
1988 is offline   Reply With Quote

Old   June 28, 2014, 01:34
Default
  #11
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
Actually, you can get comparable solution with upwind scheme too depending upon your problem case. Sometime, results obtained are comparable with the available resources. You can find many literature in which 1st order scheme give very near comparable solutions.

This is because CFD solution doesn't depend on these factors only. There are many other factors which are equally important for CFD solutions.

It's always best to monitor residuals. Generally, convergence check tell us when to stop our iterations.

I hope I tried to answer your queries.

-
Best Regards!
1988 likes this.
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
pimpleFoam: turbulence->correct(); is not executed when using residualControl hfs OpenFOAM Running, Solving & CFD 3 October 29, 2013 09:35
Delta wing cfd H_R FLUENT 16 September 19, 2013 04:41
Smoothing of IDDES delta Fuchs OpenFOAM Running, Solving & CFD 0 July 15, 2013 11:53
Differences simpleFoam vs. pimpleFoam / RASModel.H vs turbulenceModel.H uli OpenFOAM Programming & Development 7 January 26, 2013 16:01
Delta Wing Structured Grid Riaan FLUENT 3 December 31, 2004 13:03


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