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

bounding k and bounding epsilon over the limit

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

LinkBack Thread Tools Search this Thread Display Modes
Old   October 7, 2023, 12:37
Question bounding k and bounding epsilon over the limit
New Member
Join Date: Apr 2021
Posts: 3
Rep Power: 4
YinJian is on a distinguished road

I am using multiphaseInterFoam to simulate one case.The model consists of a horizontally positioned cylindrical pipe with a diameter of 115mm and a length of 5m. On the left side of the pipe is a velocity inlet with a velocity of 12 m/s, and on the right side is a wall boundary. Along the side of the pipe, there are multiple vertically oriented cylindrical small pipes with a diameter of 12mm and a length of 500mm each. The pressure outlet is located at the end of these small pipes, with a pressure of 0. For this case i am using k-epsilon turbulence model.

I have created unstructured mesh and quality of mesh seems to be right. CheckMesh report is as given below:
Checking geometry...
Overall domain bounding box (0 -0.556 -0.48451) (5 0.556 0.48451)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (2.06719e-17 4.86465e-18 1.85385e-17) OK.
Max cell openness = 4.18769e-16 OK.
Max aspect ratio = 18.761 OK.
Minimum face area = 1.70871e-12. Maximum face area = 5.02513e-05. Face area magnitudes OK.
Min volume = 7.33752e-10. Max volume = 2.49585e-07. Total volume = 0.0532077. Cell volumes OK.
Mesh non-orthogonality Max: 65.9097 average: 8.02965
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 2.57036 OK.
Coupled point location match (average 0) OK.
Mesh OK.

When I run the simulation, the values of k and epsilon suddenly become negative, and the positive values become extremely large, eventually leading to divergence.

What should I do to ensure that the simulation runs smoothly? Thanks for the help in advance.
Attached Images
File Type: png model1.png (48.0 KB, 11 views)
File Type: png mesh.png (109.2 KB, 13 views)
File Type: jpg meshcheck.jpg (69.4 KB, 10 views)
Attached Files
File Type: zip (9.8 KB, 1 views)
YinJian is offline   Reply With Quote

Old   October 8, 2023, 09:34
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,902
Rep Power: 32
hjasak will become famous soon enough
Try using upwind on convection and limited laplacian on diffusion. If this does not help, Your outlet boundary conditions are wrong.
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting:
hjasak is offline   Reply With Quote


bounding epsilon, bounding k

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
Converging Diverging Nozzle with dbnsTurbFoam Saleh Abuhanieh OpenFOAM Running, Solving & CFD 4 December 13, 2019 11:26
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 16 March 4, 2017 09:30
Bounding k and epsilon & coldEngineFoam sasanghomi OpenFOAM 1 September 13, 2013 13:12
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 1 November 25, 2008 21:21

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