Simplifying Raupach's model of shear‐stress partitioning by clarifying the relationship between β and σ and quantifying the shear‐stress probability density based on numerical simulation

Abstract

AbstractWe used computational fluid dynamics to simplify Raupach's model by clarifying the relationship between the parameter β (the ratio of element to surface drag coefficients) and the roughness‐element shape parameter σ using a k‐ω shear‐stress transport (SST) turbulence model with rigid cylindrical roughness elements. Results showed that β ranged between 102 and 242 and the parameter m (a relationship between the average and the peak surface shear stresses) ranged between 0.15 and 0.32 for partitioning of surface shear stress. β decreased with increasing σ, and the relationship followed a power law. We confirmed this relationship by comparing wind tunnel and field measurements for different types of roughness element (blocks, hemispheres, trapezoids and plants based on the results from previous studies). However, the model has some limitations due to the characteristics such as streamline and porosity of plants. The peak mean stress ratio (a) ranged between 2.09 and 3.61. We also confirmed that the probability‐density function for the surface shear stress proposed by Zou et al. (2022) described the distribution of the surface shear stress among roughness elements well, which is essential to accurately estimate surface erosion by wind.