
[Sponsors] 
April 7, 2012, 03:34 
turbulence flat plane boundary

#1 
New Member
Garrison, Liang
Join Date: Feb 2012
Posts: 15
Rep Power: 6 
Hi Foamers!
Recently I am doing a very simple case using SADDES. 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 010), 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 =1e2, 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 

April 10, 2012, 02:01 

#2 
Senior Member
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 187
Rep Power: 9 
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 

April 10, 2012, 03:09 

#3 
New Member
Garrison, Liang
Join Date: Feb 2012
Posts: 15
Rep Power: 6 
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 1e3) , 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 

April 10, 2012, 16:05 

#4 
Senior Member
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 187
Rep Power: 9 
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 

April 10, 2012, 21:50 

#5 
New Member
Garrison, Liang
Join Date: Feb 2012
Posts: 15
Rep Power: 6 
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 

April 11, 2012, 08:58 

#6 
Senior Member
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 187
Rep Power: 9 
One more hint, when you are using LES you can't see the loglaw 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 wallfunction used. If you want to capture it check them


Thread Tools  
Display Modes  


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 07: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 