SEARCH · Engineering Papers
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.
The exchange current observed at porous metal electrodes on sodium or potassium beta -alumina solid electrolytes in alkali metal vapor is quantitatively modeled with a multi-step process with good agreement with experimental results.
As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.
Electrocatalysts encapsulated by nanoscopic overlayers can control the rate of redox reactions at the outer surface of the overlayer or at the buried interface between the overlayer and the active catalyst, leading to complex behavior in the presence of two competing electrochemical reactions. This study investigated oxide encapsulated electrocatalysts (OECs) comprised of iridium (Ir) thin films coated with an ultrathin (2–10 nm thick) silicon oxide (SiO x ) or titanium oxide (TiO x ) overlayer. The performance of SiO x |Ir and TiO x |Ir thin film electrodes towards the oxygen evolution reaction (OER) and Fe(II)/Fe(III) redox reactions were evaluated. An improvement in selectivity towards the OER was observed for all OECs. Overlayer properties, namely ionic and electronic conductivity, were assessed using a combination of electroanalytical methods and molecular dynamics simulations. SiO x and TiO x overlayers were found to be permeable to H 2 O and O 2 such that the OER can occur at the MO x |Ir (M = Ti, Si) buried interface, which was further supported with molecular dynamics simulations of model SiO 2 coatings. In contrast, Fe(II)/Fe(III) redox reactions occur to the same degree with TiO x overlayers having thicknesses less than 4 nm as bare electrocatalyst, while SiO x overlayers inhibit redox reactions at all thicknesses. This observation is attributed to differences in electronic transport between the buried interface and outer overlayer surface, as measured with through-plane conductivity measurements of wetted overlayer materials. These findings reveal the influence of oxide overlayer properties on the activity and selectivity of OECs and suggest opportunities to tune these properties for a wide range of electrochemical reactions.
The electrochemical reductive dimerization of small aromatic carbonyl compounds derived from lignocellulosic biomass is a crucial C−C coupling reaction for upgrading small molecules to longchain hydrocarbons, particularly in the synthesis of drop-in sustainable aviation fuels. Although other electrochemical reduction reactions of these reactants (i.e., hydrogenation and hydrogenolysis) have undergone extensive mechanistic investigation, the understanding of dimerization remains relatively underdeveloped. Most importantly, there is a lack of understanding of the selectivity-determining step between dimerization and monomer reduction and critical factors that can affect this step. In this study, we provide a comprehensive mechanistic model to explain the competition between dimerization and monomer reduction of benzaldehyde under various conditions. Our model proposes that the selectivity between dimerization and monomer reduction depends on the competition between desorption of a ketyl radical from the electrode, necessary for dimerization, and further reduction of the ketyl radical to an alcohol on the electrode by proton-coupled electron transfer (PCET). Computationally comparing the adsorption/desorption energy and PCET activation barrier energy is challenging because conventional DFT calculations substantially underestimate the PCET kinetic barriers. In this study, we employed constrained DFTbased configuration interaction (CDFT-CI) to obtain a reliable comparison of these energies. Our mechanistic model was tested and supported by experimental results obtained with four electrodes (Cu, Pb, Bi, graphite), three pH conditions (acidic, neutral, basic), and three potentials. Our study offers a coherent mechanistic foundation that can explain how each of these conditions impacts the desorption and PCET processes and the selectivities for dimerization and alcohol production.
The performance of electrochemical devices, which play a critical role in decarbonization efforts, is often governed by proton-coupled electron transfer reactions at the electrode–electrolyte interface. These reactions are highly sensitive to the complex and dynamic microenvironment present at the electrode surface. However, characterizing this environment─particularly monitoring interfacial pH and its evolution under reaction conditions─remains challenging, necessitating the development of advanced analytical tools. Here, in this study, we introduce in-line continuous flow Raman spectroscopy (CFRS) as a spectroelectrochemical platform for quantifying interfacial pH swings generated during water-splitting. By monitoring phosphate ion speciation and controlling the hydrodynamics with a flow cell, we measure pH swings as a function of current density, flow rate, and distance from the electrode. Comparison with theoretical models reveals the impact of bulk pH, boundary layer thickness, and bubble dynamics at high current densities. Collectively, these findings establish CFRS as a platform for quantitatively investigating pH dynamics, offering critical insights for advancing electrochemical energy conversion technologies.
Three investigations conducted aboard Skylab IV and Apollo-Soyuz involved phenomena that are of interest to the biochemistry community. The formaldehyde clock reaction and the equilibrium shift reaction experiments conducted aboard Apollo Soyuz demonstrate the effect of low-g foams or air/liquid dispersions on reaction rate and chemical equilibrium. The electrodeposition reaction experiment conducted aboard Skylab IV demonstrate the effect of a low-g environment on an electrochemical displacement reaction. The implications of the three space experiments for various applications are considered.
Abstract not provided
A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.
Environmental and energy considerations demand that the Haber-Bosch process for reducing N 2 to NH 3 be replaced with electrochemical ammonia synthesis where the H atoms come from water instead of from H 2 . But a practical realization of electrochemical N 2 reduction reaction (NRR) requires the development of new generation electrocatalysts with low overpotential and high Faraday efficiency (FE). A major problem here is that the hydrogen evolution reaction (HER) competes with NRR. Herein, we consider new generation dual-site catalysts involving two different metals incorporated into a novel two-dimensional C 3 N–C 2 N heterostructure that provides a high concentration of well-defined but isolated active sites that bind two distinct metal atoms in a framework that facilitates electron transfer. We report here the mechanism and predicted kinetics as a function of applied potential for both NRR and HER for the (Fe–Ru)/C 3 N–C 2 N dual atom catalyst. These calculations employ the grand canonical potential kinetics (GCP-K) methodology to predict reaction free energies and reaction barriers as a function of applied potential. The rates are then used in a microkinetic model to predict the turn-over-frequencies (TOF) as a function of applied potential. At U = 0 V, the FE for NRR is 93%, but the current is only 2.0 mA cm –2 . The onset potential (at 10 mA cm –2 ) for ammonia on Fe–Ru/C 3 N–C 2 N is –0.22 V RHE . This leads to a calculated TOF of 434 h –1 per Fe–Ru site. In conclusion, we expect that the mechanisms for NRR and HER developed here will help lead to new generations of NRR with high TOF and FE.
Ion exchange is a powerful method to access metastable materials with advanced functionalities for energy storage applications. However, high concentrations and unfavourably large excesses of lithium are always used for synthesizing lithium cathodes from parent sodium material, and the reaction pathways remain elusive. Here, using layered oxides as model materials, we demonstrate that vacancy level and its corresponding lithium preference are critical in determining the accessible and inaccessible ion exchange pathways. Taking advantage of the strong lithium preference at the right vacancy level, we establish predictive compositional and structural evolution at extremely dilute and low excess lithium based on the phase equilibrium between Li 0.94 CoO 2 and Na 0.48 CoO 2 . Such phase separation behaviour is general in both surface reaction-limited and diffusion-limited exchange conditions and is accomplished with the charge redistribution on transition metals. Guided by this understanding, we demonstrate the synthesis of Na y CoO 2 from the parent Li x CoO 2 and the synthesis of Li 0.94 CoO 2 from Na y CoO 2 at 1-1,000 Li/Na (molar ratio) with an electrochemical assisted ion exchange method by mitigating the kinetic barriers. Our study opens new opportunities for ion exchange in predictive synthesis and separation applications. Ion exchange is a powerful method to access metastable materials for energy storage, but identifying lithium and sodium interchange in layered oxides remains challenging. Using such model materials, vacancy level and corresponding lithium preference are shown to be crucial for ion exchange pathway accessibility.
Photoelectrochemical generation of a potent organic hydride donor at silicon is demonstrated. Two different oxide-coated p-type silicon photoelectrodes reduced 1,2,3,5,6-pentamethyl-1H-benzo[d]imidazol-3-ium hexafluorophosphate, [PMBI][PF 6 ], to its corresponding imidazole hydride, PMBIH, in the presence of a proton source. Under 1 sun illumination, in acetonitrile with 1,5-diazabicyclo[4.3.0]non-5-ene (DBN) buffer, the p-Si photoelectrodes convert PMBI + to PMBIH with good Faradaic efficiencies (FEs): 78% FE at −2.3 V vs Fc +/0 for Si|TiO 2 and 83% FE at −2.6 V vs Fc +/0 for Si|SiO 2 (where Si|SiO 2 represents silicon coated with an oxide layer). Generally, the Si|TiO 2 catalyzed the reaction at milder potentials than Si|SiO 2 , but the Si|SiO 2 had better selectivity for PMBIH generation over H 2 evolution than Si|TiO 2 . In light of prior studies of these photoelectrodes, the differences in selectivity and onset potential suggest a difference in mechanism, likely an outer-sphere electron transfer (ET) mechanism at the SiO 2 surface and potentially a proton-coupled ET process at the TiO 2 surface. To help understand reaction efficiency and identify areas of improvement, a thermochemical model for understanding net hydride transfer from the semiconductor to an acceptor in solution was developed. We find that the reactions in the present system are quite downhill. The high overpotentials (even while maintaining selectivity over H 2 evolution) emphasize the need for improved catalysts. Furthermore, this approach to evaluate the thermodynamics of net hydride transfer should be broadly valuable for electrochemical and photoelectrochemical processes.
Achieving high current densities without thermal performance degradation at high temperatures is one of the main challenges for enhancing the competitiveness of photo-electrochemical energy storage systems. We describe a system that overcomes this challenge by incorporating an integrated photoelectrode with a redox flow cell, which functions as a coolant for the excess heat from the photo-absorber. We perform quantitative analyses to theoretically validate and highlight the merit of the system. Practical operation parameters, including daily temperature and redox reaction kinetics, are modeled with respect to heat and charge transfer mechanisms. Our analyses show a profound impact on the resulting solar-to-chemical efficiencies and stored power, which are 21.8% higher than that of a conventional photovoltaic-assisted energy storage system. This paves the way for reassessing the merit of photovoltaic-integrated systems, which have hitherto been underrated as renewable energy storage systems.
The oxidation of iodide in NaI-AlBr 3 , NaI-AlCl 3 , and NaI-GaCl 3 molten salts was analyzed using simulation software to extract relevant kinetic parameters. The experimental oxidation potentials were ordered AlCl 3 < AlBr 3 < GaCl 3 , with higher oxidation potentials correlating with softer Lewis acidity of the metal halide. An iodide oxidation and metal halide speciation model was developed and simulated to fit the electrochemical response, enabling determination of electrochemical charge transfer parameters and chemical equilibrium constants. NaI-AlBr 3 displayed the fastest electron transfer rates yet showed the lowest current densities. All salts revealed smaller than expected current densities, explained by equilibrium between various species, where some are not electrochemically active at the studied potentials. These equilibrium reactions are due to the various metal halide species, controlling the reactant concentration of iodide and the resultant current. We hypothesize the electrochemically active iodide species, present as a metal halide monomer (MX 3 I − ), is decreased dramatically from the expected concentration, sequestered as a more stable metal halide dimer species (M 2 X 6 I − ) with a higher oxidation potential. Traditional Tafel analysis of the experimental data supports the validity of the simulations. These results increase understanding of iodide oxidation in low-temperature Lewis acidic molten salts and inform task-specific molten salt design.
The mission of NECCES is to develop an understanding of how key electrode reactions occur in cathode materials for Li based batteries, and how they can be controlled to improve electrochemical performance, from the atomistic level to the macroscopic level through the life-time of the operating battery. This mission will allow the determination of the ultimate limits of intercalation reactions for chemical energy storage. The approach and specific research goals of NECCES are to: 1. Attain reversible multi-electron transfer in a cathode material using lithium, using the model compound ε -VOPO 4 ; 2. Close the gap between the theoretical and practical energy density for intercalation compounds, using the model compound LiNi 0.85 Co 0 .10Al 0.05 O 2 ; 3. Understand the transition from cation redox to anion redox electrochemistry; 4. Understand performance limiting transport in positive electrode structures from the local through the meso to the macroscale.
Insufficient durability of commercial catalysts is a persistent issue for PEFC applications.1 Herein, Pt on Nb-doped-TiO 2 is reported as a highly durable electrocatalyst on which the direct 4-electron reduction of oxygen to water is more facile compared to Pt/C. All performance metrics are reported comparing Pt/Nb-TiO2 with a 15% Pt loading, against a commercial Pt/C catalyst (46.5% Pt loading, Tanaka, K. K.). Nb-doped-TiO 2 with high surface area and high electronic conductivity was synthesized using the supercritical fluid method. Initially, the durability of the catalyst was characterized using accelerated stability tests (ASTs) involving 10,000 high potential cycles (DOE/FCCJ protocol) and the Pt/Nb-TiO 2 was found to retain 78% of its initial electrochemically active surface area (ECSA) compared to the 57.6 % retained by Pt/C. These observations were in excellent agreement with previous reports that the Pt particle size of Pt/C grew from 2nm to 8nm during the AST protocol along with severe corrosion and amorphization of the carbon surface.2, 3 TEM and XPS studies of the Pt-Nb-doped-TiO 2 catalyst showed that the Pt particle size grew from 3nm to 6nm and the Nb(IV) and Ti(III) in the support were oxidized to Nb(V) and Ti(IV) after the durability test. Thus, the improvement in Pt/Nb-TiO 2 ECSA retention was attributed to the lower extent of particle growth and lack of oxidative support loss upon oxidation as compared to Pt/C. The oxygen reduction reaction (ORR) performance was characterized by linear polarization using a rotating disk electrode (RDE). The electrochemical surface areas of Pt/Nb-TiO 2 and Pt/C were found to be 48m2·g -1 and 83m2·g -1 respectively, and the mass activity for the ORR at 0.9V vs. RHE were found to be 150 mA·mg-1Pt and 124 mA·mg-1Pt respectively. The improved mass activity on Pt/Nb-TiO 2 was attributed to strong metal support interaction (SMSI) between the Nb-TiO 2 support and the Pt catalyst based on the 625 meV decrease in the binding energy of the Pt4f x-ray photo-electron spectroscopy (XPS) peaks of Pt/Nb-TiO 2 compared to Pt/C. To quantify the impact of the SMSI, a kinetic model was applied to calculate the elementary reaction rate constants for the various steps of the ORR on both catalysts. The reaction rate constant (k1) for the direct 4-electron transfer pathway to produce H 2 O was significantly larger in Pt/Nb-TiO 2 as compared to Pt/C. Thus, the reduction in the electron binding energy as observed in the XPS was found to aid the facile filling of the higher energy 2p orbitals of O 2 thereby leading to improved 4-electron transfer kinetics and improved overall activity.
Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.
Solar-driven water splitting requires sufficient photovoltage to drive both water oxidation and proton reduction. Understanding the factors driving and limiting photovoltage generation is therefore crucial to optimizing photoelectrode design but has proven challenging to determine under operando conditions for photoanodes driving slow multiredox reactions such as water oxidation. In this work, operando optical spectroscopy is employed to measure the hole quasi-Fermi level (E F,p ) position in model hematite photoanodes as a function of applied bias and light intensity. The quasi-Fermi level splitting determined from these data are shown be in excellent agreement with the directly measured photovoltages, demonstrating the primarily electrochemical rather than primarily electrostatic origin of photovoltage in these photoelectrodes. E F,p pinning is observed at low light intensities and biases, indicative of hole trap states lying ∼0.2 eV above the valence band edge with a density of ∼1 nm −2 . Hole accumulation in these trap states is correlated with first order water oxidation. At higher light intensities and/or more anodic bias, E F,p becomes unpinned, assigned to saturation of these trap states, and correlated with the onset of third order water oxidation to molecular oxygen. Comparison with rate law analyses for other photoanodes indicates that such hole trap states may be a ubiquitous feature of metal oxides and suggests that materials processing strategies to suppress the density of such states would be a promising strategy to enhance photoanode performance.
Liquid water in porous media hinders the transport of reactant to catalyst layer, where the electrochemical reactions occur, which in turn, reduces performance and lifetime of a proton exchange membrane fuel cell (PEMFC). In this study, two types of gas diffusion layer (GDL) materials are used to study the two-phase saturation distribution in a PEMFC. The correlations between the effective and anisotropic transport properties and porosity of the two GDLs are solved by the pore scale model. In addition, the anisotropic liquid water permeability and the relationship between saturation pressure and capillary pressure is determined by using the Lattice-Boltzmann method. Fuel cell performance as well as liquid water distribution in the GDL using neutron radiography are obtained for validation. Under wet conditions, the cell with Freudenberg GDL performs better than that using Toray GDL, especially at high current densities. The results of the two-phase model simulation show that the peak water saturation for Toray GDL occurs in the catalyst layer and can reach 40%, while the peak water saturation for Freudenberg GDL occurs inside the GDL directly under the ribs and only reaches 25% saturation. The combined results provide key insights to enable high power density operation of a PEMFC through GDL material optimization.