Results for “plume dynamics”

DOE OSTI2026

Impact of Forest Canopy Structure on Buoyant Plume Dynamics During Wildland Fires

Heterogeneous forest canopies can generate complex turbulent structures, but in the presence of a fire plume, these interactions are not fully understood. This study investigates the influence of forest canopy heterogeneity on buoyant plume dynamics resulting from surface thermal anomalies representing wildland fires, utilizing Large Eddy Simulation (LES). The Parallelized Large-Eddy Simulation Model (PALM) was employed to simulate six canopy configurations: no canopy, homogeneous canopy, external plume-edge canopy, internal plume-edge canopy, 100 m gap canopy, and 200 m gap canopy. Each configuration was analyzed with and without a static surface heat flux patch of 5000 W ∙ m -2 , resulting in a resting buoyant plume. Simulations were conducted under three crosswind speeds: 0, 5, and 10 m ∙ s -1 . Results show that canopy structure significantly modifies plume behavior, mean flow, and turbulent kinetic energy (TKE) budgets. Plume updraft speed and tilt varied with canopy configuration and crosswind speed. Horizontal pressure gradients associated with plume-atmosphere interaction were modified based on the canopy configuration, resulting in varying crosswind speed reductions at the plume region. Strong momentum absorption was observed above the canopy for the crosswind cases, with the greatest enhancement in the gap canopies. Momentum injection from below the canopy due to the heat source was also observed, resulting in plume structure modulation based on canopy configuration. TKE was found to be the largest in the gap canopy configurations. TKE budget analysis revealed that buoyant production dominated over shear production. At the center of the heat patch, the gap canopy configurations showed enhanced buoyancy within the gap. These results improve our knowledge of fire-canopy-atmosphere interactions that can inform fire models on the impacts of canopy heterogeneity on plume dynamics and ember ejections.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Machine Learning-Driven Quantification of CO2 Plume Dynamics at Illinois Basin Decatur Project Sites Using Microseismic Data

This study utilizes machine learning to quantify CO2 plume extents by analyzing microseismic data from the Illinois Basin Decatur Project (IBDP). Leveraging a unique dataset of well logs, microseismic records, and CO2 injection metrics, this work aims to predict the temporal evolution of subsurface CO2 saturation plumes. The findings illustrate that machine learning can predict plume dynamics, revealing vertical clustering of microseismic events over distinct time periods within certain proximities to the injection well, consistent with an invasion percolation model. The buoyant CO2 plume partially trapped within sandstone intervals periodically breaches localized barriers or baffles, which act as leaky seals and impede vertical migration until buoyancy overcomes gravity and capillary forces, leading to breakthroughs along vertical zones of weakness. Between different unsupervised clustering techniques, K-Means and DBSCAN were applied and analyzed in detail, where K-means outperformed DBSCAN in this specific study by indicating the combination of the highest Silhouette Score and the lowest Davies–Bouldin Index. The predictive capability of machine learning models in quantifying CO2 saturation plume extension is significant for real-time monitoring and management of CO2 sequestration sites. The models exhibit high accuracy, validated against physical models and injection data from the IBDP, reinforcing the viability of CO2 geological sequestration as a climate change mitigation strategy and enhancing advanced tools for safe management of these operations.

Iyegbekedo, Ikponmwosa↗

DOE OSTI2024

Evaluating crown scorch predictions from a computational fluid dynamics wildland fire simulator

Abstract Background Crown scorch—the heating of live leaves, needles, and buds in the vegetative canopy to lethal temperatures without widespread combustion—is one of the most common fire effects shaping post-fire canopies. Despite the ability of computational fluid dynamic models to finely resolve fire activity and buoyant plume dynamics including heterogenous 3D distributions of forest canopy heating, these models have had only limited use in simulating fire effects and have not been used to evaluate crown scorch. Here, we demonstrate a method of evaluating crown scorch using a computational fluid dynamics model, FIRETEC, and validate this approach by simulating the experiments that were used to develop Van Wagner’s 1973 crown scorch model. Results The average scorch height prediction from FIRETEC compares well with the empirical model derived by Van Wagner, which is the most widely used empirical model for crown scorch. We further find that the 3D buoyant plume dynamics from a steady and homogeneous idealized heat source on the ground results in a spatially heterogenous crown scorch pattern reflecting complex heating dynamics that are best represented by percent scorch rather than height of scorch. Conclusions The ability of the computational fluid dynamics model to capture variation in crown scorch due to 3D buoyant plume dynamics provides direct links between forest structure, fire behavior, and fire effects that can be used by forest managers and researchers to better understand how fires result in crown damage under various environmental and management scenarios.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Informing field-scale CO 2 storage simulations with sandbox experiments: The effect of small-scale heterogeneities

Small-scale heterogeneities can significantly affect the fate of the CO 2 plume and trapping during CO 2 migration. We conducted geologic carbon storage field-scale simulations to investigate the impact of small-scale heterogeneities on plume dynamics and trapping performance. Small-scale heterogeneities have been shown to increase the amount of trapped CO 2 during buoyancy-driven flow. The trapped CO 2 saturation is validated by previous sandbox experimental work during buoyancy-driven flow in realistic heterogeneous domains and is implemented through the critical CO 2 saturation parameter (i.e., the first non-zero value in the drainage CO 2 relative permeability curve). Depending on the type and degree of heterogeneity, various critical CO 2 saturation values are exhibited. Furthermore, we investigated the effect of small-scale heterogeneities when multiple capillary pressure models are employed. This study demonstrates that an increase in critical CO 2 saturation reduces the CO 2 plume size and lateral extent, accompanying an increase in residual trapping and a decrease in solubility trapping. Lastly, we show that independent of the capillary pressure model used, an increase in critical saturation leads to similar CO 2 plume dynamics distribution and trapping performance. These results emphasize the importance of quantifying the effect of small-scale heterogeneity as they affect the large-scale behavior of the CO 2 plume.

58 GEOSCIENCES↗

DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗

DOE OSTI2024

Height-resolved emission spectroscopy and high-speed imaging of the TiAl6V4 vapor plume under laser powder bed fusion conditions

Optical emission spectroscopy is increasingly used as an in situ monitoring technique during laser powder bed fusion (LPBF) because plume emission holds elemental information not found in other in situ sensing techniques. This work explores the shape, stability, and temperature of the emission plume above the melt pool of Ti6Al4 V undergoing laser scans under LPBF-like processing conditions, using high-speed video and height-resolved spectroscopy to study the details of plume dynamics. Optical emission spectroscopy is conducted in the 480 nm to 525 nm region, where Ti emission is strong, with 0.3 mm vertical resolution above the baseplate. The Boltzmann plot method is used to determine temperature at each probed elevation, which indicates that the hottest location of the plume is occasionally elevated 0.3 mm to 0.6 mm above the scanning surface. The results show that the plume shape, stability, temperature, and spectra are highly dependent on the processing conditions. We highlight some of the complexities of optical emission spectroscopy and discuss potential challenges for implementing optical emission spectroscopy across an entire build.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2025

Factors influencing water quality in surface water and alluvial groundwaters downgradient of a reclaimed surface coal mine in the Powder River Basin of southeastern Montana, USA

Abstract Coal mining and reclamation can have a profound influence on hydrogeologic systems, with clear consequences for groundwater quality, yet their long-term influence on downgradient water quality over time following reclamation is less well documented. Geochemical trends were evaluated in water quality downgradient of a fully reclaimed landscape at the former Big Sky Mine in the Rosebud Creek watershed (southeastern Montana, USA), over a 3-year period (2020–2022), including bond release in 2022. Within 6 km downgradient from the reclaimed area, sulfate concentrations decreased from approximately 3500 to 1800 mg l −1 within the Miller Coulee alluvial aquifer. Major ions, δ 34 S SO4 values, and residence time tracers suggest that the observed decreases in sulfate concentration result from a combination of dilution by mixed-age inflows and incomplete transit of the high salinity plume from the mine boundary. Both bedrock and alluvial aquifers of the Rosebud Creek corridor contained contributions of millennia-old regional groundwater, which may serve to mitigate mine-derived high salinity waters. Rosebud Creek, which traverses the outflow zone of Miller Coulee in the study area, exhibited high sulfate concentrations during low flows and consistent downgradient increases in sulfate concentration. The possibility of plume dynamics in Miller Coulee suggests that the greatest water quality impacts may not yet have reached Rosebud Creek.

Keeshin, Skye I. (ORCID:0009000986449872)↗

DOE OSTI2026

On the formation of swelling and related flaws in laser powder bed fusion

Process monitoring in laser powder bed fusion additive manufacturing can provide insights into stochastic anomalies, melt pool and plume dynamics, and part quality. Swelling, a build anomaly where overbuilt material protrudes through the powder layer after recoating, is readily detectable in post-recoat visible light images of the powder bed. Here, this work identifies several of the underlying mechanisms driving swelling formation by analyzing the influence of processing parameters, laser scan paths, and build plate locations on the presence of swelling detected in situ. Swelling near the edge of the part and swelling in the internal region of the part are shown to correlate with different process conditions. Edge and internal swelling may be driven by different phenomena, with edge swelling predominately occurring on the edge of a part facing the laser module and correlated to clusters of near-surface voids (detected with X-ray computed tomography). A larger spot size, higher laser power, and lower scan velocity also increased the presence of edge swelling. Laser spot size and scan path influenced internal swelling, which occurred preferentially with a larger spot size and in regions with large melt pools, caused by localized heat accumulation due to non-optimal processing parameters or scan path strategies. For coupons processed with a slicer-defined maximum scan vector length, swelling seldom occurred at internal vector-stripe boundaries. These results provide a mechanistic understanding of how swelling can be linked to material flaws, insight into how some instances of swelling can be avoided, and evidence supporting the use of swelling as an in situ indicator for quality assurance and part qualification.

Anomaly↗

DOE OSTI2026

Hourly PM 2.5 Estimates across California from 2018 to 2023

This study presents a new data set of hourly PM 2.5 concentrations across California from 2018 to 2023 at a three-kilometer resolution. This data set was developed by assimilating observations from PurpleAir and the U.S. EPA Air Quality System monitors into wildfire smoke forecasts from the High-Resolution Rapid Refresh Smoke (HRRR-Smoke) model using the Gridpoint Statistical Interpolation (GSI) three-dimensional variational data assimilation framework. Archived forecasts of modeled wildfire smoke PM 2.5 from HRRR-Smoke create the background field for assimilation, which is then corrected using surface observations of total PM 2.5 . The resulting reanalysis from GSI provides an estimate of total PM 2.5 that minimizes error from both the observational and the model data. Validation results indicate strong performance, with monthly R 2 values ranging from 0.73 to 0.91 across the six-year data set, comparable to other PM 2.5 data sets. Case studies are presented for three major fire events, the 2018 Camp Fire, 2019 Kincade Fire, and 2020 Lightning Complex Fires to demonstrate the data set’s fidelity in resolving plume dynamics and local exposure patterns. Root-mean-squared error averaged over each month scales with average PM 2.5 concentrations, resulting in a low error under typical conditions but higher absolute errors during extreme smoke events. This is the first long-term, hourly PM 2.5 data set of its kind for California and enables the generation of subdaily exposure metrics, such as peak hourly concentrations, exceedance durations, and time-of-day exposure peaks. The novelty and strong validation of this data set make it a compelling resource for future studies on the impact and significance of subdaily PM 2.5 exposure.

PM2.5↗

DOE OSTI2024

Machine Learning Applications in Analyzing the Role of Shale Barriers and Baffles for CO2 Storage

This study uses machine learning to analyze microseismic data from the Illinois Basin Decatur Project (IBDP) and quantify CO₂ plume extents. By leveraging well logs, microseismic records, and CO₂ injection metrics, the research predicts subsurface CO₂ plume dynamics. Findings show vertical clustering of microseismic events near the injection well, with CO₂ periodically breaching barriers due to buoyancy. K-Means clustering performed best, achieving the highest Silhouette Score and lowest Davies-Bouldin Index. This capability is crucial for real-time monitoring and management of CO₂ sequestration sites, validated against physical models and IBDP data, reinforcing CO₂ geological sequestration's viability and enhancing management tools.

Carr, Timothy↗

DOE OSTI2024

The influence of laser energy on deuterium emission characteristics from a Zircaloy-4 plasma

Laser-produced plasma coupled with optical emission spectroscopy (OES) is a promising technique for detecting certain isotopes, with unique capabilities such as standoff and rapid detection and minimal to no sample preparation requirements. The key figure-of-merit for isotopic analysis using optical spectroscopy tools is the linewidth relative to the isotope shift. Although the isotopes of hydrogen (1H, 2H, and 3H) possess large isotopic shifts (1H–2H ≈ 180 pm, 1H–3H ≈ 240 pm), being a light element, the H transitions are susceptible to various broadening mechanisms in the plasma environment. One of the critical parameters that influence the linewidth of a transition in an LPP is the incident laser energy. In the present study, we evaluated the role of laser energy on plume expansion dynamics, deuterium emission intensity, and linewidth in a nanosecond laser-produced Zircaloy-4 plasma. The changes in 2Hα emission intensity and linewidth were investigated for varying laser fluence and time after plasma onset. Spatially resolved and spatially integrated OES were performed and compared to investigate the emission spectral features and linewidth of 2Hα. Monochromatic two-dimensional time-resolved imaging was also performed to understand the morphology of the deuterium and protium emission relative to all species in the plume. Our results showed that 1Hα and 2Hα emissions predominantly occur closer to the target. Measurements of 2Hα linewidth approached similar values at later times of plasma evolution regardless of the laser energy. The linewidths of the 2Hα transition showed insignificant differences between spatially resolved and spatially integrated measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Effects of non-equilibrium ionization and excitation on radiation absorption in plasma plumes induced by ablation of metal targets with nanosecond laser pulses

Ionization and radiation absorption in nanosecond laser-induced plasma plumes are routinely modeled using the Saha–Boltzmann equilibrium ionization model (EQM). However, the equilibrium assumption can be inaccurate during the laser pulse when non-equilibrium effects significantly impact radiation absorption. In the present work, the EQM and non-equilibrium collisional-radiative model (CRM) are compared to reveal the effect of plasma non-equilibrium on radiation absorption in non-homogeneous plumes and degree of plasma shielding. Simulations of plume expansion induced by irradiation of a copper target in 1 atm argon background gas with a 10 ns Gaussian pulse at a fluence from 8 Jcm −2 to 14 Jcm −2 are performed with a hybrid computational model that couples a lumped particle direct simulation Monte Carlo method with either CRM or EQM. The simulations show that the EQM strongly underestimates the effects of ionization and radiation absorption compared to CRM and, contrary to the CRM, predict strong ionization of the background gas. The differences between the models are explained by the qualitatively different coupling between plume expansion and dynamics of ionization and exitation processes in the CRM and EQM under conditions when the characteristic times for most radiation- and electron-induced processes are longer than the pulse duration. The CRM-based predictions are also found to agree much better with available experimental data. In conclusion, these results indicate that the model of Saha–Boltzmann equilibrium cannot be used for reliable prediction of the degree of plasma shielding in plumes induced by nanosecond laser pulses or for processing results of spectroscopic measurements at early stages of expansion of such plumes.

Physics - Plasma physics↗

DOE OSTI2026

Assessing shellfish water exposure to fecal bacteria pollution in Salish Sea: three-dimensional modeling and implications for monitoring

Fecal bacteria (FB) contamination poses significant risks to shellfish safety and management in coastal and estuarine waters. Despite extensive pollution identification and correction efforts, FB contamination in shellfish-growing areas persists in the Salish Sea, highlighting the need to identify overlooked sources and better understand FB transport from riverine and shoreline inputs to shellfish beds. To address this, a high-resolution three-dimensional hydrodynamic model coupled with FB kinetics was developed and applied to a case study site in Salish Sea—Portage Bay—to simulate freshwater plume circulation, flushing dynamics, and bacterial transport. Daily FB loading from the major freshwater inflow—Nooksack River was generated by both linear interpolation and integrating a machine learning approach (XGBoost), trained on historical hydrological and meteorological data. The model successfully reproduced both the magnitude and seasonal variation of FB concentrations in Portage Bay for the year of 2021, demonstrating that simplified FB kinetics with first-order decay due to mortality was effective in this dynamic coastal environment with short flushing time. Model results identified the Nooksack River as the dominant far-field FB source, while scenario simulations showed that near-field coastal stormwater outfalls elevated local FB levels following rainfall, particularly under low-flow conditions. The XGBoost prediction provided comparable or superior accuracy to linear interpolation, particularly during periods of missing observational data, by capturing short-term variability and event-driven loading more effectively. Integrating data-driven riverine FB inputs with mechanistic coastal numerical modeling provides a robust framework for operational forecasting of shellfish bed exposure risk and supports adaptive monitoring and management of shellfish growing areas in the Salish Sea and similar coastal systems.

Salish Sea↗

DOE OSTI2025

Expansion dynamics of femtosecond laser-induced plasmas: Influence of thermophysical plasma properties

This study investigates the expansion dynamics of femtosecond laser-induced plasmas, emphasizing the impact of plasma thermophysical properties and ambient gas composition. Through shadowgraphy experiments and multiphase computational fluid dynamics (CFD) simulations, the influence of parameters such as heat capacity, molecular weight, and thermal conductivity on plume morphology, shockwave evolution, and energy dissipation mechanisms is examined. A mixture multiphase model is implemented to capture the interaction between the plasma and the surrounding gas. Simulation results reveal that plasma expansion is strongly inertia-driven. Results show that differences in plasma properties and ambient conditions affect the shape and temperature distribution of the expanding plume. The early-stage dynamics are primarily dictated by pressure forces, whereas thermal and viscous effects play a growing role in the plume's behavior during later stages of expansion. The CFD findings show the necessity of accurate initial condition characterization, including crater geometry and plasma pressure and temperature, for reliable modeling of plasma evolution in laser ablation processes.

CFD modeling↗

DOE OSTI2025

Spatiotemporal predictions of toxic urban plumes using deep learning

Industrial accidents, chemical spills, and structural fires can release large amounts of harmful materials that disperse into urban atmospheres and impact populated areas. Computer models are typically used to predict the transport of toxic plumes by solving fluid dynamical equations. However, these models can be computationally expensive due to the need for many grid cells to simulate turbulent flow and resolve individual buildings and streets. In emergency response situations, alternative methods are needed that can run quickly and adequately capture important spatiotemporal features. Here, we present a novel deep learning model called ST-GasNet inspired by the mathematical equations that govern the behavior of plumes as they disperse through the atmosphere. ST-GasNet learns the spatiotemporal dependencies from a limited set of temporal sequences of ground-level toxic urban plumes generated by a high-resolution large eddy simulation model. On independent sequences, ST-GasNet accurately predicts the late-time spatiotemporal evolution, given the early-time behavior as an input, even when a building splits a large plume into smaller plumes. By incorporating large-scale wind boundary condition information, ST-GasNet achieves a prediction accuracy of at least 90% on test data for the entire prediction period.

Civil and Environmental Engineering↗

DOE OSTI2025

High fidelity simulations of contaminant dispersion in an urban environment with comparison to magnetic resonance imaging measurements

The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.

54 ENVIRONMENTAL SCIENCES↗