A Total Turbulent Energy Closure Model for Neutrally and Stably Stratified Atmospheric Boundary Layers

2007 ◽  
Vol 64 (11) ◽  
pp. 4113-4126 ◽  
Author(s):  
Thorsten Mauritsen ◽  
Gunilla Svensson ◽  
Sergej S. Zilitinkevich ◽  
Igor Esau ◽  
Leif Enger ◽  
...  
Keyword(s):  
Potential Temperature ◽  
Stability Limit ◽  
Turbulent Energy ◽  
Large Eddy Simulations ◽  
Turbulence Closure ◽  
Turbulent Heat ◽  
Atmospheric Conditions ◽  
Mean Flow ◽  
Wide Range ◽  
Large Eddy

Abstract This paper presents a turbulence closure for neutral and stratified atmospheric conditions. The closure is based on the concept of the total turbulent energy. The total turbulent energy is the sum of the turbulent kinetic energy and turbulent potential energy, which is proportional to the potential temperature variance. The closure uses recent observational findings to take into account the mean flow stability. These observations indicate that turbulent transfer of heat and momentum behaves differently under very stable stratification. Whereas the turbulent heat flux tends toward zero beyond a certain stability limit, the turbulent stress stays finite. The suggested scheme avoids the problem of self-correlation. The latter is an improvement over the widely used Monin–Obukhov-based closures. Numerous large-eddy simulations, including a wide range of neutral and stably stratified cases, are used to estimate likely values of two free constants. In a benchmark case the new turbulence closure performs indistinguishably from independent large-eddy simulations.

Simulation study of wake encounters with straight and deformed vortices

2016 ◽  
Vol 120 (1226) ◽  
pp. 651-674 ◽  
Author(s):  
D. Vechtel
Keyword(s):  
Simulation Study ◽  
Thermal Stratification ◽  
Parameter Variation ◽  
Decay Process ◽  
Large Eddy Simulations ◽  
Vortex Rings ◽  
Atmospheric Conditions ◽  
Wide Range ◽  
Large Eddy ◽  
The Impact

ABSTRACTA simulation study was conducted in order to investigate the influence of vortex deformation on wake encounter characteristics. Wake vortices tend to be strongly deformed during the decay process, depending on the atmospheric conditions in terms of turbulence and thermal stratification. For quantification of the influence of vortex deformation, encounters of an aircraft of the ‘Medium’ category behind a generator aircraft of the ‘Heavy’ category were simulated with straight vortices and with realistically deformed vortices derived from large-eddy simulations. All relevant parameters that influence the encounter characteristics, such as encounter angles and positions, were varied within a wide range. In order to cover all kinds of vortex deformation, encounters with different vortex ages from 16-136 seconds were simulated. Hence, all relevant phases during the vortex decay from nearly straight and wavy vortices to vortex rings were considered.The parameter variation study revealed that on average the impact on the encountering aircraft is less with deformed vortices than with straight vortices of comparable strength. Especially with vortex rings, the encountering aircraft is exposed to a much smaller impact. However, the results also show a larger aircraft response during encounters with wavy vortices just prior to vortex linking. The maximum aircraft response with wavy vortices is stronger than with straight vortices of comparable strength. Also, the strongest encounters occur under greater encounter angles with deformed vortices than with straight ones.

Investigation on hydrodynamic forces leading to coarse particle entrainment using Large-Eddy Simulations

2021 ◽  
Author(s):  
Gaston Latessa ◽  
Angela Busse ◽  
Manousos Valyrakis
Keyword(s):  
Experimental Data ◽  
Flow Field ◽  
Large Scale ◽  
Large Eddy Simulations ◽  
Flow Conditions ◽  
Mean Flow ◽  
Protection Measures ◽  
Practical Applications ◽  
Particle Entrainment ◽  
Large Eddy

<p>The prediction of particle motion in a fluid flow environment presents several challenges from the quantification of the forces exerted by the fluid onto the solids -normally with fluctuating behaviour due to turbulence- and the definition of the potential particle entrainment from these actions. An accurate description of these phenomena has many practical applications in local scour definition and to the design of protection measures.</p><p>In the present work, the actions of different flow conditions on sediment particles is investigated with the aim to translate these effects into particle entrainment identification through analytical solid dynamic equations.</p><p>Large Eddy Simulations (LES) are an increasingly practical tool that provide an accurate representation of both the mean flow field and the large-scale turbulent fluctuations. For the present case, the forces exerted by the flow are integrated over the surface of a stationary particle in the streamwise (drag) and vertical (lift) directions, together with the torques around the particle’s centre of mass. These forces are validated against experimental data under the same bed and flow conditions.</p><p>The forces are then compared against threshold values, obtained through theoretical equations of simple motions such as rolling without sliding. Thus, the frequency of entrainment is related to the different flow conditions in good agreement with results from experimental sediment entrainment research.</p><p>A thorough monitoring of the velocity flow field on several locations is carried out to determine the relationships between velocity time series at several locations around the particle and the forces acting on its surface. These results a relevant to determine ideal locations for flow investigation both in numerical and physical experiments.</p><p>Through numerical experiments, a large number of flow conditions were simulated obtaining a full set of actions over a fixed particle sitting on a smooth bed. These actions were translated into potential particle entrainment events and validated against experimental data. Future work will present the coupling of these LES models with Discrete Element Method (DEM) models to verify the entrainment phenomena entirely from a numerical perspective.</p>

scholarly journals Stability and Sensitivity Analysis of Hydrodynamic Instabilities in Industrial Swirled Injection Systems

2017 ◽  
Author(s):  
Thomas L. Kaiser ◽  
Thierry Poinsot ◽  
Kilian Oberleithner
Keyword(s):  
Stability Analysis ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Vortex Core ◽  
Large Eddy Simulations ◽  
Mode Shapes ◽  
Mean Flow ◽  
Precessing Vortex Core ◽  
Large Eddy ◽  
Good Agreement

The hydrodynamic instability in an industrial, two-staged, counter-rotative, swirled injector of highly complex geometry is under investigation. Large eddy simulations show that the complicated and strongly nonparallel flow field in the injector is superimposed by a strong precessing vortex core. Mean flow fields of large eddy simulations, validated by experimental particle image velocimetry measurements are used as input for both local and global linear stability analysis. It is shown that the origin of the instability is located at the exit plane of the primary injector. Mode shapes of both global and local linear stability analysis are compared to a dynamic mode decomposition based on large eddy simulation snapshots, showing good agreement. The estimated frequencies for the instability are in good agreement with both the experiment and the simulation. Furthermore, the adjoint mode shapes retrieved by the global approach are used to find the best location for periodic forcing in order to control the precessing vortex core.

scholarly journals Flat vs. curved rigid-lid LES computations of an open-channel confluence

2019 ◽  
Vol 21 (2) ◽  
pp. 318-334 ◽  
Author(s):  
Pedro Xavier Ramos ◽  
Laurent Schindfessel ◽  
João Pedro Pêgo ◽  
Tom De Mulder
Keyword(s):  
Free Surface ◽  
Secondary Flow ◽  
Open Channel ◽  
Pressure Field ◽  
Large Eddy Simulations ◽  
Numerical Treatment ◽  
Mean Flow ◽  
Large Eddy ◽  
Flow Case ◽  
Channel Confluence

Abstract This paper describes the application of four Large Eddy Simulations (LES) to an open-channel confluence flow, making use of a frictionless rigid-lid to treat the free-surface. Three simulations are conducted with a flat rigid-lid, at different elevations. A fourth simulation is carried out with a curved rigid-lid which is a closer approximation to the real free-surface of the flow. The curved rigid-lid is obtained from the time-averaged pressure field on the flat rigid-lid from one of the initial three simulations. The aim is to investigate the limitations of the free-surface treatment by means of a rigid-lid in the simulation of an asymmetric confluence, showing the differences that both approaches produce in terms of mean flow, secondary flow and turbulence. After validation with experimental data, the predictions are used to understand the differences between adopting a flat and a curved rigid-lid onto the confluence hydrodynamics. For the present flow case, although it was characterized by a moderately low downstream Froude number (Fr ≈ 0.37), it was found that an oversimplification of the numerical treatment of the free-surface leads to a decreased accuracy of the predictions of the secondary flow and turbulent kinetic energy.

scholarly journals The Iowa Atmospheric Observatory: Revealing the Unique Boundary Layer Characteristics of a Wind Farm

2019 ◽  
Vol 23 (2) ◽  
pp. 1-27 ◽  
Author(s):  
Eugene S. Takle ◽  
Daniel A. Rajewski ◽  
Samantha L. Purdy
Keyword(s):  
Boundary Layer ◽  
Wind Farm ◽  
Agricultural Landscape ◽  
Potential Temperature ◽  
Wind Farms ◽  
Vertical Profiles ◽  
Atmospheric Conditions ◽  
Wide Range ◽  
Water Vapor Mixing Ratio ◽  
Utility Scale

Abstract The Iowa Atmospheric Observatory was established to better understand the unique microclimate characteristics of a wind farm. The facility consists of a pair of 120-m towers identically instrumented to observe basic landscape–atmosphere interactions in a highly managed agricultural landscape. The towers, one within and one outside of a utility-scale low-density-array wind farm, are equipped to measure vertical profiles of temperature, wind, moisture, and pressure and can host specialized sensors for a wide range of environmental conditions. Tower measurements during the 2016 growing season demonstrate the ability to distinguish microclimate differences created by single or multiple turbines from natural conditions over homogeneous agricultural fields. Microclimate differences between the two towers are reported as contrasts in normalized wind speed, normalized turbulence intensity, potential temperature, and water vapor mixing ratio. Differences are analyzed according to conditions of no wind farm influence (i.e., no wake) versus wind farm influence (i.e., waked flow) with distance downwind from a single wind turbine or a large group of turbines. Differences are also determined for more specific atmospheric conditions according to thermal stratification. Results demonstrate agreement with most, but not all, currently available numerical flow-field simulations of large wind farm arrays and of individual turbines. In particular, the well-documented higher nighttime surface temperature in wind farms is examined in vertical profiles that confirm this effect to be a “suppression of cooling” rather than a warming process. A summary is provided of how the wind farm boundary layer differs from the natural boundary layer derived from concurrent measurements over the summer of 2016.

scholarly journals Three-Dimensional Turbulent Vortex Shedding From a Surface-Mounted Square Cylinder: Predictions With Large-Eddy Simulations and URANS

2014 ◽  
Vol 136 (6) ◽  
Author(s):  
B. A. Younis ◽  
A. Abrishamchi
Keyword(s):  
Experimental Data ◽  
Turbulence Model ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Large Eddy Simulations ◽  
Navier Stokes ◽  
Reynolds Stresses ◽  
Mean Flow ◽  
Navier Stokes Equations ◽  
Large Eddy

The paper reports on the prediction of the turbulent flow field around a three-dimensional, surface mounted, square-sectioned cylinder at Reynolds numbers in the range 104–105. The effects of turbulence are accounted for in two different ways: by performing large-eddy simulations (LES) with a Smagorinsky model for the subgrid-scale motions and by solving the unsteady form of the Reynolds-averaged Navier–Stokes equations (URANS) together with a turbulence model to determine the resulting Reynolds stresses. The turbulence model used is a two-equation, eddy-viscosity closure that incorporates a term designed to account for the interactions between the organized mean-flow periodicity and the random turbulent motions. Comparisons with experimental data show that the two approaches yield results that are generally comparable and in good accord with the experimental data. The main conclusion of this work is that the URANS approach, which is considerably less demanding in terms of computer resources than LES, can reliably be used for the prediction of unsteady separated flows provided that the effects of organized mean-flow unsteadiness on the turbulence are properly accounted for in the turbulence model.

Large-Eddy Simulations of Wall Bounded Turbulent Flows Using Unstructured Linear Reconstruction Techniques

2015 ◽  
Vol 137 (5) ◽  
Author(s):  
Dario Amirante ◽  
Nicholas J. Hills
Keyword(s):  
Reynolds Number ◽  
Turbulent Flows ◽  
Outer Radius ◽  
Large Eddy Simulations ◽  
Test Cases ◽  
Mean Flow ◽  
Flow Equations ◽  
Turbulent Statistics ◽  
Large Eddy ◽  
Least Squares Approach

Large-eddy simulations (LES) of wall bounded, low Mach number turbulent flows are conducted using an unstructured finite-volume solver of the compressible flow equations. The numerical method employs linear reconstructions of the primitive variables based on the least-squares approach of Barth. The standard Smagorinsky model is adopted as the subgrid term. The artificial viscosity inherent to the spatial discretization is maintained as low as possible reducing the dissipative contribution embedded in the approximate Riemann solver to the minimum necessary. Comparisons are also discussed with the results obtained using the implicit LES (ILES) procedure. Two canonical test-cases are described: a fully developed pipe flow at a bulk Reynolds number Reb = 44 × 103 based on the pipe diameter, and a confined rotor–stator flow at the rotational Reynolds number ReΩ = 4 × 105 based on the outer radius. In both cases, the mean flow and the turbulent statistics agree well with existing direct numerical simulations (DNS) or experimental data.

Large Eddy Simulations of Flow and Heat Transfer in Rotating Ribbed Duct Flows

2005 ◽  
Vol 127 (5) ◽  
pp. 486-498 ◽  
Author(s):  
Mayank Tyagi ◽  
Sumanta Acharya
Keyword(s):  
Heat Transfer ◽  
Temperature Field ◽  
Large Scale ◽  
Density Ratio ◽  
Low Frequency ◽  
Turbulent Energy ◽  
Large Eddy Simulations ◽  
Temperature Fields ◽  
Periodic Domain ◽  
Large Eddy

Large eddy simulations are performed in a periodic domain of a rotating square duct with normal rib turbulators. Both the Coriolis force as well as the centrifugal buoyancy forces are included in this study. A direct approach is presented for the unsteady calculation of the nondimensional temperature field in the periodic domain. The calculations are performed at a Reynolds number (Re) of 12,500, a rotation number (Ro) of 0.12, and an inlet coolant-to-wall density ratio Δρ/ρ of 0.13. The predicted time and space-averaged Nusselt numbers are shown to compare satisfactorily with the published experimental data. Time sequences of the vorticity components and the temperature fields are presented to understand the flow physics and the unsteady heat transfer behavior. Large scale coherent structures are seen to play an important role in the mixing and heat transfer. The temperature field appears to contain a low frequency mode that extends beyond a single inter-rib geometric module, and indicates the necessity of using at least two inter-rib modules for streamwise periodicity to be satisfied. Proper orthogonal decomposition (POD) of the flowfield indicates a low dimensionality of this system with almost 99% of turbulent energy in the first 80 POD modes.

Responses of Shallow Cumulus Convection to Large-Scale Temperature and Moisture Perturbations: A Comparison of Large-Eddy Simulations and a Convective Parameterization Based on Stochastically Entraining Parcels

2012 ◽  
Vol 69 (6) ◽  
pp. 1936-1956 ◽  
Author(s):  
Ji Nie ◽  
Zhiming Kuang
Keyword(s):  
Large Scale ◽  
Potential Temperature ◽  
Large Eddy Simulations ◽  
Cumulus Parameterization ◽  
Temperature Perturbation ◽  
Cumulus Convection ◽  
Positive Temperature ◽  
Shallow Cumulus ◽  
Scale Temperature ◽  
Large Eddy

Abstract Responses of shallow cumuli to large-scale temperature/moisture perturbations are examined through diagnostics of large-eddy simulations (LESs) of the undisturbed Barbados Oceanographic and Meteorological Experiment (BOMEX) case and a stochastic parcel model. The perturbations are added instantaneously and allowed to evolve freely afterward. The parcel model reproduces most of the changes in the LES-simulated cloudy updraft statistics in response to the perturbations. Analyses of parcel histories show that a positive temperature perturbation forms a buoyancy barrier, which preferentially eliminates parcels that start with lower equivalent potential temperature or have experienced heavy entrainment. Besides the amount of entrainment, the height at which parcels entrain is also important in determining their fate. Parcels entraining at higher altitudes are more likely to overcome the buoyancy barrier than those entraining at lower altitudes. Stochastic entrainment is key for the parcel model to reproduce the LES results. Responses to environmental moisture perturbations are quite small compared to those to temperature perturbations because changing environmental moisture is ineffective in modifying buoyancy in the BOMEX shallow cumulus case. The second part of the paper further explores the feasibility of a stochastic parcel–based cumulus parameterization. Air parcels are released from the surface layer and temperature/moisture fluxes effected by the parcels are used to calculate heating/moistening tendencies due to both cumulus convection and boundary layer turbulence. Initial results show that this conceptually simple parameterization produces realistic convective tendencies and also reproduces the LES-simulated mean and variance of cloudy updraft properties, as well as the response of convection to temperature/moisture perturbations.

Large Eddy Simulations of Flow and Heat Transfer in Rotating Ribbed Duct Flows

2004 ◽  
Author(s):  
Mayank Tyagi ◽  
Sumanta Acharya
Keyword(s):  
Heat Transfer ◽  
Temperature Field ◽  
Large Scale ◽  
Density Ratio ◽  
Low Frequency ◽  
Turbulent Energy ◽  
Large Eddy Simulations ◽  
Temperature Fields ◽  
Periodic Domain ◽  
Large Eddy

Large eddy simulations are performed in a periodic domain of a rotating square duct with normal rib turbulators. Both the Coriolis force as well as the centrifugal buoyancy force are included in this study. A direct approach is presented for the unsteady calculation of the non-dimensional temperature field in the periodic domain. The calculations are performed at a Reynolds number (Re) of 12, 500, a Rotation number (Ro) of 0.12 and an inlet coolant-to-wall density ratio (Δρ/ρ) of 0.13. The time-averaged Nusselt numbers compare satisfactorily with the data of Wagner et al. (J. Turbomachinery, Vol. 114, pp. 847–857). Time-sequences of the vorticity components and the temperature fields are presented to understand the flow physics and the unsteady heat transfer processes. Large scale coherent structures are seen to play an important role in the mixing and heat transfer. The temperature field appears to contain a low frequency mode that extends beyond a single inter-rib geometric module, and indicates the necessity of using at least two inter-rib modules for streamwise periodicity to be satisfied. Proper orthogonal decomposition (POD) of 200 snapshots indicates a low dimensionality of this system with almost 99% of turbulent energy in the first 80 POD modes.

Sign in / Sign up

Export Citation Format

Share Document