Changes between Version 69 and Version 70 of doc/tec/bc


Ignore:
Timestamp:
Sep 13, 2021 10:35:08 AM (4 years ago)
Author:
suehring
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • doc/tec/bc

    v69 v70  
    563563\end{equation*}
    564564}}}
    565 In order to allow for a dynamic adjustment of the strength of the imposed turbulence, the Reynolds stress tensor is update at regular intervals that can be adjusted individually. In other words, the turbulence generator considers for changing atmospheric stability in case of a multi-day simulation.
    566 
    567 Too date, no proper parametrization of turbulent length scales that works for all stability regimes and within the entire boundary layer is available.
    568 Hence, for the moment the integral length scales are set to
    569  {{{
    570 #!Latex
    571 \begin{equation*}
    572 l_{i,j} = 8 \cdot MIN( \Delta x_i ), \: i \in (1,2,3)\,,
    573 \end{equation*}
    574 }}}
    575 which arises from following considerations: On the one hand, from the numerical point of view the imposed perturbations should not be rapidly eliminated by the numerics.
    576 The numerical dissipation and dispersion, however, act on scales up to 8 times the grid spacing (5th order scheme, see: wiki:/doc/tec/discret]), meaning that scales < 8 times the grid spacing are rapidly dispersed and dissipated due to numerical errors. In order to trigger further turbulence development within the model domain, the imposed correlated turbulence should be on scales larger than the numerically-affected scales.
    577 On the other hand, however, imposing too large length scales could trigger structures that exist throughout the entire model domain, particularly under stable conditions, which in turn could bias the simulation results.
    578 Hence, as a compromise, length scales are set to 8 times the minimum grid spacing.
    579 
    580 Note, for ''z''>''z,,i,,'' the components of the stress tensor, length- and timescales are set to zero so that effectively no synthetic turbulence is imposed above the boundary-layer height (also saving computational costs).  //
    581  
    582 '''Please note, method 2 currently undergoes extensively testing.'''
     565In order to allow for a dynamic adjustment of the strength of the imposed turbulence, the Reynolds stress tensor is update at regular intervals that can be adjusted individually. In other words, the turbulence generator considers for changing atmospheric stability in case of a multi-day simulation. For further information about this approach we refer to [https://doi.org/10.5194/gmd-14-5435-2021 Kadasch et al. (2021)].
    583566
    584567At this point we emphasize that using the turbulence generator from [#xie2008 Xie and Castro (2008)] only generates turbulence which is correlated in space and time but not necessarily generate realistic turbulent structures. Large coherent structures like e.g. hexagonal pattern as typically observed in a convective boundary layer, however, cannot be generated by this method.
    585568Further, we want to add that turbulence is only added to the three wind components. No perturbations are added to the subgrid-scale turbulent-kinetic energy and potential temperature.
    586569
    587 If switched on, the turbulence generator imposed turbulent fluctuations on all lateral boundaries with Dirichlet boundary conditions for the velocity components. For example, if the offline nesting is switched on, where all four lateral boundaries are non-cyclic, the turbulence generator applied at all lateral boundaries, even though a lateral boundary is also an outflow boundary layer.
    588 
    589570== References ==
    590571* [=#holtslag] '''Holtslag AAM, Bruin HARD.''' 1988. Applied modelling of the night-time surface energy balance over land. J. Appl. Meteorol. 27: 689–704.
     
    600581* [=#lund1998] '''Lund TS, Wu X, Squires KD.''' 1998. Generation of turbulent inflow data for spatially-developing boundary layer simulations. J. Comput. Phys. 140: 233–258.
    601582
     583
     584* [=#kadasch2021] '''Kadasch E, Sühring M, Gronemeier T, Raasch S.''' 2021. Mesoscale nesting interface of the PALM model system 6.0. Geosci. Model Dev., 14: 5435–5465.
     585
    602586* [=#kataoka2002] '''Kataoka H, Mizuno M.''' 2002. Numerical flow computation around aerolastic 3d square
    603587  cylinder using inflow turbulence. Wind Struct. 5: 379–392.