Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive flow”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Electrified vapour deposition at ultrahigh temperature and atmospheric pressure for nanomaterials synthesis

Vapour-phase synthesis methods have shown promise for the scalable synthesis of nanomaterials and coatings. However, the vaporization of different precursors for the synthesis of a broad nanomaterial space, particularly at atmospheric pressure, while maintaining compositional and structural control of the final product is challenging. Here we report the generation of an ultrahigh-temperature atomic vapour at atmospheric pressure based on electrified heating, for the growth of multi-elemental nanomaterials and thin films. This process relies on a reactor design whereby solid-state precursors are vaporized within a semi-confined space beneath an electrified heater that can reach ~3,000 K. The proximity of the heater rapidly breaks down the bonds of metal salt precursors and decomposes them into an atomic vapour that expands into a high-temperature (>2,000 K), highly reactive and high-flux vapour (10 21 –10 22 atoms per cm 2 per second) that travels upwards in a directional flow. When mixed with entrained ambient gases, the highly reactive atomic species rapidly nucleate and grow into the desired final products, including alloys, oxides, sulfides and thin films, which can be deposited on a low-temperature substrate. This EVD approach can synthesize a broad range of functional nanomaterials at atmospheric pressure, including single-phase multi-elemental nanomaterials formed under thermodynamically non-equilibrium conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reactive Transport Modeling of Anthropogenic Carbon Mineralization in Stacked Columbia River Basalt Reservoirs

Numerical simulation of CO2 storage in basalts and related reactive lithologies requires modeling complex, coupled hydrologic and chemical processes, including multi-phase flow and transport, partitioning of CO2 into the aqueous phase, and chemical interactions with aqueous fluids and rock minerals. We conducted reactive transport simulations of the Wallula pilot-scale CO2 injection into the flow tops of the Grande Ronde Basalt using our PNNL STOMP-CO2 simulator with the ECKEChem reactive module. Our mineralization simulation of the ~1,000 tons of injected CO2 into the interflow zones was based on the hydrologic transport model we previously developed. For this work, the simulations considered geochemical reactions involving the basalt components, precipitates, formation brine, and injected CO2. In our benchmark case, carbonate minerals precipitated, resulting in ~20% of the CO2 being mineralized in 10 years. Increasing the reaction rate of a single primary mineral phase (clinopyroxene) by an order of magnitude resulted in a carbon mineralization reaction extent of ~90% over the same time interval. Based on these initial sensitivity analysis results, it is clear that a thorough understanding of primary mineral dissolution rates is required for accurately predicting long-term fate and transport of injected CO2 into basalt formations. Our reactive transport numerical simulations will be key components of commercial-scale CO2 storage operation permitting, de-risking, and optimization in mafic and ultramafic reservoirs.

Cao, Ruoshi↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

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

Heterogeneity of neutrophils in arterial hypertension

Highlights: • Normal-density and low-density neutrophils from patients with hypertension influx more Na{sup +} compared to normotensive individuals. • Normal-density neutrophils from patients with hypertension generate low levels of reactive oxygen species. • Normal-density and low-density neutrophils from patients with hypertension exhibit lower levels of apoptosis. • Normal-density neutrophils from patients with hypertension maintain an activated and degranulated phenotype. • Low-density neutrophils from patients with hypertension show atypical mitochondrial morphology. Cellular heterogeneity and diversity are recognized to contribute to the functions of neutrophils under homeostatic and pathological conditions. We previously suggested that the chronic inflammatory responses associated with hypertension (HTN) are related to the participation of different subpopulations of neutrophils. Two populations of neutrophils can be obtained by density gradient centrifugation: normal-density neutrophils (NDN) and low-density neutrophils (LDN). However, the lack of standardized functional protocols has limited phenotypic characterization and functional comparisons of LDN and NDN. Based on their capability to incorporate Na{sup +}, maturity and activation stage, we characterized NDN and LDN in blood samples from ten patients with HTN and ten healthy individuals (HI) using flow cytometry. We compared the levels of reactive oxygen species (ROS), generation of neutrophil extracellular traps (NETs) and levels of apoptosis in NDN and LDN. In general, the NDN and LDN subpopulations from patients with HTN exhibited higher levels of sodium influx and ROS, and lower levels of apoptosis than the corresponding NDN and LDN subsets from HI. Transmission electron microscopy revealed NDN and LDN from patients with HTN exhibited alterations to mitochondrial morphology and fewer cytoplasmic granules than the corresponding HI subpopulations. Our results indicate both the NDN and LDN subpopulations enhance the effects of inflammation that contribute to the pathophysiology of HTN. Further detailed studies are required to characterize the events during ontogeny of the myeloid lineage that result in the diverse phenotypic characteristics of each subpopulation of LDN and NDN.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Tracking the phase changes in micelle-based NiGa nanocatalysts for methanol synthesis under activation and working conditions

The hydrogenation of CO 2 into high energy density fuels such as methanol, where the required H 2 is obtained from renewable sources, is of utmost importance for a sustainable society. In recent years, NiGa alloys have attracted attention as promising catalyst material systems for the hydrogenation of CO 2 into methanol at ambient pressures. They thus represent an energy-saving alternative to the Cu-based catalysts employed in today's catalytic industry that require high pressures for the CO 2 hydrogenation. However, the underlying reaction mechanisms for the NiGa system are still under debate. One of the challenges here is to unravel the evolution and coexistence of the different species in the heterogeneous NiGa catalyst system under activation and reaction conditions. To shed light on their evolution under activation in H 2 and their catalytic roles under CO 2 hydrogenation working conditions on well-defined Ni 3 Ga 1 and Ni 5 Ga 3 nanoparticle (NP) catalysts, we employed a multi-probe approach in this study. It included advanced machine learning-based analysis of operando X-ray absorption spectroscopy data combined with operando powder X-ray diffraction and near ambient pressure X-ray photoelectron spectroscopy measurements, as well as reactivity studies using bed-packed mass flow reactors. In addition, we employed atomic force microscopy and scanning transmission electron microscopy for structural characterization. Under H 2 activation at 1 bar total pressure, we concluded the formation of metallic Ni, starting for Ni 3 Ga 1 at 300 °C, and for Ni 5 Ga 3 at 400 °C. At higher temperatures, the formation of NiGa alloys follows. The α'-Ni 3 Ga 1 alloy phase is predominantly formed for the Ni 3 Ga 1 NPs, while the coexistence of α'-Ni 3 Ga 1 , δ-Ni 5 Ga 3 and Ga 2 O 3 phases is observed for the Ni 5 Ga 3 NPs after the H 2 activation. The formation of the Ga 2 O 3 phase also results in the presence of excess metallic Ni. Under CO 2 hydrogenation reaction conditions, Ga partially oxidizes again to form a Ga 2 O 3 -rich particle shell for both NP compositions, yet, to a larger extent for the Ni 3 Ga 1 NPs, which, in turn, feature a higher amount of excess Ni. We reveal that metallic Ni is responsible for the high selectivity of the Ni 3 Ga 1 NPs towards the production of methane in our catalytic tests. In contrast, the Ni 5 Ga 3 NPs display a strong selectivity toward methanol production (>92%), more than one order of magnitude higher than that for the Ni 3 Ga 1 NPs, which we ascribe to the presence of the δ-Ni 5 Ga 3 phase.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗
DOE OSTI2022

The importance of freeze–thaw cycles for lateral tracer transport in ice-wedge polygons

Abstract. A significant portion of the Arctic coastal plain is classified as polygonal tundra and plays a vital role in soil carbon cycling. Recent research suggests that lateral transport of dissolved carbon could exceed vertical carbon releases to the atmosphere. However, the details of lateral subsurface flow in polygonal tundra have not been well studied. We incorporated a subsurface transport process into an existing state-of-the-art hydrothermal model. The model captures the physical effects of freeze–thaw cycles on lateral flow in polygonal tundra. The new modeling capability enables non-reactive tracer movement within subsurface. We utilized this new capability to investigate the impact of freeze–thaw cycles on lateral flow in the polygonal tundra. Our study indicates the important role of freeze–thaw cycles and the freeze-up effect in lateral tracer transport, suggesting that dissolved species could be transported from the middle of the polygon to the sides within a couple of thaw seasons. Introducing lateral carbon transport into the climate models could substantially reduce the uncertainty associated with the impact of thawing permafrost.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Automatic Time Step Control to Resolve Hydromechanically Driven Fault Reactivation, Spontaneous Nucleation, and Seismic Arrest

Abstract A physical understanding of the progression from flow‐driven (quasi‐static) poromechanical deformation to dynamic fault rupture is critical to the resilient operations of several engineering systems. These processes are bridged by a progression from fault reactivation to the spontaneous nucleation of unstable sliding. Toward addressing this challenge, novel automatic time step size control methods are developed to enable accurate and efficient simulation of these dynamics and transitions from the first principles. The controllers combine local models for discretization error and Coulomb failure conditions to automatically adjust the time step size across several orders of magnitude. The methods do not require additional empirical or theoretical input and can resolve the pre‐rupture, interseismic, and seismic periods to the allowed accuracy. The computational results reveal that the proposed methods automatically capture the onset of reactivation and nucleation for homogeneous and heterogeneous fields. Hydrodynamic and structural heterogeneity lead to disparate critical nucleation sizes compared to those predicted by theory. The results highlight its potential in predicting induced seismicity in realistic subsurface engineering systems and at practical scales.

Environmental Sciences & Ecology↗
DOE OSTI2025

A quantitative risk assessment framework for fault reactivation in underground hydrogen storage: Coupled simulation and deep learning approach

Underground hydrogen storage (UHS) is emerging as a critical solution for large-scale energy storage. However, like all subsurface fluid injection activities, UHS poses the risk of injection-induced fault reactivation. Accurate risk assessment is essential to ensuring the safety and efficiency of UHS operations. This study presents the development of deep-learning surrogate models for fault reactivation prediction in UHS, trained on a comprehensive database of fully coupled fluid flow-geomechanics simulations. Our findings reveal that analytical models often yield unreliable estimates, with errors up to 54% in the allowable injection pressure, potentially leading to a 40% reduction in UHS operational capacity. The developed surrogate models were incorporated into a quantitative risk assessment (QRA) framework, enabling probabilistic evaluation of fault reactivation risk while accounting for uncertainties in the input variables. Site-specific features, such as horizontal stress gradients, fault’s dip and strike angles, and operational parameters like bottom-hole injection pressure and well-fault distance, were identified as the primary drivers of fault reactivation across various stress regimes. Whereas other hydraulic, geological, and poroelastic reservoir properties were found to have a secondary impact. Notably, we observed that the risk of fault reactivation for a critically oriented fault with a static friction coefficient greater than 0.55 remains below 10% in a normal faulting stress regime. However, the risk significantly increases as the stress regime transitions from normal to strike-slip and ultimately to reverse faulting conditions. These findings underscore the importance of rigorous site characterization and comprehensive QRA evaluations to optimize UHS performance and minimize geomechanical risks.

25 ENERGY STORAGE↗
DOE OSTI2022

The temperatures of ejecta transporting in vacuum and gases

In this work, we measure continuous thermal radiance from evolving clouds of liquid metal fragments ejected into vacuum, nonreactive, and reactive gas. We implement a model for the thermalization of the ejecta and gas and use this to constrain the absolute temperature of the ejecta cloud. This model enables further analyses of ejecta thermal behavior under a variety of conditions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

In situ characterization of material extrusion printing by near-infrared spectroscopy

Material extrusion printing of reactive resins and inks present a unique challenge due to the time-dependent nature of the rheological and chemical properties they possess. As a result, careful print optimization or process control is important to obtain consistent, high quality prints via additive manufacturing. Here, we present the design and use of a near-infrared (NIR) flow through cell for in situ chemical monitoring of reactive resins during printing. Differences between in situ and off-line benchtop measurements are presented and highlight the need for in-line monitoring capability. Additionally, in-line extrusion force monitoring and off-line post inspection using machine vision is demonstrated. By combining NIR and extrusion force monitoring, it is possible to follow cure reaction kinetics and viscosity changes during printing. When combined with machine vision, the ability to automatically identify and quantify print artifacts can be incorporated on the printing line to enable real-time, artificial intelligence-assisted quality control of both process and product. Together, these techniques form the building blocks of an optimized closed-loop process control strategy when complex reactive inks must be used to produce printed hardware.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Coupled Transport, Reactivity, and Mechanics in Fractured Shale Caprocks

Abstract Shales are low‐permeability caprocks that confine fluid, such as CO 2 , nuclear waste, and hydrogen, in storage formations. Stress‐induced fractures in shale caprocks provide pathways for fluid to leak and potentially contaminate fresh water aquifers. Fractured shales are also increasingly considered as resources for CO 2 sequestration, enhanced geothermal, and unconventional energy recovery. Injecting reactive fluids into shales introduces chemical disequilibrium, causing an onset of a series of dissolution, precipitation, and fines mobilization mechanisms. The reactions have rapid kinetics and significant impact on porosity and permeability; consequently, flow and storage properties of caprocks. While previous research has explored the separate effects of these reactions, this study aims to uncover their simultaneous occurrence and collective influence. This study unveils these highly coupled transport and reactivity mechanisms by tracking and visualizing the reaction‐induced alterations in the matrix, microcracks, and fractures of shales over time. We conducted brine injection experiments sequentially at pH 4 and 2 in a naturally fractured Wolfcamp shale sample while simultaneously imaging the dynamic processes using X‐ray computed tomography (CT). CT images are validated by finer resolution images obtained using micro‐CT and scanning electron microscopy. We also tracked the sample permeability and fluid chemistry using brine permeability and inductively coupled plasma mass spectrometry, respectively. Findings show that fluid primarily flowed through fractures, dissolving reactive minerals and mobilizing fines on fracture surfaces. Dissolution of fracture asperities under confining stress resulted in the closing of fractures. Clogging in narrow fracture pathways, caused by fines accumulation, diverted fluid flow into matrix pores.

58 GEOSCIENCES↗
DOE OSTI2023

Fuel Flow Modeling under Thermal Fluidics Considerations within Operational Domains of Molten Salt Reactors

Molten Salt Reactors (MSRs) offer significant versatility supported by their design flexibility and safety advantages due to high temperatures, design customization capabilities, strong negative reactivity feedback effects, and salt chemistry characteristics. Different control schedules and different mass flow rates in the MSR primary systems impact neutronics and thermal hydraulics coupling. The result is two-fold. First, it leads to adaptability in various applications when the customization of a MSR unit is desired. Second, it calls for the quantification of dynamics and safety characteristics to assess the operational domain and make sure reactor stability and inherent safety characteristics are maintained. This paper provides an in-depth transient analysis of the operational implications of MSRs having heat-generating fuel and heat-transporting salt mixed forming liquid fuel salts contained within primary systems. A coupling model was developed to simulate fuel flow in molten salt reactors. A zero-dimensional reactor kinetics model was used with a one-dimensional heat transfer model in the reactor core. The temperature reactivity feedbacks of the fuel salt resulting from the thermal-hydraulics model were used in the reactor kinetics model to complete the coupled code. The model was applied to simulate the Molten Salt Breeder Reactor (MSBR) response at steady state and transients. The developed model was shown to be a suitable tool for the dynamic analysis of molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

LaNi x Fe 1–x O 3 as flexible oxygen or carbon carriers for tunable syngas production and CO 2 utilization

The current study reports LaFe 1–x Ni x O 3–δ redox catalysts as flexible oxygen or carbon carriers for CO2 utilization and tunable production of syngas at relatively low temperatures (~700 °C), in the context of a hybrid redox process. Specifically, perovskite-structured LaFe 1–x Ni x O 3–δ with seven different compositions (x = 0.4–1) were prepared and investigated. Cyclic experiments under alternating methane and CO 2 flows indicated that all the samples exhibited favorable reactive performance: CH 4 and CO 2 conversions varied between 85% and 98% and 70–88%, respectively. While H 2 /CO ratio from Fe-rich redox catalysts was ~2.3:1 in the methane conversion step, Ni-rich catalysts produced a concentrated (~ 93.7 vol%) hydrogen stream via methane cracking. The flexibility of LaFe 1–x Ni x O 3–δ to produce syngas (or hydrogen) with tunable compositions was found to be governed by the iron/nickel (Fe/Ni) ratio. Redox catalysts with higher Fe contents act as a lattice oxygen carrier via chemical looping partial oxidation (CLPOx) of methane whereas those with higher Ni contents function as a carbon carrier via chemical looping methane cracking (CLMC) scheme. XRD analysis and temperature-programmed reactions revealed that both types of catalysts involve the formation of La 2 O 3 and Ni 0 /Ni-Fe phases under the methane environment. The ability to re-incorporate La 2 O 3 and Ni/Fe into a perovskite structure gives rise to oxygen-carrying capacity whereas stable Ni 0 or Ni/Fe phases would catalyze methane cracking without lattice oxygen exchange in the reaction cycles. Here, temperature programmed oxidation and Raman spectroscopy indicated the presence of graphitic and amorphous carbon species, which were effectively gasified by CO 2 to produce concentrated CO. Stability tests over LaFe 0.5 Ni 0.5 O 3 and LaNiO 3 revealed that the redox performance was stable over a span of 50 cycles.

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

Investigation of structural effects of aromatic compounds on sooting tendency with mechanistic insight into ethylphenol isomers

Small aromatic molecules with oxygen-containing functional groups are a promising class of fuel additives, as they can be readily sourced from depolymerized lignin. These oxygenated aromatic compounds (OACs) show a lower sooting tendency than aromatic hydrocarbons, but OACs having alkyl groups such as ethylphenol show a higher sooting tendency than other OACs such as phenol and anisole, despite the oxygen moiety. In this study, we investigate the relationship between chemical structure and soot precursor formation to explain observed differences in the sooting tendency of OACs and to gain insight into how alkyl or oxygenated substituents on the aromatic ring affect soot precursor formation. The weakest bond for 15 aromatic compounds was identified and cleavage of these bonds was shown to generate either benzyl or phenoxy radicals. A linear relationship between standard enthalpy of formation (ΔH f o ) of these radicals and the yield sooting index (YSI) was found, and thus ΔH f ° can be applied as a metric to estimate YSIs of various aromatic compounds; higher ΔH f o of a radical indicates an increase in the radical reactivity and leads to more soot precursor formation. Flow reactor experiments were performed for 2-ethylphenol and 3-ethylphenol to elucidate how ortho and meta substitution effects the sooting tendency. Soot precursors were identified from the experiment and their formation pathways were investigated computationally. 2-ethylphenol produces more oxygenated products than 3-ethylphenol since the ortho position has increased resonance stabilization of radical intermediates, which leads to lower YSI. Overall, these results further inform the selection of potential biomass-derived fuel blendstocks that have favorable sooting tendencies.

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

A Geo‐Structurally Based Correction Factor for Apparent Dissolution Rates in Fractured Media

Abstract Field measurements of apparent geochemical weathering reaction rates in subsurface fractured porous media are known to deviate from laboratory measurements by multiple orders of magnitude. To date, there is no geologically based explanation for this discrepancy that can be used to predict reaction rates in field systems. Proposed correction factors are typically based on ad hoc characterizations related to geochemical kinetic models. Through a series of high‐fidelity reactive transport simulations of mineral dissolution within explicit 3D discrete fracture networks, we are able to link the geo‐structural attributes with reactive transport observations. We develop a correction factor to linear transition state theory for the prediction of the apparent dissolution rate based on measurable geological properties. The modified rate law shows excellent agreement with numerical simulations, indicating that geological structure could be a primary reason for the discrepancy between laboratory and field observations of apparent dissolution rates in fractured media.

58 GEOSCIENCES↗