CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

KEpsilon model incorrect result

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 26, 2022, 05:37
Default KEpsilon model incorrect result
  #1
New Member
 
Join Date: Apr 2022
Posts: 9
Rep Power: 3
vronti is on a distinguished road
Hello
I am trying to simulate a realizableKE model.
Geometry: 0.006x0.009mm rectangle.
BC: top bottom - walls; sides - symmetry; left - inlet; right - outlet

I calculated the velocity for the Reynolds numbers 500, 2000, 5000, and 9000 to be [0.04959, 0.19837, 0.4959, 0.8926] m/s with the hydraulic diameter being the distance between the top and bottom walls. (I consider this case as a parallel plate as I am using symmetry on either side)

The simulation ran perfectly fine for RE500 and RE2000 but the 5000 Reynolds number is causing problems. I will attach the case setup for reference.

Though, I defined no-slip condition for top bottom walls velocity is unevenly spread on the walls which I have no clue why.
Screenshot 2022-08-23 115955.jpg

Simulation converged to 1e-4
Code:
Time = 3000

smoothSolver:  Solving for Ux, Initial residual = 0.000166362, Final residual = 1.08171e-06, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 9.70676e-07, Final residual = 9.70676e-07, No Iterations 0
smoothSolver:  Solving for Uz, Initial residual = 1.86604e-05, Final residual = 1.41014e-07, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.000271834, Final residual = 1.46852e-06, No Iterations 2
GAMG:  Solving for p, Initial residual = 1.46851e-06, Final residual = 4.37859e-07, No Iterations 1
GAMG:  Solving for p, Initial residual = 4.3786e-07, Final residual = 4.3786e-07, No Iterations 0
GAMG:  Solving for p, Initial residual = 4.3786e-07, Final residual = 4.3786e-07, No Iterations 0
time step continuity errors : sum local = 2.11366e-07, global = 6.86507e-09, cumulative = 0.000122676
smoothSolver:  Solving for epsilon, Initial residual = 3.00512e-06, Final residual = 4.96645e-08, No Iterations 2
bounding epsilon, min: -0.00358805 max: 3.25773 average: 0.12346
smoothSolver:  Solving for k, Initial residual = 5.60576e-07, Final residual = 5.60576e-07, No Iterations 0
ExecutionTime = 9182.55 s  ClockTime = 9221 s
If anyone has any idea why this is happening please help

Best regards,

Case: flowDevelop_5000.zip
vronti is offline   Reply With Quote

Reply

Tags
kepsilon, openfoam, turbulance

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
Use of k-epsilon and k-omega Models Jade M Main CFD Forum 40 January 27, 2023 07:18
Table bounds warnings at: END OF TIME STEP CFXer CFX 4 July 16, 2020 23:44
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 06:20
The density cross the interface is incorrect in the KEpsilon model in interFoam? duan OpenFOAM Programming & Development 0 August 14, 2015 11:40
problem with wallFunctions (kEpsilon model) dinonettis OpenFOAM 12 April 22, 2010 21:37


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