Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase”

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

High-fidelity pool boiling simulations on multiple nucleation sites using interface capturing method

Boiling has proved to be one of the most efficient means for heat transfer and is a very important phenomenon during severe accident scenarios in light water reactors. High-fidelity pool boiling simulations can provide a numerical database for improving mechanistic boiling models by allowing for specific evaluation of interactions among bubbles. Previously published pool boiling simulations investigated two nucleation sites in which bubble growth at one site suppressed nucleation at the other site. Based on previous study results, more complicated interface-capturing simulations on pool boiling were conducted using PHASTA code with locally refined unstructured mesh. First, different boundary conditions (BCs) were assessed to support robustness and reproducibility of the boiling model. Then, a scale study was conducted at a larger domain with nine nucleation sites where either nine or four nucleation sites are activated. Involving more nucleation sites increased the complexity of bubble interactions from surrounding sites. Finally, bubble departure behavior influenced by wall heat flux was investigated. In this work, when heat flux was increased, the order of bubble departure changed, but diagonal bubbles always departed after one another. The departure time interval between the first and second bubble reduced as heat flux increased. The corresponding frequency was almost linearly proportional to the heat flux. In addition, bubble departure behavior was found to be greatly influenced by the nucleation site pattern. Multiple nucleation sites resulted in superimposed inhibitive effects from surrounding sites to each bubble, which extensively delayed the departure. This new observation was not discussed in previously published works. The work presented here provides new insight on the fundamental understanding of boiling phenomena, contributes to the development of a 3D multiphase computational fluid dynamics (M-CFD) model, and provides a more comprehensive database for data-driven pool boiling studies.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗
DOE OSTI2021

Multiplicity scaling of light nuclei production in relativistic heavy-ion collisions

Using the nucleon coalescence model based on kinetic freeze-out nucleons from the 3D MUSIC+UrQMD and the 2D VISHNU hybrid model with a crossover equation of state, we study the multiplicity dependence of deuteron ( d ) and triton ( t ) production from central to peripheral Au+Au collisions at s NN = 7.7, 14.5, 19.6, 27, 39, 62.4 and 200 GeV and Pb+Pb at s NN = 2.76 TeV, respectively. It is found that the ratio N t N p / N d 2 of the proton yield N p , deuteron yield N d and triton yield N t exhibits a scaling behavior in its multiplicity dependence, i.e., decreasing monotonically with increasing charged-particle multiplicity. A similar multiplicity scaling of this ratio is also found in the nucleon coalescence calculation based on kinetic freeze-out nucleons from a multiphase transport (AMPT) model. The scaling behavior of N t N p / N d 2 can be naturally explained by the interplay between the sizes of light nuclei and the nucleon emission source. We further argue that the multiplicity scaling of N t N p / N d 2 can be used to validate the production mechanism of light nuclei, and the collision energy dependence of this yield ratio can further serve as a baseline in the search for the QCD critical point in relativistic heavy-ion collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Impact of globally spin-aligned vector mesons on the search for the chiral magnetic effect in heavy-ion collisions

In high-energy heavy-ion collisions, the chiral magnetic effect (CME) is predicted to arise from the interplay between the chirality imbalance of quarks in the nuclear medium and the intense magnetic field, and will cause a charge separation along the magnetic field direction. While the search for the CME is still ongoing in experiments at Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC), the CME-sensitive observables need to be scrutinized to exclude the non-CME contributions. In this work, we examine the influence of globally spin-aligned ρ mesons on the γ 112 correlator, the R Ψ 2 (ΔS) correlator, and the signed balance functions, via a toy model and a multiphase transport model (AMPT). The global spin alignment of vector mesons could originate from non-CME mechanisms in heavy-ion collisions, and is characterized by the 00-component of the spin density matrix, ρ 00 . Here, we find that the CME observables show similar dependence on ρ 00 , and could receive a positive (negative) contribution from ρ-decay pions, if the ρ 00 of ρ mesons is larger (smaller) than 1/3. Since pions are the most abundant particles in such collisions, the ρ 00 measurements for ρ mesons become crucial to the interpretation of the CME data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗
DOE OSTI2021

SiOC coatings on yttria stabilized zirconia microspheres using a fluidized bed coating process

In this study, defect-free SiOC coatings were prepared on yttria stabilized zirconia (YSZ) microspheres by a fluidized bed coating process. Effects of the rheological properties of the coating solution on the coating process were elucidated. An impact regime diagram was constructed, which demonstrated that the coating mechanisms were collision/impact. During the fluidized bed coating, longer spouted time resulted in wider dispersion, longer residence time, and more circulatory motion of particles; the fluid distributed more uniformly throughout the column, as demonstrated in our Multiphase Flow with Interface eXchange (MFiX) simulations. Two-step pyrolysis in Ar achieved complete coating layers, which were comprised of SiOC, SiO 2 , SiC, and graphite. The two stage mass loss during the pyrolysis corresponded to simultaneous reactions due to depolymerization and hydrocarbon loss from 400 to 600°C. Carbon cluster size in the pyrolyzed samples was calculated to be 25 ± 2 Å. In conclusion, this work provides a new method for producing SiOC coatings on micron spheres, with nuclear TRISO fuel particles as the application.

42 ENGINEERING↗
DOE OSTI2025

Shining light on nanoscale ‘vine-on-stick’ eutectic structures in the Al-Ce-Ni system

Multiphase eutectics often comprise entangled solid phases with nanoscale periodicity, making it difficult to unravel their 3D connectivity using conventional 2D techniques. Among such systems, the three-phase eutectic Al-Al 11 Ce 3 -Al 3 Ni stands out for its ultrafine (∼100 nm) interphase spacing and promising creep resistance, yet its microstructure remains relatively unexplored despite its potential for high-temperature applications. Here, we use scanning hard X-ray microscopy with an unprecedented ∼10 nm pixel size to resolve its 3D morphology. Reconstructions reveal a novel “vine-on-stick” motif, wherein Al 11 Ce 3 wraps around Al 3 Ni pillars. Analysis of phase tortuosities confirms that Al 11 Ce 3 exhibits more convoluted morphologies than Al 3 Ni. No orientation relationship was observed between the two intermetallics, suggesting that growth is controlled by local solute gradients rather than epitaxy. Fragmentation of intermetallics along their longitudinal axis suggests a Rayleigh-type breakup mechanism in solid state. The “vine-on-stick” pattern may generalize to other alloy systems with variable interfacial anisotropies and low mutual solubilities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A distributed voltage inference framework for cyber-physical attacks detection and localization in active distribution grids

The transition to active distribution grids with real-time monitoring and control depends on the proliferation of advanced communication networks and devices. This paradigm shift towards a cyber-physical architecture also introduces new vulnerabilities for adversaries to exploit and launch sophisticated cyber-physical attacks targeting grid observability. Current research highlights the challenges in distinguishing attacks on voltage phasor or nodal injection measurements and isolating multi-source attack locations in a multiphase distribution grid. The attack detection and localization methods in literature face accuracy issues, applications across diverse attack scenarios, or scalability limits. Here, to bridge these gaps, this paper proposes a distributed Voltage Inference framework for real-time detection and localization of cyber-physical attacks, addressing scalability, adaptability, and accuracy challenges in state-of-the-art methods. The proposed methodology leverages the distributed nature of the Voltage Inference framework through a two-step process of prediction and correction, together with a tractable graph partitioning approach, providing a reliable solution to identify compromised measurement sources and facilitate isolation. Extensive testing on IEEE 13 and 123-node distribution feeders underscores the algorithm’s efficacy, enhancing the security and resilience of active distribution grids against evolving cyber threats. Additionally, Hardware-in-the-Loop (HIL) implementation validates the proposed strategy’s practical applicability in real-world scenarios.

active distribution grids↗
DOE OSTI2022

Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales

Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.

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

Coupled hydro-thermal flow and radionuclide transport driven by spatial variation of heat-generating radioactive wastes in shale formations

Deep geologic disposal of multiple nuclear waste packages with various heat sources can induce nonuniform hydro-thermal behaviors in the near-field of the repository, consequently influencing the long-term radionuclide transport in the far-field once waste form breach initiates. Here, this study looks into three cases with variation in the spatial order of six groups of heat sources (10th, 50th, 75th, 90th, 95th, and 99th percentiles of heat outputs generated from 1,981 as-loaded dual-purpose canisters in the field site) in a shale-hosted repository with respect to the uni-directional groundwater flow (from west to east): (1) cooler waste packages from west to east, (2) hotter waste packages from west to east, and (3) hottest waste packages in the middle of the repository. Our field-scale PFLOTRAN simulation represents heat-driven multiphysics coupled mechanisms, including multiphase flow, heat transfer, and chemical/radioactive transport, and also, calculates the onset of waste form breach based on temperature-dependent canister vitality. The results from this sensitivity study will quantify the short- (less than 1 × 10 3 years) and long-term (up to 1 × 10 6 years) impacts of sporadic heat pulses from waste package on the spatio-temporal perturbation in hydro-thermal flow quantities and the rate of radionuclide transport in both near- and far-field of the repository system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Phase Stability and Kinetics of Topotactic Dual Ca 2+ –Na + Ion Electrochemistry in NaSICON NaV 2 (PO 4 ) 3

Recent reports of reversible calcium plating and stripping have rekindled interest in the development of Ca-ion batteries (CIBs) as next-generation energy storage devices. This technology has the potential to overcome the limitations of conventional Li-ion batteries, but CIBs are plagued by a paucity of suitable cathode materials. To date, NaSICON-structured NaV 2 (PO 4 ) 3 has been demonstrated as a successful cathode candidate, exhibiting reversible (de)intercalation of 0.6 mol Ca 2+ along with stable cycling performance. However, a complex multiphase mixture forms on discharge so the Ca-ion charge storage mechanism in the NaSICON framework is poorly understood. Here in this work, we report on an investigation of the structure and/or Na + /Ca 2+ environment(s) of a variety of chemically prepared NaSICON Ca x Na y V 2 (PO 4 ) 3 phases which were characterized using synchrotron XRD, SEM-EDS, 23 Na NMR, and TEM. Highly calciated CaV 2 (PO 4 ) 3 , Ca 1.5 V 2 (PO 4 ) 3 , and CaNaV 2 (PO 4 ) 3 phases can be prepared at high temperature, but -unlike Ca 0.6 NaV 2 (PO 4 ) 3 -these materials are electrochemically inactive. To better understand the fundamental factors impacting successful Ca 2+ electrochemistry in this system, DFT was employed to examine the Ca x Na y V 2 (PO 4 ) 3 phase diagram and Ca 2+ diffusion mechanism. Theoretical insights show that phase separation into Na-rich and Ca-rich phases is a reason for the capacity limitation and demonstrate that Na + ions in the host materials assist the migration of neighboring Ca 2+ ions, enabling reversible electrochemistry in Ca x Na y V 2 (PO 4 ) 3 . This investigation of fundamental principles affecting reversible Ca 2+ (de)intercalation in Ca x Na y V 2 (PO 4 ) 3 allows for the development of design principles to enable the discovery of a variety of successful cathodes for CIBs.

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

Resonant Soft X-ray Scattering Reveals Hierarchical Structure in a Multicomponent Vapor-Deposited Glass

Multiphase vapor-deposited glasses are an important class of materials for organic electronics, particularly organic photovoltaics and thermoelectrics. These blends are frequently regarded as molecular alloys and there have been few studies of their structure at nanometer scales. Here, in this work, we show that a codeposited system of TPD (N,N'-bis(3-methylphenyl)-N,N'-diphenylbenzidine) and DO37 (disperse orange 37), two small molecule glass-formers, separates into amorphous, compositionally distinct phases with a domain size and spacing ca. 10s of nanometers that depends on substrate temperature during deposition. Domains rich in one of the two components become larger and more pure at higher deposition temperatures. We use resonant soft X-ray scattering (RSoXS) complemented with atomic force microscopy (AFM) and photoinduced force microscopy to measure the phase separation, topography, and purity of the deposited films. A forward-simulation approach to RSoXS analysis, the National Institute of Standards and Technology RSoXS Simulation Suite (NIST RSoXS simulation suite), is used with models developed from AFM images to evaluate the energy dependence of scattering across multiple length scales and interpret the RSoXS with respect to structure within the films. We find that the RSoXS is sensitive to a length scale of phase separation buried within the film that is consistent with the surface composition profile, and correlates to the topography to an extent that depends on substrate temperature. We demonstrate that vacuum scattering, which is often ignored in RSoXS analysis, contributes significantly to the features and energy dependence of the RSoXS pattern, and then illustrate how to properly account for vacuum scattering to analyze films with significant roughness. We then use this analysis framework to understand structure development mechanisms that occur during vapor deposition of a TPD-DO37 codeposited glass with results that outline paths to tune morphology in multicomponent materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dara: Automated Multiple-Hypothesis Phase Identification and Refinement from Powder X-ray Diffraction

Powder X-ray diffraction (XRD) is a foundational technique for characterizing crystalline materials. However, the reliable interpretation of XRD patterns, particularly in multiphase systems, remains a manual and expertise-demanding task. As a characterization method that only provides structural information, multiple reference phases can often be fit to a single pattern, leading to potential misinterpretation when alternative solutions are overlooked. To ease humans’ efforts and address the challenge, we introduce Dara (data-driven automated Rietveld analysis), a framework designed to automate the robust identification and refinement of multiple phases from powder XRD data. Dara performs an exhaustive tree search over all plausible phase combinations within a given chemical space and validates each hypothesis using the BGMN Rietveld refinement routine. Key features include structural database filtering, automatic clustering of isostructural phases during tree expansion, and peak-matching-based scoring to identify promising phases for refinement. When ambiguity exists, Dara generates multiple hypothesis which can then be decided between by human experts or with further characterization tools. By enhancing the reliability and accuracy of phase identification, Dara enables scalable analysis of realistic complex XRD patterns and provides a foundation for integration into multimodal characterization workflows, moving toward fully self-driving materials discovery.

Biological databases↗
DOE OSTI2020

Mediated Fuel Cells: Soluble Redox Mediators and their Applications to Electrochemical Reduction of O2 and Oxidation of H2, Alcohols, Biomass, and Complex Fuels

Mediated fuel cells are electrochemical devices that produce power in a manner similar to that of conventional proton exchange membrane fuel cells (PEMFCs). They differ from PEMFCs in their use of redox mediators dissolved in liquid electrolyte to conduct oxidation of the fuel or reduction of the oxidant, typically O2, in bulk solution. The mediators transport electrons (and often protons) between the electrode and the catalysts or chemical reagents in solution. This strategy can help overcome many of the challenges associated with conventional fuel cells, including managing complex multi-phase reactions (as in O2 reduction) or the use of challenging or heterogeneous fuels, such as hydrocarbons, polyols and biomass. Mediators are also commonly used in enzymatic fuel cells, where direct electron transfer from the electrode to the enzymatic active site can be slow. This review provides a comprehensive survey of historical and recent mediated fuel cell efforts, including applications using chemical and enzymatic catalysts. Our research on in this field has been exploring a number of different topics and has been supported different funding sources: the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences (molecular mediators for catalytic oxygen reduction); the Great Lakes Bioenergy Research Center, DOE Office of Science DE-SC0018409 (research on biomass-based fuels; U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Catalysis Program, DE-FG02-05ER15690 (copper-catalyzed oxidation reactions); and the NIH NIGMS, R01 GM100143 and R35 GM134929 (use of molecular mediators for organic chemical synthesis).

Anson, Colin W.↗
DOE OSTI2021

Bridging Scales in Bioenergy and Catalysis: A Review of Mesoscale Modeling Applications, Methods, and Future Directions

Between the molecular and reactor scales, which are familiar to the chemical engineering community, lies an intermediate regime, here termed the “mesoscale,” where transport phenomena and reaction kinetics compete on similar time scales. Bioenergy and catalytic processes offer particularly important examples of mesoscale phenomena owing to their multiphase nature and the complex, highly variable porosity characteristic of biomass and many structured catalysts. In this review, we overview applications and methods central to mesoscale modeling as they apply to reaction engineering of biomass conversion and catalytic processing. A brief historical perspective is offered to put recent advances in context. Applications of mesoscale modeling are described, and several specific examples from biomass pyrolysis and catalytic upgrading of bioderived intermediates are highlighted. Methods including reduced order modeling, finite element and finite volume approaches, geometry construction and import, and visualization of simulation results are described; in each category, recent advances, current limitations, and areas for future development are presented. Owing to improved access to high-performance computational resources, advances in algorithm development, and sustained interest in reaction engineering to sustainably meet societal needs, we conclude that a significant upsurge in mesoscale modeling capabilities is on the horizon that will accelerate design, deployment, and optimization of new bioenergy and catalytic technologies.

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

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗
DOE OSTI2024

3D Quantification of Pore Networks and Anthropogenic Carbon Mineralization in Stacked Basalt Reservoirs

Basalt formations are promising candidates for the geologic storage of anthropogenic CO 2 due to their storage capacity, porosity, permeability, and reactive geochemical trapping ability. The Wallula Basalt Carbon Storage Pilot Project demonstrated that supercritical CO 2 injected into >800 m deep Columbia River Basalt Group stacked reservoir flow tops mineralizes to ankerite-siderite-aragonite on month-year timescales, with 60% of the 977 metric tonnes of CO 2 converted within two years. The potential impacts of mineral precipitation and consequent changes on rock porosity, pore structure, pore size, and pore size distributions have likely been underestimated hitherto. Herein, we address these knowledge gaps using X-ray Micro Computed Tomography (XMT) to evaluate the pore network architecture of sidewall cores recovered two years after CO 2 injection. In the present study, we performed a detailed quantitative analysis of the CO 2 -reacted basalt cores by XMT imaging. Reconstructed 3D images were analyzed to study the distribution and volumetric details of porosity and carbonate nodules in the cores along with the various other phases, providing insight into paragenesis and carbonate growth mechanisms, including mineralogic/chemical zonation. Finally, these findings are being used to parametrize multiphase reactive transport models to predict the fate and transport of subsurface CO 2 , enabling scale-up to commercial-scale geologic carbon storage in basalts and other reactive mafic-ultramafic formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeling Novel Aqueous Particle and Cloud Chemistry Processes of Biomass Burning Phenols and Their Potential to Form Secondary Organic Aerosols

Phenols emitted from biomass burning contribute significantly to secondary organic aerosol (SOA) formation through the partitioning of semivolatile products formed from gas-phase chemistry and multiphase chemistry in aerosol liquid water and clouds. The aqueous-phase SOA (aqSOA) formed via hydroxyl radical (˙OH), singlet molecular oxygen ( 1 O 2 *), and triplet excited states of organic compounds ( 3 C*), which oxidize dissolved phenols in the aqueous phase, might play a significant role in the evolution of organic aerosol (OA). However, a quantitative and predictive understanding of aqSOA has been challenging. Here, we develop a stand-alone box model to investigate the formation of SOA from gas-phase ˙OH chemistry and aqSOA formed by the dissolution of phenols followed by their aqueous-phase reactions with ˙OH, 1 O 2 *, and 3 C* in cloud droplets and aerosol liquid water. Here, we investigate four phenolic compounds, i.e., phenol, guaiacol, syringol, and guaiacyl acetone (GA), which represent some of the key potential sources of aqSOA from biomass burning in clouds. For the same initial precursor organic gas that dissolves in aerosol/cloud liquid water and subsequently reacts with aqueous phase oxidants, we predict that the aqSOA formation potential (defined as aqSOA formed per unit dissolved organic gas concentration) of these phenols is higher than that of isoprene-epoxydiol (IEPOX), a well-known aqSOA precursor. Cloud droplets can dissolve a broader range of soluble phenols compared to aqueous aerosols, since the liquid water contents of aerosols are orders of magnitude smaller than cloud droplets. Further, our simulations suggest that highly soluble and reactive multifunctional phenols like GA would predominantly undergo cloud chemistry within cloud layers, while gas-phase chemistry is likely to be more important for less soluble phenols. But in the absence of clouds, the condensation of low-volatility products from gas-phase oxidation followed by their reversible partitioning to organic aerosols dominates SOA formation, while the SOA formed through aqueous aerosol chemistry increases with relative humidity (RH), approaching 40% of the sum of gas and aqueous aerosol chemistry at 95% RH for GA. Our model developments of biomass-burning phenols and their aqueous chemistry can be readily implemented in regional and global atmospheric chemistry models to investigate the aqueous aerosol and cloud chemistry of biomass-burning organic gases in the atmosphere.

09 BIOMASS FUELS↗