SEARCH · Engineering Papers

Results for “Irreversible processes”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Filtration Performance of Simulated 200 West Area Waste Feeds

This report describes the scaled experimental system and approach used to examine dead-end filtration performance of representative 200W waste feeds. The scaled system, which was originally designed and assembled to test Tank Side Cesium Removal (TSCR) system performance with higher-than-expected solid loadings in 2021 (Schonewill et al. 2021), was repurposed to conduct the current experiments at ~1/145 of full scale (based on throughput). Six experimental runs were conducted with five different 200W waste feed simulants: three using a DEF module scaled for TSCR and three using a DEF module scaled for the 200W process modules (based on the current design for the Advanced Modular Pretreatment System). Each experiment was run continuously for multiple days with an operating approach prototypic of the full-scale system. Staff performing the experimental runs monitored performance, obtained data from calibrated process instruments, and collected samples for observation and analysis. The measured data are presented with a focus on assessing DEF performance – specifically, the filters’ differential pressure response to the five waste simulants, frequency and efficacy of backwashing, and baseline recovery between experimental runs; data related to ion exchange column performance are also discussed in cases where the opportunity arose. The experimental campaign demonstrated that the DEFs satisfied their primary function of protecting the ion exchange column from solid intrusion for all the representative simulants used. The filters readily handled solids loadings of =500 ppm (and even greater), especially the modules scaled to the 200W process modules. Adjustments to the processing flow rate and reductions in feed temperature were observed to affect the rate of differential pressure increase on the filters, but neither adversely affected the ability of the DEFs to perform their primary function. Backflushing reliably recovered filter performance in all runs, although it did not prevent irreversible fouling for one simulant. The run that exhibited irreversible fouling established that both the quantity and the nature of the solids being filtered need to be considered when projecting filter performance.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Coupled Thermo-Hydrological-Mechanical-Chemical Behavior of Anisotropic Granite for Geologic Disposal of High-Level Radioactive Waste: A Core-Scale Laboratory Investigation

The coupled thermo-hydrological-mechanical-chemical (THMC) behavior of rock within an Excavation Damaged Zone (EDZ) is critical for the safety and long-term performance of a geological repository for high-level radioactive wastes. While many laboratory experiments have been conducted to investigate EDZ rocks, the flow and deformation characteristics resulting from anisotropic rock textures and microcrack distribution under triaxial loading and elevated temperatures remain poorly understood. Particularly, cracks at various scales serve as fast paths for fluid flow and solute transport and present as focal points of mechanical weakness, which complicate the coupled THMC processes in anisotropic EDZ rocks and challenge modeling predictions. Here, in this study, a series of core-scale experiments was conducted on three granite samples under repository-relevant conditions. These rock samples were obtained from the Grimsel Underground Research Laboratory (URL), featured by anisotropic minerals (represented by bedding layers) and microcrack distributions and coarse cm-scale grain sizes. During the experiments, samples were subjected to an elevated temperature at 90 °C and different triaxial loading conditions either by radial (normal to bedding layers) or axial (parallel to bedding layers) compaction. Water was injected into the samples, and the rock permeability evolutions and effluent water chemistry were monitored closely. For intact samples, thermal expansion of minerals at 90 °C resulted in a large, 75% irreversible permeability reduction and rock strengthening under radial compaction, while thermal impact was limited to a 15% permeability reduction under axial compaction. In contrast, for a sample containing open cracks, the growth of fractures during the experiment resulted in an abrupt permeability increase and fast failure at 90 °C. The effluent water chemistry indicates much more considerable mineral dissolution from large shear sliding than that in rocks dominated by mechanical compaction. These results helped better understand the coupled THMC processes in anisotropic rocks containing cracks, evaluate the behaviors of EDZ rocks, and predict the long-term evolution of EDZ for the performance of the repositories.

Coupled THMC processes↗

DOE OSTI2025

Short-Term Inhibition and Long-Term Enhancement of Irreversible Trace Metal Binding to Goethite in Multi-Metal Systems

Iron (oxyhydr)oxide minerals are important sorbents of trace metals in aquatic environments. Adsorption–desorption hysteresis has been documented for trace metals after aging, suggesting that they become incorporated over time. We previously found that ion size controls the extent of irreversible metal binding of single metals to goethite. In this study, we evaluate how the presence of multiple metals influences irreversible binding. Mixed Ni–Zn and Ni–Cd solutions were aged with goethite at pH 7 for 2 days, 30 days, and 60 days, after which isotope exchange experiments were performed to assess the lability of the adsorbed metal pools. After 2 days of aging, Ni and Cd in mixed-metal systems were more labile than in single-metal systems, indicating that competitive adsorption may partially block irreversible binding on short time scales. After 60 days, all three metals had larger irreversibly bound fractions than in single-metal systems. X-ray absorption near-edge structure (XANES) spectroscopy indicates that irreversibly bound fractions correlate with incorporation into goethite. This study demonstrates that pools of coadsorbed metals are bioaccessible on a time scale of minutes. Cooperative effects promote the structural sequestration of even large trace metals, causing reactive transport models to likely overpredict metal mobility. Long-term metal entrapment processes hinder micronutrient availability, impact contaminant sequestration and critical mineral recovery, and yield complex pathways through which metal solubilization may be enhanced or hindered during redox cycling.

58 GEOSCIENCES↗

DOE OSTI2021

Non-equilibrium insertion of lithium ions into graphite

Graphite has been regarded as the most important anode material for currently used lithium-ion batteries due to its two-dimensional (2D) nature hosting ionic intercalations. However, the kinetic insertion of Li ions is still not well known microscopically. In this work, we investigate the real-time intercalation process of Li ions using in situ transmission electron microscopy. We observe the lithium insertion process at the atomic scale, in which the graphite layers undergo expansion, forming wrinkles and finally inhomogeneous cracks as the Li ions accumulate, different from the proposed models. Leveraging on theoretical simulations, Li-ion migration driven by an external electrical field is suggested to be induced into the irreversible wrinkled structures. This non-equilibrium behavior that occur in lithium-ion batteries can be more pronounced at a high charging rate, which will practically degrade the capacity of graphite. Furthermore, this work unveils the reaction scenario of the non-equilibrium Li-ion insertion, which benefits the understanding of the performance of graphite-based energy-storage devices.

25 ENERGY STORAGE↗

DOE OSTI2021

Visualizing plant cell wall changes proves the superiority of hydrochloric acid over sulfuric acid catalyzed γ-valerolactone pretreatment

Thermochemical pretreatment is one of the key steps to process lignocellulosic biomass for production of biofuels and biomaterials. γ-valerolactone (GVL) has been used as a green and renewable solvent to efficiently dissolve lignin under acidic conditions and enhance subsequent enzyme digestibility. Further improvement is still required to lower the pretreatment temperature in order to reduce sugar degradation and irreversible lignin condensation, as well as the capital cost. In this study, we compared the use of HCl and H 2 SO 4 as catalysts during GVL pretreatment, and found the performance of GVL-HCl at 100 °C was comparable to that of GVL-H 2 SO 4 at 120 °C in terms of xylan and lignin removal and enzyme digestibility. We further monitored the lignin removal and cellulose accessibility in the cell walls at different pretreatment time points by imaging the changes in lignin auto-fluorescence and in CtCBM3-GFP binding with confocal laser scanning microscopy (CLSM), respectively. We found GVL-HCl pretreatment at a relative low temperature (100 °C) could rapidly remove lignin in the compound middle lamella and cell corner areas, which concurred with the increase of CtCBM3-GFP binding in these areas. Real-time imaging of cell wall degradation by cellulases further revealed that the secondary cell walls could be digested from both cell lumen and CML sides, and eventually fully deconstructed within 24 h. Furthermore, our results provide new insight into the effects of chloride anions on plant cell wall structure during GVL pretreatment, and offer potential routes to further optimize and enhance the efficiency of GVL-based pretreatment.

09 BIOMASS FUELS↗

DOE OSTI2024

On Gibbs Equilibrium and Hillert Nonequilibrium Thermodynamics

During his time at Royal Institute of Technology (Kungliga Tekniska högskolan) in Sweden, the present author learned nonequilibrium thermodynamics from Mats Hillert. The key concepts are the separation of internal and external variables of a system and the definitions of potentials and molar quantities. In equilibrium thermodynamics derived by Gibbs, the internal variables are not independent and can be fully evaluated from given external variables. While irreversible thermodynamics led by Onsager focuses on internal variables though often mixed with external variables. Hillert integrated them together by first emphasizing their differences and then examining their connections. His philosophy was reflected by the title of his book “Phase Equilibria, Phase Diagrams and Phase Transformations” that puts equilibrium, nonequilibrium, and internal processes on equal footing. Here, in the present paper honoring Hillert, the present author reflects his experiences with Hillert and his work in last 40 years and expresses his gratitude for all the wisdom and support from him in terms of “Hillert nonequilibrium thermodynamics” and discusses some recent topics that the present author has been working on.

36 MATERIALS SCIENCE↗

DOE OSTI2022

The Role of Oxygen in Lithiation and Solid Electrolyte Interphase Formation Processes in Silicon-Based Anodes

Silicon oxides (SiO x ) have been considered as promising alternatives to pure Si in high energy anodes in lithium-ion batteries (LIBs) due to their improved cycling stability. However, their fundamental lithiation mechanism has not yet been systematically investigated, and potential collateral downsides remain unclear. In this work, we report on the role of oxygen in lithiation/delithiation and solid electrolyte interphase (SEI) formation processes in SiO x thin film model electrodes with different oxygen contents. Here, we show that the SiO x anodes with higher oxygen content experience smaller volume change and form a thinner and more stable SEI, both of which are beneficial for cycling stability. However, these SiO x anodes also show an irreversible lithiation at around 0.7 V attributed to the reduction of Si oxides, leading to lower first cycle coulombic efficiency that is undesirable for practical applications. Overall, these results offer a balanced perspective on the advantages and disadvantages that oxygen brings to Si-based anodes in LIBs.

25 ENERGY STORAGE↗

DOE OSTI2026

Cooperative Agreement To Analyze variabiLity, change and predictabilitY in the earth SysTem (CATALYST)

CATALYST proposes to perform foundational coordinated research in a team-oriented collaborative effort aimed at advancing a robust understanding of modes of Earth system variability and change using models, observations and process studies. The proposed research will address the DOE/BER mission by exploring the limits to predictability, identifying fundamental underlying mechanisms, quantifying interactions among modes of variability, and discovering tipping points in the Earth system to understand the current and future impacts of these phenomena on regional and global climate. Four fundamental gaps are identified in our knowledge of the Earth system: 1) What are the limits to predictability on various timescales? 2) What are the interactions among modes of Earth system variability? 3) How may modes of Earth system variability change in response to changes in external forcing, and what are the tipping points involved with those changes? 4) How are high impact events connected to modes of Earth system variability and how may they change in the future? Related to those gaps in our knowledge, we formulate four research objectives to address those gaps using a combination of Earth system models (ESMs) and machine learning (ML) methods. Research Objective 1 (RO1) addresses the first gap above and proposes to understand modes of variability and their limits of predictability on subseasonal to decadal timescales using ESMs and ML. Research Objective 2 (RO2) addresses the second gap and proposes to use a hierarchy of models to understand relevant processes and feedbacks related to how modes of variability interact with each other. Research Objective 3 (RO3) is designed to study the third gap and proposes to examine the role of external forcings in changes of modes of Earth system variability and their interactions, and the likelihood and predictability of tipping points and irreversible changes. Research Objective 4 (RO4) will address the fourth gap and proposes to use high resolution ESMs, regionally refined models (RRMs), and ML methods to investigate the relationships between high impact events (e.g. flash droughts and precipitation extremes, atmospheric rivers (ARs), tropical cyclones (TCs), storm surge/sea level rise), the synoptic systems that produce them, and their changes related to modes of Earth system variability. The research will involve the use of the Community Earth System Model (CESM), Energy Exascale Earth System Model (E3SM), CMIP multi-model data sets, a hierarchy of simpler models, and numerous observational data sets. In the course of the proposed research, CATALYST will contribute to metrics and diagnostics that will be integrated in Coordinated Model Evaluation Capabilities (CMEC), particularly with regards to the Quasi-biennial Oscillation (QBO) and its interactions with the Madden-Julian Oscillation (MJO), high atmospheric pressure blocking, and new precipitation metrics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Process for generating high purity synthesis gas hydrogen from heavy oil or hydrocarbons

The present invention provides a steam reforming process for heavy oil or hydrocarbons using a circulating fluidized bed reactor, the process having a reforming step and a regeneration step, wherein the reforming step and the regeneration step comprise a fluidized reactor containing a fluidizable nickel-containing reforming catalyst and produce hydrogen as a product of the reforming bed. The invention produces high purity hydrogen in the synthesis gas product stream and avoids irreversible fouling on the catalyst.

Srinivas, Girish↗

DOE OSTI2023

Novel method to control explosive shock sensitivity: Formulations, experiments and modeling

The microstructure of a heterogeneous high explosive (HE) affects its shock initiation sensitivity and detonation performance. Alteration of the void content and/or void structure (i.e., bulk heating or mechanical damage) therefore changes the shock initiation behavior. Controlling and predicting the change in shock sensitivity after an HE has undergone microstructural changes addresses an important and challenging goal for the design and understanding of novel energetic material formulations. In this study, we aimed to develop HE systems to precisely tune shock sensitivity by a thermal treatment prior to use. Specifically, we incorporated a small fraction (1 wt% or less) of thermally expandable microspheres (TEMs) during the formulation process of various plastic-bonded explosives (PBX). TEMs typically consist of a thermoplastic acrylonitrile shell (10-50 µm diameter), which encapsulates an inert low boiling hydrocarbon. Upon heating, the TEMs expand as the shell softens while the hydrocarbon gasifies, increasing the internal pressure and expanding the particle by as much as 120 vol% (irreversibly). Here, we present our progress on shock sensitivity comparisons of HE formulations doped with TEMs after heating and expansion. Experiments, using Composition C-4 thermally cycled to 120°C, showed an increase in shock sensitivity. Additionally, mesoscale modeling revealed that the TEM itself does not act as a hotspot, but instead has a secondary effect on the response of nearby voids. We conclude that TEMs indeed do provide a method to tune shock sensitivity. Further, our results emphasize the non-uniformity of shock waves at the mesoscale, and that upstream defects influence the reactivity of downstream defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Atomic resolution tracking of nerve-agent simulant decomposition and host metal–organic framework response in real space

Gas capture and sequestration are valuable properties of metal–organic frameworks (MOFs) driving tremendous interest in their use as filtration materials for chemical warfare agents. Recently, the Zr-based MOF UiO-67 was shown to effectively adsorb and decompose the nerve-agent simulant, dimethyl methylphosphonate (DMMP). Understanding mechanisms of MOF-agent interaction is challenging due to the need to distinguish between the roles of the MOF framework and its particular sites for the activation and sequestration process. Here, we demonstrate the quantitative tracking of both framework and binding component structures using in situ X-ray total scattering measurements of UiO-67 under DMMP exposure, pair distribution function analysis, and theoretical calculations. The sorption and desorption of DMMP within the pores, association with linker-deficient Zr6 cores, and decomposition to irreversibly bound methyl methylphosphonate were directly observed and analyzed with atomic resolution.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Cross-flow separation characteristics and piloting of graphene oxide nanofiltration membrane sheets and tubes for kraft black liquor concentration

Dewatering of weak black liquor (WBL) in the kraft cycle by evaporation is highly energy intensive. Membranes are an attractive alternative for energy-efficient dewatering, but existing commercial polymeric or ceramic membranes are either degraded in BL or have high capital costs. Our recent works have demonstrated the engineering of graphene oxide (GO) nanofiltration membranes, their stability and promising performance in BL conditions, and preliminary scale-up into sheets and tubes. Here, we describe in detail the separation characteristics of GO membrane sheets and tubes under real BL conditions and crossflow operation. Recycle-mode piloting of a GO tubular membrane showed average “production flux” of 16 L/m2/h (LMH) and high rejections of lignin (98.3%), total solids (66%), and total organic carbon (83%), with no signs of irreversible fouling identified. A corresponding GO sheet membrane produced an average flux of ~25 LMH and maintained high lignin rejection of ~97% during a slipstream pilot at a kraft mill site using WBL with ~16 wt% total solids (TS). Finally, we piloted a Dow/DuPont XUS1808 polyamide composite reverse osmosis (RO) membrane for last-mile processing of the GO nanofiltration membrane permeate. The RO membrane showed a steady state flux of 19 LMH at 65 bar and produced ~0.02 wt% TS water product, which is highly suitable for reuse in pulp washing operations in the kraft process. The results have strong positive implications for the industrial application of GO membranes in BL concentration and other related applications.

Materials Science↗

DOE OSTI2022

Designing Sites in Heterogeneous Catalysis: Are We Reaching Selectivities Competitive With Those of Homogeneous Catalysts?

A critical review of different prominent nanotechnologies adapted to catalysis is provided, with focus on how they contribute to the improvement of selectivity in heterogeneous catalysis. Ways to modify catalytic sites range from the use of the reversible or irreversible adsorption of molecular modifiers to the immobilization or tethering of homogeneous catalysts and the development of well-defined catalytic sites on solid surfaces. The latter covers methods for the dispersion of single-atom sites within solid supports as well as the use of complex nanostructures, and it includes the post-modification of materials via processes such as silylation and atomic layer deposition. All these methodologies exhibit both advantages and limitations, but all offer new avenues for the design of catalysts for specific applications. Because of the high cost of most nanotechnologies and the fact that the resulting materials may exhibit limited thermal or chemical stability, they may be best aimed at improving the selective synthesis of high value-added chemicals, to be incorporated in organic synthesis schemes, but other applications are being explored as well to address problems in energy production, for instance, and to design greener chemical processes. Finally, the details of each of these approaches are discussed, and representative examples are provided. We conclude with some general remarks on the future of this field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Mixing Matters: Nanoscale Heterogeneity and Stability in Metal Halide Perovskite Solar Cells

The structural stability of the metal halide perovskite (MHP) absorber material is crucial for the long-term solar cell stability in this thin-film photovoltaic technology. Here, we use mixed A-site FA0.83Cs0.17PbI3 to demonstrate that nanoscale compositional heterogeneity can serve as initiation sites for more macroscale, irreversible phase segregation, which causes device performance degradation. Probing compositional heterogeneity on length scales that has not been detected with conventional characterization techniques, we analyze the tetragonal to cubic phase transition behavior to indirectly determine the level of nanoscale compositional heterogeneity in the initial films. Further, we show that the thermal annealing conditions of the MHP layer during film processing influence the initial nanoscale compositional heterogeneity, and changing these processing conditions can be used to improve the device performance stability. The insights into structural degradation mechanisms initiated by nanoscale compositional heterogeneity and the proposed mitigation strategies will help guide the way toward long-term stable MHP solar cells.

14 SOLAR ENERGY↗

DOE OSTI2025

Observations of Subduction, Downward Heat Flux and Dense Filament Collapse in the Northern Gulf of Mexico

Submesoscale processes are important contributors to the global heat budget and generally support upward heat transport through restratification. However, in salinity‐stratified regions, such as the northern Gulf of Mexico with its influx of freshwater from the Mississippi‐Atchafalaya river system, temperature can act like a passive tracer and submesoscale processes can contribute to downward heat transport. Oceanic heat content is a factor in many environmental risks the region faces, for example, hurricane intensification, and marine heatwaves. During the 2022 field campaign of the Submesoscales Under Near‐Resonant Inertial Shear Experiment, a sampling plan was developed to study such submesoscale processes in high resolution. Over 31 hr, four assets (two research ships and two remotely controlled boats) drove in parallel across a dense filament, capturing its evolution in time and space. The observations show that surface waters, warmed by daytime solar radiation, were subducted and that the associated overturning circulation transported heat below the surface layer where it was later irreversibly mixed away. The estimated downward heat flux was as strong as the concurrent net air‐sea heat flux into the ocean. The filament was then observed to rapidly collapse which we attribute to boundary layer turbulence and the breakdown of geostrophic balance. The collapsing fronts display behaviors indicative of gravity currents. These observations highlight how in salinity‐stratified regions, frontal dynamics can be associated with downward heat flux and how the submesoscale can play an important role in the oceanic heat budget.

58 GEOSCIENCES↗

DOE OSTI2025

Coarse-grained molecular dynamics simulation of solvent-dependent cellulose nanofiber interactions

Associations between cellulose are important both in biofuel production and in the use of cellulose for biomaterials. Cellulose nanofibers (CNFs) are sustainable, strong, light-weight alternatives to traditional materials in manufacturing, but are challenging to obtain due to irreversible aggregation in solution during preparative fibrillation. Therefore, it is imperative to understand the underlying factors driving aggregation with a view to designing solvents that can effectively compete with interfiber interactions, hence reducing aggregation. Molecular dynamics (MD) simulation at atomic detail can provide useful information on local interactions. However, the length and timescales accessible are too short to fully capture association processes. Here, we provide a method for accessing the longer length and timescales required using coarse-grained (CG) MD simulations with a MARTINI force field to calculate the interaction behavior of CNFs in three selected solvents: NaOH-urea-water, acetone, and neat water. The CG results are consistent with our prior all-atom MD and with previous experimental results. While acetone is found not to be an effective solvent, urea and ionic moieties in NaOH-urea-water not only solvate the fibrils but also improve the confinement of water molecules around them as shown by the solvent residence times and mean-square displacements. Overall, the presence of urea and ions reduces the likelihood of aggregation in multi-CNF systems relative to neat water irrespective of whether the hydrophobic or hydrophilic CNF surfaces are interacting. In conclusion, the CG method shows clear promise for selecting potential high-performance solvents for experimental prioritization in bioenergy and biomaterials research in a relatively fast manner as well as for understanding the aggregation and rheological behavior of CNF-solvent systems.

aggregation↗

DOE OSTI2024

Ameliorating the sodium storage performance of hard carbon anode through rational modulation of binder

Hard carbon anodes have emerged as promising candidates for sodium-ion batteries due to their inherent advantages. Nevertheless, the surface imperfections in these materials often culminate in irreversible electrolyte consumption, fostering the development of a heterogeneous and fragile solid electrolyte interface (SEI), thereby compromising the initial Coulombic efficiency (ICE). Here, drawing inspiration from the catalytic potential of C=O (carbonyl) bonds in directing preferential salt reduction, we introduce a novel strategy that leverages the modulation of the binder, a long-term overlooked pivotal components in the electrode process. Specifically, Polymethyl methacrylate (PMMA), abundant in C=O groups, is partially substituted for PVDF, ensuring robust adhesion of the electrode material to the current collector while preserving superior mechanical properties. The accurate combination of two binders with delightful compatibility in the state-of-art electrode process, can promote a uniform formation of the SEI on the hard carbon surface enriched in inorganic components, which can ensure long-term interfacial stability and suppresses excessive solvent decomposition and facilitates Na + transfer at the interface. Consequently, the initial Coulombic efficiency of the hard carbon anode with 70 %PMMA binder achieves 86 %, with prominent cycling stability (88 % capacity retention over 500 cycles) at a high current density of 1.2 A g −1 . When paired with high loading cathodes to assemble the pouch cell, it also demonstrates stable operational scenarios.

25 ENERGY STORAGE↗

DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗