CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

P Divergence for k-Epsilon Model (Simple Foam)

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

LinkBack Thread Tools Search this Thread Display Modes
Old   January 6, 2021, 21:05
Unhappy P Divergence for k-Epsilon Model (Simple Foam)
New Member
almondymarvin's Avatar
Join Date: Jan 2021
Location: Calgary, AB
Posts: 1
Rep Power: 0
almondymarvin is on a distinguished road

I am trying to simulate water flow inside an emitter (drip irrigation) using RANS and LES (Samagorin–Lilly), I have constructed an structured quad mesh in ANSA ( it is a very coarse mesh of 25600 volume elements.

I am just trying to perform RANS as first step using k-Epsilon model. I reused the pitzDaily simpleFoam tutorial to try and run it using my mesh, I changed everything that was needed (PIDs match to those ones in the tutorial). But for some reason my simulation explodes and diverges, specially p.

I left the initial conditions as they are, this is a preliminary simulation, will change them later.

Here is the last iteration:

Time = 5

smoothSolver: Solving for Ux, Initial residual = 0.579153, Final residual = 0.0324501, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.530541, Final residual = 0.0388022, No Iterations 2
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0258436, No Iterations 2
time step continuity errors : sum local = 9.49971e+40, global = -4.51565e+40, cumulative = -4.51565e+40
smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.0233028, No Iterations 2
bounding epsilon, min: -9.60072e+58 max: 2.86832e+60 average: 1.98086e+57
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0699013, No Iterations 1
bounding k, min: -5.66487e+50 max: 1.4514e+52 average: 3.68843e+48
ExecutionTime = 0.47 s ClockTime = 0 s

Time = 6

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.086456, No Iterations 8
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0864664, No Iterations 8
*Code Stops*

What could be the possible causes of my divergence?
I have a suspicion it could be my mesh but not 100% sure, I tried another mesh that was for another school assignment in the same file and it did converge! I am very confused at the moment

Any comment would be helpful, thanks!
almondymarvin is offline   Reply With Quote


cfd, k-epsilon, openfoam, rans

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
[CGNS] CGNS converters available mbeaudoin OpenFOAM Meshing & Mesh Conversion 137 December 14, 2018 05:20
error with reactingFoam BakedAlmonds OpenFOAM Running, Solving & CFD 4 June 22, 2016 03:21
simpleFoam parallel AndrewMortimer OpenFOAM Running, Solving & CFD 12 August 7, 2015 19:45
Overflow Error in Multiphase Modelling with Two Continuous Fluids ashtonJ CFX 6 August 11, 2014 15:32
[blockMesh] Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Meshing & Mesh Conversion 10 April 2, 2007 15:00

All times are GMT -4. The time now is 17:56.