Changes between Version 1 and Version 2 of doc/tec/canopy

Timestamp:

Oct 16, 2017 10:17:11 AM (7 years ago)

Author:

kanani

Comment:

--

doc/tec/canopy

@@ -1 @@
-A description of the embedded plant canopy model will follow.
+= 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., [#raupach1996 Raupach et al., 1996;] [#dupont2009 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 [#shaw1992 Shaw and Schumann (1992)] and [#watanabe2004 Watanabe (2004)] by adding the term ''C'',,ui,, to the momentum equation (see Eq. one in Sect [wiki:/doc/tec/gov governing equations]:
+{{{
+#!Latex
+\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, ''C'',,ui,, represents the net resolved-scale dynamical effect of the canopy, averaged over the respective grid volume. ''c'',,d,, is the canopy drag coefficient with typical values around ''0.2'' (e.g., [#cescatti2004 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., [#yi2008 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. [wiki:/doc/tec/sgs turbulence closure]:
+{{{
+#!Latex
+\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 [#shaw1992 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., [#shaw2003 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 ([#cassiani2008 Cassiani et al., 2008]; [#finnigan2009 Finnigan et al., 2009]; [#dupont2009 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 [wiki:/doc/tec/gov governing equations]:
+{{{
+#!Latex
+\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 [#brown1966 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):
+{{{
+#!Latex
+\begin{align*}
+  &
+  Q_{\theta}(z)=Q_{\theta}(z_\mathrm{c})\,\exp\left(-\eta\,\text{LAI}\right),
+\end{align*}
+}}}
+with
+{{{
+#!Latex
+\begin{align*}
+   \text{LAI} = \int_{z}^{z_\mathrm{c}}\;{\text{LAD}}\,\mathrm{d}z
+\end{align*}
+}}}
+where ''Q'',,θ,,''(z'',,c,,'')'' is the prescribed heat flux at the top of the canopy layer ''z'',,c,, 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 [wiki:/doc/tec/gov governing equations]):
+{{{
+#!Latex
+\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., [#watanabe2004 Watanabe, 2004]).
+
+== Recent applications ==
+PALM simulations with the embedded canopy model were recently performed by [#kanani2014 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 [#kanani2014b 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 [#kanani2014c Kanani et al. (2014a], [#kanani2014d b]).
+
+[[Image(09.png,600px,border=1)]]
+
+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 km^3^).
+
+== References ==
+
+* [=#raupach1996] '''Raupach MR, Finnigan JJ, Brunet Y.''' 1996. Coherent eddies and turbulence in vegetation canopies: the mixing-layer analogy. Bound.-Lay. Meteorol. 78: 351–382.
+
+* [=#dupont2009] '''Dupont S, Brunet Y.''' 2009. Coherent structures in canopy edge flow: a large-eddy simulation study. J. Fluid Mech. 630: 93–128.
+
+* [=#shaw1992] '''Shaw RH, Schumann U.''' 1992. Large-eddy simulation of turbulent flow above and within a forest. Bound.-Lay. Meteorol. 61: 47–64.
+
+* [=#watanabe2004] '''Watanabe T.''' 2004. Large-eddy simulation of coherent turbulence structures associated with scalar ramps over plant canopies. Bound.-Lay. Meteorol. 112: 207–341.
+
+* [=#cescatti2004] '''Cescatti A, Marcolla B.''' 2004. Drag coefficient and turbulence intensity in conifer canopies. Agr. Forest Meteorol. 121: 197–206.
+
+* [=#yi2008] '''Yi C.''' 2008. Momentum transfer within canopies. J. Appl. Meteorol. 47: 262–275.
+
+* [=#shaw2003] '''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.
+
+* [=#cassiani2008] '''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.
+
+* [=#finnigan2009] '''Finnigan JJ, Shaw RH, Patton EG.''' 2009. Turbulence structure above a vegetation canopy. J. Fluid Mech. 637: 387–424.
+
+* [=#brown1966] '''Brown KW, Covey W.''' 1966.  The energy-budget evaluation of the micro-meteorological transfer process within a cornfield. Agr. Meteorol. 3: 73–96.
+
+* [=#kanani2014] '''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.
+
+* [=#kanani2014b] '''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. [http://dx.doi.org/10.1007/s10546-014-9986-3 doi].
+
+* [=#kanani2014c] '''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. [http://dx.doi.org/10.5446/14311 doi].
+
+* [=#kanani2014d] '''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. [http://dx.doi.org/10.5446/14368 doi].
+
+
+
+