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

help with boundary conditions

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 11, 2012, 13:57
Default help with boundary conditions
  #1
New Member
 
Join Date: Jun 2009
Posts: 27
Rep Power: 7
dshawul is on a distinguished road
I am trying to do ABL simulation using OF on an open area. It is a simple case where a uniform velocity is applied at the inlet and a uniform pressure at the outlet. I am getting high values of pressure at the bottom of the inlet (about 400pa as shown in the figure). The more I refine the mesh the larger its value is, so I am wondering if my boundary conditions are correct or not ? I would really really appreciate your help in solving this problem. I have attached the case file with the boundary conditions I am using.
Thank you
Attached Images
File Type: jpg test1.jpg (15.8 KB, 26 views)
Attached Files
File Type: zip test.zip (8.8 KB, 6 views)
dshawul is offline   Reply With Quote

Old   January 12, 2012, 06:15
Default Questions & clues about your case
  #2
Member
 
Frederic Collonval
Join Date: Apr 2009
Location: Technische Universitaet Munich - Lehrstuhl fuer Thermodynamik
Posts: 53
Rep Power: 7
fcollonv is on a distinguished road
Hello Daniel,

Your test case looks very similar to a turbulent flat plate case (I don't know what ABL means - abbreviation are not recommended...). If so, here is a known test case that could be useful: http://www.grc.nasa.gov/WWW/wind/valid/fpturb/fpturb.html
You will see that for incompressible flow, the first point of the plate is a singularity (hence your high value there that get higher but more localized when the mesh is finer). It's therefore better to put the inlet ahead of the interesting section (see nasa test case)

I have also some questions:
* why are you using a 3D mesh, when your case is obviously 2D ??
* why are you setting the pressure to 1 at the outlet? The value of the pressure in incompressible case is relative, so the reference is usually 0 (as set in your system/fvSolution::SIMPLE). And consequently the outlet is usually set to that reference.

Hope that it will help you,

Frederic
__________________
Frederic Collonval
Technische Universität München
Thermodynamics Dpt.
fcollonv is offline   Reply With Quote

Old   January 12, 2012, 09:31
Default
  #3
New Member
 
Join Date: Jun 2009
Posts: 27
Rep Power: 7
dshawul is on a distinguished road
Dear Frederic
Thank you very much! That was really helpful. ABL = atmospheric boundary layer flow. I didn't know that p would be singular at bottom of the inlet. As in the flat plate case, I have an internal boundary layer developed in my test case too. I want to avoid that by applying a log-law inlet velocity profile. In which case the profile should come out intact at the other end too. I wanted to try uniform profile before doing that which caused this "problem".

I have corrected the value of p to 0 at the outlet and will be using 2D mesh for my next tests.

I have other questions if you don't mind:

My y+ values are very big and I can not bring it down by using a fine mesh close to the wall because of another constraint that yp > Ks. Ks is the equivalent sand grain roughness for a rough wall as used in nutRoughWall of OF. I will be using rough walls later which could have ks values as high as 15m. What does OF do when y+~=20000 , does it still apply log-law ?

thanks again
Daniel

Last edited by dshawul; January 12, 2012 at 10:38.
dshawul is offline   Reply With Quote

Old   January 16, 2012, 05:15
Default Law of the wall and big y+
  #4
Member
 
Frederic Collonval
Join Date: Apr 2009
Location: Technische Universitaet Munich - Lehrstuhl fuer Thermodynamik
Posts: 53
Rep Power: 7
fcollonv is on a distinguished road
Hey

Quote:
Originally Posted by dshawul View Post
My y+ values are very big and I can not bring it down by using a fine mesh close to the wall because of another constraint that yp > Ks. Ks is the equivalent sand grain roughness for a rough wall as used in nutRoughWall of OF. I will be using rough walls later which could have ks values as high as 15m. What does OF do when y+~=20000 , does it still apply log-law ?
As far as I know, it will apply the log-law for any y+ bigger than the y+_laminar. But it's up to you if you want to implement something different or not using it...

Cheers,

Fred
__________________
Frederic Collonval
Technische Universität München
Thermodynamics Dpt.
fcollonv 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
symmetry boundary conditions in cfx lost.identity CFX 41 May 22, 2013 07:21
Impinging Jet Boundary Conditions Anindya Main CFD Forum 24 January 11, 2012 13:40
OpenFOAM Variable Velocity Boundary Conditions NickolasPl OpenFOAM Programming & Development 2 May 19, 2011 05:37
boundary conditions and mesh exporting vaina74 Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 May 27, 2010 09:38
Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15


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