streamWiseVel_powerLaw Subroutine

private subroutine streamWiseVel_powerLaw(velSW, velTau, distToBnd, viscKine, nElems)

This routines computes streamWise velocity component from friction velocity and distance to boundary using two equation Werner and Wengle profile. Ref to following paper for Schmitt three layer equations. Haussmann, M. et al. (2019) ‘Large-eddy simulation coupled with wall models for turbulent channel flows at high Reynolds numbers with a lattice Boltzmann method — Application to Coriolis mass flowmeter’, Computers & Mathematics with Applications. Elsevier Ltd, 78(10), pp. 3285–3302.


real(kind=rk), intent(out) :: velSW(:)

Stream-wise velocity component from wall model

real(kind=rk), intent(in) :: velTau(:)

Friction velocity computd from wall model

real(kind=rk), intent(in) :: distToBnd(:)

Distance to the boundary in the discrete normai direction

real(kind=rk), intent(in) :: viscKine(:)

Kinematic viscosity

integer, intent(in) :: nElems

Number of elements in input and output arrays



integer, private :: iElem
real(kind=rk), private :: yPlus