|
[Sponsors] |
May 25, 2016, 13:08 |
atmBoundaryLayerInletVelocity
|
#1 |
Member
Almond Wong
Join Date: May 2016
Posts: 68
Rep Power: 10 |
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 |
|
May 25, 2016, 13:32 |
|
#2 |
Member
Almond Wong
Join Date: May 2016
Posts: 68
Rep Power: 10 |
If someone could provide some example, it will be awesome!!!
Last edited by BakedAlmonds; May 25, 2016 at 17:14. |
|
May 26, 2016, 12:35 |
|
#3 |
Senior Member
Paulo Vatavuk
Join Date: Mar 2009
Location: Campinas, Brasil
Posts: 200
Rep Power: 18 |
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 |
|
May 31, 2016, 09:29 |
|
#4 |
Member
Almond Wong
Join Date: May 2016
Posts: 68
Rep Power: 10 |
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? |
|
March 15, 2017, 09:21 |
|
#5 |
New Member
Nathan
Join Date: Mar 2017
Posts: 2
Rep Power: 0 |
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 |
|
|
|