|
[Sponsors] |
January 11, 2012, 13:57 |
help with boundary conditions
|
#1 |
New Member
Join Date: Jun 2009
Posts: 27
Rep Power: 16 |
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 |
|
January 12, 2012, 06: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: 17 |
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, 09:31 |
|
#3 |
New Member
Join Date: Jun 2009
Posts: 27
Rep Power: 16 |
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. |
|
January 16, 2012, 05: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: 17 |
Hey
Quote:
Cheers, Fred
__________________
Frederic Collonval Technische Universität München Thermodynamics Dpt. |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Impinging Jet Boundary Conditions | Anindya | Main CFD Forum | 25 | February 27, 2016 12:58 |
symmetry boundary conditions in cfx | lost.identity | CFX | 41 | May 22, 2013 07:21 |
OpenFOAM Variable Velocity Boundary Conditions | NickolasPl | OpenFOAM Programming & Development | 2 | May 19, 2011 05:37 |
[Netgen] boundary conditions and mesh exporting | vaina74 | OpenFOAM Meshing & Mesh Conversion | 2 | May 27, 2010 09:38 |
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues | michele | OpenFOAM Meshing & Mesh Conversion | 2 | July 15, 2005 04:15 |