CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Verification & Validation

Jump in epsilon values near the wall : low re k-epsilon model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 25, 2013, 07:09
Default Jump in epsilon values near the wall : low re k-epsilon model
  #1
Member
 
Malik
Join Date: Dec 2012
Location: Austin, USA
Posts: 52
Rep Power: 4
malaboss is on a distinguished road
Hi FOAMers !
I am now studying an airfoil profile which small chord makes it look like a cylinder.
The thickness of the airfoil is about several cm and the velocity at the inlet is 10m/s.
I considered my fluid as incompressible due to low velocities.
k at the inlet is 1.5 ( 10% of turbulence)
epsilon at the inlet is (19,2, I choosed the highest mixing length that I considered to be about half of the airfoil thickness)

First, I tried to use a (high reynolds) k epsilon turbulence model with a wall function. It gave me some weird as epsilon near the wall was really high : several tens of thousands.
I saw in the following tutorial /tutorials/incompressible/simpleFoam/pitzDaily similar values.
Is it OK to have such great values near the wall for epsilon ?


Then, I decided to use a low Re kepsilon model launderSharma, and by initializing k an epsilon near to zero close to the wall (as it is done in the incompressible/boundaryFoam/boundaryLaunderSharma tutorial) I find a high jump for my espilon values as you can see in the pictures :
the first one is the espilonfield and the second one is the plotting of epsilon values in the radial direction.
I already checked y+ values and i the average is lower than 1.

I guess I should not have that result. Do you know why i get this and how to fix it ?
If you need it, I can provide more information about my dictionnaries.

Thanks to all for your help
Attached Images
File Type: jpg epsilonFieldWall_19.2.jpg (29.4 KB, 62 views)
File Type: jpg epsilonJumpWall_19.2.jpg (28.5 KB, 50 views)
malaboss is offline   Reply With Quote

Old   February 1, 2013, 17:36
Default
  #2
Member
 
Join Date: Jun 2012
Posts: 55
Rep Power: 5
maHein is on a distinguished road
High values for dissipation near the wall are normal. Just have a look at the validation simulation of a flow ofer a flat plate performed by NASA http://www.grc.nasa.gov/WWW/wind/val.../fpturb01.html. You could just convert your data to dimensionless epsilon+ and compare it to the validation results.
maHein is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
Natural convection in a closed domain STILL NEEDING help! Yr0gErG FLUENT 3 June 12, 2013 02:12
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug unoder OpenFOAM Installation 11 January 30, 2008 21:30
Values of k & epsilon in the moving wall ghlee Main CFD Forum 1 August 16, 1999 12:42
Wall functions Abhijit Tilak Main CFD Forum 6 February 5, 1999 02:16


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