# help with boundary conditions

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

January 11, 2012, 14:57
help with boundary conditions
#1
New Member

Join Date: Jun 2009
Posts: 27
Rep Power: 8
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
 test1.jpg (15.8 KB, 27 views)
Attached Files
 test.zip (8.8 KB, 6 views)

 January 12, 2012, 07:15 Questions & clues about your case #2 Member   Frederic Collonval Join Date: Apr 2009 Location: Technische Universitaet Munich - Lehrstuhl fuer Thermodynamik Posts: 53 Rep Power: 8 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.

 January 12, 2012, 10:31 #3 New Member   Join Date: Jun 2009 Posts: 27 Rep Power: 8 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 11:38.

January 16, 2012, 06:15
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: 8
Hey

Quote:
 Originally Posted by dshawul 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post lost.identity CFX 41 May 22, 2013 07:21 Anindya Main CFD Forum 24 January 11, 2012 14:40 NickolasPl OpenFOAM Programming & Development 2 May 19, 2011 05:37 vaina74 Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 May 27, 2010 09:38 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 20:30.