Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Burn”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Synopsis of SURF and SURFplus

Developed by Menikoff and Shaw, the SURF reactive burn model builds on the Ignition and Growth concept by incorporating the lead shock pressure directly into the volumetric hot-spot burn rate. The plus-extension augments the model with a late time surface burn rate due to carbon clustering. Together, SURF and SURFplus allow for robust modeling of both conventional and insensitive high explosives. Practically, the SURF burn models require shock detection and the advection of the lead shock pressure as an additional material field. Calibrated parameters for given equations of state and thermodynamic closure are dependent on initial temperature and density. An implementation in Python is given.

42 ENGINEERING↗
DOE OSTI2026

Lighting a Star on Earth—Again and Again

At Lawrence Livermore National Laboratory’s (LLNL) National Ignition Facility (NIF), scientists are learning what it takes to keep a fusion reaction burning. Their latest experiments, presented this week at the American Physical Society’s Division of Plasma Physics annual meeting, build on a series of successful ignition shots to reveal how capsule design and laser precision can sustain record-setting performance. At NIF, researchers are consistently producing plasmas dominated by self-heating from fusion reactions, achieving target gain greater than two — meaning the fusion output energy exceeds twice the energy used to initiate the reaction. Using the world’s most energetic laser to create these conditions, LLNL scientists can explore the detailed physics of ignition and burn that underpin the Laboratory’s stockpile modernization mission. “It’s a bit like lighting a campfire with a match and a single breath,” said LLNL physicist David Schlossberg, who is co-leading the work. “If you have too little fuel, it burns out quickly; too much, and it won’t light at all. We’ve found the balance between capsule mass and drive strength that starts the burn and keeps it going efficiently.” NIF is the only laboratory in the world to have achieved fusion ignition, a historic breakthrough first realized in December 2022. Since then, the facility has repeated ignition ten times, most recently in October 2025, demonstrating the robustness of the experimental platform and the reproducibility of ignition physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Mass dependence of overshooting beneath convective envelopes

State of the art stellar evolution codes [Paxton et al., 2010, Demarque et al., 2004, Weiss and Schlattl, 2008, Siess et al., 2013, Christensen-Dalsgaard, 2008] evoke a diffusive process to model convective overshooting. The diffusion coefficient [Freytag et al., 1996, Pratt et al., 2017] used in the model can be set to change with the classification of the convective zone as non-burning, H-burning, He-burning, or metal-burning. It can also be set to change at defined evolutionary points, such as the bottom of the asymptotic giant branch, or during the third dredge up [Herwig, 2000, Lugaro et al., 2003]. Aside from these abrupt changes, the diffusion coefficient is typically locked to a percentage of the pressure scale height measured at the convective boundary. However, there is no theoretical reason for convective overshooting, or indeed other convective properties, to change in the same way that the pressure scale height changes as a star evolves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Deflagration to Detonation Transition Update: XDDT Code Modularization

A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Assessment of a reconfiguration of the InterSpread Plus US national FMD model as a potential tool to analyze a foot-and-mouth disease outbreak on a single large cattle feedlot in the United States

An incursion of foot-and-mouth disease (FMD) into the United States remains a concern of high importance and would have devastating socioeconomic impacts to the livestock and associated industries. This highly transmissible and infectious disease poses continual risk for introduction into the United States (US), due to the legal and illegal global movement of people, animals, and animal products. While stamping out has been shown to effectively control FMD, depopulation of large cattle feedlots (>50,000 head) presents a number of challenges for responders due to the resources required to depopulate and dispose of large numbers of animals in a timely and effective manner. However, evaluating alternative strategies for FMD control on large feedlots requires a detailed within-farm modeling approach, which can account for the unique structure of these operations. To address this, we developed a single feedlot, within-farm spread model using a novel configuration within the InterSpread Plus (ISP) framework. As proof of concept we designed six scenarios: (i) depopulation - the complete depopulation of the feedlot, (ii) burn-through – a managed “burn-through” where the virus is allowed to spread through the feedlot and only movement restriction and biosecurity are implemented, (iii) firebreak-NV – targeted depopulation of infected pens and adjacent pens without vaccination; (iv) firebreak - targeted depopulation of infected pens and adjacent pens with vaccination of remaining pens; (v) harvest-NV - selective harvest of pens where a 100% movement restriction is applied for 28-30 days, then pens are set for selection to be sent to slaughter, while allowing a controlled “burn-through” without vaccination; and (vi) harvest - selective harvest of pens with vaccination. Overall, the burn-through scenario (ii) had the shortest epidemic duration (31d (30, 33)) median (25th, 75th percentiles), while the firebreak scenario (iv) had the longest (47d (38,55)). Additionally, we found that scenarios implementing depopulation delayed the peak day of infection and reduced the total number of pens infected compared to non-depopulation scenarios.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Long-range transported continental aerosol in the eastern North Atlantic: three multiday event regimes influence cloud condensation nuclei

Abstract. The eastern North Atlantic (ENA) is a region dominated by pristine marine environment and subtropical marine boundary layer clouds. Under unperturbed atmospheric conditions, the regional aerosol regime in the ENA varies seasonally due to different seasonal surface-ocean biogenic emissions, removal processes, and meteorological regimes. However, during periods when the marine boundary layer aerosol in the ENA is impacted by particles transported from continental sources, aerosol properties within the marine boundary layer change significantly, affecting the concentration of cloud condensation nuclei (CCN). Here, we investigate the impact of long-range transported continental aerosol on the regional aerosol regime in the ENA using data collected at the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility on Graciosa Island in 2017 during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign. We develop an algorithm that integrates number concentrations of particles with optical particle dry diameter (Dp) between 100 and 1000 nm, single scattering albedo, and black carbon concentration to identify multiday events (with duration >24 consecutive hours) of long-range continental aerosol transport in the ENA. In 2017, we detected nine multiday events of long-range transported particles that correspond to ∼ 7.5 % of the year. For each event, we perform HYSPLIT 10 d backward trajectories analysis, and we evaluate CALIPSO aerosol products to assess, respectively, the origins and compositions of aerosol particles arriving at the ENA site. Subsequently, we group the events into three categories, (1) mixture of dust and marine aerosols, (2) mixture of marine and polluted continental aerosols from industrialized areas, and (3) biomass burning aerosol from North America and Canada, and we evaluate their influence on aerosol population and cloud condensation nuclei in terms of potential activation fraction and concentrations at supersaturation of 0.1 % and 0.2 %. The arrival of plumes dominated by the mixture of dust and marine aerosol in the ENA in the winter caused significant increases in baseline Ntot. Simultaneously, the baseline particle size modes and CCN potential activation fraction remained almost unvaried, while cloud condensation nuclei concentrations increased proportionally to Ntot. Events dominated by a mixture of marine and polluted continental aerosols in spring, fall, and winter led to a statistically significant increase in baseline Ntot, a shift towards larger particular sizes, a higher CCN potential activation fractions, and cloud condensation nuclei concentrations of >170 % and up to 240 % higher than during baseline regime. Finally, the transported aerosol plumes characterized by elevated concentration of biomass burning aerosol from continental wildfires detected in the summertime did not statistically contribute to increase baseline aerosol particle concentrations in the ENA. However, particle diameters were larger than under baseline conditions, and CCN potential activation fractions were >75 % higher. Consequentially, cloud concentration nuclei concentrations increased by ∼ 115 % during the period affected by the biomass burning events. Our results suggest that, through the year, multiday events of long-range continental aerosol transport periodically affect the ENA and represent a significant source of CCN in the marine boundary layer. Based on our analysis, in 2017, the multiday aerosol plume transport dominated by a mixture of dust and marine aerosol, a mixture of marine and polluted continental aerosols, and biomass burning aerosols caused increases in the NCCN baseline regime of, respectively, 6.6 %, 8 %, and 7.4 % at SS 0.1 % (and, respectively, 6.5 %, 8.2 %, and 7.3 % at SS 0.2 %) in the ENA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Landsat near-infrared (NIR) band and ELM-FATES sensitivity to forest disturbances and regrowth in the Central Amazon

Forest disturbance and regrowth are key processes in forest dynamics, but detailed information on these processes is difficult to obtain in remote forests such as the Amazon. We used chronosequences of Landsat satellite imagery (Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus) to determine the sensitivity of surface reflectance from all spectral bands to windthrow, clear-cut, and clear-cut and burned (cut + burn) and their successional pathways of forest regrowth in the Central Amazon. We also assessed whether the forest demography model Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM-FATES, accurately represents the changes for windthrow and clear-cut. The results show that all spectral bands from the Landsat satellites were sensitive to the disturbances but after 3 to 6 years only the near-infrared (NIR) band had significant changes associated with the successional pathways of forest regrowth for all the disturbances considered. In general, the NIR values decreased immediately after disturbance, increased to maximum values with the establishment of pioneers and early successional tree species, and then decreased slowly and almost linearly to pre-disturbance conditions with the dynamics of forest succession. Statistical methods predict that NIR values will return to pre-disturbance values in about 39, 36, and 56 years for windthrow, clear-cut, and cut + burn disturbances, respectively. The NIR band captured the observed, and different, successional pathways of forest regrowth after windthrow, clear-cut, and cut + burn. Consistent with inferences from the NIR observations, ELM-FATES predicted higher peaks of biomass and stem density after clear-cuts than after windthrows. ELM-FATES also predicted recovery of forest structure and canopy coverage back to pre-disturbance conditions in 38 years after windthrows and 41 years after clear-cut. The similarity of ELM-FATES predictions of regrowth patterns after windthrow and clear-cut to those of the NIR results suggests the NIR band can be used to benchmark forest regrowth in ecosystem models. Our results show the potential of Landsat imagery data for mapping forest regrowth from different types of disturbances, benchmarking, and the improvement of forest regrowth models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ecosystem leaf area, gross primary production, and evapotranspiration responses to wildfire in the Columbia River basin

Wildfires impact vegetation mortality and productivity and are increasing in intensity, frequency, and spatial area in the western United States. The rates of vegetation recovery after fires play a major role in the reestablishment of biomass and ecosystem functioning (e.g., structure, resilience, and productivity), but such recovery rates are poorly understood. Here we use remotely sensed data products from the Moderate Resolution Imaging Spectroradiometer (MODIS) to quantify the resistance and resilience of leaf area index (LAI), gross primary production (GPP), and evapotranspiration (ET) to 138 wildfires of various burn severity across the Columbia River basin (CRB) of the Pacific Northwest in 2015. Increasing burn severity caused lower resistance and resilience for all three variables. Resistance and resilience are highest in grasslands, intermediate in savanna, and lowest in needleleaf evergreen forests, consistent with the adaptation of these vegetation types to fire. LAI has consistently lower resistance and resilience than GPP and ET, which is consistent with physical and physiological mechanisms that compensate for reduced LAI. Resilience is influenced by precipitation, vapor pressure deficit (VPD), and burn severity across all three vegetation types; however, burn severity plays a more minor role in grasslands. Increasing wildfire severity will reduce the resistance and resilience and lengthen the recovery time of vegetation structure and fluxes with climate change, with significant consequences for the provision of ecosystem functioning and implications for model predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of Global Fire Simulations in CMIP6 Earth System Models

Fire is the primary form of terrestrial ecosystem disturbance on a global scale and an important Earth system process. Most Earth system models (ESMs) have incorporated fire modeling, with 19 of them submitting model outputs of fire-related variables to the Coupled Model Intercomparison Project Phase 6 (CMIP6). This study provides the first comprehensive evaluation of CMIP6 historical fire simulations by comparing them with multiple satellite-based products and charcoal-based historical reconstructions. Our results show that most CMIP6 models simulate the present-day global burned area and fire carbon emissions within the range of satellite-based products. They also capture the major features of observed spatial patterns and seasonal cycles, the relationship of fires with precipitation and population density, and the influence of the El Niño–Southern Oscillation (ENSO) on the interannual variability of tropical fires. Regional fire carbon emissions simulated by the CMIP6 models from 1850 to 2010 generally align with the charcoal-based reconstructions, although there are regional mismatches, such as in southern South America and eastern temperate North America prior to the 1910s and in temperate North America, eastern boreal North America, Europe, and boreal Asia since the 1980s. The CMIP6 simulations have addressed three critical issues identified in CMIP5: (1) the simulated global burned area being less than half of that of the observations, (2) the failure to reproduce the high burned area fraction observed in Africa, and (3) the weak fire seasonal variability. Furthermore, the CMIP6 models exhibit improved accuracy in capturing the observed relationship between fires and both climatic and socioeconomic drivers and better align with the historical long-term trends indicated by charcoal-based reconstructions in most regions worldwide. However, the CMIP6 models still fail to reproduce the decline in global burned area and fire carbon emissions observed over the past 2 decades, mainly attributed to an underestimation of anthropogenic fire suppression, and the spring peak in fires in the Northern Hemisphere midlatitudes, mainly due to an underestimation of crop fires. In addition, the model underestimates the fire sensitivity to wet–dry conditions, indicating the need to improve fuel wet-ness estimation. Based on these findings, we present specific guidance for fire scheme development and suggest a postprocessing methodology for using CMIP6 multi-model outputs to generate reliable fire projection products.

Wildfire, Earth system models↗
DOE OSTI2022

Combustion characteristics and detailed simulations of surrogates for a Tier II gasoline certification fuel

An experimental and numerical study of combustion of a gasoline certification fuel (‘indolene’), and four (S4) and five (S5) component surrogates for it, is reported for the configurations of an isolated droplet burning with near spherical symmetry in the standard atmosphere, and a single cylinder engine designed for advanced compression ignition of pre-vaporized fuel. The intent was to compare performance of the surrogate for these different combustion configurations and to assess the broader applicability of the kinetic mechanism and property database for the simulations. A kinetic mechanism comprised of 297 species and 16797 reactions was used in the simulations that included soot chemistry and accounted for unsteady transport, liquid diffusion inside the droplet, radiative heat transfer, and variable properties. The droplet data showed a clear preference for the S5 surrogate in terms of burning rate. The simulations showed generally very good agreement with measured droplet, flame, and soot shell diameters. Measurements of combustion timing, in-cylinder pressure, and mass-averaged gas temperature were also well predicted with a slight preference for the S5 surrogate. Preferential vaporization was not evidenced from the evolution of droplet diameter but was clearly revealed in simulations of the evolution of mixture fractions inside the droplets. The influence of initial droplet diameter (Do) on droplet burning was strong, with S5 burning rates decreasing with increasing Do due to increasing radiation losses from the flame. Flame extinction was predicted for Do =3.0 mm as a radiative loss mechanism but not predicted for smaller Do for the conditions of the simulations.

Guo, Songtao↗
DOE OSTI2021

Combustion studies of n-butyl acetate synthesized by a new biological process and comparisons with neat n-butyl acetate

Blending petroleum fuels with biofuels is a common approach for stemming the depletion of crude oil while also mitigating the impact of their combustion on the environment. It has recently been considered that n-butyl acetate (BA, C6H12O2, boiling point of 399K) could be a viable biofuel additive to diesel fuel. This paper reports the combustion dynamics of BA droplets in the absence of external convection. Two grades of BA were examined: one synthesized by a new process that uses a solventogenic Clostridium strain through an extractive fermentation process using n-hexadecane as the extractant (SBA); the other commercially available as a high-purity (99.9%) 'neat' BA grade produced by Fischer esterification (NBA). The mass concentration of by-products from the synthesis process amounted to approximately 6% of the total composition and included n-butanol, n-hexadecane, iso-propyl alcohol and ethyl acetate. The platform used to study the combustion of commercial and synthesized BA was an isolated droplet burning with spherical symmetry. Initial droplet diameters were fixed at 0.6 mm. Experiments were carried out in a drop tower to promote one-dimensional transport dynamics, in the standard atmosphere, and with ignition accomplished by spark discharge. Differences in droplet burning rates and flame structures were undetectable between NBA and SBA, despite the differences in burning among the SBA components individually. The results presented show that the new synthesis process for butyl acetate yields a sustainable alternative to conventional methodologies with burning characteristics that are identical to NBA in stagnant gas transport fields.

Wang, Yujie↗
DOE OSTI2022

Resilience in soil bacterial communities of the boreal forest from one to five years after wildfire across a severity gradient

Wildfires can represent a major disturbance to ecosystems, including soil microbial communities belowground. Furthermore, fire regimes are changing in many parts of the world, altering and often increasing fire severity, frequency, and size. The boreal forest and taiga plains ecoregions of northern Canada are characterized by naturally-occurring stand-replacing wildfires on a 40–350 year basis. We previously studied the effects of wildfire on soil microbial communities one year post-fire across 40 sites, spanning a range of burn severity. Here, we return to the same sites five years post-fire to test a series of hypotheses about the effects of fire on bacterial community composition. We ask questions on two themes: which factors control bacterial community composition during post-fire recovery, and how does the importance of different fire-responsive traits change during post-fire recovery? We find the following: Five years post-fire, vegetation community, moisture regime, pH, total carbon, texture, and burned/unburned all remained significant predictors of bacterial community composition with similar predictive value (R 2 ). Bacterial communities became more similar to unburned sites five years post-fire, across the range of severity, suggesting resilience, while general structure of co-occurrence networks remained similar one and five years post-fire. Fast growth potential, as estimated using predicted 16S rRNA copy numbers, was no longer significantly correlated with burn severity five years post-fire, indicating the importance of this trait for structuring bacterial community composition may be limited to relatively short timescales. Many taxa that were enriched in burned sites one year post-fire remained enriched five years post-fire, although the degree to which they were enriched generally decreased. Specific taxa of interest from the genera Massilia, Blastococcus, and Arthrobacter all remained significantly enriched, suggesting that they may have traits that allow them to continue to flourish in the post-fire environment, such as tolerance to increased pH or ability to degrade pyrogenic organic matter. This hypothesis-based work expands our understanding of the post-fire recovery of soil bacterial communities and raises new hypotheses to test in future studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Modeling Volatility-Based Aerosol Phase State Predictions in the Amazon Rainforest

Organic aerosol (OA) is a complex matrix of various constituents—fresh (primary organic aerosols—POA) and aged via oxidation (secondary organic aerosols—SOA), generated from biogenic, anthropogenic, and biomass burning sources. The viscosity of OA can be critical in influencing new particle formation, reactive uptake processes that impact evaporation-growth kinetics, and the lifetime of particles in the atmosphere. We report this work utilizes a well-defined relationship between volatility and viscosity for pure compounds, which we incorporated within the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to simulate the phase state and viscosity of bulk OA during the dry-to-wet transition season (September–October) in the Amazon rainforest during 2014. Our simulations indicate spatial and temporal heterogeneity in aerosol phase state often not captured by global-scale models. We show the strong role of water associated with organic aerosol ($w_s$) as the dominant factor that can be used to quantitatively estimate OA viscosity. Analysis of WRF-Chem simulations across the entire atmospheric column indicates a strong inverse log-linear relationship between $w_s$ and OA viscosity with a correlation coefficient approaching 1, in the background and biomass burning-influenced conditions. At high altitudes where relative humidity (RH) and temperatures are low, our simulations indicate that OA exists in a semisolid-/solid-like phase state, consistent with previous studies. OA hygroscopicity is strongly correlated (ca. –0.8) with OA viscosity at RH ca. 30–50%, but this RH range is found mostly at low OA concentrations and the middle troposphere (ca. 6–10 km altitudes) in our simulated domain. OA hygroscopicity is uncorrelated with viscosity at higher-RH (near surface) and lower-RH (upper troposphere) regimes. At the urban site near surface, where day–night differences in RH are significant, RH is found to drive the phase state. At the background forested site near surface, where day–night RH differences are small, biomass burning-influenced OA is semisolid and a significant OA associated with background conditions is liquid-like. Simulations indicate a long tail of OA viscosity frequency distributions extending in the semisolid/solid regimes over background biogenic-influenced conditions due to the role of low-volatility OA components such as monoterpene oxidation products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Enhancing fire emissions inventories for acute health effects studies: integrating high spatial and temporal resolution data

Daily fire progression information is crucial for public health studies that examine the relationship between population-level smoke exposures and subsequent health events. Issues with remote sensing used in fire emissions inventories (FEI) lead to the possibility of missed exposures that impact the results of acute health effects studies. This paper provides a method for improving an FEI dataset with readily available information to create a more robust dataset with daily fire progression. High temporal and spatial resolution burned area information from two FEI products are combined into a single dataset, and a linear regression model fills gaps in daily fire progression. The combined dataset provides up to 71% more PM 2.5 emissions, 69% more burned area, and 367% more fire days per year than using a single source of burned area information. The FEI combination method results in improved FEI information with no gaps in daily fire emissions estimates. The combined dataset provides a functional improvement to FEI data that can be achieved with currently available data.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2024

Validation of prediction capability of operating space for plasma initiation in MAST-U

DYON is a plasma initiation modelling code that solves the differential equation system of the full circuit equations (plasma current, active coil currents and eddy currents in full passive structures) and 0D global energy and particle balance equations (Kim 2022 Nucl. Fusion 62 126012). In order to test the capability of the full electromagnetic plasma initiation model to predict individual discharges in experiments and thus the operating space in the device, a dedicated experimental database was built in MAST-U by scanning the prefilled gas pressure p 0 and the induced loop voltage V loop . In the experimental operating space of p 0 and V loop the lower and the upper limits of p 0 are determined by the plasma breakdown failure and the plasma burn-through failure, respectively. The lower limit of V loop is determined by the plasma burn-through failure. By directly reading the control room data used in each discharge (i.e. currents in the solenoid, poloidal field coils, and toroidal field coils, p 0 , and gas puffing rate), the full electromagnetic DYON consistently predicted the failed breakdown, failed burn-through, and successful plasma initiation discharges in the experimental database, demonstrating its capability to predict the operating space for inductive plasma initiation. The Paschen curve calculated with the effective connection length in MAST-U indicates a much higher p 0 required for plasma breakdown than the experimental data, indicating that individual field line evaluation is necessary to calculate the quantitative requirements for Townsend breakdown. The demonstration in this paper shows that the full electromagnetic DYON could be a useful simulation tool to assess the feasibility of inductive plasma initiation and to optimise operating scenarios in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Data and scripts associated with the manuscript evaluating the hydrologic responses of the Pacific Northwest watersheds to wildfires (v2)

This data package is associated with the publication “Evaluating Post-fire Watershed Response to Varying Burn Severity and Precipitation Regimes Using Fully-distributed and Integrated Hydrologic Models” submitted to Journal of Hydrology (Li et al. 2025). In this study, we employed the Advanced Terrestrial Simulator (ATS), an integrated watershed model that couples surface flow, subsurface flow, and canopy biophysical processes, to investigate post-fire hydrologic responses in a few selected watersheds with varying burn severity.The data package contains the required input data (meteorological forcing, Leaf Area Index, wildfire burn severities, etc.) to run the model, configuration files, the Jupyter notebooks in Python to pre-process and post-process data, the figures in the manuscript, and the modeling output files. The variables include watershed-averaged evapotranspiration, watershed-averaged surface/subsurface/canopy water content, and river discharge at watershed outlet.The data package contains a file-level metadata that lists and describes all the files contained in the data package (ATS_flmd.csv), a data dictionary file that defines columns headers across all csv files contained in the data package (ATS_dd.csv), a data package level readme file (the current file), and four zipped folders.The ‘data’ folder provides data needed to run the model in .h5, .i2s, .xyz, .shp, and .exo formats. The sub-folders are for each data types. The ‘model’ folder provides input files (.xml format) and essential model outputs. Each sub-folder provides the files from each simulated watershed. The ‘notebooks’ folder provides the Jupyter notebooks (.ipynb format) for pre- and post- processing model files, and for producing the figures in the manuscript. The ‘figures’ folder provides the figures associated with manuscript in .pdf and .png formats.The ‘model’ folder and the ‘data’ folder have been split into 5GB-large pieces using the Linux command ‘split -b 5120m model.zip model.zip.’ and ‘split -b 5120m data.zip data.zip.’, respectively. They can be merged back using the Linux command ‘cat model.zip.* > model.zip’ and ‘cat data.zip.* > data.zip’, respectively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Post-fire time series of sensor and geochemistry sample data from surface water, groundwater, precipitation, soil, and vegetation across Oak Creek watershed, Washington

This dataset supports a broader study examining wildfire impacts on hydrologic connectivity across 5 sites within the Oak Creek watershed and the resulting biogeochemical impacts. Stream sites were selected using the Advanced Terrestrial Simulator (ATS) hydrologic model to identify locations with varying groundwater contributions and hydrologic responses across different burn severity scenarios. The Retreat Fire burned from July 23 to August 2 in 2024, affecting the five study sites at varying burn severities. Each site is equipped with YSI EXO2 sondes logging sub-hourly throughout the year, and grab samples are collected approximately every six weeks. YSI sondes are used to measure temporally resolved proxies for groundwater inputs (specific conductivity) and organic matter (fluorescent dissolved organic matter; fDOM) along with basic water quality and depth. Grab samples of surface water, groundwater, and precipitation are analyzed for water stable isotopes and conductivity to understand endmembers for hydrologic mixing Grab samples of surface water, groundwater, soil water, and litter/vegetation/soil leachates are analyzed for organic matter composition measured by Fourier-Transform Ion Cyclotron Resonance Mass Spectrometry (FTICR-MS) to understand organic matter dynamics. Game camera photos are provided in a separate data package available at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3018598. Future versions of this dataset will include time series data from YSI EXO2 sondes (fDOM, dissolved oxygen, temperature, depth, specific conductance, turbidity, pH), BaroTROLL sensors (air temperature and barometric pressure), rain gauges (precipitation), and data from the soil and vegetation samples. Because this study is ongoing, this data package will be updated regularly to include newly collected data and the additional data types. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) a folder of field photos; (2) a folder of surface water sample data; (3) a folder of raw Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) data; (4) a data checks report; (5) file-level metadata; (6) data dictionary; (7) field metadata; (8) readme; (9) international generic sample number (IGSN) mapping file; and (10) field protocols. The sample data subfolder contains (1) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data and averages; (2) total dissolved nitrogen data and averages; (3) stable water isotopes and averages; (4) methods codes; (5) FTICR-MS methods; and (15) a subfolder of 9.4 Tesla (9.4T) FTICR-MS data. This folder contains the processed data and three subfolders, one containing the .xml files, one containing the water CoreMS output files, and the other containing instructions and scripts for processing the files in CoreMS (https://github.com/EMSL-Computing/CoreMS). All files are .csv, .pdf, .R, .xml, .d, .html, .Rmd, .py, .cal, .json, .jpg, .jpeg, .png, .mov, or .mp4.

Biogeochemistry↗
DOE OSTI2026

Data and scripts associated with a manuscript analyzing ELM-FATES parameter sensitivity under pre-fire and postfire scenarios using machine learning

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “Fire Severity-Dependent Shifts in Vegetation Parameter Sensitivity: A Pre- and Post-Fire Analysis Using ELM-FATES and Explainable AI” submitted to Journal of Advances in Modeling Earth Systems (Zahura et al. 2026). The study examines vegetation physiological parameters controlling pre-fire and post-fire vegetation dynamics. To support this analysis, 73 vegetation parameters in Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (Fisher et al., 2018) , which is coupled with E3SM (Energy Exascale Earth System Model) land model (ELM, ELM-FATES), were perturbed using a Sobol sequence to generate 1,024 ensemble members for two plant functional types: needleleaf evergreen extratropical trees (NEET) and C3 grass. Simulations were conducted for the pre-fire period (2016) and post-fire period (2018–2023). Burn severity was represented by modifying the Nesterov index in FATES to 75,000, 150,000, and 300,000 for low, moderate, and high severity, respectively. A no-fire scenario was also included. Simulations were performed for 16 grid cells in the American River Watershed across different burn severities and plant functional types. XGBoost (eXtreme Gradient Boosting) models were trained using the parameter ensembles and ELM-FATES-simulated outputs, including leaf area index (LAI), gross primary productivity (GPP), aboveground biomass, vegetation evaporation, transpiration, and soil evaporation. Models were trained separately for each year and burn severity, followed by SHAP (SHapley Additive exPlanations) analysis to identify changes in dominant parameters after fire disturbance. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package contains the ELM-FATES simulation data. The scripts and data related to the analysis will be added later. The inputs and outputs from ELM-FATES are inside the “FATES” folder. “FATES_domain_surface” contains the domain and surface netcdfs that were used to run ELM-FATES in the study area. “FATES_parameters” contains the 1024 ensembles that were generated using Sobol sequence. “FATES_outputs” folder contains ELM-FATES simulated variables. All files are .csv and .nc (NetCDF).

Aboveground biomass↗