skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Lagrangian subgrid-scale modeling applied to evolving firebrand particle transport
The transport and deposition of firebrand particles is an important fire spread mechanism in wildland fires. These particles can be transported by wind over large distances and can ignite secondary fires upon landing. The transport of firebrands by wind is a complex, multiscale process that is largely controlled by interactions between the firebrand particles and the atmospheric wind. To account for the complex temporal evolution of atmospheric turbulence over large scales, the use of large-eddy simulation (LES) techniques is necessary. However, filtering of subgrid-scale (SGS) turbulence in LES hinders the accuracy of particle transport models. In this work, we employ a Lagrangian SGS model in an LES framework to investigate the effects of small-scale turbulence on the transport of mass- and size-changing firebrand particles. The impact of SGS turbulence was analyzed by comparing landing and trajectory statistics for firebrand and regular (fixed size and mass) particles under different Stokes numbers. It was found that the presence of SGS turbulence modifies the particle transport behavior, which is characterized by smaller spanwise dispersions but larger travel distances along the streamwise direction compared with particles under no SGS turbulence. As expected, the enhanced velocity field produced by the SGS model has larger influence on the statistics of firebrand particles compared with regular particles due to the time-evolving reduction in particle mass and size induced by pyrolysis.  more » « less
Award ID(s):
2043103
PAR ID:
10397584
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Stanford University
Date Published:
Journal Name:
Center for Turbulence Research Proceedings of the Summer Program 2022
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Langmuir turbulence consists of Langmuir circulation (LC) generated at the surface of rivers, lakes, bays, and oceans by the interaction between the wind-driven shear and surface gravity waves. In homogeneous shallow water, LC can extend to the bottom of the water column and interact with the bottom boundary layer. Large-eddy simulation (LES) of LC in shallow water was performed with the finite volume method and various forms of subgrid-scale (SGS) model characterized by different near-wall treatments of the SGS eddy viscosity. The wave forcing relative to wind forcing in the LES was set following the field measurements of full-depth LC during the presence of LC engulfing a water column 15 m in depth in the coastal ocean, reported in the literature. It is found that the SGS model can greatly impact the structure of LC in the lower half of the water column. Results are evaluated in terms of (1) the Langmuir turbulence velocity statistics and (2) the lateral (crosswind) length scale and overall cell structure of LC. LES with an eddy viscosity with velocity scale in terms of S and Ω (where S is the norm of the strain rate tensor and Ω is the norm of the vorticity tensor) and a Van Driest wall damping function (referred to as the S-Omega model) is found to provide best agreement with pseudo-spectral LES in terms of the lateral length scale and overall cell structure of LC. Two other SGS models, namely the dynamic Smagorinsky model and the wall-adapting local-eddy viscosity model are found to provide less agreement with pseudo- spectral LES, for example, as they lead to less coherent bottom convergence of the cells and weaker associ ated upward transport of slow downwind moving fluid. Finally, LES with the S-Omega SGS model is also found to lead to good agreement with physical measurements of LC in the coastal ocean in terms of Langmuir turbulence decay during periods of surface heating 
    more » « less
  2. Spot fires pose a major risk and add to the already complex physics, which makes fire spread so hard to predict, especially in the wildland urban interface. Firebrands can not only cross fuel breaks and thwart other suppression efforts but also directly damage infrastructure and block evacuation routes. Transport models and computational fluid dynamics tools often make simplifications when predicting spot fire risk, but there is a relative lack of experimental data to validate such parameterizations. To this end, we present a field experiment performed at the University of California Berkeley Blodgett Research Forest in California where we recorded the flame and firebrands emanating from a nighttime hand-drawn pile fire using high-frequency imaging. We used image-processing to characterize the fire intensity and turbulence as well as particle tracking velocimetry to measure ejected firebrand kinematics as they are lofted by the plume. We further collected embers that settled around the fire at varying distances and measured their size, shape, density, and settling distributions. We also examine existing physics-based time-averaged models of firebrand lofting and note discrepancies between such models, often used due to their speed and simplicity, and our experimental observations. Finally, we discuss some implications our observations could have on future modeling efforts by considering the time-dependent fire dynamics, intermittency in the plume turbulence, and in the firebrand generation rate. To the best of our knowledge, these are the first in situ observations of firebrand generation and lofting from representative fuels, addressing a major source of data gap and uncertainty in the wildland fire literature. 
    more » « less
  3. Abstract We present a new ensemble of 36 numerical experiments aimed at comprehensively gauging the sensitivity of nested large-eddy simulations (LES) driven by large-scale dynamics. Specifically, we explore 36 multiscale configurations of the Weather Research and Forecasting (WRF) Model to simulate the boundary layer flow over the complex topography at the Perdigão field site, with five nested domains discretized at horizontal resolutions ranging from 11.25 km to 30 m. Each ensemble member has a unique combination of the following input factors: (i) large-scale initial and boundary conditions, (ii) subgrid turbulence modeling in thegray zoneof turbulence, (iii) subgrid-scale (SGS) models in LES, and (iv) topography and land-cover datasets. We probe their relative importance for LES calculations of velocity, temperature, and moisture fields. Variance decomposition analysis unravels large sensitivities to topography and land-use datasets and very weak sensitivity to the LES SGS model. Discrepancies within ensemble members can be as large as 2.5 m s−1for the time-averaged near-surface wind speed on the ridge and as large as 10 m s−1without time averaging. At specific time points, a large fraction of this sensitivity can be explained by the different turbulence models in the gray zone domains. We implement a horizontal momentum and moisture budget routine in WRF to further elucidate the mechanisms behind the observed sensitivity, paving the way for an increased understanding of the tangible effects of the gray zone of turbulence problem. Significance StatementSeveral science and engineering applications, including wind turbine siting and operations, weather prediction, and downscaling of climate projections, call for high-resolution numerical simulations of the lowest part of the atmosphere. Recent studies have highlighted that such high-resolution simulations, coupled with large-scale models, are challenging and require several important assumptions. With a new set of numerical experiments, we evaluate and compare the significance of different assumptions and outstanding challenges in multiscale modeling (i.e., coupling large-scale models and high-resolution atmospheric simulations). The ultimate goal of this analysis is to put each individual assumption into the wider perspective of a realistic problem and quantify its relative importance compared to other important modeling choices. 
    more » « less
  4. 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
  5. Abstract In large‐eddy simulations, subgrid‐scale (SGS) processes are parameterized as a function of filtered grid‐scale variables. First‐order, algebraic SGS models are based on the eddy‐viscosity assumption, which does not always hold for turbulence. Here we apply supervised deep neural networks (DNNs) to learn SGS stresses from a set of neighboring coarse‐grained velocity from direct numerical simulations of the convective boundary layer at friction Reynolds numbersReτup to 1243 without invoking the eddy‐viscosity assumption. The DNN model was found to produce higher correlation between SGS stresses compared to the Smagorinsky model and the Smagorinsky‐Bardina mixed model in the surface and mixed layers and can be applied to different grid resolutions and various stability conditions ranging from near neutral to very unstable. The DNN model can capture key statistics of turbulence ina posteriori(online) tests when applied to large‐eddy simulations of the atmospheric boundary layer. 
    more » « less