scholarly journals Near-Surface Intensification of Tornado Vortices

2007 ◽  
Vol 64 (7) ◽  
pp. 2176-2194 ◽  
Author(s):  
D. C. Lewellen ◽  
W. S. Lewellen

Abstract An idealized analytical model and numerical large-eddy simulations are used to explore fluid-dynamic mechanisms by which tornadoes may be intensified near the surface relative to conditions aloft. The analytical model generalizes a simple model of Barcilon and Fiedler and Rotunno for a steady supercritical end-wall vortex to more general vortex corner flows, angular momentum distributions, and time dependence. The model illustrates the role played by the corner flow swirl ratio in determining corner flow structure and intensification; predicts an intensification of near-surface swirl velocities relative to conditions aloft of Iυ ∼ 2 for supercritical end-wall vortices in agreement with earlier analytical, numerical, and laboratory results; and suggests how larger intensification factors might be achieved in some more general corner flows. Examples of the latter are presented using large-eddy simulations. By tuning the lateral inflow boundary conditions near the surface, quasi-steady vortices exhibiting nested inner and outer corner flows and Iυ ∼ 4 are produced. More significantly, these features can be produced without fine tuning, along with an additional doubling (or more) of the intensification, in a broad class of unsteady evolutions producing a dynamic corner flow collapse. These scenarios, triggered purely by changes in the far-field near-surface flow, provide an attractive mechanism for naturally achieving an intense near-surface vortex from a much larger-scale less-intense swirling flow. It is argued that, applied on different scales, this may sometimes play a role in tornadogenesis and/or tornado variability. This phenomenon of corner flow collapse is considered further in a companion paper.

Author(s):  
Daniel P. Stern ◽  
George H. Bryan ◽  
Chia-Ying Lee ◽  
James D. Doyle

AbstractRecent studies have shown that extreme wind gusts are ubiquitous within the eyewall of intense tropical cyclones (TCs). These gusts pose a substantial hazard to human life and property, but both the short-term (i.e., during the passage of a single TC) and long-term (over many years) risk of encountering such a gust at a given location is poorly understood. Here, simulated tower data from large-eddy simulations of idealized TCs in a quiescent (i.e., no mean flow or vertical wind shear) environment are used to estimate these risks for the offshore region of the United States. For both a category 5 and category 3 TC, there is a radial region where nearly all simulated towers experience near-surface (the lowest 200 m) 3-s gusts exceeding 70 m s−1 within a 10-minute period; on average, these towers respectively sample peak 3-s gusts of 110 and 80 m s−1. Analysis of an observational dropsonde database supports the idealized simulations, and indicates that offshore structures (such as wind turbines) in the eyewall of a major hurricane are likely to encounter damaging wind speeds. This result is then incorporated into an estimate of the long-term risk, using analyses of the return period for major hurricanes from both a best-track database and a statistical-dynamical model forced by reanalysis. For much of the nearshore region of the Gulf of Mexico and southeastern US coasts, this analysis yields an estimate of a 30-60% probability of any given point experiencing at least one 70 m s−1 gust within a 30-year period.


2009 ◽  
Vol 48 (6) ◽  
pp. 1161-1180 ◽  
Author(s):  
Francis L. Ludwig ◽  
Fotini Katopodes Chow ◽  
Robert L. Street

Abstract This paper demonstrates the importance of high-quality subfilter-scale turbulence models in large-eddy simulations by evaluating the resolved-scale flow features that result from various closure models. The Advanced Regional Prediction System (ARPS) model was used to simulate neutral flow over a 1.2-km square, flat, rough surface with seven subfilter turbulence models [Smagorinsky, turbulent kinetic energy (TKE)-1.5, and five dynamic reconstruction combinations]. These turbulence models were previously compared with similarity theory. Here, the differences are evaluated using mean velocity statistics and the spatial structure of the flow field. Streamwise velocity averages generally differ among models by less than 0.5 m s−1, but those differences are often significant at a 95% confidence level. Flow features vary considerably among models. As measured by spatial correlation, resolved flow features grow larger and less elongated with height for a given model and resolution. The largest differences are between dynamic models that allow energy backscatter from small to large scales and the simple eddy-viscosity closures. At low altitudes, the linear extent of Smagorinsky and TKE-1.5 structures exceeds those of dynamic models, but the relationship reverses at higher altitudes. Ejection, sweep, and upward momentum flux features differ among models and from observed neutral atmospheric flows, especially for Smagorinsky and TKE-1.5 coarse-grid simulations. Near-surface isopleths separating upward fluxes from downward are shortest for the Smagorinsky and TKE-1.5 coarse-grid simulations, indicating less convoluted turbulent interfaces; at higher altitudes they are longest. Large-eddy simulation (LES) is a powerful simulation tool, but choices of grid resolution and subfilter model can affect results significantly. Physically realistic dynamic mixed models, such as those presented here, are essential when using LES to study atmospheric processes such as transport and dispersion—in particular at coarse resolutions.


2007 ◽  
Vol 124 (3) ◽  
pp. 405-424 ◽  
Author(s):  
Inanc Senocak ◽  
Andrew S. Ackerman ◽  
Michael P. Kirkpatrick ◽  
David E. Stevens ◽  
Nagi N. Mansour

1992 ◽  
Vol 242 ◽  
pp. 51-78 ◽  
Author(s):  
P. J. Mason ◽  
D. J. Thomson

The ability of a large-eddy simulation to represent the large-scale motions in the interior of a turbulent flow is well established. However, concerns remain for the behaviour close to rigid surfaces where, with the exception of low-Reynolds-number flows, the large-eddy description must be matched to some description of the flow in which all except the larger-scale ‘inactive’ motions are averaged. The performance of large-eddy simulations in this near-surface region is investigated and it is pointed out that in previous simulations the mean velocity profile in the matching region has not had a logarithmic form. A number of new simulations are conducted with the Smagorinsky (1963) subgrid model. These also show departures from the logarithmic profile and suggest that it may not be possible to eliminate the error by adjustments of the subgrid lengthscale. An obvious defect of the Smagorinsky model is its failure to represent stochastic subgrid stress variations. It is shown that inclusion of these variations leads to a marked improvement in the near-wall flow simulation. The constant of proportionality between the magnitude of the fluctuations in stress and the Smagorinsky stresses has been empirically determined to give an accurate logarithmic flow profile. This value provides an energy backscatter rate slightly larger than the dissipation rate and equal to idealized theoretical predictions (Chasnov 1991).


2017 ◽  
Vol 145 (8) ◽  
pp. 3025-3048 ◽  
Author(s):  
Jeffrey D. Mirocha ◽  
Katherine A. Lundquist

To facilitate multiscale simulation using the Weather Research and Forecasting Model, vertical mesh refinement for one-way concurrent nested simulation was recently introduced. Grid refinement in the vertical dimension removes the requirement of different grid aspect ratios on the bounding versus the nested domain, such that results from refinement are in the horizontal directions only, and thereby can also reduce numerical errors on the bounding domain arising from large aspect ratios in the presence of complex terrain. Herein, the impacts of vertical grid refinement on the evolving downstream flow in nested large-eddy simulations are evaluated in relation to other model configuration choices, including turbulence subfilter-scale (SFS) stress models, mesh configuration, and alternative methods for calculating several near-surface flow parameters. Although vertical nesting requires coarsening of the vertical grid on the bounding domain, leading to a smaller range of resolved turbulence scales in the nest’s lateral boundary conditions, parameter values within the nested domains are generally only minimally impacted, relative to nesting using the same vertical grid on each domain. Two dynamic SFS models examined herein generally improved the simulated mean wind speed, turbulence kinetic energy, stresses and spectra, on both domains, and accelerated equilibration rates within nested domains, relative to two constant coefficient models. A new method of extrapolating horizontal velocity components to near-surface locations at nested domain lateral boundaries, and a correction to the calculation of deformation elements near the surface, are each shown to slightly alter the mean parameter values, yet only minimally impact equilibration rates within the nested domain.


Author(s):  
Björn Maronga ◽  
Dan Li

AbstractWe revisit the longstanding problem of grid sensitivity, i.e., the lack of grid convergence in large-eddy simulations (LES) of the stable boundary layer. We use a comprehensive set of LES of the well-known Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study 1 (GABLS1) case with varying grid spacings between 12.5 m and 1 m to investigate several physical processes and numerical features that are possible causes of grid sensitivity. Our results demonstrate that there are two resolution regimes in which grid sensitivity manifests differently. We find that changing the numerical advection schemes and the subgrid-scale models alters the simulation results, but the options tested do not fully address the grid-sensitivity issue. Moreover, sensitivity runs suggest that the surface boundary condition and the interaction of the surface with the near-surface flow, as well as the mixing with the free atmosphere, are unlikely to be the causes of the observed grid sensitivity. One interesting finding is that the grid sensitivity in the fine grid-spacing regime (grid spacings $$\le 2\,\mathrm{m}$$ ≤ 2 m ) is closely related to the reduction in the energy content of large-scale turbulence, leading to less turbulence kinetic energy and hence lower boundary-layer heights. The present work demonstrates that there is still an urgent need to address this grid-sensitivity issue in order to perform reliable LES of the stable boundary layer.


Atmosphere ◽  
2020 ◽  
Vol 11 (4) ◽  
pp. 434 ◽  
Author(s):  
Seung-Bu Park ◽  
Jong-Jin Baik

The decay of the Convective Boundary Layer (CBL) is studied using large-eddy simulations of free and advective CBLs, in which surface heat supply is suddenly cut off. After the cutoff, coherent convective circulations last about one convective time scale and then fade away. In the mixed layer, the decay time scale increases with height, indicating that nonlocal eddies decay slower than near-surface local eddies. The slower decay of turbulence in the middle of CBL than near-surface turbulence is reconfirmed from the analysis of pattern correlations of perturbations of vertical velocity. Perturbations of potential temperature and scalar concentration decay faster and slower than vertical velocity perturbations, respectively. A downward propagation of negative heat flux and its oscillation are found and a quadrant analysis reveals that warmer air sinking events are responsible for the downward propagation. The fourth quadrant events seem to be induced by demixing of air parcels, entrained from above the CBL. The advective CBL simulation with geostrophic wind illustrates that near-surface eddies are mechanically generated and they decelerate flow from the bottom up in the CBL/residual layer. The two-dimensional spectra show the height- and scale-dependent characteristics of decaying convective turbulence again in the free and advective boundary layer simulations.


Author(s):  
Aaron Wang ◽  
Ying Pan ◽  
Paul M. Markowski

AbstractThis work explores the influence of Weighted Essentially Non-Oscillatory (WENO) schemes on Cloud Model 1 (CM1) large-eddy simulations (LES) of a quasi-steady, horizontally homogeneous, fully developed, neutral atmospheric boundary layer (ABL). An advantage of applying WENO schemes to scalar advection in compressible models is the elimination of acoustic waves and associated oscillations of domain-total vertical velocity. Applying WENO schemes to momentum advection in addition to scalar advection yields no further advantage, but has an adverse effect on resolved turbulence within LES. As a tool designed to reduce numerically generated spurious oscillations, WENO schemes also suppress physically realistic instability development in turbulence-resolving simulations. Thus, applying WENO schemes to momentum advection reduces vortex stretching, suppresses the energy cascade, reduces shear-production of resolved Reynolds stress, and eventually amplifies the differences between the surface-layer mean wind profiles in the LES and the mean wind profiles expected in accordance with the filtered law of the wall (LOTW). The role of WENO schemes in adversely influencing surface-layer turbulence has inspired a concept of anti-WENO (AWENO) schemes to enhance instability development in regions where energy-containing turbulent motions are inadequately resolved by LES grids. The success in reproducing the filtered LOTW via AWENO schemes suggests that improving advection schemes is a critical component toward faithfully simulating near-surface turbulence and dealing with other "Terra Incognita" problems.


Sign in / Sign up

Export Citation Format

Share Document