Three-dimensional simulations and wind-tunnel experiments on airflow over isolated forest stands

Prediction of windthrow risk to individual or groups of retained trees in harvested stands requires an improved understanding of canopy airflow dynamics. Large-eddy simulations were used to simulate wind-tunnel experiments in two and three dimensions to compare with observations for model validation and to address parameter space considerations for the design of subsequent retention pattern experiments. The three-dimensional simulations were similar to the observed wind-tunnel data for the statistical profiles for but there were greater differences in skewness and kurtosis. These results were obtained using a common leaf-area drag formulation without either skin friction or speed dependent drag that enables scaling with U ₀ (ambient wind speed) and h (height of the canopy). This scaling results in a single non-dimensional parameter h/h _c where h _c (x, y, z) is the momentum range resulting from the canopy drag. The validity of the model scaling was tested using two-dimensional simulations. The irrotational component of the flow (potential flow) was found to be important when defining vertical domain limitations and has significant implications for time dependent flow (i.e. turbulent conditions) when considering retention pattern design. The sudden onset of drag associated with the isolated stand presents some unexpected challenges. The horizontal scales of the shearing instabilities were simulated in two dimensions and found to range between 2h for early times to 7h for later times. The early-time horizontal scales are in the range of logical retention pattern scales and as such need to be taken into account as part of the parameter space, i.e. a range of retention pattern lengths need consideration.