CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Surface Roughness

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 6, 2011, 05:16
Default Surface Roughness
  #1
New Member
 
Jochem
Join Date: May 2010
Posts: 28
Rep Power: 16
Jochem is on a distinguished road
Hello,

I am trying to simulate windflow over a flat plane. In this case I am trying to have the same windprofile at the inlet and the outlet.

For the inlet I am using atmBoundaryinletVelocity. I use a certain surface roughness(Z_0) for the inlet and after some trial and error this works. Then I use nutkRoughWallfunction (from turbineSiting tutorial) to apply I certain roughness to the flat plane. If I then use the same roughness as for the inlet the wind profile should be constant over the hole Case, I guess?

On the forum I found how to fill in the parameters for nutkRoughWallfunction:

Ks = 20*Z_0 and Cs = 0.327 (always)

I tried this for different roughness but I never get the same wind profile on the inlet and the outlet. Can someone tell me why? In the next figure you can see the difference between the windspeed.

Best regards,

Jochem
Attached Images
File Type: jpg Windprofile.jpg (53.3 KB, 63 views)
Jochem is offline   Reply With Quote

Old   July 6, 2011, 06:33
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 22
Bernhard is on a distinguished road
What kind of turbulence modelling are you using?
How does your mesh look like?
What is the boundary condition at the outlet?
What does this profile look like half-way the domain?
Is the inlet profile identical to the outlet profile for smooth walls? Or is it really the roughness that introduces the problem?
Bernhard is offline   Reply With Quote

Old   July 6, 2011, 07:17
Default
  #3
New Member
 
Jochem
Join Date: May 2010
Posts: 28
Rep Power: 16
Jochem is on a distinguished road
Hello Bernard,

Thanks for the quick reply.

I am using the kepsilon turbelence model with the simpleFoam solver.

For now i am just using blockMesh and apllying roughness to the ground patch. If I am able to apply different roughness, I will use snappyHexMesh. It seemed better to start with the most simple case.

The boundary at the outlet for U is inletoutlet for p is uniformfixedValue.

I also tried to plot half way the domain, but actually the wind profile starts changing from the begin of the Case, so I dont think It has something to do with a wrong boundary condition for the outlet.

I havent try to set the surface for smooth walls, I will try to do this now.

I hope you can see something is wrong with these answers, I am struggling with this problem for some time.

Regards,

Jochem
Jochem is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
nutRoughWallFunction and Surface Roughness Documented AlanR OpenFOAM Programming & Development 64 January 6, 2024 08:13
[Gmsh] Error : Self intersecting surface mesh, computing intersections & Error : Impossible velan OpenFOAM Meshing & Mesh Conversion 3 October 22, 2015 12:05
[Gmsh] Problem with Gmsh nishant_hull OpenFOAM Meshing & Mesh Conversion 23 August 5, 2015 03:09
Surface roughness alefem FLOW-3D 6 June 3, 2010 13:49
CFX4.3 -build analysis form Chie Min CFX 5 July 13, 2001 00:19


All times are GMT -4. The time now is 07:41.