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

turbulence flat plane boundary

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 7, 2012, 04:34
Default turbulence flat plane boundary
  #1
New Member
 
Garrison, Liang
Join Date: Feb 2012
Posts: 15
Rep Power: 14
garrison is on a distinguished road
Hi Foamers!
Recently I am doing a very simple case using S-A-DDES.
The computational domain is 12m long, of which the left side is inlet, the right side is outlet. Inlet velocity is 1m/s, I choose water instead of air so the Reynold number is 1.2e7.
I intend to get a classic turbulent boundary layer velocity profile like this one


but I failed, I got one like this (pls see my attachment)
It seems I can capture the laminar zone of the BL(Yplus from 0-10), but fail to simulate the log zone, and also, the thickness of the BL should be like 0.17m (from the turbulent BL theory),


but the one I get is far less than 0.17m, it's about 0.01m. You can find in my attachment that the velocity in BL grows to the freestream velocity at about yPlus =80. but theoretically, yPlus=80 should be the log zone, right?

So, I don't know what's wrong with it? My mesh is like this: first point yPlus is below 1, with about 50 points in the BL (the thickness I get from the above formula, 0.17m).
delta t =1e-2, the CFL number is less then 1.

Seems like my boundary layer has not fully developed to a turbulent one. but I got my Reynold number base on the length of the plane (Rex) 1.2e7 which I think is large enough. And I run the case for about 500 seconds, with the inlet velocity 1m/s, the plane length 12m, it should flow through the domain for like 40 times.

Can anyone give me any hint?
Thanks in advance

Garry
Attached Images
File Type: jpg BL.jpg (16.4 KB, 49 views)
garrison is offline   Reply With Quote

Old   April 10, 2012, 03:01
Default
  #2
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 273
Rep Power: 18
anishtain4 is on a distinguished road
Dear Garrison

first of all you should scale your y axis too, it's u+ not u/uref, as you see your sublayer shape is not the same either,
second; tell us a little about your solver, your model and if any wall function you are using, and boundary conditions
PS: you mentioned fully developed, there isn't such a thing on external flows
anishtain4 is offline   Reply With Quote

Old   April 10, 2012, 04:09
Default
  #3
New Member
 
Garrison, Liang
Join Date: Feb 2012
Posts: 15
Rep Power: 14
garrison is on a distinguished road
Thank you for you reply Anishtain!

Actually I did rescale the y axis,
u+=ux/u*
u*=sqrt(tau/density)
tau=du/dy*mu
du/dy=( u[i+1] - u[i] ) / ( y[i+1] - y[i] )

Am I right?

But I didn't change the x scale to a log one. The new attachment is of a log x axis.

here's my BC:
inlet:
U = 1m/s
p = zeroGradient
nuTilda = 0

outlet:
U = zeroGradient
p = 0
nuTilda = zeroGradient

top:
U = zeroGradient
p = zeroGradient
nuTilda = zeroGradient

bottomwall:
U = 0
p = zeroGradient
nuTilda = 0

I am confused about the inlet and the outlet conditions. At the beginning of the plate, the BL shout be a laminar one and develop to a turbulent BL with the growth of x (at about Rex= 1e6 ) (please correct me if I'm wrong), according to this I should set the inlet nuTilda to zero.
am I right?

I tried to set the nuTilda to a larger value (like 1e-3) , then the thickness of the BL did grow a lot, but the velocity profile just keep the same shape, when changing the different nuTilda, I even got negtive Ux

Or should I set the inlet U with a turbulent shape?

And I use pisoFoam, LES(DES) - SADDES turbulent model, without any wall function, the yPlus is less than 1. The LES/DES method should be used to a 3D mesh, but I just have 1 grid in spanwise direction with cyclic BC (not empt), can that be the reason of my problem?

If anyone could tell me how to set the BCs for nuTilda or U or P, it will really help me lot!!

PS. in the attachment below, red points represents x=y ( sublayer ? ), and black points is the u profile I got.

Best Regards
Attached Images
File Type: jpg a.jpg (20.1 KB, 36 views)
garrison is offline   Reply With Quote

Old   April 10, 2012, 17:05
Default
  #4
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 273
Rep Power: 18
anishtain4 is on a distinguished road
when you set the inlet boundary condition of turbulent to zero it means you have laminar inlet, which in your case is probable to lead to a laminar boundary layer till the end, since you are not adding any perturbation yourself. there even exist some experiments that flow has not been turbulent for orders larger than 1e5.
There are some methods to determine the k or epsilon in RANS models, but about LES read some papers, I'm not aware of it.
by LES modeling you should have enough grid to capture the kolmogrov scales, read some introductory notes on turbulence about this, 2d turbulent does not exist
anishtain4 is offline   Reply With Quote

Old   April 10, 2012, 22:50
Default
  #5
New Member
 
Garrison, Liang
Join Date: Feb 2012
Posts: 15
Rep Power: 14
garrison is on a distinguished road
Thank you Anishtain,
Then next I'll try with some finer grids and more grids in spanwise to make it a real 3D case.
Best Regards
garrison is offline   Reply With Quote

Old   April 11, 2012, 09:58
Default
  #6
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 273
Rep Power: 18
anishtain4 is on a distinguished road
One more hint, when you are using LES you can't see the log-law profile or sublayer as it is unsteady and turbulence, these profiles are averaged ones which would be captured with RANS models when there is no wall-function used. If you want to capture it check them
anishtain4 is offline   Reply With Quote

Reply


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
Turbulent Boundary Layer on a Flat Plate Hoshang Garda FLUENT 1 November 27, 2013 11:24
2D airfoil simulation, turbulence, boundary layer, unstructured grid louis.langouche@gmail.com FLUENT 0 March 3, 2010 09:19
transient simulation of a rotating rectangle icesniffer CFX 1 August 8, 2009 08:25
Wall turbulence, viscosity, boundary layer Patrick Godon Main CFD Forum 1 November 5, 2003 16:39
The Boundary Condition about the Flat Plate boing Main CFD Forum 1 January 6, 2002 17:53


All times are GMT -4. The time now is 01:33.