fricVel_reichardt_fixedPoint Subroutine

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

This routine computes friction velocity from reichardt wall model profile using fixed-point iteration method

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(inout) :: velTau(:)

Friction velocity computed from wall model. it is inout to provide velTau from previous timestep as initial velTau for fixed-point or Newton iteration solver

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

Stream-wise velocity component from which friction velocity is computed

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


Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: error
integer, private :: iElem
integer, private :: iter
real(kind=rk), private :: velTau_old
real(kind=rk), private :: velTau_new
real(kind=rk), private :: yPlus
real(kind=rk), private :: yPlus_fac