Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction Kinetics”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Enhancing sorption kinetics by oriented and single crystalline array-structured ZSM-5 film on monoliths

Abstract To enhance the reaction kinetics without sacrificing activity in porous materials, one potential solution is to utilize the anisotropic distribution of pores and channels besides enriching active centers at the reactive surfaces. Herein, by designing a unique distribution of oriented pores and single crystalline array structures in the presence of abundant acid sites as demonstrated in the ZSM-5 nanorod arrays grown on monoliths, both enhanced dynamics and improved capacity are exhibited simultaneously in propene capture at low temperature within a short duration. Meanwhile, the ZSM-5 array also helps mitigate the long-chain HCs and coking formation due to the enhanced diffusion of reactants in and reaction products out of the array structures. Further integrating the ZSM-5 array with Co 3 O 4 nanoarray enables comprehensive propene removal throughout a wider temperature range. The array structured film design could offer energy-efficient solutions to overcome both sorption and reaction kinetic restrictions in various solid porous materials for various energy and chemical transformation applications.

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

Catalytic Strategy for Conversion of Triacetic Acid Lactone to Potassium Sorbate

Here this study shows a new route to produce potassium sorbate (KS) from triacetic acid lactone (TAL), which is a chemical platform that can be biologically synthesized from natural sources. Sorbic acid and its potassium salt (KS) are widely used as preservatives in foods and pharmaceuticals. Three steps are used to produce KS from TAL: 1) hydrogenation of TAL into 4-hydroxy- 6-methyltetrahydro-2-pyrone (HMP), 2) dehydration of HMP to parasorbic acid (PSA), and 3) ring-opening and hydrolysis of PSA to KS. TAL can be fully hydrogenated over Ni/SiO2 to give near quantitative yields of HMP. A three-step reaction kinetics model was developed for dehydration of HMP into PSA. This model was used to show that the highest PSA yield occurs at low temperatures. An experimental PSA yield of 84.2% with respect to TAL was obtained which agreed with the prediction of the reaction kinetics model. KOH was used as a coreactant for the ring-opening hydrolysis of PSA to produce >99.9% yield of KS from PSA. Tetrahydrofuran (THF) was used to purify the TAL derived-KS (TAL-KS). The TAL-KS had a KS purity of 95.5%. The overall yield of TAL-KS with respect to TAL was calculated to be 77.3%. TAL-KS produced in this study had similar antimicrobial activities as commercial KS.

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

Thermite and intermetallic projectiles examined experimentally in air and inert gas environments

Intermetallic (aluminum and zirconium) and thermite (aluminum and molybdenum trioxide) projectiles were launched using a high velocity impact ignition testing system. The experiments were designed to simulate reactivity in high (argon) and low (air) altitude environments. The projectiles were launched into a chamber that included a steel target plate for projectile penetration before impacting a rear witness plate. The chamber was semi-sealed and instrumented for quasi-static pressure data. The results provide an understanding of energy release from the projectile materials and of the environmental influence on performance. The transient pressure traces provide insight into reaction kinetics. A bifurcation in transient pressure rise was an indication of a shift in reaction kinetics from the inherent reactive material to metal oxidation with the environment. The bifurcation was delayed by about 0.15 ms for the intermetallic relative to the thermite, evidence that the thermite reaction proceeded faster upon impact than the intermetallic. The two-step process (impact ignition of the reactive material followed by metal oxidation) was shown to produce higher energy conversion efficiencies than projectiles composed of pure fuel (i.e., aluminum) reported previously. Both reactive materials showed energy conversion efficiencies greater than 30% (for air) and 50% (for argon), and an explanation of underestimated efficiency and energy losses is provided. These results have implications for advancing formulations for ballistic applications. Structural reactive materials can be used to modify the effective reactivity of metal-containing formulations in varied atmospheric environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electrolyte-Assisted Hydrogen Cycling in Lithium and Sodium Alanates at Low Pressures and Temperatures

An investigation of electrolyte-assisted hydrogen storage reactions in complex aluminum hydrides (LiAlH4 and NaAlH4) reveals significantly reduced reaction times for hydrogen desorption and uptake in the presence of an electrolyte. LiAlH4 evolves ~7.8 wt% H2 over ~3 h in the presence of a Li-KBH4 eutectic at 130 °C compared to ~25 h for the same material without the electrolyte. Similarly, NaAlH4 exhibits 4.8 wt% H2 evolution over ~4 h in the presence of a diglyme electrolyte at 150 °C compared to 4.4 wt% in ~15 h for the same material without the electrolyte. These reduced reaction times are composed of two effects, an increase in reaction rates and a change in the reaction kinetics. While typical solid state dehydrogenation reactions exhibit kinetics with rates that continuously decrease with the extent of reaction, we find that the addition of an electrolyte results in rates that are relatively constant over the full desorption window. Fitting the kinetics to an Avrami-Erofe’ev model supports these observations. The desorption rate coefficients increase in the presence of an electrolyte, suggesting an increase in the velocities of the reactant-product interfaces. In addition, including an electrolyte increases the growth parameters, primarily for the second desorption steps, resulting in the observed relatively constant reaction rates. Similar effects occur upon hydrogen uptake in NaH/Al where the presence of an electrolyte enables hydrogenation under more practical low temperature (75 °C) and pressure (50 bar H2) conditions.

25 ENERGY STORAGE↗
DOE OSTI2021

Towards ternary binders involving limestone additions — A review

The review summarises literature to examine the transition from portland limestone cement system to composite ternary binder systems involving limestone. Interest in limestone addition as an ideal component in multicomponent binder systems has surged as evident from the large volume of literature published in the recent past. A ternary blended system, with co-substitution of limestone, has the potential to complement the reaction of the supplementary cementitious materials (SCMs). The direct addition of limestone powder helps to attain higher substitution levels of portland cement clinker, improve early age properties, and supplement SCM's reactivity. However, the dilution of hydrates could hamper the long-term benefits. In this review, the interaction of fine limestone is classified and elaborated under two broad umbrellas: physical and chemical interactions. The physical interactions can manifest in three ways, namely, filler action, shearing action and improved packing, which alters reaction rate and extent at early ages. The chemical interactions also modify the reaction kinetics and phase assemblage due to nucleation of C-S-H on calcite surfaces, preservation of the ettringite phase and formation of carboaluminates. Two different forms of carboaluminate hydrates — hemicarboaluminate and mono-carboaluminate can be present in the hydrate matrix depending on the balance between carbonate ions and aluminates in the pore solution. Several factors such as replacement level, particle size, choice of SCM, its reactivity and reactive aluminates content, sulphate levels, curing temperature, and duration of curing can control the carboaluminate formation, reaction degree of SCMs and chemical interaction from limestone additions. A combination of physical and chemical effects makes fine limestone a potential material for co-substitution with aluminosilicate based SCMs, mainly fly ash, slag, and calcined clay. In this review, the factors affecting limestone-SCM composites are summarised based on a detailed literature survey. The effects of SCM-limestone cement composites on hydration kinetics, reaction chemistry, the reactivity of SCMs, the stability of hydrated phases, and contribution to the physical structure development and macroscopic properties by evaluating hydration and mechanical properties are discussed. The importance of AFm (Al{sub 2}O{sub 3}–Fe{sub 2}O{sub 3}-mono) phases in various deterioration mechanisms in concrete and their influence on performance characteristics in different exposure environment is critically appraised.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A self-adjusting platinum surface for acetone hydrogenation

We show that platinum displays a self-adjusting surface that is active for the hydrogenation of acetone over a wide range of reaction conditions. Reaction kinetics measurements under steady-state and transient conditions at temperatures near 350 K, electronic structure calculations employing density-functional theory, and microkinetic modeling were employed to study this behavior over supported platinum catalysts. The importance of surface coverage effects was highlighted by evaluating the transient response of isopropanol formation following either removal of the reactant ketone from the feed, or its substitution with a similarly structured species. The extent to which adsorbed intermediates that lead to the formation of isopropanol were removed from the catalytic surface was observed to be higher following ketone substitution in comparison to its removal, indicating that surface species leading to isopropanol become more strongly adsorbed on the surface as the coverage decreases during the desorption experiment. This phenomenon occurs as a result of adsorbate–adsorbate repulsive interactions on the catalyst surface which adjust with respect to the reaction conditions. Reaction kinetics parameters obtained experimentally were in agreement with those predicted by microkinetic modeling when the binding energies, activation energies, and entropies of adsorbed species and transition states were expressed as a function of surface coverage of the most abundant surface intermediate (MASI, C 3 H 6 OH*). It is important that these effects of surface coverage be incorporated dynamically in the microkinetic model (e.g., using the Bragg–Williams approximation) to describe the experimental data over a wide range of acetone partial pressures.

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

An improved detailed chemical kinetic model for C3-C4 linear and iso-alcohols and their blends with gasoline at engine-relevant conditions

Propanol and butanol isomers have received significant research attention as promising fuel additives or neat biofuels. Robust chemical kinetic models are needed that can provide accurate and efficient predictions of combustion performance across a wide range of engine relevant conditions. Herein, this study seeks to improve the understanding of ignition and combustion behavior of pure C3-C4 linear and iso-alcohols, and their blends with gasoline at engine-relevant conditions. In this work, a kinetic model with improved thermochemistry and reaction kinetics was developed based on recent theoretical calculations of H-atom abstraction and peroxy radical reaction rates. Kinetic model validations are reported, and the current model reproduces the ignition delay times of the C3 and C4 alcohols well. Variations in reactivity over a wide range of temperatures and other operating conditions are also well predicted by the current model. Recent ignition delay time measurements from a rapid compression machine of neat iso-propanol and iso-butanol [Cheng et al., Proc. Combust Inst. (2020)] and blends with a research grade gasoline [Goldsborough et al., Proc. Combust Inst. (2020)] at elevated pressure (20–40 bar) and intermediate temperatures (780–950 K) were used to demonstrate the accuracy of the current kinetic model at conditions relevant to boosted spark-ignition engines. The effects of alcohol blending with gasoline on the autoignition behavior are discussed. The current model captures the suppression of reactivity in the low-temperature and negative-temperature-coefficient (NTC) region when either isopropanol and isobutanol are added to a research grade gasoline. Sensitivity and reaction flux analysis were performed to provide insights into the relevant fuel chemistry of the C3-C4 alcohols.

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

All-Solid-State Lithium–Sulfur Batteries Enhanced by Redox Mediators

We report redox mediators (RMs) play a vital role in some liquid electrolyte-based electrochemical energy storage systems. However, the concept of redox mediator in solid-state batteries remains unexplored. Here, we selected a group of RM candidates and investigated their behaviors and roles in all-solid-state lithium-sulfur batteries (ASSLSBs). The soluble-type quinone-based RM (AQT) shows the most favorable redox potential and the best redox reversibility that functions well for lithium sulfide (Li 2 S) oxidation in solid polymer electrolytes. Accordingly, Li 2 S cathodes with AQT RMs present a significantly reduced energy barrier (average oxidation potential of 2.4 V) during initial charging at 0.1 C at 60 °C and the following discharge capacity of 1133 mAh g s -1 . Using operando sulfur K-edge X-ray absorption spectroscopy, we directly tracked the sulfur speciation in ASSLSBs and proved that the solid-polysulfide-solid reaction of Li 2 S cathodes with RMs facilitated Li 2 S oxidation. In contrast, for bare Li 2 S cathodes, the solid-solid Li 2 S-sulfur direct conversion in the first charge cycle results in a high energy barrier for activation (charge to ~4 V) and low sulfur utilization. The Li 2 S@AQT cell demonstrates superior cycling stability (average Coulombic efficiency 98.9% for 150 cycles) and rate capability owing to the effective AQT-enhanced Li-S reaction kinetics. This work reveals the evolution of sulfur species in ASSLSBs and realizes the fast Li-S reaction kinetics by designing an effective sulfur speciation pathway

25 ENERGY STORAGE↗
DOE OSTI2022

Investigating Particle Size‐Dependent Redox Kinetics and Charge Distribution in Disordered Rocksalt Cathodes

Abstract Understanding how various redox activities evolve and distribute in disordered rocksalt oxides (DRX) can advance insights into manipulating materials properties for achieving stable, high‐energy batteries. Herein, the authors present how the reaction kinetics and spatial distribution of redox activities are governed by the particle size of DRX materials. The size‐dependent electrochemical performance is attributed to the distinct cationic and anionic reaction kinetics at different sizes, which can be tailored to achieve optimal capacity and stability. Overall, the local charged domains in DRX particles display random heterogeneity caused by the isotropic delithiation pathways. Owing to the kinetic limitation, the micron‐sized particles exhibit a holistic “core‐shell” charge distribution, whereas sub‐micron particles show more uniform redox reactions throughout the particles and ensembles. Sub‐micron DRX particles exhibit increasing anionic redox activities yet inferior cycling stability. In summary, engineering particle size can effectively modulate how cationic and anionic redox activities evolve and distribute in DRX materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Low Catalyst Loading Enhances Charge Accumulation for Photoelectrochemical Water Splitting

Abstract Solar water oxidation is a critical step in artificial photosynthesis. Successful completion of the process requires four holes and releases four protons. It depends on the consecutive accumulation of charges at the active site. While recent research has shown an obvious dependence of the reaction kinetics on the hole concentrations on the surface of heterogeneous (photo)electrodes, little is known about how the catalyst density impacts the reaction rate. Using atomically dispersed Ir catalysts on hematite, we report a study on how the interplay between the catalyst density and the surface hole concentration influences the reaction kinetics. At low photon flux, where surface hole concentrations are low, faster charge transfer was observed on photoelectrodes with low catalyst density compared to high catalyst density; at high photon flux and high applied potentials, where surface hole concentrations are moderate or high, slower surface charge recombination was afforded by low‐density catalysts. The results support that charge transfer between the light absorber and the catalyst is reversible; they reveal the unexpected benefits of low‐density catalyst loading in facilitating forward charge transfer for desired chemical reactions. It is implied that for practical solar water splitting devices, a suitable catalyst loading is important for maximized performance.

Liu, Tianying↗
DOE OSTI2023

Low Catalyst Loading Enhances Charge Accumulation for Photoelectrochemical Water Splitting

Abstract Solar water oxidation is a critical step in artificial photosynthesis. Successful completion of the process requires four holes and releases four protons. It depends on the consecutive accumulation of charges at the active site. While recent research has shown an obvious dependence of the reaction kinetics on the hole concentrations on the surface of heterogeneous (photo)electrodes, little is known about how the catalyst density impacts the reaction rate. Using atomically dispersed Ir catalysts on hematite, we report a study on how the interplay between the catalyst density and the surface hole concentration influences the reaction kinetics. At low photon flux, where surface hole concentrations are low, faster charge transfer was observed on photoelectrodes with low catalyst density compared to high catalyst density; at high photon flux and high applied potentials, where surface hole concentrations are moderate or high, slower surface charge recombination was afforded by low‐density catalysts. The results support that charge transfer between the light absorber and the catalyst is reversible; they reveal the unexpected benefits of low‐density catalyst loading in facilitating forward charge transfer for desired chemical reactions. It is implied that for practical solar water splitting devices, a suitable catalyst loading is important for maximized performance.

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

Fundamental Studies of Bifunctional Catalysts for Tandem Reactions

During the past few years of our funding from DOE Basic Energy Sciences, we have focused on the catalytic conversion of biomass-derived 5-hydroxymethyl furfural (HMF) as a platform molecule for sustainable synthesis of value-added chemicals. The combination of the diene and aldehyde functionalities in HMF enables catalytic production of acetalized HMF derivatives with diol or epoxy reactants to allow reversible synthesis of norcantharimide derivatives upon Diels-Alder reaction with maleimides. Whereas the electron-withdrawing nature of aldehyde group inhibits the reactivity of diene group in HMF to undergo Diels-Alder reactions with dienophiles, by converting the aldehyde group in HMF to an electron-donating hydroxyl group, HMF can append maleimide-based chemicals through Diels-Alder coupling. Acetalization of HMF not only produces an electron-rich diene for Diels-Alder reaction, but it also alters the reversibility of acetal formation and thereby allows for the controlled release of chemicals by hydrolysis. Therefore, we synthesized various norcantharimide derivatives from biomass-derived HMF by acetalization over an acid catalyst, Amberlyst-15, followed by Diels-Alder reaction with maleimide. The norcantharimides release the starting materials by retro Diels-Alder reaction that is triggered by acetal hydrolysis under acidic (≤pH 3) conditions. In subsequent work, we explored the synthesis and properties of functional polyurethanes and polyesters with tunable properties from biomass-derived (HMF)-Acetone-HMF (HAH) monomers. HAH can be selectively hydrogenated over Cu and Ru catalysts to produce partially-hydrogenated HAH (PHAH) and fully-hydrogenated HAH (FHAH). HAH contains functional groups, including hydroxyl, furan, enone, and ketone functionalities, that can be exploited to further tune the polymer properties. The π-electron conjugation between the enone and furan groups in HAH monomer stabilizes the furan group to prevent the uncontrolled degradation of furans during reactions, such as etherification. Therefore, HAH is a renewable monomer that can be used to synthesize polymers, having high molecular weight, symmetric functionalities, an inexpensive production price, and being derived from renewable resources without competition with food resources. We demonstrated the synthesis of functional polyurethanes and a renewable polyester from HAH-derived monomers, and 4,4’-methylenebis(phenyl isocyanate) (MDI). We showed that it is essential to be able to selective hydrogenate C=C double bonds in HAH to be able to achieve effective Diels Alder reactions of the furan ring. Therefore, we studied the hydrogenation at temperatures from 313 – 393 K of HAH over Pd, Ru, and Cu based catalysts. HAH was selectively hydrogenated to produce partially-hydrogenated monomers (PHAH) over Cu and Ru catalysts and to fully-hydrogenated HAH monomers (FHAH) over the Ru catalyst. Pd based catalysts yielded a mixture of partially and fully hydrogenated monomers. Reaction kinetics models were employed to quantify the kinetic behavior for hydrogenation over Ru, Cu, and Pd catalysts. A 5-step pathway exhibited over Pd and Ru catalysts consists of both series and parallel reaction steps, where HAH is both converted to fully hydrogenated products sequentially via series reactions of partially hydrogenated intermediates, as well as converted directly in parallel reactions to form the fully hydrogenated products. In contrast, a 3-step pathway over the Cu catalyst consists only of the consecutive reaction steps, where the final product was formed via series reactions of intermediate products. Additionally, reaction over the Cu catalyst did not hydrogenate the furan rings of the HAH molecule and yielded a different final product than those hydrogenation over Pd and Ru catalysts. Using the results from our detailed reaction kinetics studies, we developed universal conditions for the maximum production of the various hydrogenated products for both batch and plug flow systems. In subsequent work, we employed reaction schemes to optimize the yields of products from selective hydrogenations of HAH in isopropanol solvent and in the presence of liquid water. Reaction schemes consisting of 7, 9, and 11 steps were examined to describe the rates of formation of the observed products and reaction intermediates for hydrogenation of HAH over Ru and Pd catalysts, and a 3-step scheme was studied over Cu catalysts. Rate constants and activation energies were calculated using these reaction schemes, and we then apply these reaction schemes to explore the effects of water addition on the hydrogenation pathways. The effects of water addition to isopropanol (IPA) solvents on the hydrogenation of HAH were markedly different over Pd, Ru, and Cu catalysts. Over the Pd catalyst, the addition of water to IPA increased hydrogenation rates and promoted hydrogenation of furan rings. The addition of water to IPA yielded significant carbon losses over the Ru catalyst, and slowed hydrogenation steps over Cu, while significantly inhibiting hydrogenation of the ketone group. This behavior opened routes toward increased production rates of PHAH=O (a partially hydrogenated form of HAH containing a C=O bond), a product in which the diene groups of the furan rings were not hydrogenated. Addition of water also allowed increased feed concentrations of HAH that were previously not possible in pure IPA solvents.

09 BIOMASS FUELS↗
DOE OSTI2025

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗
DOE OSTI2026

Kinetics Modeling and Reactor Design Study of Glucose-to-Terpenes Cell-Free Conversion

Cell-free systems offer many advantages over traditional biological conversion by eliminating biological growth constraints. It also offers easy manipulation and finetuning of the reaction conditions for each individual enzyme. The conversion of cellulosic glucose to Limonene, a terpene, is a promising pathway for producing fuels and chemicals. Recent advances in developing cell-free systems focuses on bench scale optimization of terpene yield and to demonstrate its feasibility towards commercialization [1,2]. There is significant knowledge gap regarding reaction kinetics of these cell-free systems to further study how it will perform at larger scale. We present here, our studies on reaction kinetics and reactor design implications of cell-free glucose to Limonene conversion to facilitate the further development and commercialization of this process. We developed a novel kinetic model based on the metabolic-network structure of the cell-free system with multi-substrate reversible Michaelis-Menten rate law. To estimate kinetic parameters for this system of rate equations, we employed Bayesian optimization to perform global search with the assistance of gaussian processes to balance exploration and exploitation. The model parameters estimated showed good results compared with experimental data. The estimated parameters were used to perform sensitivity analysis. We found that Hexokinase is one of the most critical enzymes that affect the conversion of the glucose. We also observed that abundance of co-factors is also critical to the conversion of glucose to limonene. We investigated packed bed reactors with enzymes immobilized on the surface of particles to convert glucose stream into Limonene for larger scale production. The reactor design such as particle size, enzyme loading, and flow rate are found to be critical for improving yields. [1] Dudley, Q.M., Nash, C.J. and Jewett, M.C., 2019. Synthetic Biology, 4(1), p.ysz003. [2] Korman, T.P., Opgenorth, P.H. and Bowie, J.U., 2017. Nature communications, 8(1), p.15526.

09 BIOMASS FUELS↗
DOE OSTI2020

Dendrite Growth Morphology Modeling in Liquid and Solid Electrolytes

The main goal of this project is to develop a multi-scale modeling approach that connects micron-scale phase-field models and atomic-scale density functional theory (DFT)-based simulations via parameter- and relationship-passing in order to predict Li-metal dendrite morphology evolution, in both liquid and solid electrolytes. The key hypothesis of the DFT-informed phase-field multiscale modeling approach is that it can capture the electrochemical-mechanical driving forces and incorporate the roles of nano-meter-thin solid electrolyte interphase (SEI) in liquid electrolytes as well as of the microstructures of micro-meter-thick solid electrolytes (SEs) for all-solid-state batteries. In this project, we have formulated and implemented phase-field models to incorporate the electrochemical driving forces in liquid electrolytes and then incorporate mechanical driving forces to simulate dendrite growth in solid electrolytes with resolved microstructures. We have implemented two treatments for the SEI: an explicit model to include the microstructure of the SE or SEI in the phase field model and an implicit model to simulate the impact of nano-meter thick SEI in liquid electrolytes by varying the electrode/electrolyte interfacial properties. The key interfacial properties, including the electronic and ionic transport properties, the charge transfer reaction kinetics, and mechanical properties, were computed by DFT-based calculations. At the DFT-based model, one key advancement is to directly predict the charge transfer reaction kinetics at a complex Li/SEI/electrolyte interface by linking DFT with density functional tight binding (DFTB) calculations. As the main accomplishments, we have demonstrated two successful predictions in both solid electrolyte and liquid electrolyte based on this multiscale approach. The predicted intergranular Li dendrite growth in LLZO revealed the importance of trapped electrons at internal interfaces in the microstructure of LLZO. The predicted electroplating morphology of mossy Li and faceted Mg agreed well with experiments. The insights provided by the multiscale model and the model enabled electrolyte and SEI design will accelerate the development of Li-metal electrode for high energy density batteries, that meet DOE’s target on cell density (>350 Wh/kg) and cost below $100/kWhuse for EV applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Analytical modeling for redox flow battery design

Deeper market penetration of redox flow batteries requires optimization of the cell performance. Though important for optimization, detailed analytical solutions have not been developed for electrolyte flow, mass and charge transport, and reaction kinetics within redox flow batteries. To this end, here we present analytical solutions to active species concentration and over-potential based on advection-diffusion transport for ions and Bulter-Volmer model for interface reaction kinetics. The solutions were validated with results from a finite element model. These solutions were then applied to investigate the relationship between over-potential and state of charge, current density, reaction rate constant, flow velocity, diffusivity, total active species concentration, and electrode structure. Explicit formulas were identified for minimum activation over-potential and limiting current density as well as their dependence on electrolyte properties, operation conditions, and electrode structure. With our new mathematical formulas, this work provides a theoretical framework for flow battery design.

25 ENERGY STORAGE↗
DOE OSTI2020

Flexible and Modular Simultaneous Modeling of Flow and Reactive Transport in Rivers and Hyporheic Zones

Investigations of coupled multiphysics processes in rivers and hyporheic zones have extensively used numerical models. Most existing models use a sequential, one-way coupling between the surface and subsurface domains. Such one-way coupling potentially introduces error. To overcome this, a fully coupled model, hyporheicFoam, was developed using the open-source computational platform OpenFOAM. It captures the coupled flow and multicomponent reactive transport processes within both surface and subsurface domains and across their interface. The coupling between two domains is implemented by mapping conservative flux boundary conditions at the interface through an iterative algorithm. Reactive transport is enabled by specifying a reaction network. To start, we have implemented reaction kinetics following the double Monod-type model with inhibition. The model capability is illustrated through modeling of both conservative and reactive hyporheic flow and transport through dune bedforms. With the novel coupled model, it is now possible to quantify reactions wherein the reactants and products are constantly exchanging between domains and have feedbacks. hyporheicFoam can simulate large, three-dimensional cases owing to the computational flexibility and power offered by the code structure and parallel design of OpenFOAM.

58 GEOSCIENCES↗
DOE OSTI2025

Selective Hydrogenation of Furfural Acetone over a Cu Catalyst: Combined Theoretical and Experimental Study

The selective hydrogenation of biomass-derived hydroxymethyl furfural (HMF)-acetone-HMF (HAH) presents an alternative route to producing highly functionalized polyesters and polyurethanes. While HAH undergoes furan ring hydrogenation over Pd, Ru, and Ni catalysts, furan ring hydrogenation is not observed over Cu. Herein, we combine reaction kinetics experiments and density functional theory calculations to elucidate the selective hydrogenation behavior of HAH over Cu catalysts. We identified furfural acetone (FAc) as a suitable surrogate for modeling HAH hydrogenation over Cu surfaces and performed reaction kinetics experiments between temperatures of 313–393 K and a H 2 partial pressure of 55 bar. Similar to the behavior of HAH, hydrogenation of FAc follows a consecutive two-step hydrogenation pathway over Cu and does not undergo furan ring hydrogenation. The apparent activation energy barriers for hydrogenation of the aliphatic double bond (0.58 eV) and carbonyl (0.43 eV) of FAc measured in a continuous flow reactor setup are consistent with those reported in prior studies for batch HAH hydrogenation. Reaction orders with respect to each reactant, including H 2 and FAc, were determined to be nearly one. Density functional theory (DFT; GGA-PBE-D3) calculations on Cu(111) showed that the hydrogenation of the aliphatic double bond of FAc is more facile than the hydrogenation of the furan ring, which displays weak interactions with the Cu surface. We determined an apparent activation energy barrier for FAc hydrogenation (0.58 eV) that agreed with our DFT predictions (highest barrier for FAc hydrogenation of 0.57 eV). Finally, our DFT calculations further show that weak interactions between the furan ring and Cu surface are responsible for the selective hydrogenation behavior.

Biomass↗