Consistent inflow boundary conditions for modelling the neutral equilibrium atmospheric boundary layer for the SST k-ω model

2017 ◽  
Vol 24 (5) ◽  
pp. 465-480 ◽  
Author(s):  
Yi Yang ◽  
Zhuangning Xie ◽  
Ming Gu
2020 ◽  
Author(s):  
Eckhard Kadasch ◽  
Matthias Sühring ◽  
Tobias Gronemeier ◽  
Siegfried Raasch

Abstract. In this paper, we present a newly developed mesoscale nesting interface for the PALM model system 6.0, which enables PALM to simulate the atmospheric boundary layer under spatially heterogeneous and non-stationary synoptic conditions. The implemented nesting interface, which is currently tailored to the mesoscale model COSMO, consists of two major parts: (i) the preprocessor INIFOR, which provides initial and time-dependent boundary conditions from mesoscale model output and (ii) PALM's internal routines for reading the provided forcing data and superimposing synthetic turbulence to accelerate the transition to a fully developed turbulent atmospheric boundary layer. We describe in detail the conversion between the sets of prognostic variables, transformations between model coordinate systems, as well as data interpolation onto PALM's grid, which are carried out by INIFOR. Furthermore, we describe PALM's internal usage of the provided forcing data, which besides the temporal interpolation of boundary conditions and removal of any residual divergence includes the generation of stability-dependent synthetic turbulence at the inflow boundaries in order to accelerate the transition from the turbulence-free mesoscale solution to a resolved turbulent flow. We demonstrate and evaluate the nesting interface by means of a semi-idealized benchmark case. We carried out a large-eddy simulation (LES) of an evolving convective boundary layer on a clear-sky spring day. Besides verifying that changes in the inflow conditions enter into and successively propagate through the PALM domain, we focus our analysis on the effectiveness of the synthetic turbulence generation. By analysing various turbulence statistics, we show that the inflow in the present case is fully adjusted after having propagated for about 1.5 eddy turn-over times downstream, which corresponds well to other state-of-the-art methods for turbulence generation. Furthermore, we observe that numerical artefacts in the form of under-resolved convective structures in the mesoscale model enter the PALM domain, biasing the location of the turbulent up- and downdrafts in the LES. With these findings presented, we aim to verify the mesoscale nesting approach implemented in PALM, point out specific shortcomings, and build a baseline for future improvements and developments.


2014 ◽  
Vol 142 (9) ◽  
pp. 3418-3424 ◽  
Author(s):  
M. Inoue ◽  
G. Matheou ◽  
J. Teixeira

An arrangement of a large-eddy simulation (LES) is described that facilitates a spatially developing thermally stratified atmospheric boundary layer (ABL). When the inflow and outflow boundary conditions are specified, the LES of stably stratified ABL turns out to be challenging because spurious reflections of waves at the boundary accumulate inside the domain. To tackle this problem, a fringe method with an auxiliary LES running concurrently is applied to enforce upstream/downstream boundary conditions. An artificial forcing term is applied within a fringe region located at the beginning of the main LES domain in order to ensure statistically stationary inflow boundary conditions. The auxiliary LES, which is horizontally homogeneous in a doubly periodic domain, is used to determine the inflow condition of the main LES domain. The present scheme is used to provide an Eulerian perspective of the stratocumulus to shallow cumulus cloud (Sc–Cu) transition, one of the key cloud regimes over the subtropical ocean. In this study, the transition is triggered by increasing the sea surface temperature (SST) and the LES runs until a statistically steady evolution of the Sc–Cu transition is achieved. The flow statistics are compared with those from a recycling-type method and it is found that the fringe method is more suitable for the current applications.


2021 ◽  
Vol 14 (9) ◽  
pp. 5435-5465
Author(s):  
Eckhard Kadasch ◽  
Matthias Sühring ◽  
Tobias Gronemeier ◽  
Siegfried Raasch

Abstract. In this paper, we present a newly developed mesoscale nesting interface for the PALM model system 6.0, which enables PALM to simulate the atmospheric boundary layer under spatially heterogeneous and non-stationary synoptic conditions. The implemented nesting interface, which is currently tailored to the mesoscale model COSMO, consists of two major parts: (i) the preprocessor INIFOR (initialization and forcing), which provides initial and time-dependent boundary conditions from mesoscale model output, and (ii) PALM's internal routines for reading the provided forcing data and superimposing synthetic turbulence to accelerate the transition to a fully developed turbulent atmospheric boundary layer. We describe in detail the conversion between the sets of prognostic variables, transformations between model coordinate systems, as well as data interpolation onto PALM's grid, which are carried out by INIFOR. Furthermore, we describe PALM's internal usage of the provided forcing data, which, besides the temporal interpolation of boundary conditions and removal of any residual divergence, includes the generation of stability-dependent synthetic turbulence at the inflow boundaries in order to accelerate the transition from the turbulence-free mesoscale solution to a resolved turbulent flow. We demonstrate and evaluate the nesting interface by means of a semi-idealized benchmark case. We carried out a large-eddy simulation (LES) of an evolving convective boundary layer on a clear-sky spring day. Besides verifying that changes in the inflow conditions enter into and successively propagate through the PALM domain, we focus our analysis on the effectiveness of the synthetic turbulence generation. By analysing various turbulence statistics, we show that the inflow in the present case is fully adjusted after having propagated for about two to three eddy-turnover times downstream, which corresponds well to other state-of-the-art methods for turbulence generation. Furthermore, we observe that numerical artefacts in the form of grid-scale convective structures in the mesoscale model enter the PALM domain, biasing the location of the turbulent up- and downdrafts in the LES. With these findings presented, we aim to verify the mesoscale nesting approach implemented in PALM, point out specific shortcomings, and build a baseline for future improvements and developments.


2021 ◽  
Author(s):  
Benedikt Seitzer ◽  
Bernd Leitl ◽  
Frank Harms

<p>Large-eddy simulations are increasingly used for studying the atmospheric boundary layer. With increasing computational resources even obstacle-resolving Large-eddy simulations became possible and will be used in urban climate studies more frequently. In these applications, grid sizes are in the order of a few meters. Whereas major urban structures can be resolved in general, details like aerodynamically rough surface structures can not be resolved explicitly. Based on the original fields of application, boundary conditions in Large-eddy simulations were initially formulated for surfaces of homogeneous roughness and for wall-distances much larger than the roughness sublayer height (Hultmark et al., 2013). The height of the roughness sublayer depends on the size of small-scale obstacles present on the surface exposed to the flow (Raupach et al., 1991). Typically, boundary conditions are evaluated between the surface and the first grid level. Thus, grid resolution in obstacle-resolved Large-Eddy simulations should also be a question of scales and therefore has to be chosen carefully (Basu and Lacser, 2017; Maronga et al., 2020). <br />In several wind tunnel experiments presented here, we measured the near-wall influence of differently scaled and shaped objects on a flow and its turbulence characteristics. Experimental setups were replicated numerically using the PALM model (Maronga et al. 2019). In a first, more generic experiment, the flow over horizontally homogeneous surfaces of different roughness was investigated. In a second experiment, the spatial separation of the turbulence scales was investigated in a more complex flow case. These experiments lead to considerations on model grid sizes in urban type Large-eddy simulations. The limitations of interpreting simulation results within the urban canopy layer are highlighted. There is an urgent need to reconsider how near-wall results of urban large-eddy simulations are generated and interpreted in the context of practical applications like flow and transport modelling in urban canopies. <br /><br /><em><strong>References</strong></em><br /><em>Basu, S. and Lacser, A. (2017). A Cautionary Note on the Use of Monin–Obukhov Similarity Theory in Very High-Resolution Large-Eddy Simulations. Boundary-Layer Meteorol, 163(2):351–355.</em></p> <p><em>Hultmark, M., Calaf, M., and Parlange, M. B. (2013). A new wall shearstress model for atmospheric boundary layer simulations. J Atmos Sci,70(11):3460–3470.</em></p> <p><em>Maronga, B., et al. (2020). Overview of the PALM model system 6.0. Geosci Model Dev Discussions, 06(June):1–63.</em></p> <p><em>Maronga, B., Knigge, C., and Raasch, S. (2020). An Improved Surface Boundary Condition for Large-Eddy Simulations Based on Monin–Obukhov Similarity Theory: Evaluation and Consequences forGrid Convergence in Neutral and Stable Conditions. Boundary-Layer Meteorol, 174(2):297–325.</em></p> <p><em>Raupach, M. R., Antonia, R. A., and Rajagopalan, S. (1991). Rough-wall turbulent boundary layers. Appl Mech Rev, 44(1):1–25</em></p>


2012 ◽  
Vol 204-208 ◽  
pp. 4490-4494
Author(s):  
Xiang Long Yang ◽  
Zhong Wei Huang ◽  
Lei Yang

The Influence of six kinds of combinations of different boundary conditions on the maintenance of homogeneity of atmospheric boundary layer in computational wind engineering was investigated. The inlet condition was applied by either prescribing velocity and turbulent quantities or using a periodic boundary. The top boundary condition includes symmetry, prescribing velocity and turbulent quantities, constant static pressure and applying driving shear stress. Numerical results show that the stream-wise velocity and the turbulent dissipation rate are almost immune to boundary conditions, but the turbulent kinetic energy is affected by boundary conditions dramatically. Best maintenance can be obtained by using periodic boundary condition at inlet and outlet and applying driving shear stress on the top of the domain.


2016 ◽  
Vol 820 ◽  
pp. 351-358
Author(s):  
Andrea Vranková ◽  
Milan Palko

The aim of the paper was to define the input options over the boundary layer, as the entrance boundary conditions for simulation in ANSYS. The boundary layer is designed for use in external aerodynamics of buildings (part of the urban structure) for selected sites occurring in the territory of the Slovak Republic.


Sign in / Sign up

Export Citation Format

Share Document