scholarly journals Anabatic Flow along a Uniformly Heated Slope Studied through Large-Eddy Simulation

Atmosphere ◽  
2021 ◽  
Vol 12 (7) ◽  
pp. 850
Author(s):  
Carlo Cintolesi ◽  
Dario Di Santo ◽  
Francesco Barbano ◽  
Silvana Di Sabatino

Anabatic flows are common phenomena in the presence of sloping terrains, which significantly affect the dynamics and the exchange of mass and momentum in the low-atmosphere. Despite this, very few studies in the literature have tackled this topic. The present contribution addresses this gap by utilising high-resolved large-eddy simulations for investigating an anabatic flow in a simplified configuration, commonly used in laboratory experiments. The purpose is to analyse the complex thermo-fluid dynamics and the turbulent structures arising from the anabatic flow near the slope. In such a flow, three main dynamic layers are identified and reported: the conductive layer close to the surface, the convective layer where the most energetic motion develops, and the outer region, which is almost unperturbed. The analysis of instantaneous fields reveals the presence of thermal plumes, which are stable turbulent structures enhancing vertical transport and mixing of momentum and temperature. Such structures are generated by thermal instabilities in the conductive layer that trigger the rise of the plumes above them. Their evolution along the slope is described, identifying three regions responsible for the plumes generation, stabilisation, and merging. To the best of the authors’ knowledge, this is the first numerical experiment describing the along-slope behaviour of the thermal plumes in the convective layer.

Energies ◽  
2021 ◽  
Vol 14 (11) ◽  
pp. 3078
Author(s):  
Carlo Cintolesi ◽  
Francesco Barbano ◽  
Silvana Di Sabatino

Thermal convective flows are common phenomena in real urban canyons and strongly affect the mechanisms of pollutant removal from the canyon. The present contribution aims at investigating the complex interaction between inertial and thermal forces within the canyon, including the impacts on turbulent features and pollutant removal mechanisms. Large-eddy simulations reproduce infinitely long square canyons having isothermal and differently heated facades. A scalar source on the street mimics the pollutant released by traffic. The presence of heated facades triggers convective flows which generate an interaction region around the canyon-ambient interface, characterised by highly energetic turbulent fluxes and an increase of momentum and mass exchange. The presence of this region of high mixing facilitates the pollutant removal across the interface and decreases the urban canopy drag. The heating-up of upwind facade determines favourable convection that strengthens the primary internal vortex and decreases the pollutant concentration of the whole canyon by 49% compare to the isothermal case. The heating-up of the downwind facade produces adverse convection counteracting the wind-induced motion. Consequently, the primary vortex is less energetic and confined in the upper-canyon area, while a region of almost zero velocity and high pollution concentration (40% more than the isothermal case) appears at the pedestrian level. Finally, numerical analyses allow a definition of a local Richardson number based on in-canyon quantities only and a new formulation is proposed to characterise the thermo-dynamics regimes.


2016 ◽  
Vol 126 ◽  
pp. 375-381 ◽  
Author(s):  
Zun Cai ◽  
Xiao Liu ◽  
Cheng Gong ◽  
Mingbo Sun ◽  
Zhenguo Wang ◽  
...  

Author(s):  
Matthew J. Churchfield ◽  
Ye Li ◽  
Patrick J. Moriarty

This paper presents our initial work in performing large-eddy simulations of tidal turbine array flows. First, a horizontally periodic precursor simulation is performed to create turbulent flow data. Then those data are used as inflow into a tidal turbine array two rows deep and infinitely wide. The turbines are modelled using rotating actuator lines, and the finite-volume method is used to solve the governing equations. In studying the wakes created by the turbines, we observed that the vertical shear of the inflow combined with wake rotation causes lateral wake asymmetry. Also, various turbine configurations are simulated, and the total power production relative to isolated turbines is examined. We found that staggering consecutive rows of turbines in the simulated configurations allows the greatest efficiency using the least downstream row spacing. Counter-rotating consecutive downstream turbines in a non-staggered array shows a small benefit. This work has identified areas for improvement. For example, using a larger precursor domain would better capture elongated turbulent structures, and including salinity and temperature equations would account for density stratification and its effect on turbulence. Additionally, the wall shear stress modelling could be improved, and more array configurations could be examined.


2015 ◽  
Vol 17 (1) ◽  
pp. 245-255 ◽  
Author(s):  
Matteo Colli ◽  
Luca G. Lanza ◽  
Roy Rasmussen ◽  
Julie M. Thériault

Abstract The use of windshields to reduce the impact of wind on snow measurements is common. This paper investigates the catching performance of shielded and unshielded gauges using numerical simulations. In Part II, the role of the windshield and gauge aerodynamics, as well as the varying flow field due to the turbulence generated by the shield–gauge configuration, in reducing the catch efficiency is investigated. This builds on the computational fluid dynamics results obtained in Part I, where the airflow patterns in the proximity of an unshielded and single Alter shielded Geonor T-200B gauge are obtained using both time-independent [Reynolds-averaged Navier–Stokes (RANS)] and time-dependent [large-eddy simulation (LES)] approaches. A Lagrangian trajectory model is used to track different types of snowflakes (wet and dry snow) and to assess the variation of the resulting gauge catching performance with the wind speed. The collection efficiency obtained with the LES approach is generally lower than the one obtained with the RANS approach. This is because of the impact of the LES-resolved turbulence above the gauge orifice rim. The comparison between the collection efficiency values obtained in case of shielded and unshielded gauge validates the choice of installing a single Alter shield in a windy environment. However, time-dependent simulations show that the propagating turbulent structures produced by the aerodynamic response of the upwind single Alter blades have an impact on the collection efficiency. Comparison with field observations provides the validation background for the model results.


2020 ◽  
Author(s):  
Ji Sung Na ◽  
Taekyun Kim ◽  
Emilia Kyung Jin ◽  
Seung-Tae Yoon ◽  
Won Sang Lee ◽  
...  

Abstract. The role of the refreezing effect in the ice shelf–ocean boundary layer (IOBL) flow with a super-cooled, plume beneath the ice shelf is investigated using the large-eddy simulation. To reveal the detailed physical processes and characteristics of the IOBL flow, a model domain is initialized and forced by in situ observations and a comparison is made between two simulations, one with the refreezing effect and one without. The simulated velocity, potential temperature, and salinity field are validated with in situ observations performed in Terra Nova Bay in the western Ross Sea in 2016/2017, confirming that the vertical structures in the simulation results agree well with observations. In particular, it is evident that, when the refreezing effect is considered, the IOBL flow can be more realistically resolved, especially upward advection from the sub-ice shelf plume and the ice front eddy. Beneath the ice shelf, two district regions (the inner and outer regions) are identified based on flow characteristics and the refreezing pattern. In the inner region, stratification and stable conditions are observed with negative momentum flux and low refreezing rates. Meanwhile, in the outer region, high shear impact and unstable conditions with a heat flux of −9 to −52 W m−2 are observed, demonstrating the high refreezing rate and the entrainment of super-cooled water from the sub-ice shelf plume. A total of 94 % of the refreezing events occur in the outer region, with a maximum refreezing rate of 1.86 m yr−1 at the ice front.


2013 ◽  
Vol 70 (6) ◽  
pp. 1513-1527 ◽  
Author(s):  
Jing Huang ◽  
Elie Bou-Zeid

Abstract This study seeks to quantitatively and qualitatively understand how stability affects transport in the continuously turbulent stably stratified atmospheric boundary layer, based on a suite of large-eddy simulations. The test cases are based on the one adopted by the Global Energy and Water Cycle Experiment (GEWEX) Atmospheric Boundary Layer Study (GABLS) project, but with a largely expanded stability range where the gradient Richardson number (Rig) reaches up to around 1. The analysis is mainly focused on understanding the modification of turbulent structures and dynamics with increasing stability in order to improve the modeling of the stable atmospheric boundary layer in weather and climate models, a topic addressed in Part II of this work. It is found that at quasi equilibrium, an increase in stability results in stronger vertical gradients of the mean temperature, a lowered low-level jet, a decrease in vertical momentum transport, an increase in vertical buoyancy flux, and a shallower boundary layer. Analysis of coherent turbulent structures using two-point autocorrelation reveals that the autocorrelation of the streamwise velocity is horizontally anisotropic while the autocorrelation of the vertical velocity is relatively isotropic in the horizontal plane and its integral length scale decreases as stability increases. The effects of stability on the overall turbulent kinetic energy (TKE) and its budget terms are also investigated, and it is shown that the authors' large-eddy simulation results are in good agreement with previous experimental findings across varied stabilities. Finally, Nieuwstadt's local-scaling theory is reexamined and it is concluded that the height z is not a relevant scaling parameter and should be replaced by a constant length scale away from the surface, indicating that the z-less range starts lower than previously assumed.


2002 ◽  
Vol 473 ◽  
pp. 83-101 ◽  
Author(s):  
DOUGLAS G. DOMMERMUTH ◽  
JAMES W. ROTTMAN ◽  
GEORGE E. INNIS ◽  
EVGENY A. NOVIKOV

We present some preliminary results from using large-eddy simulation to compute the late wake of a sphere towed at constant speed through a non-stratified and a uniformly stratified fluid. The wake is computed in each case for two values of the Reynolds number: Re = 104, which is comparable to that used in laboratory experiments, and Re = 105. An important aspect of the simulation is the use of an iterative procedure to relax the initial turbulence field so that the normal and shear turbulent stresses are properly correlated and the turbulent production and dissipation are in equilibrium. For the lower Reynolds number our results compare well with existing laboratory experimental results. For the higher Reynolds number we find that even though the turbulence is more developed and the wake contains finer structure, most of the similarity properties of the wake are unchanged compared with those observed at the lower Reynolds number.


Sign in / Sign up

Export Citation Format

Share Document