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

trying why my runs diverge - simple case I built

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

Like Tree1Likes
  • 1 Post By Tobermory

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 28, 2023, 15:55
Default trying why my runs diverge - simple case I built
  #1
New Member
 
Join Date: Oct 2023
Posts: 5
Rep Power: 2
chemotkbilo is on a distinguished road
So I have the following case that I built, the mesh lis ok (used checkMesh), and developed the the required constants. However, when I run simpleFoam or foamRun, it starts to solve up to certain point and then it starts to produce nan.
For example, see below

HTML Code:
[CODE]SmoothSolver:  Solving for Ux, Initial residual = 0.91362581, Final residual = 9.1691167e-06, No Iterations 36
smoothSolver:  Solving for Uy, Initial residual = 0.97344182, Final residual = 7.8785232e-06, No Iterations 34
smoothSolver:  Solving for Uz, Initial residual = 0.99999992, Final residual = 8.4295095e-06, No Iterations 37
GAMG:  Solving for p, Initial residual = 1, Final residual = nan, No Iterations 1000
time step continuity errors : sum local = nan, global = nan
GAMG:  Solving for p, Initial residual = nan, Final residual = nan, No Iterations 1000
time step continuity errors : sum local = nan, global = nan, cumulative = nan
smoothSolver:  Solving for epsilon, Initial residual = nan, Final residual = nan, No Iterations 1000
smoothSolver:  Solving for k, Initial residual = nan, Final residual = nan, No Iterations 1000
ExecutionTime = 15.673765 s  ClockTime = 15 s

Courant Number mean: nan max: nan
Time = 0.0469628s

smoothSolver:  Solving for Ux, Initial residual = nan, Final residual = nan, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = nan, Final residual = nan, No Iterations 1000
smoothSolver:  Solving for Uz, Initial residual = nan, Final residual = nan, No Iterations 1000
[/CODE]
My files are provided below - any help would be appreciated.
Attached Files
File Type: zip pitzDai2.zip (24.1 KB, 4 views)
chemotkbilo is offline   Reply With Quote

Old   October 29, 2023, 14:08
Default
  #2
Senior Member
 
Join Date: Apr 2020
Location: UK
Posts: 640
Rep Power: 13
Tobermory will become famous soon enough
Looks like you are trying to run as a transient simulation. Try changing over to steady state:

- in fvSchemes, change ddtSchemes default from Euler to steadyState
- in fvSolution, comment out the PIMPLE section and add the following:
Code:
SIMPLE
{
    nNonOrthogonalCorrectors 0;
    consistent      yes;

    residualControl
    {
        p               1e-2;
        U               1e-3;
        "(k|epsilon)" 1e-3;
    }
}

relaxationFactors
{
    fields
    {
        p               0.3;
    }
    equations
    {
        U               0.7;
        k               0.7;
        epsilon         0.7;
    }
}
- in controlDict, change deltaT to 1

Now run with simpleFoam, and this should converge/behave nicely. Good luck!
chemotkbilo likes this.
Tobermory is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
FSI. simple pipe case crashes Svensen ANSYS 0 November 19, 2018 13:59
Dambreak case can't be simulated by the simple coupled level set and vof (sclsvof) saatt OpenFOAM Running, Solving & CFD 1 January 21, 2018 18:40
Simple test case for k-epilon model mathletic Main CFD Forum 2 September 28, 2012 04:04
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 13:24
Need help on simple CFD case. (using CFD-ACE+) Sean Main CFD Forum 1 September 30, 2005 11:05


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