Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “species dynamics”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Do We Already have Samples of CERES H Chondrite Haliites and the CERES-HEBE Link

We investigate the hypothesis that halite grains in the brecciated H chondrites Zag and Monahans originate from Ceres. Evidence includes mineralogy of the halites consistent with formation on a large, carbonaceous, aqueously active body close to the H chondrite parent body >4 Ga ago. Evidence also includes orbital simularities between 1 Ceres and the purported H chondrite parent body (HPB) 6 Hebe, possibly facilitating a gentle transfer between the bodies. Halite grains in the Monahans and Zag Hchondrites are exogenous to the H chondrite parent body and were transported to the HPB >4 Ga ago. Examination of minerals and carboanceous materials entrained within the halites shows that the halite parent body (HaPB) is consistent with a carbonaceous body [1]. It is probably a large body due to the variety of entrained carbonaceous materials which probably accreted from multiple sources. The halite grains contain intact, HaPB-origin, ancient fluid inclusions indicating that transfer between the HaPB and the HPB was a gentle process resulting in a ΔT of <25 degC. Ejection from the HaPB may have been via cryovolcanic processes similar to those on modern-day Enceladus, which have been interpreted to include halite from spectroscopic observations. The ΔTmax to preserve the brinebearing halite restricts the impact velocity to the HPB at less than ≈350-700 m/s, depending upon the fraction of kinetic energy used heat the sample. [2-6]. Therefore the HaPB and HPB must have shared nearby orbits at the time of the HaPB-HPB transfer. Evidence presented elsewhere indicates asteroid 6 Hebe is a favored candidate for the HPB based on reflectance spectrum similarity with H chondrites and dynamical arguments [7,8]. The modern orbits of Ceres and Hebe are reasonably similar, with aphelion/perihelion of Ceres and Hebe of 2.99/2.55 and 2.91/1.94 AU, respectively. Initial calculations indicate an approximate mean infall velocity of 1.20 to 1.38 km/s. While higher than 350- 700 m/s, the orbits may have been more favorable >4 Ga ago. Additional dynamical factors need to be investigated. A combination of factors suggests Ceres as the HaPB. It is a carbonaceous body with suggestions of past aqueous activity [9], which is consistent with the mineral species found in H chondrite halites. Ceres is also a large body capable of accreting the range of carbonaceous materials noted [5]. It is relatively near to purported HPB Hebe, which is required to preserve halite fluid inclusions. The above evidence defines a hypothesized scenario featuring ejection of halite grains from Ceres onto Hebe. Halite was then entrained in H-chondrite near-surface breccias and ejected from Hebe for transport to Earth.

Fries, Marc D.↗
DOE OSTI2021

Structure and Stability of Pristine and Carboxylate-Covered Anatase TiO 2 (001) in Aqueous Environment

The interactions of carboxylic acids and water with TiO 2 surfaces are important in applications ranging from solar cells to biomedical devices. Here we focus on the aqueous interface with the minority (001) surface of anatase TiO 2 (A-001) and the behavior of adsorbed formic and acetic acid monolayers at this interface. We investigated the structure and stability of the pristine reconstructed and formic/acetic acid covered A-001 surfaces in contact with water using density functional theory (DFT) calculations and ab initio molecular dynamics (AIMD) simulations. The (1 × 4) reconstruction of the pristine surface is found to be stable in aqueous environment, within the time scale of our simulation. Here, carboxylic acids adsorb in deprotonated bidentate (BD) form on A-001, with the dissociated proton transferred to a surface oxygen to form a bridging hydroxyl. Of the two possible configurations, BD bridging and BD straddling, of the adsorbed species, the latter is found to rapidly transform to a monodentate structure during our simulations. Further investigation of mixed acetate–formate monolayers on A-001 in water indicates that also BD bridging species can become unstable at the boundaries between formate and acetate-covered regions, transforming to a monodentate form that does not prevent water adsorption on the TiO 2 surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1980

Physical and dynamical conditions in interstellar clouds

Techniques for measuring the local and global parameters that describe molecular clouds are presented and illustrated by selected results. The local parameters which characterize a given location in a cloud are gas kinetic temperature, dust temperature, total density, magnetic field strength, abundance of species i, volume cooling and heating rates for gas and dust, and the thermal and turbulent velocities. Global properties which characterize the cloud as a whole include the size expressed as a length or area, the orientation and shape, the mass, and the integrated heating and cooling rates. If the cloud as a whole is collapsing or rotating, then the collapse velocity and the rotational velocity or angular momentum are global properties. An intermediate class of parameters which shares features of local and global properties includes the column density of species i, the total column density, and the average density.

Evans, N. J., II↗
DOE OSTI2024

A transient site balance model for atomic layer etching

We present a transient site balance model of plasma-assisted atomic layer etching of silicon (Si) with alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ). Molecular dynamics (MD) simulation results are used to provide parameters for the model. The model couples the dynamics of a top monolayer surface region ('top layer') and a perfectly mixed subsurface region ('mixed layer'). The differential equations describing the rates of change of the Cl coverage in the two layers are transient mass balances. Model predictions include Cl coverages and rates of etching of various species from the surface as a function of Cl 2 or Ar + fluence. The simplified phenomenological model reproduces the MD simulation results well over a range of conditions. Comparing model predictions directly to experimental optical emission spectroscopy data, as reported in a previous paper, provides further evidence of the accuracy of the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Lithium degradation in lithium–sulfur batteries: insights into inventory depletion and interphasial evolution with cycling

The promise of high energy density lithium–sulfur batteries with long cycle life is currently tempered by the rapid degradation of lithium-metal anodes with cycling. An in-depth understanding of the dynamical behavior in liquid electrolytes, including the mechanisms underlying depletion of lithium inventory and evolution of lithium interphases, is crucial to make Li–S batteries a reality. We use here an anode-free full cell configuration, pairing a Li 2 S cathode with a bare nickel current collector with no lithium metal on it, to quantitatively estimate the lithium inventory loss per cycle. Lithium inventory loss is shown to be the main factor limiting the overall cyclability of Li–S batteries. Time-of-flight secondary ion mass spectrometry measurements on the deposited lithium reveal the presence of substantial metallic lithium even after most of the active lithium inventory has been depleted. The trapped metallic lithium is rendered electrochemically inactive by the growth of a resistive electrolyte decomposition interphase on the lithium surface. The bulk of the deposited lithium is shown to be composed of various fully reduced interphasial components, including several hydrogen-containing species that show a substantial reduction in intensity with cycling. This indicates considerable gas evolution and is also correlated with the loss of lithium inventory. The use of an anode-free full cell configuration provides a framework for accurate assessment of the dynamics of lithium inventory depletion and characterization of the accompanying interphasial evolution with cycling. The insights gained will prove invaluable to the development of strategies for extending the cycle life of energy-dense Li–S batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Coupled Thermal-Hydraulic Analysis and Species Mass Transport in a Versatile Experimental Salt Irradiation Loop (VESIL) using COBRA-TF

With the resurgence of interest in Molten Salt Reactor (MSR) technology among commercial and governmental agencies, there is a need for new experiments and new modelling capabilities to characterize the unique phenomena present in this fluid fuel system. A Versatile Experimental Salt Irradiation Loop (VESIL) is currently under investigation at Idaho National Laboratory. It is designed to be a dynamic fuel-bearing salt loop placed in the Advanced Test Reactor. One of the key phenomena this proposed experiment plans to elucidate is fission product speciation in the fuel-salt and the subsequent effects this can have on fuel-salt properties, fission product transport and interactions, source term generation, corrosion and redox potential control. Specifically, noble gases (Xe & Kr) will bubble out to a plenum or off-gas system and noble metals (Mo, Tc, Te etc.) will precipitate and deposit in specific zones in the loop. This work extends the mass transfer and species interaction models in the general species transport capability of COBRA-TF (CTF) and applies these models to give a preliminary estimation of fission product behavior in the proposed VESIL design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2009

Infrared Heterodyne Spectroscopy and its Unique Application to Planetary Studies

Since the early 1970's the infrared heterodyne technique has evolved into a powerful tool for the study of molecular constituents, temperatures, and dynamics in planetary atmospheres. Its extremely high spectral resolution (Lambda/(Delta)Lambda/>10(exp 6)) and highly accurate frequency measurement (to 1 part in 10(exp 8)) enabled the detection of nonthermal/natural lasing phenomena on Mars and Venus; direct measurements of winds on Venus, Mars, and Titan; study of mid-infrared aurorae on Jupiter; direct measurement of species abundances on Mars (ozone, isotopic CO2), hydrocarbons on Jupiter, Saturn., Neptune, and Titan, and stratospheric composition in the Earth's stratosphere (O3, CIO, N2O, CO2 ....). Fully resolved emission and absorption line shapes measured by this method enabled the unambiguous retrieval of molecular abundances and local temperatures and thermal structure in regions not probed by other techniques. The mesosphere of Mars and thermosphere of Venus are uniquely probed by infrared heterodyne spectroscopy. Results of these studies tested and constrained photochemical and dynamical theoretical models describing the phenomena measured. The infrared heterodyne technique will be described. Highlights in its evolution to today's instrumentation and resultant discoveries will be presented, including work at Goddard Space Flight Center and the University of Koln. Resultant work will include studies supporting NASA and ESA space missions and collaborations between instrumental and theoretical groups.

Kostiuk, Theodore↗
DOE OSTI2024

High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore

Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.

hydrogen plasma↗
NASA NTRS2001

Concentration Field of Reactants and Products Species in a Reacting Vortex Ring

The proposed paper will present experimental and numerical results on the concentration fields of both reactants and products species in a reacting vortex ring that is generated from the interaction between a diffusion flame and a laminar vortex ring. Flame-vortex interactions are canonical configurations used to study the underlying processes occurring in complicated turbulent reacting flows. This type of configuration contains many of the fundamental aspects of the coupling between fluid dynamics and combustion that could be investigated with more controllable conditions than are possible under direct investigations of turbulent flames. The current configuration has been studied experimentally by Chen and Dahm and Chen et al. under microgravity conditions, and by Park and Shin, and You et al. under normal gravity conditions. This configuration is similar to that used in the analyses of Karagozian and Manda of their 2-D vortex pair in which both fuel and entrained oxidizer are present. The vortex ring used in this study is generated by issuing methane into an air environment through the exit of an axisymmetric nozzle. The experiments were conducted under microgravity conditions in order to remove the undesirable effects of buoyancy that can affect both the flame structure and ring dynamics resulting in possibly asymmetric and nonrepeatable interactions. The experimental technique of diode laser wavelength modulation spectroscopy (WMS) is used to measure concentration fields of reactants, CH4 and O2, products, H2O, CO2, OH, and temperature fields which can be inferred from either line pairs of O2 or OH lines. This technique has been investigated previously by Silver and Bomse et al. This is the first time that the technique has been applied to reacting vortex rings under microgravity conditions. The effect of ring circulation and fuel volume on the species concentration fields will be investigated. The experimental results will be compared to the current numerical results, and used to validate the numerical studies. In addition, the existence of burned cores during the interactions will be determined, and the increase in reactant consumption with increased ring circulation will be examined. Numerical studies were also conducted by solving the Navier-Stokes and mixture fraction equations with the assumptions of unity Lewis and Schmidt numbers. Equilibrium chemistry and flamelet libraries were used to obtain the temperature and species mass fraction fields. The numerical results will serve as guidelines in conducting the experimental studies. Ring circulation and fuel volume effects on the interactions and species concentration fields will be investigated and compared to experimental results.

Chen, Shin-Juh↗
DOE OSTI2023

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

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

Revealing Au13 as Elementary Clusters During the Early Formation of Au Nanocrystals

Understanding the formation mechanism of nanocrystals in solution is very fundamental to the development of materials science. For a metal nanocrystal, the cluster mediated formation mechanism is still poorly understood. Particularly, identifying what types of clusters are dominant and how they evolve into a nanocrystal in the early nucleation stage remains a great challenge. Here, using liquid-cell transmission electron microscopy, we directly observe the formation of ultra-small Au clusters (~0.84 nm) in the presence of PAA-Na. These clusters, which correspond to the size of Au13 cluster, coalesce to form nanocrystals. Our molecular dynamics simulations suggest Au13 in an aqueous environment have higher stability when compared to other cluster sizes and provide atomistic details of growth by cluster coalescence. Collectively, our demonstration of Au13 as the dominant species with elaboration of their coalescence kinetics sheds light on nonclassical nanocrystal formation mechanisms and offers useful guidelines for designing innovative pathways for synthesis of nanomaterials.

Au13 clusters, cluster coalescence, liquid-cell TE↗
DOE OSTI2023

Robust Interfacial Effect in Multi-interface Environment through Hybrid Reconstruction Chemistry for Enhanced Energy Storage

Electrochemical-oxidation-driven reconstruction has emerged as an efficient approach for developing advanced materials, but the reconstructed microstructure still faces challenges including inferior conductivity, unsatisfying intrinsic activity, and active-species dissolution. In this report we present hybrid reconstruction chemistry that synergistically couples electrochemical oxidation with electrochemical polymerization (EOEP) to overcome these constraints. During the EOEP process, the metal hydroxides undergo rapid reconstruction and dynamically couple with polypyrrole (PPy), resulting in an interface-enriched microenvironment. We observe that the interaction between PPy and the reconstructed metal center (i.e., Mn > Ni, Co) is strongly correlated. Theoretical calculation results demonstrate that the strong interaction between Mn sites and PPy breaks the intrinsic limitation of MnO 2 , rendering MnO 2 with a metallic property for fast charge transfer and enhancing the ion-adsorption dynamics. Operando Raman measurement confirms the promise of EOEP-treated Mn(OH) 2 (E-MO/PPy) to stably work under a 1.2 V potential window. The tailored E-MO/PPy exhibits a high capacitance of 296 F g –1 at a large current density of 100 A g –1 . Here our strategy presents breakthroughs in upgrading the electrochemical reconstruction technique, which enables both activity and kinetics engineering of electrode materials for better performance in energy-related fields.

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

Accumulation of Large Ion Populations with High Ion Densities and Effects Due to Space Charge in Traveling Wave-Based Structures for Lossless Ion Manipulations (SLIM) IMS-MS

The accumulation of very large ion populations in traveling wave (TW)-based Structures for Lossless ion Manipulations (SLIM) has been studied to better understand aspects of ‘in-SLIM’ ion accumulation, and particularly its use in conjunction with ion mobility spectrometry (IMS). A linear SLIM ion path was implemented that had a ‘gate’ for blocking and accumulating ions for arbitrary time periods. Removing the gate potential caused ions to exit and the spatial distributions of accumulated ions were examined. The ion populations for a set of peptides were observed to increase with increased accumulation times until space change effects became significant, after which the peptide precursor ion populations decreased due to growing space charge-related ion activation, reactions, and losses. Ions closest to the gate were preferentially activated, and increasingly so with added storage times, TW amplitude, and the size of the total ion population. For analytical applications, we found that relatively weak TWs in the accumulation/storage region should be utilized to prevent or minimize ion losses and/or ion heating effects that can also lead to fragmentation. Our results further indicate that an accumulation region close to the SLIM entrance is advantageous for both speeding accumulation, minimizing ion heating, and avoiding ion population profiles that result in IMS peak tailing. Importantly, space charge-driven separations were observed for large populations of accumulated species and attributed to the opposing effects of space charge and the TW. In these separations, ion species form distributions or peaks, sometimes moving against the TW, and ordered in the SLIM based on their mobilities. Such separations may offer utility for ion pre-fractionation of ions and increasing the dynamic range measurements, increasing the resolving power of IMS separations by decreasing peak widths for accumulated ion populations, and other purposes benefitting from separations of extremely large ion populations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Exploring the Effects of Geopolitical Shifts on Global Wildlife Trade

Abstract International wildlife trade is a major driver of species extinction and biological invasions. Anticipating environmental risks requires inferences about trade patterns, which are shaped by geopolitics. Although the future cannot be predicted, scenarios can help deal with the uncertainty of future geopolitical dynamics. We propose a framework for generating and analyzing scenarios based on four geopolitical storylines, distinguished by combinations of international trade barrier strength and domestic law enforcement degree across countries supplying and demanding wildlife. We then use historical data on bird trade to classify countries into geopolitical profiles and confirm that trade barriers and law enforcement allow predicting bird trade patterns, supporting our scenarios’ plausibility and enabling projections for future global bird trade. Our framework can be used to examine the consequences of geopolitical changes for wildlife trade and to advise policy and legislation. Reducing demand for wildlife and ameliorating global inequality are key for curbing trade related risks.

Ribeiro, Joana (ORCID:0000000237555785)↗
NASA NTRS1980

Analysis of solid propellant combustion in a closed vessel including secondary reaction

A theory for combustion of solid propellants in a closed vessel is presented allowing for residual exothermic chemical reaction in the bulk of the gas in the vessel. Particular attention is given to propellants exhibiting thick gaseous flame zones such as nitrocellulose, double-base and nitramine propellants. For these, the reaction at high pressures is assumed to involve mainly the oxidation of residual hydrocarbons by NO. It is shown that the direct dynamic coupling between the exothermicity, the molecular weight reduction and the changing pressure can influence the dp/dt-p traces obtained, in a manner not directly related to mass burning rate of the solid. Energy and species conservation equations are derived for the bulk of the vessel in differential form; the system is solved numerically. The results show the effect of extended chemical reaction upon measurable combustion characteristics such as dp/dt-p and burn rate pressure exponent, demonstrating its potential importance in interpretation of closed vessel firing data, depending on the pace of the residual gas phase reactions.

Benreuven, M.↗
NASA NTRS1999

The Interactions of a Flame and Its Self-Induced Boundary Layer

The interaction of a laminar flame with its self-generated boundary layer in a rectangular channel was numerically simulated using the two-dimensional, reacting, Navier-Stokes equations. A two species chemistry model was implemented which simulates the stoichiometric reaction of acetylene and air. Calculations were performed to investigate the effects of altering the boundary condition of the wall temperature, the Lewis number, the dynamic viscosity, and the ignition method. The purpose of this study was to examine the fundamental physics of the formation of the boundary layer and the interaction of the flame as it propagates into the boundary layer that its own motion has created.

Ott, James D.↗
DOE OSTI2025

Molecular level insight into non-bilayer structure formation in thylakoid membranes: a molecular dynamics study

In oxygenic photosynthetic organisms, the light reactions are performed by protein complexes embedded in the lipid bilayer of thylakoid membranes (TMs). The organization of the bulk lipid molecules into bilayer structures provide optimal conditions for the build-up of the proton motive force (pmf) and its utilization for ATP synthesis. However, the lipid composition of TMs is dominated by the non-bilayer lipid species monogalactosyl diacylglycerol (MGDG), and functional plant TMs, besides the bilayer, contain large amounts of non-bilayer lipid phases. Bulk lipids have been shown to be associated with lumenal, stromal-side and marginal-region proteins and proposed to play roles in the self-assembly and photoprotection of the photosynthetic machinery. Furthermore, it has recently been pointed out that the generation and utilization of pmf for ATP synthesis according to the ‘protet’ or protonic charge transfer model Kell (Biochim Biophys Acta Bioenerg 1865(4):149504, 2024), requires high MGDG content Garab (Physiol Plant 177(2):e70230, 2025). In this study, to gain better insight into the structural and functional roles of MGDG, we employed all atom and coarse-grained molecular dynamics simulations to explore how temperature, hydration levels and varying MGDG concentrations affect the structural and dynamic properties of bilayer membranes constituted of plant thylakoid lipids. Our findings reveal that MGDG promotes increased membrane fluidity and dynamic fluctuations in membrane thickness. MGDG-rich stacked bilayers spontaneously formed inverted hexagonal phases; these transitions were enhanced at low hydration levels and at elevated but physiologically relevant temperatures. It can thus be inferred that MGDG plays important roles in heat and drought stress mechanisms.

59 BASIC BIOLOGICAL SCIENCES↗
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↗