SEARCH · Search NASA
Results for “scaling reactions”
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.
Improved Solid-State Reaction Method for Scaled-Up Synthesis of Ceramic Proton-Conducting Electrolyte Materials
Protonic ceramic electrochemical cells (PCECs) represent promising technologies in the production of clean electricity, decarbonized hydrogen, chemicals, and fuels at intermediate temperatures. One of the challenges in commercializing PCECs is to produce the electrolyte materials on a large scale. The conventional solid-state reaction (SSR) method suffers from tedious synthesis procedures and low phase purity of the products due to the formation of unwanted secondary phases. Herein, we report an improved SSR (i-SSR) method for kilogram-scale production of high phase-purity electrolyte material BaZr 0.4 Ce 0.4 Y 0.1 Yb 0.1 O 3–δ (BZCYYb4411). In this method, the ball-milled precursor powders are pelletized prior to calcination, which effectively reduces the length of the diffusion paths between the components during perovskite phase formation. In this study, the synthesis procedure and calcination temperature are carefully optimized for efficient and repeatable production based on the powder crystallization behavior. A combined technoeconomic analysis and life cycle assessment modeling suggest that the i-SSR method could reduce the total production cost and greenhouse gas emissions by up to 19% and 39%, respectively, compared to the conventional SSR method. The high quality of the synthesized electrolyte material is corroborated by the excellent electrical conductivity and electrochemical performance of the fabricated PCEC cells.
High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor
Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.
Dissolution of periodic arrays of grains: Upscaling of pore-scale simulations with fast reactions
Upscaling methods are frequently used to derive transport equations at the macroscopic scale from more fundamental equations formulated at the pore scale. These methods typically give a suitable structure for the macroscopic equations and can also provide explicit expressions for the constitutive parameters, such as permeability and dispersion coefficients. Introducing chemical reactions complicates upscaling in at least two important ways. First, the interplay between chemical reactions and transport processes introduces a new length scale, which can be much smaller than the convective or dispersive length scales. A small reactive length scale breaks one of the key assumptions in upscaling; that there is a significant separation in length between the pore-scale and macro-scale processes. The second complication is that if reactions take place at mineral surfaces (dissolution or precipitation) then the pore space itself is evolving in time. In this paper we suggest ways in which these difficulties can be approached, based on analysis of pore-scale simulation data. First, we noticed that the concentration field in successive unit cells has an almost identical spatial variation, with a single scaling factor for each unit cell that is proportional to the incoming reactant flux. Using pore-scale simulations to determine the mass transfer coefficient in a few unit cells, we can calculate the concentration field in the whole domain, even when dissolution is entirely transport limited. Second, we have noticed a time-dependent mapping of the grain shapes from different unit cells. From these observations, we can deduce constitutive relations where the only time-varying parameter is the porosity. Here, we show that a model based on these ideas can quantitatively account for the pore-scale simulation data.
Development of Post-Consumer Textile Waste Degradation Protocol and Related Downstream Transformations: Cooperative Research and Development (Final Report)
Tereform, Inc. (Tereform) is developing molecular deconstruction processes to transform waste materials into chemical building blocks. During the CRADA, Tereform deconstructed real-world post-consumer substrates into chemical monomers, validated their performance on laboratory scales, demonstrated feasibility of the degradation products in downstream transformations, and successfully scaled the reaction to kilogram-scale. These results were used to develop and refine technoeconomic analysis and lifecycle assessments to evaluate the economic feasibility and environmental impacts of the process.
Prediction of local concentration fields in porous media with chemical reaction using a multi scale convolutional neural network
The study of solute transport in porous media is of interest in many chemical engineering systems. Some example applications include packed bed catalytic reactors, filtration devices, and batteries. The pore scale modeling of these systems is time consuming and may require large computing resources, for this reason computational fluid dynamics (CFD) simulations are not practical if a large number of simulations is required, like in multiscale modeling, where a model at a large scale calls for pore scale simulations. It has been shown that neural networks can be trained with a dataset of flow simulations and then predict fields orders of magnitude faster, and with less computational resources, in new domains. However, it is crucial to provide the neural network with an effective description of the domain and the undergoing operating conditions to be able to train models that generalize accurately in unseen samples. Therefore, research is needed to employ neural networks in new complex systems. The appropriate training of a network for predicting coupled flow and solute transport processes is an outstanding problem due to the complex interplay between geometry and operating conditions. In this work, we train a multi scale convolutional neural network (MSNet) with a diverse dataset of simulations of transport and chemical reaction in porous media to predict the local concentration fields in images of porous media. Our dataset contains a wide diversity of sphere pack arrangements under different operating conditions (Péclet and Reynolds numbers). Further, we train a robust model by employing different input descriptors that represent the medium and the different operating conditions of each system. Our trained model is able to provide nearly instantaneous predictions, compared to around twenty hours of the CFD workflow, with less than 3.5% error on new geometries and transport conditions. Thus the model could be easily integrated in a multiscale workflow where fast response is needed.
Revisiting the Effects of Seed Liquid Water Content on Isoprene Secondary Organic Aerosol Growth
Aerosol liquid water content (LWC) has been shown to enhance partitioning of water-soluble oxygenated organic vapors into the aerosol particles and facilitate aqueous chemistry, impact SOA formation, composition, and physicochemical properties. Here, this study examines the effect of aerosol liquid water content (LWC) on isoprene SOA formation using effloresced and deliquesced ammonium sulfate ((NH 4 ) 2 SO 4 ) and sodium chloride (NaCl) seed particles in an environmental chamber with photo-oxidation extending over hours under moderate relative humidity (50%-60%). Our findings indicate that isoprene SOA mass yield is not significantly influenced (p value >> 0.05) by the phase state or the chemical identity of the seed particles within a 3–4-hour reaction time scale. Detailed molecular analysis using advanced high resolution mass spectrometry indicates that while both effloresced and deliquesced seed types produced similar molecular formula distributions, deliquesced (NH 4 ) 2 SO 4 slightly favored the production of less volatile compounds. Depth-profiling analysis by single particle mass spectrometry confirmed that isoprene SOA predominantly located on the surfaces of particles with core-shell morphologies, with partial dissolution of the effloresced inorganic salt observed even below the deliquescence relative humidity (DRH) of salt. This finding suggests that assuming dry salt core particles do not interact with SOA coatings, based on salt deliquescence points, may overlook important aqueous chemistry. These results emphasize the critical role of reaction time scale and mixing state on SOA formation and call for further research to understand LWC effects on SOA from different precursors.
Dynamics of disordered intermediates during the two-electron alkaline MnO 2 conversion reaction for grid-scale batteries
Battery technologies beyond Li-ion are likely needed for extensive integration of grid-scale storage. The rechargeable Zn-MnO 2 chemistry has the potential for high sustainability, high safety, and low cost, using Earth-abundant basis materials. In an alkaline electrolyte, the MnO 2 cathode can cycle reversibly if modified by including a Bi-containing additive, although the cycling mechanism remains mostly unknown. This work presents an account of the intermediate species involved in the electrochemical transformation from layered δ-MnO 2 to Mn(OH) 2 and back. During charge, a disordered intermediate with a structure resembling layered β-MnOOH exists stably for an extended period, corresponding to a regime known to have unexpected electrochemical activity of Bi. During discharge, β-MnOOH exists only briefly and is never the majority material, revealing that the cycling mechanism is asymmetric. In conclusion, these findings represent a significant advance in mechanistic knowledge and can enable engineering to develop the system for commercial use.
Catalytic Resonance Theory: The Catalytic Mechanics of Programmable Ratchets
Catalytic reaction networks of multiple elementary steps operating under dynamic conditions via a programmed input oscillation are difficult to interpret and optimize due to reaction system complexity. To understand these dynamic systems, individual elementary catalytic reactions oscillating between catalyst states were evaluated to identify their three fundamental characteristics that define their ability to promote reactions away from equilibrium. First, elementary catalytic reactions exhibit directionality to promote reactions forward or backward from equilibrium as determined by a ratchet directionality metric comprised of the input oscillation duty cycle and the reaction rate constants. Second, catalytic ratchets are defined by the catalyst state of strong or weak binding that permit reactants to proceed through the transition state. Third, elementary catalytic ratchets exhibit a cutoff frequency which defines the transition in applied frequency for which the catalytic ratchet functions to promote chemistry away from equilibrium. All three ratchet characteristics are calculated from chemical reaction parameters including rate constants derived from linear scaling parameters, reaction conditions, and catalyst state. The characteristics of the reaction network’s constituent elementary catalytic reactions provided an interpretation of complex reaction networks and a method of predicting the behavior of dynamic surface chemistry on oscillating catalysts.
Generalized Brønsted‐Evans‐Polanyi Relationships for Reactions on Metal Surfaces from Machine Learning
Abstract Brønsted‐Evans‐Polanyi (BEP) relationships, i. e., a linear scaling between reaction and activation energies, lie at the core of computational design of heterogeneous catalysts. However, BEPs are not general and often require reparameterization for each class of reactions. Here we construct generalized BEPs (gBEPs), which can predict activation energies for a diverse dataset of reactions of C, O, N and H containing molecules on metal surfaces. In a first step we develop a set of descriptors based on scaling relationships that can capture the change in chemical identity of reactants during the reaction. Subsequently, we use the reaction energy, these descriptors and a single descriptor for the surface structure to parameterize machine learning based regression approaches for the prediction of activation energies. The best approach we developed shows a Mean Absolute Error (MAE) of 0.11 eV for the training set (80 % of the data set) and 0.23 eV for the test set (20 % of the data set). The methodology presented here allows to calculate activation energies within fractions of seconds on a typical personal computer and due to its generality, accuracy and simplicity in application it might prove to be useful in transition metal catalyst design.
Dynamics of the gold–silicon eutectic reaction studied at limited length scales using in situ $\mathrm{TEM}$ and $\mathrm{STEM}$
We report the dynamics of the gold–silicon eutectic reaction in limited dimensions were studied using in situ transmission electron microscopy and scanning transmission electron microscopy heating experiments. The phase transformation, viewed in both plan-view and cross-section of the film, occurs through a complex combination of dislocation and grain boundary motion and diffusion of silicon along gold grain boundaries, which results in a dramatic change in the microstructure of the film. The conversion observed in cross-section shows that the eutectic mixture forms at the Au–Si interface and proceeds into the Au film at a discontinuous growth rate. This complex process can lead to a variety of microstructures depending on sample geometry, heating temperature, and the ratio of gold to silicon which was found to have the largest impact on the eutectic microstructure. The eutectic morphology varied from dendrites to hollow rectangular structures to Au–Si eutectic agglomerates with increasing silicon to gold ratio.
Metal-Free, Oxidative Condensation of Heteroaryl- and Aryl Carbaldehydes towards Functionalized Carbonitriles
Functionalization of pyridines with the carbonitrile group for advanced cyclization reactions routinely incorporates metals, toxic oxidants, and cyanide reagents, often in linear synthetic sequences. Here, motivated by the strategic need of a heteroaryl carbonitrile synthon for further elaboration, we disclose an efficient, metal-free, one-pot reaction sequence that uses ammonium acetate as the nitrogen source and bis-acetoxyiodobenzene as the oxidant with low loadings. The developed method tolerates a diverse substrate scope including heteroaryl, aryl, and alkyl substrates with similar efficiency and expands the practitioner's reaction. Method development, optimization, substrate scope, preliminary mechanistic observations, competition experiments, and a scale up reaction are delineated herein.
Impacts on mechanical strength of chemical reactions induced by hydrous supercritical CO2 in Boise Sandstone
Geomechanics experiments were used to assess mechanical alteration of Boise Sandstone promoted by reactions with supercritical carbon dioxide (scCO 2 ) and water vapor. During geologic carbon storage, scCO 2 is injected into subsurface reservoirs, forming buoyant plumes. At brine-plume interfaces, scCO 2 can dissolve into native brines, and water from brines can partition into scCO 2 , forming hydrous scCO 2 . This study investigates the effect of hydrous scCO 2 on the strength of Boise Sandstone. Samples are first exposed to recirculating hydrous scCO 2 for 24 h at 70 °C and 13.8 MPa scCO 2 pressure. Samples are reacted with scCO 2 with added water contents up to 500 mL. After scCO 2 exposure, samples are deformed at room temperature under confining pressures of 3.4, 6.9, and 10.3 MPa. The results demonstrate that hydrous scCO 2 induces chemical reactions in Boise Sandstone, with ions migrating from the solid into the hydrous scCO 2 phase. At the longer time-scales, these reactions could lead to mechanical weakening in the samples; however, on the scale of our experiments, the strength changes are within sample variability. Because the solubility of water in scCO 2 is extremely low (0.008 mol H 2 O per 1 mol CO 2 ), the mineral dissolution of Boise Sandstone was under 0.002 wt.%. Additionally, mineral grains and pore throats in Boise Sandstone are cemented with quartz, which is not susceptible to dissolution at these conditions. Our results indicate that humidity in scCO 2 plumes is unlikely to sustain chemical reactions and induce long term strength changes in quartz cemented sandstones due to resistant mineralogies and low water solubility.
ChemComp: Compiling and Computing with Chemical Reaction Networks
The exponential growth in computing demands driven by scientific computing, data analytics, and artificial intelligence is pushing conventional CMOS-based high-performance computing systems to their physical and energy efficiency limits. As we approach the era of post-exascale computing, disruptive approaches are necessary to overcome these barriers and achieve substantial gains in energy efficiency. Analog and hybrid digital-analog computing systems have emerged as promising alternatives, offering the potential for orders-of-magnitude improvements in efficiency. Among these, biochemical computing stands out as a novel paradigm capable of leveraging the natural efficiency of chemical reactions, which have shown promise in solving optimization problems by converging to steady states. By scaling up reaction networks or reaction vessel sizes, biochemical systems present an opportunity to meet the high-performance demands of modern computing tasks. Despite their promise, significant theoretical and practical challenges remain, particularly in formulating and mapping computational problems to chemical reaction networks (CRNs) and designing viable biochemical computing devices. This paper addresses these challenges by introducing new ideas to ChemComp, a compilation and emulation framework for chemical computation. This work describes the mechanisms through which solutions to ordinary differential equations (ODEs) that can be represented as CRN systems can be achieved. Furthermore, we explain the design principles of an ODE dialect implemented as a multi-level intermediate representation (MLIR) compiler extension that will be coupled with existing infrastructure. We demonstrate the potential of our framework through a case study emulating a simplified chemical reservoir computing device. This work establishes foundational tools and methodologies necessary to harness the computational power of chemistry, paving the way for the development of energy-efficient, high-performance computing systems tailored to contemporary and future computational needs.
Simplified Models of the Bicarbonate Buffer for Scaled Simulations of CO 2 Electrolyzers
Bicarbonate electrolytes are used in a range of chemical processes; however, resolved simulation of these electrolytes is difficult, as disparate reaction time scales lead to numerical stiffness and the formation of fine boundary layers. Based on several physically motivated approximations, we reduce the full set of chemical reactions within a bicarbonate electrolyte to a simpler subset, eliminating the numerical stiffness. We supported this simplification via a two-variable singular perturbation expansion and demonstrated that under neutral conditions (6 < pH < 10) typically found within a fluid-fed CO 2 electrolyzer, the simplifications lead to negligible error. Here, we also discuss two alternative simplifications, one valid at high pH and another valid at arbitrary pH. These simplifications reduce the condition number of the matrices resulting from spatiotemporal discretization by up to 10 orders of magnitude and enable three-dimensional (3D) simulation of CO 2 electrolyzers containing carbonate solutions.
Fracture Intensity Impacts on Reaction Front Propagation and Mineral Weathering in Three‐Dimensional Fractured Media
Abstract Studying reaction front propagation in heterogeneous natural settings is challenging, but numerical simulations can provide insight into the varying spatial and temporal scales of reaction front propagation. Here, the impact of increasing fracture intensity on mineral dissolution rates, and the extent of reaction front propagation is investigated using reactive transport simulations in upscaled discrete fracture network domains with varied fracture intensity. Domain‐averaged dissolution rates vary less than 0.5 log units regardless of the fracture intensity, but the spatial distribution of reactions is controlled by the location and number of dead‐end fractures and the number of connected flowpaths through the domain. Higher fracture intensities lead to more weathering in the domain because of more available mineral for water‐rock interactions. We find that reaction fronts propagate through the primary flowpaths in the first 10,000 years of the simulation for a 10‐m length domain, then propagate into secondary flowpaths and dead‐end fractures between 10,000 and 100,000 years, and finally into the matrix over timescales of hundreds of thousands of years. The domain‐averaged reaction rates decrease through time corresponding to a transition from dissolution in advection‐dominated, fast‐flowing pathways, to dissolution in transport‐limited zones of disconnected fractures and matrix. Matrix dissolution, or dissolution under transport‐limited conditions, is the dominant process at late times in these simulations. The results of these simulations recreate the observed paradox found in nature where highly fractured hillslopes tend to be more weathered but have slower weathering rates, while hillslopes with fewer fractures, are less weathered but have higher dissolution rates.
Mineral reaction kinetics constrain the length scale of rock matrix diffusion
Mass transport by aqueous fluids is a dynamic process in shallow crustal systems, redistributing nutrients as well as contaminants. Rock matrix diffusion into fractures (void space) within crystalline rock has been postulated to play an important role in the transient storage of solutes. The reacted volume of host rock involved, however, will be controlled by fluid-rock reactions. Here we present the results of a study which focusses on defining the length scale over which rock matrix diffusion operates within crystalline rock over timescales that are relevant to safety assessment of radioactive and other long-lived wastes. Through detailed chemical and structural analysis of natural specimens sampled at depth from an active system (Toki Granite, Japan), we show that, contrary to commonly proposed models, the length scale of rock matrix diffusion may be extremely small, on the order of centimetres, even over timescales of millions of years. This implies that in many cases the importance of rock matrix diffusion will be minimal. Additional analyses of a contrasting crystalline rock system (Carnmenellis Granite, UK) corroborate these results.