An official website of the United States government
Abstract Mesoscale‐to‐microscale coupling is an important tool for conducting turbulence‐resolving multiscale simulations of realistic atmospheric flows, which are crucial for applications ranging from wind energy to wildfire spread studies. Different techniques are used to facilitate the development of realistic turbulence in the large‐eddy simulation (LES) domain while minimizing computational cost. Here, we explore the impact of a simple and computationally efficient Stochastic Cell Perturbation method using momentum perturbation (SCPM‐M) to accelerate turbulence generation in boundary‐coupled LES simulations using the Weather Research and Forecasting model. We simulate a convective boundary layer (CBL) to characterize the production and dissipation of turbulent kinetic energy (TKE) and the variation of TKE budget terms. Furthermore, we evaluate the impact of applying momentum perturbations of three magnitudes below, up to, and above the CBL on the TKE budget terms. Momentum perturbations greatly reduce the fetch associated with turbulence generation. When applied to half the vertical extent of the boundary layer, momentum perturbations produce an adequate amount of turbulence. However, when applied above the CBL, additional structures are generated at the top of the CBL, near the inversion layer. The magnitudes of the TKE budgets produced by SCPM‐M when applied at varying heights and with different perturbation amplitudes are always higher near the surface and inversion layer than those produced by No‐SCPM, as are their contributions to the TKE. This study provides a better understanding of how SCPM‐M reduces computational costs and how different budget terms contribute to TKE in a boundary‐coupled LES simulation.
more »
« less
The Influence of WENO Schemes on Large-Eddy Simulations of a Neutral Atmospheric Boundary Layer
Abstract This 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.
more »
« less
- Award ID(s):
- 1821885
- PAR ID:
- 10353440
- Date Published:
- Journal Name:
- Journal of the Atmospheric Sciences
- ISSN:
- 0022-4928
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Hurricanes are among the costliest natural disasters in the United States and regularly inflict severe damage on urban infrastructure. Accurate forecasts are therefore essential for preparedness and limiting these extreme events' economic toll. Numerical weather‑prediction (NWP) models—such as the Weather Research and Forecasting (WRF) system—are powerful forecasting tools. However, some of their physical parameterizations were neither designed for nor tested with real hurricanes. This thesis addresses that gap by evaluating two key parameterizations in WRF: (i) subgrid‑scale (SGS) turbulence schemes and (ii) surface‑roughness and urban canopy treatments. The first part of the study investigates how SGS eddy‑viscosity choices affect hurricane intensity, turbulence, and wind profiles. Large‑eddy simulations (LES) of five major hurricanes were run with a 1.5‑order, three‑dimensional turbulent‑kinetic‑energy (TKE) SGS scheme. Each storm was simulated under three eddy‑viscosity settings— default, halved, and doubled—yielding 15 cases. A parallel set of 10 cases employed an alternative nonlinear backscatter and anisotropy (NBA) SGS scheme. Two idealized LES runs and one fine-grid (~80 m) nested simulation brought the total to 33. Reducing SGS stress intensified storms by raising boundary‑layer wind speeds and lowering the altitude of peak winds, improving surface‑wind forecasts by ~9 % and minimum sea‑level pressure by ~29 % relative to the default setting. These results reveal that standard SGS models are overly dissipative because they overlook the rotational suppression of turbulence, underscoring the need for SGS schemes tailored to hurricane dynamics. The second part assesses how aerodynamic roughness length (z0) and urban‑canopy schemes shape near‑surface winds over cities. For four land‑falling hurricanes affecting Houston and New Orleans, increasing z0 in the Single‑Layer Urban Canopy Model (SLUCM) reduced modeled wind speeds and cut mean absolute error (MAE) by ~20 %, whereas decreasing z0 introduced large positive biases. Additional experiments compared three urban options—Bulk (no‑urban), SLUCM, and the multi‑layer Building Energy Model (BEM). The Bulk scheme delivered the most accurate surface‑wind forecasts in every nested domain, while SLUCM slightly outperformed BEM in the limited vertical‑profile data. These findings highlight the need to recalibrate urban schemes and surface‑drag parameters when applying WRF to hurricane‑force winds.more » « less
-
Extensive experimental evidence highlights that scalar turbulence exhibits anomalous diffusion and stronger intermittency levels at small scales compared to that in fluid turbulence. This renders the corresponding subgrid-scale dynamics modeling for scalar turbulence a greater challenge to date. We develop a new large eddy simulation (LES) paradigm for efficiently and dynamically nonlocal LES modeling of the scalar turbulence. To this end, we formulate the underlying nonlocal model starting from the filtered Boltzmann kinetic transport equation, where the divergence of subgrid-scale scalar fluxes emerges as a fractional-order Laplacian term in the filtered advection–diffusion model, coding the corresponding superdiffusive nature of scalar turbulence. Subsequently, we develop a robust data-driven algorithm for estimation of the fractional (noninteger) Laplacian exponent, where we, on the fly, calculate the corresponding model coefficient employing a new dynamic procedure. Our a priori tests show that our new dynamically nonlocal LES paradigm provides better agreement with the ground-truth filtered direct numerical simulation data in comparison to the conventional static and dynamic Prandtl–Smagorinsky models. Moreover, in order to analyze the numerical stability and assessing the model's performance, we carry out comprehensive a posteriori tests. They unanimously illustrate that our new model considerably outperforms other existing functional models, correctly predicting the backscattering phenomena and, at the same time, providing higher correlations at small-to-large filter sizes. We conclude that our proposed nonlocal subgrid-scale model for scalar turbulence is amenable for coarse LES and very large eddy simulation frameworks even with strong anisotropies, applicable to environmental applications.more » « less
-
Abstract Mesoscale climate models provide indispensable tools to understand land‐atmosphere interactions over urban regions. However, uncertainties in urban canopy parameters (UCPs) and parameterization schemes lead to degraded representation of the drag effect in complex built terrains. In particular, for the widely applied single‐layer urban canopy model (SLUCM) coupled with the Weather Research and Forecasting (WRF) model, near‐surface horizontal wind speed is known to be overestimated systematically. In this study, idealized large eddy simulations (LES) and WRF‐SLUCM simulations are conducted to study the separate effect of UCPs and aerodynamic parameterization on atmospheric boundary layer processes and rainfall variabilities in Phoenix, Arizona. For LES that explicitly resolves surface geometry, significant differences between three‐dimensional (3D) versus two‐dimensional (2D) representation of urban morphology are found in the surface layer and above. When surface drag is parameterized following SLUCM, surface morphologies have little impacts on the mean momentum transfer. WRF‐SLUCM simulation results, incorporated with 3D urban morphology data, indicate that simply refining the frontal area index will reduce the surface drag, which further amplifies the systematic positive bias of SLUCM in predicting horizontal wind speed. Replacing the drag parameterization in SLUCM by LES‐based aerodynamic parameters has evident impacts on near‐surface wind speed. The impact of urban roughness representation becomes the most evident during rainfall periods, due to the important role of surface drag in dictating moisture convergence. Our study underlines that apart from intensive efforts in obtaining detailed UCPs, it is also critical to enhance the urban momentum exchange parameterization schemes.more » « less
-
Abstract The impacts of spurious numerical salinity mixing on the larger‐scale flow and tracer fields are characterized using idealized simulations. The idealized model is motivated by realistic simulations of the Texas‐Louisiana shelf and features oscillatory near‐inertial wind forcing. can exceed the physical mixing from the turbulence closure in frontal zones and within the mixed layer. This suggests that simulated mixing processes in frontal zones are driven largely by . Near‐inertial alongshore wind stress amplitude is varied to identify a base case that maximizes the ratio of to in simulations with no prescribed horizontal mixing. We then test the sensitivity of the base case with three tracer advection schemes (MPDATA, U3HC4, and HSIMT) and conduct ensemble runs with perturbed bathymetry. Instability growth is evaluated using the volume‐integrated eddy kinetic energy and available potential energy . While all schemes have similar total mixing, the HSIMT simulations have over double the volume‐integrated and 20% less relative to other schemes, which suppresses the release of and reduces the by roughly 25%. This results in reduced isohaline variability and steeper isopycnals, evidence that enhanced suppresses instability growth. Differences in and between the MPDATA and U3HC4 simulations are marginal. However, the U3HC4 simulations have 25% more . Experiments with variable horizontal viscosity and diffusivity coefficients show that small amounts of prescribed horizontal mixing improve the representation of the ocean state for all advection schemes by reducing the and increasing the .more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1175/JAS-D-21-0033.1
