CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CD-adapco

atmospheric roughness height

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 27, 2002, 17:30
Default atmospheric roughness height
  #1
Alain BASTIDE
Guest
 
Posts: n/a
Hi,

I work with atmospheric boundary layer.

wind profile of inlet boundary layer is : u=u_tau/kappa*ln(z/z_0) where roughness height z_0=0.07 meter (flat plane)

Turbulence intensity profile : I_t=1/ln(z/z_0)

I use a law of wall : u_plus=1/kappa*ln(z+/z+_0)

In StarCD Methodology page 6-6, i found that y_p (height of the center of cell) should be 4 times larger than the roughness height.

My questions :

1. I'm looking for publications on roughness height and CFD (witch grid can i take), could you help me?

2. Witch y+ range can i give?

3. Is this modelisation good?

4. Witch lenght scale can i take (i just know that L_30m=100m)?

PS : I try a "power law" for the wind profil and i obtain bad results because i use low wind and power low is for strong wind

Thanks
  Reply With Quote

Old   December 27, 2002, 21:51
Default Re: atmospheric roughness height
  #2
CJ Tune
Guest
 
Posts: n/a
For high-Re k-e models, the boundary layer has already been assumed a shape with the 'law of the wall' function. If you want to predict a resultant b/l, I think you should be calculating it directly with a low-Re k-e model.
  Reply With Quote

Old   December 28, 2002, 03:49
Default Re: atmospheric roughness height
  #3
Alain BASTIDE
Guest
 
Posts: n/a
Hi,

My Reynolds number is near 10^5. I don't know that a low Re k-e model is good.

  Reply With Quote

Old   December 29, 2002, 22:52
Default Re: atmospheric roughness height
  #4
CJ Tune
Guest
 
Posts: n/a
I've copied the following from www.adapco-online.com:

------------------------------------------------------- "The models for the RANS equation typically break down near the walls, where there is a very thin boundary layer with large velocity gradients. At high Re, this layer is so thin that it is difficult to get a sufficient number of fluid cells into this zone. Even if this is done, the forms of the k-e equation suitable for the interior do not work very well at the wall, as the flow transitions to a laminar sublayer.

STAR-CD offers three alternatives to deal with this situation. The usual procedure is to use an algebraic "Law of the wall" model in the first cell at the boundary. The boundary layer is assumed to be contained within this cell, and the shear stress at the wall, and the boundary conditions for the K and e equations are computed from the log law. For this to be applicable, the cell next to the wall must be big enough; the value of Y+ at the cell centroid should be greatar than 30. This is rarely a problem, except in the neighbourhood of separation points and reattachment regions. It can be checked easily by plotting the Y+ at the boundary in PROSTAR, after the solution.

At the next level of complexity are the two layer models. In these models, a low Reynolds number representation of the turbulence equation is used in the viscosity influenced near wall region. Since transport equations must be solved in this region, there must be a sufficient grid density in the direction normal to the wall. Typically about 10 cells are placed in the region of Y+ < 30. This usually results in a large cell count, increasing the computational requirements. There is also a (one time) penality where STAR-CD has the determine the wall distances of the viscous layer cells.

In these two layer models, the flow in the wall layer is treated by one of two ways; a one equation model where near-wall flow is simulated by a low Reynolds number model consisting of a transport equation for k and an algebraic prescription for the turbulence length scale; or a mixing length model where an algebraic expression is used for the turbulent viscosity. In either case, the solution is matched to that of the standard k- e equations at the edge of the viscosity-influenced region.

The third alternative for dealing with wall effects is a Low Reynolds Number K-e model. The meshing requirements near the wall are similar to the two layer models. This model uses the standard transport equation for K, but augments the equation for e with an extra term which depends on the distance from the nearest wall. Thus in principle, this model is valid all the way down to the solid wall. The startup computational cost is higher, because the wall distances need to be computed for all the cells." -------------------------------------------------------

With the Hi-Re models, you'll need to adjust the near-wall cell layer (to be within the recommended 30-100 Y+ value) before you get a nice answer. Hi-Re turb models will usually give bad results when accurate calculation of the near-wall flow is required, eg. conjugate heat transfer, which needs an accurate velocity profile to in turn get a good thermal b/l. ICEM even has a commercial product called 'Y+' that is supposed to automatically adjust your near wall cells according to your previous sim results.

OR: you can just TREAT it as if you were calculating for a Low-Re flow with the Low-Re k-e model. Problem now is that you've got to mesh it really fine at the walls. Doesn't say anywhere that it is wrong to use it if your (mean) flow Reynolds number is in the turbulent regime (in fact, they are the same, except for the a slight mod in the e equation).

You could try the two-layer models too.
  Reply With Quote

Old   December 30, 2002, 06:12
Default Re: atmospheric roughness height
  #5
Anton Lyaskin
Guest
 
Posts: n/a
Hi,

by the way, what profile for turbulence length scale do you use?
  Reply With Quote

Old   December 30, 2002, 06:26
Default Re: atmospheric roughness height
  #6
Alain BASTIDE
Guest
 
Posts: n/a
Hi,

I use a linear profile for turbulence length scale but i don't know if that is a good modelisation

  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
how to choose the roughness height icestar FLUENT 10 August 31, 2012 09:00
Roughness Height Bregazzi FLUENT 1 February 19, 2009 07:03
Urgent! 2*Roughness Height >grid height rodi FLUENT 2 January 15, 2008 08:56
atmospheric roughness height Alain BASTIDE Main CFD Forum 0 December 27, 2002 17:32
ROUGHNESS HEIGHT Sumon CFX 3 June 12, 2001 09:38


All times are GMT -4. The time now is 16:26.