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

Turbulent boundary layer on a flat plate

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 16, 2009, 10:40
Default Hi everyone, I am running s
  #1
New Member
 
Sébastien Bocquet
Join Date: Mar 2009
Posts: 21
Rep Power: 8
seb62 is on a distinguished road
Hi everyone,

I am running some 2D calculations over a flat plate, and comparing the openFOAM results with flu3m (Navier-Stokes CFD code from ONERA), Fluent and experimental/theoretical data. I am using 2 meshes.

Here is a description of the case :

Application : simpleFoam
Kinematic viscosity nu=1.57e-5 mē/s
Plate length L=1m
Re (end of plate) = 4.4e6
divScheme : Gauss upwind
Turbulence model : Spalart Allmaras
First mesh : 0.15 < Y+ (first node) < 0.2
Second mesh (identical to first mesh with 35 first cells close to the plate removed) :
20 < Y+ (first node) < 23

Boundary conditions :

Inlet : inlet
Top : pressureOutlet
Exit : extrapolatedOutlet
Plate : wallFunctions


U(inlet)=69.177m/s
nut(inlet)=1.5e-5 mē/s
nuTilda(inlet)=1.5e-5 mē/s
p(inlet)=p(top)=pref=0

Initial conditions = boundary conditions at inlet


For the first mesh :

Plots of U+=f(y+) :

At x=0.12 :




At x=0.94 :




For the second mesh :

Plots of U+=f(y+) and mut/mu=f(y+) :

At x=0.12 :



At x=0.94 :



Comments :

The mesh with y+ around 20 gives satisfying results especially on the velocity profile at x=0.94 and the local friction Cf. At x=0.12, we have an offset on the U+ velocity profile in the log region and a smaller turbulent viscosity compared to flu3m and Fluent

The finer mesh (y+ around 0.15) induces a significant offset on the U+ profile in the log region, with a smaller turbulent viscosity than the other calculations and an underestimated local friction (2.5e-3 compared to 3e-3 for the theoretical/experimental results at the end of the plate).

If someone encountered a similar problem...

Thanks in advance
seb62 is offline   Reply With Quote

Old   January 17, 2009, 04:30
Default Hi Sebastien, I'm working on
  #2
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
Hi Sebastien,
I'm working on a similar case, a RANS over a flat plate at Mach 0.06 and Re 1.3x10^5.
If you choose a wall-function treatment, you have to put your first cell at y+ around 30, if you want to use a y+ less than 1, you have to use a lowRe model, so if find some lowRe model suitable for a flat plate you can use it!
Ivan
ivan_cozza 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
Turbulent Boundary Layer on a Flat Plate Hoshang Garda FLUENT 1 November 27, 2013 11:24
Turbulent boundary layer on a flat plate seb62 OpenFOAM Running, Solving & CFD 0 January 16, 2009 10:33
Turbulent boundary layer on a flat plate seb62 OpenFOAM Running, Solving & CFD 0 January 16, 2009 10:30
Turbulent Bounadry layer on a ZPG flat plate Eldude Main CFD Forum 2 October 16, 2008 02:44
boundary layer of a flat-plate dwi s Main CFD Forum 1 November 9, 2000 02:45


All times are GMT -4. The time now is 23:54.