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

P Divergence for k-Epsilon Model (Simple Foam)

Register Blogs Community New Posts Updated Threads Search

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

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 (https://imgur.com/a/4tp2tA5) 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

Reply

Tags
cfd, k-epsilon, openfoam, rans


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 11:59.