Turbulence structure in a diabatically heated forest canopy composed of fractal Pythagoras trees

Abstract

We investigate the turbulent flow through a heterogeneous forest canopy by high-resolution numerical modeling. For this purpose, a novel approach to model individual trees is implemented in our large-eddy simulation (LES). A group of sixteen fractal Pythagoras trees is placed in the computational domain and the tree elements are numerically treated as immersed boundaries. Our objective is to resolve the multiscale flow response starting at the diameter of individual tree elements up to the depth of the atmospheric surface layer. A reference run, conducted for the forest flow under neutral thermal stratification, produces physically meaningful turbulence statistics. Our numerical results agree quantitatively with data obtained from former field-scale LESs and wind tunnel experiments. Furthermore, the numerical simulations resolve vortex shedding behind individual branches and trunks as well as the integral response of the turbulent flow through the heterogeneous forest canopy. A focus is the investigation of the turbulence structure of the flow under stable thermal stratification and in response to the heating of the fractal tree crowns. For the stratified flows, statistical quantities, e.g. turbulent kinetic energy and vorticity, are presented and the turbulent exchange processes of momentum and heat are considered in detail. The onset and formation of coherent structures such as elevated shear layers above the diabatically heated forest canopy are analyzed. For the stably stratified flow, temperature ramps above the forest canopy were simulated in agreement with previous observations. Thermally driven vortices with a typical diameter of the canopy height were simulated when the tree crowns were diabatically heated. The impact of the coherent flow structures on the heat flux is investigated. Keywords Atmospheric boundary layer • Forest canopy • Fractal trees • Heat • Turbulence IntroductionFlow through forest canopies transports momentum, heat, aerosol particles and trace gases as carbon dioxide, oxygen and water vapor. The associated atmospheric exchange processes between the forest and the atmosphere occur on a broad range of scales. Photosynthesis and transpiration on the leaves [32] are examples representing Communicated by R. Klein.

Vegetation in Drylands: Effects on Wind Flow and Aeolian Sediment Transport

“…To get around this issue, flow properties can be calculated using much coarser meshes, through large eddy simulations (LES) [52,[165][166][167]. Dupont et al [52] coupled a LES model with a physically based saltation model to resolve turbulent wind flow over heterogeneous vegetated landscapes composed of trees and shrubs.…”

On the Interpretation of Gravity Wave Measurements by Ground-Based Lidars

“…In general, EULAG owes its versatility to a unique design that combines a rigorous theoretical formulation in generalized curvilinear coordinates ( [40]) with non-oscillatory forward-in-time (NFT) differencing for fluids built on MPDATA, which is based on the convexity of upwind advection ( [34,38]) and a robust, exact-projection-type elliptic Krylov solver [34]. The flow solver has been applied to a wider range of scales simulating various problems, like boundary layer flow through complex geometries [41,42], baroclinic instability [43], gravity waves [44][45][46] or even solar convection [47]. Recently, the solver has been advanced to unstructured meshes; see [48].…”

A modified Darcy’s law

“…EULAG solves the governing equations of motion either in an Eulerian or in a Lagrangian form. Here, we solve the Boussinesq-approximated governing equations in a form as described by Schröttle and Dörnbrack [11]. The immersed boundary (IMB) method, used for the representation of the obstacles in the model, has already been tested successfully against wind tunnel experiments for the flow past a complex building in a neutral boundary layer, e.g., [13].…”

“…Numerically, the porous medium can be treated by immersed boundaries, e.g., [7,11,13]. To implement immersed boundaries, the momentum equations contain additional external forcing terms.…”