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

atmBoundaryLayerInletVelocity

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 3 Post By vatavuk

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 25, 2016, 13:08
Default atmBoundaryLayerInletVelocity
  #1
Member
 
Almond Wong
Join Date: May 2016
Posts: 68
Rep Power: 10
BakedAlmonds is on a distinguished road
Hi, can some explain the terms use for this file?

While looking into the turbineSiting tutorial:
the inlet is:
inlet
{
type atmBoundaryLayerInletVelocity;
#include "include/ABLConditions"
}

and the subsequent ABLConditions file has this:
Uref 10.0;
Zref 20;
zDir (0 0 1);
flowDir (1 0 0);
z0 uniform 0.1;
zGround uniform 935.0;
value $internalField;

I would like to know what does the terms in the ABLConditions file represent.

Many thanks in advance
BakedAlmonds is offline   Reply With Quote

Old   May 25, 2016, 13:32
Default
  #2
Member
 
Almond Wong
Join Date: May 2016
Posts: 68
Rep Power: 10
BakedAlmonds is on a distinguished road
If someone could provide some example, it will be awesome!!!

Last edited by BakedAlmonds; May 25, 2016 at 17:14.
BakedAlmonds is offline   Reply With Quote

Old   May 26, 2016, 12:35
Default
  #3
Senior Member
 
Paulo Vatavuk
Join Date: Mar 2009
Location: Campinas, Brasil
Posts: 200
Rep Power: 18
vatavuk is on a distinguished road
Hi BakedAlmonds,

In the C++ documentation you may look at the code for atmBoundaryLayer.H, there you will find the description of the variables.
This boundary condition uses the logaritmic profile for the wind inlet velocity:

U = (Ustar/kappa) ln( (z - zGround + z0)/z0 )

where:
kappa is the von Karman constant (=0.41)
z is the vertical coordinate
zGround is the z position of the ground
z0 is the terrain rugosity

Ustar can be calculated imposing that at height Zref you have velocity Uref this results in

Ustar=Uref*kappa/( ln( (Zref+z_0)/z_0 )

Concluding
Uref is the velocity a height Zref
Zref is the height (z-zGround) where the velocity is Uref
zDir is the vector that indicates the vertical direction
flowDir is the vector that indicates the wind direction
z0 is the terrain rugosity
zGround is the z coordinate of the ground

I hope this helps.

Paulo
aow, alanaj5 and Hamzeh193 like this.
vatavuk is offline   Reply With Quote

Old   May 31, 2016, 09:29
Default
  #4
Member
 
Almond Wong
Join Date: May 2016
Posts: 68
Rep Power: 10
BakedAlmonds is on a distinguished road
Hi Paulo, thanks for the reply. I had an attempt to try out the atm function with the following as my boundary conditions:

ABLConditions:
Uref 10.0;
Zref 500;
zDir (0 0 1);
flowDir (1 0 0);
z0 uniform 0.1;
zGround uniform 0.0;
value $internalField;

I set my top and sides patches to symmetry.

Then I realise two weird occurrences. Hope someone could assist me on this.

1) I tried to verify the inlet velocity profile so I added probes to check.

My inlet patch are bounded by vertices :
( -1100 , -1100 , 0)
( -1100 , 1100 , 0)
( -1100 , 1100 ,1100)
( -1100 , -1100 , 1100)

then i place my probes as
( -1100 0 10 )
.
.
.
.
( -1100 0 500 )

I realise that the velocity at Zref the very first time step is less than the Uref. Shouldn't it be very close if not more than that specified?

2) The velocity at the walls act like a viscous BL. Shouldn't it act as if its frictionless?
BakedAlmonds is offline   Reply With Quote

Old   March 15, 2017, 09:21
Default
  #5
New Member
 
Nathan
Join Date: Mar 2017
Posts: 2
Rep Power: 0
Nathan77 is on a distinguished road
Hello,
I am a French student, and I try to integrate a logarithmic wind profile in a simulation highlighting the eccentricity around a building. Could you help me by explaining how to use the atmboundarylayerinletvelocity function? How can we integrate it in place of the uniform wind field that we currently have? Would you have any examples on which we could rely?

Thanks in advance
Nathan
Nathan77 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



All times are GMT -4. The time now is 18:31.