Home


About us
PALM
PALM-4U
Installation
Release notes
Documentation
Getting Help
Gallery
Developers area

Version 2 (modified by kanani, 7 years ago) (diff)

--

Canopy model

The embedded plant canopy model allows for studying the turbulent flow inside and above vegetation canopy. It is well-known that vegetation canopy effects on the surface--atmosphere exchange of momentum, energy and mass can be rather complex and can significantly modify the structure of the ABL, particularly in its lower part (e.g., Raupach et al., 1996; Dupont and Brunet, 2009). It is thus not possible to describe such processes by means of the roughness length and surface fluxes of sensible and latent heat. The canopy model in PALM accounts for the vertically extended drag, release of heat, plant evaporation and leaf-air-interactions that are functions of height within the canopy layer.

Dynamical effects of the plant canopy are based on the assumption that the canopy acts as a~sink for momentum due to form (pressure) and viscous drag forces. This sink for momentum is modeled following Shaw and Schumann (1992) and Watanabe (2004) by adding the term Cui to the momentum equation (see Eq. one in Sect governing equations:

\begin{align*} & \frac{\partial u_i}{\partial t} = \ldots - \underbrace{c_{\mathrm{d}}\,\text{LAD}\,\sqrt{u_i^2} u_{i}}_{C_{u_{i}}}\;. \end{align*}

Here, Cui represents the net resolved-scale dynamical effect of the canopy, averaged over the respective grid volume. cd is the canopy drag coefficient with typical values around 0.2 (e.g., Cescatti and Marcolla, 2004), and LAD is the leaf area density (available leaf area per unit volume). As an example, LAD is rather constant with height within crop fields, whereas it is often very heterogeneous in forests, where most of the leaf area is concentrated in the trees' crown space (e.g., Yi, 2008).

The effect of the canopy on the SGS turbulence is considered by adding a similar sink term to the prognostic equation for SGS-TKE (see Sect. turbulence closure:

\begin{align*} & \frac{\partial e}{\partial t} = \ldots - \underbrace{2 c_{\mathrm{d}} \text{LAD} \sqrt{u_i^2}\,e}_{C_e}\;. \end{align*}

This approach was suggested by Shaw and Schumann (1992) and is based on the assumption that SGS-TKE is dissipated by the canopy due to the rapid dissipation of wake turbulence in the lee of plant elements. This rapid break-down of turbulence is also known as the spectral shortcut (e.g., Shaw and Patton, 2003). This type of canopy model has been successfully applied by various authors to study turbulent flows inside and above homogeneous as well as heterogeneous canopies such as forest edges (Cassiani et al., 2008; Finnigan et al., 2009; Dupont and Brunet, 2009, among others).

In case of incoming solar radiation the plant canopy acts as a source for heat. It is assumed that this warming of the foliage by solar radiation results in a warming of the surrounding air. This process is considered by adding a source term Cθ to the prognostic equation for θ (see Eq. three in Sect governing equations:

\begin{align*} & \frac{\partial \theta_\mathrm{l}}{\partial t} = \ldots + \underbrace{\frac{\partial Q_{\theta}}{\partial z}}_{C_\theta}\;. \end{align*}

In order to account for the fact that solar radiation can penetrate different layers of the canopy, based on the leaf area, an exponential decay function for the upward vertical kinematic heat flux Qθ after Brown and Covey (1966)is used. Qθ is derived at each height inside the canopy by means of the downward cumulative leaf area index (LAI):

\begin{align*} & Q_{\theta}(z)=Q_{\theta}(z_\mathrm{c})\,\exp\left(-\eta\,\text{LAI}\right), \end{align*}

with

\begin{align*} \text{LAI} = \int_{z}^{z_\mathrm{c}}\;{\text{LAD}}\,\mathrm{d}z \end{align*}

where Qθ(zc) is the prescribed heat flux at the top of the canopy layer zc and η is the extinction coefficient set to 0.6. Additionally, contributions by sinks/sources for q and s are considered in the canopy model by adding additional terms Cφ to the scalar transport equations (see Eqs. four and five in Sect governing equations):

\begin{align*} & \frac{\partial \varphi}{\partial t} = \ldots - \underbrace{c_\varphi {\text{LAD}} \sqrt{u_i^2} \left(\varphi-\varphi_{\mathrm{c},0}\right)}_{C_\varphi}\;, \end{align*}

where φ ∈ {q, s} and cφ is a user-defined scalar exchange coefficient. φc,0 and φ are the scalar concentrations at a leaf surface and in the surrounding air volume, respectively. This approach is based on the assumption that the scalar sink/source strength depends on the concentration gradient between the leaf surface and the surrounding air (e.g., Watanabe, 2004).

Recent applications

PALM simulations with the embedded canopy model were recently performed by Kanani et al. (2014c) to study the flow adjustment downstream of a transition from an unforested (clearing) to a forested surface patch. In this study the LES results were validated against multidimensional field and wind-tunnel data. In the high-resolution follow-up study of Kanani-Sühring and Raasch (2015), a detailed analysis of the turbulent scalar transport within the canopy layer was successfully performed for the first time by means of LES. Figure 11 shows exemplarily the flow at a forest edge, where an internal boundary layer developed above the forest due to the extended drag of the canopy. See also associated animations in Kanani et al. (2014a, b).

Figure 11: Snapshot of the absolute value of the 3-D rotation vector of the velocity field above a forest canopy downstream of a grassland-to-forest transition (forest volume marked by green isosurface). Pink and yellow colors illustrate strong and weak turbulence, respectively. A neutrally-stratified open-channel flow was simulated with the mean flow direction from left to right, i.e. perpendicular to the windward forest edge, using an equidistant grid spacing of 3 m. The figure shows only a subregion of the simulation domain (2 x 1 x 0.4 km3).

References

  • Raupach MR, Finnigan JJ, Brunet Y. 1996. Coherent eddies and turbulence in vegetation canopies: the mixing-layer analogy. Bound.-Lay. Meteorol. 78: 351–382.
  • Dupont S, Brunet Y. 2009. Coherent structures in canopy edge flow: a large-eddy simulation study. J. Fluid Mech. 630: 93–128.
  • Shaw RH, Schumann U. 1992. Large-eddy simulation of turbulent flow above and within a forest. Bound.-Lay. Meteorol. 61: 47–64.
  • Watanabe T. 2004. Large-eddy simulation of coherent turbulence structures associated with scalar ramps over plant canopies. Bound.-Lay. Meteorol. 112: 207–341.
  • Cescatti A, Marcolla B. 2004. Drag coefficient and turbulence intensity in conifer canopies. Agr. Forest Meteorol. 121: 197–206.
  • Yi C. 2008. Momentum transfer within canopies. J. Appl. Meteorol. 47: 262–275.
  • Shaw RH, Patton EG. 2003. Canopy element influences on resolved- and subgrid-scale energy within a large-eddy simulation. Agr. Forest Meteorol. 115: 5–17.
  • Cassiani M, Katul GG, Albertson JD. 2008. The effects of canopy leaf area index on airflow across forest edges: large-eddy simulation and analytical results. Bound.-Lay. Meteorol. 126: 433–460.
  • Finnigan JJ, Shaw RH, Patton EG. 2009. Turbulence structure above a vegetation canopy. J. Fluid Mech. 637: 387–424.
  • Brown KW, Covey W. 1966. The energy-budget evaluation of the micro-meteorological transfer process within a cornfield. Agr. Meteorol. 3: 73–96.
  • Kanani F, Träumner K, Ruck B, Raasch S. 2014c. What determines the differences found in forest edge flow between physical models and atmospheric measurements? -- an LES study. Meteorol. Z. 23: 33–49.
  • Kanani-Sühring F, Raasch S. 2015. Spatial variability of scalar concentrations and fluxes downstream of a clearing-to-forest transition: a large-eddy simulation study. Bound.-Lay. Meteorol. 155: 1–27. ​doi.
  • Kanani F, Maronga B, Knoop H, Raasch S. 2014a. Large-eddy simulation of a forest-edge flow -- adjustment of a turbulent flow to the changing surface conditions at a clearing-to-forest transition. Computer animation. ​doi.
  • Kanani F, Maronga B, Knoop H, Raasch S. 2014b. Large-eddy simulation of the scalar transport in a forest-edge flow -- spatial variability of the scalar distribution and the scalar transport downstream of a clearing-to-forest transition. Computer animation. ​doi.

Attachments (2)

Download all attachments as: .zip

                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |