Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distribution factor”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The role of the turbulence driving mode for the initial mass function

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating, and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ∼1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Influence of Ink Formulation and Drying Conditions on Ionomer Distribution in High-Performance Roll-to-Roll-Coated Gas-Diffusion Electrodes

To enable mass production of fuel cell membrane electrode assemblies (MEAs) catalyst layers production will require continuous roll-to-roll (R2R) coating processes. Gas diffusion electrodes (GDEs) are advantageous for mass production because the catalyst layer can be directely coated on the microporous layer of the gas diffusion media without the need for a decal-transfer process. It is known that the water-to-alcohol ratio in the catalyst ink influences the interactions of the ionomer with the catalyst leading to different distributions of ionomer in spray-coated catalyst layers. It is also known that during drying of colloidal mixtures, like fuel cell inks, factors such as drying rate, particle size, and agglomeration influence how the materials distribute themselves throughout the thickness of the dired film. Thus far there have only been limited studies to understand how process conditions such as ink formulation and drying temperature influence the distribution of ionomer and catalyst coated using scalable methods. This understanding is especially important for GDEs since it is known that having a sufficient amount of ionomer at the catalyst layer-membrane interface is critical for high performance. In this study we have focused on determining how the ratio of water to 1-propanol in the catalyst ink ink and drying temperature influence the distribution of ionomer throughout the thickness of the catalyst layer. Using a combination of Kelvin probe and x-ray photoelectron spectroscopy we show that an ionomer-rich surface is promoted by a higher drying rate and a water-rich catalyst ink. In contrast, a 1-propanol catalyst ink leads to a lower concentration of ionomer on the top surface. Using x-ray computed tomography, we are able to characterize the ionomer distribution throughout the thickness of the layer. We find that, in addition to promoting an ionomer-rich top surface, water-rich inks lead to a more homogenous distribution of ionomer, whereas a 1-propanol-rich ink leads to a more irregular distribution. It is found that MEA performance is improved by selecting conditions and ink formulations that promote ionomer enrichment at the top surface to facilitate a good interface with the membrane. MEAs prepared with a 75 wt% water catalyst ink with a 0.9 I/C have equivalent performance to spray-coated GDEs. Critically, these R2R-coated GDEs do not need an additional ionomer overlayer like the spray-coated GDEs do, reducing the number of processing steps in a manufacturing setting. This work shows that with the appropriate selection of materials, ink formulation, and processing conditions gas-diffusion electrodes are a viable pathway for fuel cell manufacturing.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗
DOE OSTI2025

Using convolutional neural networks to accelerate three-dimensional coherent synchrotron radiation computations

Calculating the effects of coherent synchrotron radiation (CSR) is one of the most computationally expensive tasks in accelerator physics. Here, we use convolutional neural networks (CNNs), along with a latent conditional diffusion (LCD) model, trained on physics-based simulations to speed up calculations. Specifically, we produce the 3D CSR wakefields generated by electron bunches in circular orbit in the steady-state condition. Two datasets are used for training and testing the models: wakefields generated by three-dimensional Gaussian electron distributions and wakefields from a sum of up to 25 three-dimensional Gaussian distributions. The CNNs are able to accurately produce the 3D wakefields ∼250–1000 times faster than the numerical calculations, while the LCD achieves a gain of a factor of ∼34. We also test the extrapolation and out-of-distribution generalization ability of the models. They generalize well on distributions with larger spreads than what they were trained on but struggle with smaller spreads.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Quantifying Drivers of Methane Hydrobiogeochemistry in a Tidal River Floodplain System

The influence of coastal ecosystems on global greenhouse gas (GHG) budgets and their response to increasing inundation and salinization remains poorly constrained. In this study, we have integrated an uncertainty quantification (UQ) and ensemble machine learning (ML) framework to identify and rank the most influential processes, properties, and conditions controlling methane behavior in a freshwater floodplain responding to recently restored seawater inundation. Our unique multivariate, multiyear, and multi-site dataset comprises tidal creek and floodplain porewater observations encompassing water level, salinity, pH, temperature, dissolved oxygen (DO), dissolved organic carbon (DOC), total dissolved nitrogen (TDN), partial pressure of carbon dioxide (pCO 2 ), nitrous oxide (pN 2 O), methane (pCH 4 ), and the stable isotopic composition of methane (δ 13 CH 4 ). Additionally, we incorporated topographical data, soil porosity, hydraulic conductivity, and water retention parameters for UQ analysis using a previously developed 3D variably saturated flow and transport floodplain model for a physical mechanistic understanding of factors influencing groundwater levels and salinity and, therefore, CH 4 . Principal component analysis revealed that groundwater level and salinity are the most significant predictors of overall biogeochemical variability. The ensemble ML models and UQ analyses identified DO, water level, salinity, and temperature as the most influential factors for porewater methane levels and indicated that approximately 80% of the total variability in hourly water levels and around 60% of the total variability in hourly salinity can be explained by permeability, creek water level, and two van Genuchten water retention function parameters: the air-entry suction parameter α and the pore size distribution parameter m. These findings provide insights on the physicochemical factors in methane behavior in coastal ecosystems and their representation in local- to global-scale Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Radionuclide surrogate aerosolization, resuspension and suppression in hazardous situations

Radioactive fire following nuclear accidents causes resuspension of radioactive contamination that can spread over large areas. Although many papers focus on particle deposition and resuspension, a limited number of studies model its suppression. Even less data is available for the comparison of emissions across particle size, surface, and mode of resuspension. In this study, cerium dioxide surrogate for internal hazards such as radionuclide particle matter was aerosolized into a large chamber and subsequently resuspended after sedimentation. Resuspension was induced by common modes of human movement including walking, wind, and driving on common surfaces including polyvinyl chloride (PVC), artificial grass (turf), and concrete in a 5.8 m long, 3.7 m wide and 3 m high tent. The resuspension was experimentally measured for two particle sizes (1 and 10 μm) at three different heights (0.3 m, 1.2 m, and 2.4 m) and compared across size, surface, and mode of resuspension. Particle size distribution was measured using an aerodynamic particle sizer and the resuspended particles were collected with the wetted wall cyclone aerosol collector and quantitated using inductively coupled plasma mass spectrometry. The testing was repeated with the F-500 fire retardant application to study its effect on resuspension. In this study, artificial grass flooring was found to generate the most resuspension, with walking on turf for the 1 μm particle testing yielding the highest factor overall. Concrete had the lowest resuspension, showing variations in the size distribution of aerosol as a function of height from the source resuspension factors. The application of fire retardant was found to significantly reduce surrogate resuspension, regardless of resuspension method or testing surface. Here, the results provide impactful information regarding the different modes of resuspension including human transport (walking and driving) and wind for the resuspension of radionuclide surrogate particles deposited onto common surfaces and the effect of fire retardant application on their resuspension, leading to higher levels of safety for nuclear energy plants and facilities.

42 ENGINEERING↗
DOE OSTI2024

Precision Measurements of the Neutron Magnetic Form Factor to High Momentum Transfer using Durand’s Method

Protons and neutrons, collectively known as nucleons, along with electrons, constitute the fundamental building blocks of the visible universe. Understanding their internal structure is crucial for addressing key scientific questions about our origin and existence. Elastic electron-nucleon scattering provides insights into the spatial distributions of charge and current within nucleons through their electromagnetic form factors. Accurate knowledge of these form factors over a broad range of Q2, the squared four-momentum transfer in the scattering process, reveals details about the nucleon's internal structure. However, high-Q2 data of the nucleon electromagnetic form factor is scarce due to the challenges associated with such measurements. This thesis reports preliminary results from high-precision measurements of the neutron magnetic form factor (GMn) to unprecedented Q2 using Durand's method, also known as the "ratio" method. Systematic errors are greatly reduced by ext

Datta, Provakar↗
DOE OSTI2022

DSO+T: Valuation Methodology and Economic Metrics (DSO+T Study: Volume 4)

This report summarizes a rigorous valuation analysis methodology used by the Distribution System Operator with Transactive (DSO+T) study to estimate the financial benefits and costs of adopting Transactive Energy coordination of distributed energy resources for key stakeholders (for example distribution system operators and customers). This was achieved by modeling the value exchanges between stakeholders and determining the annualized costs and revenues experienced by stakeholders, enabling the evaluation of overall impact on stakeholder’s annualized cash flow. Extensive work was conducted developing methods to estimate the operating costs of distribution system operators at a level of granularity that would allow the financial impact of implementing a transactive energy approach to be estimated. This work included developing parametric models for labor and software costs, distribution system capital and maintenance costs, growth rates, and factors to determine annualized costs of capital investments. Simulation results were used to calculate wholesale energy costs and revenues from retail sales. Valuation analysis methods were also developed for other stakeholders including customers, the transmission system operator, independent system operator, and generators. The resulting capability allows a complete mapping of the flow of financial value between stakeholders that can be directly integrated with the results of demand flexibility simulations. Example results are provided for a business-as-usual case and compared to a transactive energy case as well as to actual cost data.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

High-volume tunable resonator for axion searches above 7 GHz

In this research, we present results from an experimental demonstration of a tunable thin-shell axion haloscope whose geometry decouples its overall volume from its resonant frequency, thereby evading the steep sensitivity degradation at high frequencies. An aluminum 2.6-l (41⁢λ 3 ) prototype, which tunes from 7.1 to 8.0 GHz, was fabricated and characterized at room temperature. An axion-sensitive, straightforwardly tunable TM 010 mode is clearly identified with a room-temperature quality factor, Q, of approximately 5000. The on-resonance E-field distribution is mapped and found to agree with numerical calculations. Anticipating future cryogenic operation, we develop an alignment protocol relying only on rf measurements of the cavity, maintaining a form factor of 0.57 across the full tuning range. These measurements demonstrate the feasibility of cavity-based haloscopes with operating volume V $\gg$ λ 3 . We discuss plans for future development and the parameters required for a thin-shell haloscope exploring the postinflationary axion parameter space (approximately 4 to 30 GHz) at Dine-Fischler-Srednicki-Zhitnitsky sensitivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

β + Gamow-Teller Strengths from Unstable 14 O via the (d, 2 He) Reaction in Inverse Kinematics

For the first time, the (d, 2 He) reaction was successfully used in inverse kinematics to extract the Gamow-Teller transition strength in the β + direction from an unstable nucleus. The new technique was made possible by the use of an active-target time-projection chamber and a magnetic spectrometer, and opens a path to addressing a range of scientific challenges, including in astrophysics and neutrino physics. Here, in this Letter, the nucleus studied was 14 O, and the Gamow-Teller transition strength to 14 N was extracted up to an excitation energy of 22 MeV. The data were compared to shell-model and state-of-the-art coupled-cluster calculations. Shell-model calculations reproduce the measured Gamow-Teller strength distribution up to about 15 MeV reasonably well, after the application of a phenomenological quenching factor. In a significant step forward to better understand this quenching, the coupled-cluster calculation reproduces the full strength distribution well without such quenching, owing to the large model space, the inclusion of strong correlations, and the coupling of the weak interaction to two nucleons through two-body currents.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Precision Measurements of the Neutron Magnetic Form Factor to High Momentum Transfer using Durand's Method

Protons and neutrons, collectively known as nucleons, along with electrons, constitute the funda- mental building blocks of the visible universe. Understanding their internal structure is crucial for addressing key scientific questions about our origin and existence. Elastic electron-nucleon scatter- ing provides insights into the spatial distributions of charge and current within nucleons through their electromagnetic form factors. Accurate knowledge of these form factors over a broad range of Q2, the squared four-momentum transfer in the scattering process, reveals details about the nucleon’s internal structure. However, high-Q2 data of the nucleon electromagnetic form factor is scarce due to the challenges associated with such measurements. This thesis reports preliminary results from high-precision measurements of the neutron magnetic form factor (Gn M ) to unprecedented Q2 using Durand’s method, also known as the “ratio” method. Systematic errors are greatly reduced by extracting Gn M from the ratio of neutron-coincident (D(e, e'n)) to proton-coincident (D(e, e'p)) quasi-elastic electron scattering from deuteron. The scattered electrons were detected in the BigBite spectrometer, which features multiple Gas Elec- tron Multiplier (GEM) layers with large active area for high-precision tracking at very high rates. Simultaneous nucleon detection was performed by the Super BigBite spectrometer, which utilizes a dipole magnet with large solid angle acceptance at forward angles and a novel hadron calorimeter with very high and comparable detection efficiencies for both protons and neutrons. This setup could handle very high luminosity, making high-Q2 measurements feasible. Data were collected at five Q2 points: 3, 4.5, 7.4, 9.9, and 13.6 (GeV/c)2. Preliminary results are reported for all, with the lowest two Q2 points in good agreement with existing world data, while the higher points significantly extend the Q2 range in which Gn M is known accurately. The precision of the highest Q2 point is expected to remain unmatched for years to come.

Datta, Provakar↗
DOE OSTI2021

PV Fleet Performance Data Initiative: Performance Index-Based Analysis

In this analysis, we report on 250 PV systems throughout the US, comprising 157 MWdc of system capacity and over 10,000 monthly PI values based on high-frequency (subhourly) energy data and satellite-based resource data. The distribution of PI values is analyzed, and multiple causes of underperformance are assessed, including first-year startup issues, snowfall and inverter downtime. An initial distribution of raw monthly PI values was collected with mean measured / modeled performance of PI = 0.935. After correcting for the three identifiable loss factors mentioned above, an average monthly performance of PI = 0.994 was obtained, with a distribution closely following a Gumbel Extreme Value distribution. In particular, inverter availability was found to contribute a system energy loss of 2.3% on average, except in the first six months of operation when availability losses are closer to 8%. Other startup issues beyond inverter downtime such as partial string outage contributed an additional 5% underperformance in the first year of operation across the fleet. Winter performance was also found to be 5%-10% below summer performance on average, likely due to snowfall. A simple linear relationship was found between snow loss and monthly snowfall accumulation in cm, indicating between 6% - 40% loss in months with nonzero snowfall, depending on climate.

14 SOLAR ENERGY↗
DOE OSTI2022

A novel formulation of the unintegrated gluon distribution for DIS

In this talk, we revisit inclusive DIS in the small x limit and derive a new factorization formula that accounts for leading powers in both Bjorken and Regge limits. In this semi-classical description, we obtain a new unintegrated gluon distribution which encompasses both the dipole operator and the gluon Parton Distribution Function with an explicit dependence on the longitudinal momentum fraction x.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Modification of Quark-Gluon Distributions in Nuclei by Correlated Nucleon Pairs

We extend the QCD Parton Model analysis using a factorized nuclear structure model incorporating individual nucleons and pairs of correlated nucleons. Our analysis of high-energy data from lepton deep-inelastic scattering, Drell-Yan, and W and Z boson production simultaneously extracts the universal effective distribution of quarks and gluons inside correlated nucleon pairs, and their nucleus-specific fractions. Such successful extraction of these universal distributions marks a significant advance in our understanding of nuclear structure properties connecting nucleon- and parton-level quantities. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Gluon splitting at small x: a unified derivation for the JIMWLK, DGLAP and CSS equations

We revisit the calculation of the next-to-leading order (NLO) corrections to dijet production in electron-ion collisions at small x. We focus on the back-to-back configuration where the relative transverse momentum P ⊥ of the measured jets is much larger than both their momentum imbalance K ⊥ and the target saturation momentum Q s (x, A). In this regime, we present for the first time a complete calculation of the real NLO corrections, associated with gluon emissions outside the jet cones, at leading power in 1/P ⊥ . Our result exhibits TMD factorisation, with the same hard factor as at tree-level and a NLO correction to the Weiszäcker-Williams (WW) gluon transverse momentum dependent (TMD) distribution which involves four Wilson-line operators. By studying different kinematical regimes for K ⊥ and for the radiated gluon, we recover all the quantum evolutions that were previously identified for this process at NLO: the B-JIMWLK high-energy evolution and the CSS evolution of the gluon WW TMD, and the DGLAP evolution of the gluon PDF. When both K ⊥ and the transverse momentum transferred by the target are large compared to Q s , all the Wilson-line operators boil down to the unintegrated gluon distribution and our NLO result for the gluon TMD can be used to isolate the transverse-momentum dependent gluon splitting function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Distributed Magnetic Field and Temperature Monitoring for Superconducting Radio Frequency Cavities

The overall objective of the proposed Phase I program was to design, construct and demonstrate a fiber optic sensing system capable of providing temperature and magnetic field measurements with an enhanced spatial resolution that can be implemented over a large surface area (cryomodules) to survey superconducting radio frequency cavities and magnets. A magnetic field sensor capable of detecting fluxes on the order of 1 μT is required to detect the distribution of trapped flux on the cavity surface. A unique distributed magnetic field sensor was successfully designed and constructed to demonstrate the detection of magnetic fluxes less than 500 nT. The sensor leveraged the ultra-high sensitivity of Sentek’s picoDAS to measure the magnetostriction induced vibrations in a commercially available Metglas 2605 SC ribbon that was in physical contact with sensing fiber. Static magnetic fields were detected by applying an alternating current a copper wire proximate to the Metglas 2605SC ribbon to create an AC bias magnetic field. In an alternative approach, an AC bias magnetic field was applied to a special magnetic field sensing fiber with Metglas 2605SC cladding successfully detect a magnetic field flux of a 3 μT. Exhaustive testing was performed to characterize the dependency of sensor response on the direction of the applied magnetic field. Although the special sensing fiber based magnetic field sensor did not exhibit an observable dependence on the direction of the magnetic field, the Metglas 2605SC ribbon sensor exhibited a clear directional dependence. A wide variety of polymer materials were evaluated to enhance the temperature response of an FBG based sensor at cryogenic temperatures. The processing and performance challenges provided the motivation to develop a new simple cryogenic temperature sensor that uses a commercially available fiber optic splice protector. The EVA hot melt tube that becomes adhered to the optical fiber and the polyolefin outer tube that shrinks upon heating in the fusion splicer heater provide the high thermal expansion coefficient necessary to impart a significant strain on the FBG when exposed to cryogenic temperatures. The temperature sensitivity (Δ𝜆𝐵𝑟𝑎𝑔𝑔~ 62 𝑝𝑝𝑚/℃) of the FBG-based sensor was on par with the best reported to date. The simple design, use of readily available cost-effective materials, and well-established processing techniques lends this approach to the creation of hundreds to thousands of temperature sensors on one single optical fiber length. The inherently small form factor also allows for co-location with the distributed magnetic field sensor. In preparation for field testing of the prototype sensing system at the Jefferson Labs in potential Phase II program, several different cable designs were evaluated to package the sensors. The preliminary successful demonstration of fully functional sensing cables provides the foundation for subsequent development efforts to advance the Technology Readiness Level of the technology. The technical feasibility of the proposed approach was successfully demonstrated in this Phase I effort.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Integrated Land Suitability Assessment for Depots Siting in a Sustainable Biomass Supply Chain

A sustainable biomass supply chain would require not only an effective and fluid transportation system with a reduced carbon footprint and costs, but also good soil characteristics ensuring durable biomass feedstock presence. Unlike existing approaches that fail to account for ecological factors, this work integrates ecological as well as economic factors for developing sustainable supply chain development. For feedstock to be sustainably supplied, it necessitates adequate environmental conditions, which need to be captured in supply chain analysis. Using geospatial data and heuristics, we present an integrated framework that models biomass production suitability, capturing the economic aspect via transportation network analysis and the environmental aspect via ecological indicators. Production suitability is estimated using scores, considering both ecological factors and road transportation networks. These factors include land cover/crop rotation, slope, soil properties (productivity, soil texture, and erodibility factor) and water availability. This scoring determines the spatial distribution of depots with priority to fields scoring the highest. Two methods for depot selection are presented using graph theory and a clustering algorithm to benefit from contextualized insights from both and potentially gain a more comprehensive understanding of biomass supply chain designs. Graph theory, via the clustering coefficient, helps determine dense areas in the network and indicate the most appropriate location for a depot. Clustering algorithm, via K-means, helps form clusters and determine the depot location at the center of these clusters. An application of this innovative concept is performed on a case study in the US South Atlantic, in the Piedmont region, determining distance traveled and depot locations, with implications on supply chain design. The findings from this study show that a more decentralized depot-based supply chain design with 3depots, obtained using the graph theory method, can be more economical and environmentally friendly compared to a design obtained from the clustering algorithm method with 2 depots. In the former, the distance from fields to depots totals 801,031,476 miles, while in the latter, it adds up to 1,037,606,072 miles, which represents about 30% more distance covered for feedstock transportation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Microscale spatial distribution and soil organic matter persistence in top and subsoil

The spatial distribution of organic substrates and microscale soil heterogeneity significantly influence organic matter (OM) persistence as constraints on OM accessibility to microorganisms. However, it is unclear how changes in OM spatial heterogeneity driven by factors such as soil depth affect the relative importance of substrate spatial distribution on OM persistence. This work evaluated the decomposition and persistence of 13 C and 15 N labeled water-extractable OM inputs over 50 days as either hotspot (i.e., pelleted in 1–2 mm-size pieces) or distributed (i.e., added as OM < 0.07 μm suspended in water) forms in topsoil (0–0.2 m) and subsoil (0.8–0.9 m) samples of an Andisol. We observed greater persistence of added C in the subsoil with distributed OM inputs relative to hotspot OM, indicated by a 17% reduction in cumulative mineralization of the added C and a 10% higher conversion to mineral-associated OM. A lower substrate availability potentially reduced mineralization due to OM dispersion throughout the soil. NanoSIMS (nanoscale secondary ion mass spectrometry) analysis identified organo-mineral associations on cross-sectioned aggregate interiors in the subsoil. On the other hand, in the topsoil, we did not observe significant differences in the persistence of OM, suggesting that the large amounts of particulate OM already present in the soil outweighed the influence of added OM spatial distribution. Here, we demonstrated under laboratory conditions that the spatial distribution of fresh OM input alone significantly affected the decomposition and persistence of OM inputs in the subsoil. On the other hand, spatial distribution seems to play a lower role in topsoils rich in particulate OM. The divergence in the influence of OM spatial distribution between the top and subsoil is likely driven by differences in soil mineralogy and OM composition.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Predicting Emission Source Terms in a Reduced-Order Fire Spread Model—Part 1: Particulate Emissions

A simple, easy-to-evaluate, surrogate model was developed for predicting the particle emission source term in wildfire simulations. In creating this model, we conceptualized wildfire as a series of flamelets, and using this concept of flamelets, we developed a one-dimensional model to represent the structure of these flamelets which then could be used to simulate the evolution of a single flamelet. A previously developed soot model was executed within this flamelet simulation which could produce a particle size distribution. Executing this flamelet simulation 1200 times with varying conditions created a data set of emitted particle size distributions to which simple rational equations could be tuned to predict a particle emission factor, mean particle size, and standard deviation of particle sizes. These surrogate models (the rational equation) were implemented into a reduced-order fire spread model, QUIC-Fire. Using QUIC-Fire, an ensemble of simulations were executed for grassland fires, southeast U.S. conifer forests, and western mountain conifer forests. Resulting emission factors from this ensemble were compared against field data for these fire classes with promising results. Also shown is a predicted averaged resulting particle size distribution with the bulk of particles produced to be on the order of 1 μm in size.

54 ENVIRONMENTAL SCIENCES↗