Boundary condition specifies a atmospheric boundary layer (ABL) velocity inlet profile given the friction velocity value, flow direction n and direction of the parabolic coordinate z. More...
#include <incompressibleRASModels/atmBoundaryLayerInletVelocityFvPatchVectorField.H>
Boundary condition specifies a atmospheric boundary layer (ABL) velocity inlet profile given the friction velocity value, flow direction n and direction of the parabolic coordinate z.
U = (Ustar/K) ln((z - zGround + z0)/z0) where: Ustar is the frictional velocity K is karman's constant z0 is the surface roughness lenght z is the verical coordinate zGround is the minumum coordinate value in z direction. and: Ustar = K Uref/ln((Zref + z0)/z0) where: Uref is the reference velocity at Zref Zref is the reference height.
Reference: D.M. Hargreaves and N.G. Wright "On the use of the k-epsilon model in commercial CFD software to model the neutral atmospheric boundary layer" Journal of Wind Engineering and Industrial Aerodynamics 95(2007) 355-369.
NOTE: D.M. Hargreaves and N.G. Wright recommend Gamma epsilon in the k-epsilon model should be changed from 1.3 to 1.11 for consistency. The roughness height (Er) is given by Er = 20 z0 following the same reference