Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic parameters”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A wide range experimental study and further development of a kinetic model describing propane oxidation

This work presents an experimental and kinetic modeling study of propane oxidation. Ignition delay times of propane were measured in a high-pressure shock tube and in rapid compression machines in the temperature range 689 - 1700 K at equivalence ratios of 0.5, 1.0 and 2.0 in 'air', for a wide range of pressures from 20 to 90 bar. CO and H 2 O mole fraction profiles for propane oxidation were measured in a shock tube behind reflected shock waves in the temperature range 1370-1840 K at equivalence ratios of 0.5, 1.0 and 2.0 and at a pressure of approximately 1.3 atm. Moreover, propane oxidation was studied using a jet-stirred reactor coupled to a synchrotron vacuum ultraviolet photoionization mass spectrometer at low temperatures in the range 565 - 690 K and at a pressure of 1 atm. Further, this wide range of experimental datasets for propane oxidation was used to reoptimize and update our previous kinetic mechanisms, AramcoMech3.0 and NUIGMech1.1. In the current mechanism, NUIGMech1.3, the thermo-chemical parameters of all species relevant to low-temperature propane oxidation chemistry, including propyl-peroxyl, hydroperoxyl-propyl, hydroperoxyl-propyl-peroxyl, and carbonyl-hydroperoxide radicals, are updated based on newly calculated values at the CCSD(T)-F12/TZ-F12//B2PLYPD3/TZ///B2PLYP-D3/TZ level of theory. The improvements made in the thermochemical values and in the kinetic parameters for the low-temperature propane oxidation reactions in NUIGMech1.3 result in better model agreement with the new IDTs and speciation data, including carbon monoxide, formaldehyde, propene, acetaldehyde and various minor products such as ethylene, acetic acid, acrolein as well as various hydroperoxide and cyclic ether species.

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

Chemical Modeling of the Capture of Radioiodine Species

Iodine is a common fission product resulting from the transmutation of uranium fuel in nuclear power reactors. The short-lived radioactive isotope, 131 I, with a half-life of around eight days, would be a particular health concern if released into the environment. Capture materials, such as activated carbon, are used routinely in nuclear power plants and medical isotope production facilities to remove this hazard during normal operations to reduce emissions to acceptable levels. Iodine, unlike other common fission products such as noble gases, is challenging because of the wide range of molecular forms and oxidation states (from +VII in IF 7 to -I in iodide) that it displays, which can interconvert both in the gas phase and in the solution phase for example. The work in this project is focused on three computational tasks. In the first, we consider the binding of iodine-containing molecules of varying oxidation states to the major component of activated carbon – pristine graphene-like sheets. The second task focused on molecular binding and possible chemical transformations resulting from the reaction of molecular iodine with chemical functional groups on the edge of the graphene sheets or surface defects in the sheet itself. The final task was associated with radiolytic breakdown of 131 I to 131 Xe in iodine containing molecules and investigated the resulting molecular fragmentation process both in the gas phase and on the graphene-like surface. For each of the tasks, the intent was to learn about the thermodynamics and kinetics of the chemical transformations and use the calculations to provide parameters for input into kinetic models describing these processes for interpretation of experimental data in future work.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Laser-to-proton conversion efficiency studies for proton fast ignition

We assess the conversion efficiency from intense picosecond laser pulses to multi-MeV ion beams for a wide range of laser and target parameters, using 2D kinetic particle-in-cell simulations. Scalings are addressed in a quasi-one-dimensional geometry, leaving out beam divergence. Then, we study the conversion efficiency into a narrow spatial band along the laser axis for flat and hemispherical targets in large-scale 2D simulations. Combining these findings allows us to calculate the energy requirements for ignition of a compressed ICF target with an intense proton beam in a fast-ignition scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

High-Temperature Reaction Kinetics of the e aq – and HO 2 • Radicals with Iron(II) Ions in Aqueous Solutions

Pulsed electron radiolysis was used to determine the chemical reaction kinetics and Arrhenius parameters for iron(II) reactions in aqueous solutions under irradiation. The second-order Fe 2+ reactions with the hydrated electron (e aq – ) and the perhydroxyl radical (HO 2 • ), arising from water radiolysis, were measured to high temperatures using custom-built flow-through cells with a multichannel optical detection system. The reaction with the HO 2 • radical was found to proceed via the formation of a metal-ion adduct species, Fe 2+ –HO 2 • . Additionally, the adduct’s molar extinction coefficient and its first-order decay rate coefficients are reported.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

CaMn 0.9 Ti 0.1 O 3 based redox catalysts for chemical looping – Oxidative dehydrogenation of ethane: Effects of Na 2 MoO 4 promoter and degree of reduction on the reaction kinetics

Reduction kinetics and stability of 20 wt% Na 2 MoO 4 -promoted CaMn 0.9 Ti 0.1 O 3 were investigated for its applications in Chemical Looping – Oxidative Dehydrogenation (CL-ODH) of ethane, a potential alternative for ethylene production with higher efficiency and lower emissions. Here, the present work reports a kinetics model and parameters for a Na 2 MoO 4 -promoted, Ti-doped CaMnO 3 (CaMn 0.9 Ti 0.1 O 3 ) redox catalyst under H 2 and C 2 H 4 . A first-order reaction model provides the best fit for the reduction of Na 2 MoO 4 /CaMn 0.9 Ti 0.1 O 3 under H 2 , while the C 2 H 4 reduction is well described by an Avrami–Erofe’ev model. The activation energy for C 2 H 4 oxidation is approximately three times higher than that for H 2 conversion, showing that the activation of C 2 H 4 is significantly more difficult on the surface of the redox catalyst. The reduction rate of Na 2 MoO 4 /CaMn 0.9 Ti 0.1 O 3 under H 2 at 750 °C is more than two orders of magnitude greater than that under C 2 H 4 , while the reduction rate of unpromoted CaMn 0.9 Ti 0.1 O 3 is comparable under H 2 and C 2 H 4 , showing that the addition of Na 2 MoO 4 effectively suppresses C 2 H 4 combustion relative to H 2 oxidation. The kinetics results for Na 2 MoO 4 /CaMn 0.9 Ti 0.1 O 3 confirm its excellent selectivity towards hydrogen combustion, making it a promising candidate under CL-ODH. Additionally, the stability of the CaMn 0.9 Ti 0.1 O 3 @ Na 2 MoO 4 core-shell structure, which was the underlying reason for the excellent selectivity, was examined under both shallow and deep reductions. It was determined that deep reduction of the redox catalyst, e.g. higher than 80% solid conversion, would lead to loss of sodium and hence to decreased selectivity for hydrogen combustion. In contrast, the core-shell structure was well-maintained, exhibiting excellent performance after 50 redox cycles when deep reduction of the redox catalyst was avoided. This study offers a basis for both the CL-ODH reactor design and redox catalyst optimizations.

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

Mitigation of Molten Salt Corrosion

The desire to increase the peak operational temperature of Concentrated Solar Power (CSP) plants to lower the levelized cost of CSP-derived electricity, has led to consideration of an earth-abundant, low-cost, thermally-stable, molten MgCl 2 -KCl-NaCl (40-40-20 mole%) salt for high-temperature thermal energy storage and heat transfer. The reliable, long-term containment of such a molten salt is challenging, however, due to the corrosion of cost-effective, high-temperature structural metal alloys in this salt (particularly if this salt is contaminated with oxygen-bearing species). This project has been focused on: i) development of equipment for evaluating corrosion kinetics in molten chlorides under well-controlled thermal, fluid dynamic, and atmospheric conditions at CSP-relevant temperatures, ii) use of such equipment to evaluate the influence of forced convection on corrosion kinetics (to determine the role of liquid phase diffusion on such corrosion), and iii) identifying a corrosion-mitigation strategy that can be industrially applied to achieve a projected annual corrosion ≤30 μm at CSP-relevant temperatures. The overall goals of this project have been: i) to understand the rate-limiting step(s) and associated kinetic laws governing corrosion in molten chlorides at 750°C, ii) to determine key operational parameters affecting the corrosion kinetics in molten salts, and iii) to demonstrate the tailoring of such parameters to appreciably reduce such corrosion.

14 SOLAR ENERGY↗
DOE OSTI2026

Mechanistic Insights for Plasma-Catalytic CO 2 Reduction over TiO 2 in a Dielectric Barrier Discharge Reactor

Reaction kinetics experiments coupled with phenomenological kinetic modeling and parameter estimation are used to elicit insights into the mechanism and active sites for the plasma-catalytic dissociation of CO 2 on TiO 2 . Experimental and model insights showed that gas-phase reactions contribute at least two-thirds of the overall product formation at explored conditions; weak temperature dependence, strong sensitivity to specific energy input (SEI), apparent first order in CO 2 , and positive influence of cofed argon (Ar) and oxygen (O 2 ) for the gas-phase contributions all suggest that expected plasma reaction steps such as electron-impact and high-energy collisions are the dominant modes for CO 2 dissociation. The Arrhenius-like expression for gas contributions resulted in a preexponential of 4.40 × 10 –3 s –1 , an E SEI,g of 7.90 × 10 –4 mol/kJ, and an E a,g of 1.00 × 10 –3 J/mol. For surface contributions, the small apparent barrier of 16.3 kJ/mol, relatively weaker dependence on SEI, first-order dependence on CO 2 , and insensitivity to cofed Ar and O 2 all point to CO 2 dissociation on TiO 2 surface facets without vacancies and aided by plasma (leading to vibrationally excited CO 2 and/or a reactive surface with significant surface charge accumulation). The Arrhenius-like expression resulted in a preexponential of 7.81 × 10 –2 s –1 , an E SEI,s of 1.90 × 10 –3 mol/kJ, and an E a,s of 1.63 × 10 4 J/mol. The derived kinetic model further enabled a systematic evaluation of the effect of inputs (plasma power, flow rate, CO 2 inlet concentration, and temperature) to identify process trends and optimal operating conditions.

catalyst↗
DOE OSTI2023

Towards a heat- and mass-balanced kinetic model of TATB decomposition

We report TATB (1,3,5-triamino-2,4,6-trinitrobenzene) was thermally degraded by two small-scale analytical methods – simultaneous differential scanning calorimetry and thermogravimetric analysis (SDT) and a hot-stage microscope with Fourier Transform Infrared (FTIR) analysis capabilities. SDT used ramped heating, isothermal soaking, and thermal pretreatment at various conditions. The heat flow and mass loss were monitored during various treatment conditions to derive chemical decomposition kinetics and Arrhenius parameters. FTIR experiments used isothermal heating, and changes were monitored spectroscopically. Solid samples generated at specific conditions were collected from both methods and were analyzed by DMSO extraction followed by chemical speciation by optical and mass spectrometric methods. Characterization provided the following reaction insights: TATB decreases in a sigmoidal pattern in isothermally heated samples. Other soluble products gradually increase in concentration and then abruptly decline in concentration during the second exotherm, such as diamino-dinitro-benzofurazan and amino-nitro-benzodifurazan. FTIR showed gradual changes in the amino and nitro functionality, shifting positions and decreasing intensity for the first 40 min. Then the solid gradually appeared more like an amorphous C with N incorporated, similar to previous studies on thermally degraded TATB-type materials. Extracted residues (DMSO-soluble components removed) examined by FTIR showed an abrupt change in chemical composition between 40- and 45-min isothermal treatment, indicating early forming solids are different than later forming residues. A reliable mass- and energy-balanced global reaction network must include at least two autocatalytic reactions, either in parallel or series, and at least one must have an explicit initiation reaction having a low activation energy.

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

Observing and Modeling the Sequential Pairwise Reactions that Drive Solid‐State Ceramic Synthesis

Abstract Solid‐state synthesis from powder precursors is the primary processing route to advanced multicomponent ceramic materials. Designing reaction conditions and precursors for ceramic synthesis can be a laborious, trial‐and‐error process, as heterogeneous mixtures of precursors often evolve through a complicated series of reaction intermediates. Here, ab initio thermodynamics is used to model which pair of precursors has the most reactive interface, enabling the understanding and anticipation of which non‐equilibrium intermediates form in the early stages of a solid‐state reaction. In situ X‐ray diffraction and in situ electron microscopy are then used to observe how these initial intermediates influence phase evolution in the synthesis of the classic high‐temperature superconductor YBa 2 Cu 3 O 6+ x (YBCO). The model developed herein rationalizes how the replacement of the traditional BaCO 3 precursor with BaO 2 redirects phase evolution through a low‐temperature eutectic melt, facilitating the formation of YBCO in 30 min instead of 12+ h. Precursor selection plays an important role in tuning the thermodynamics of interfacial reactions and emerges as an important design parameter in planning kinetically favorable synthesis pathways to complex ceramic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Study of Intermittent Turbulence in Stable Arctic Boundary Layers

Abstract Polar boundary layers are difficult to model due to the existence of intermittent turbulence within stable layers. Here we present a case study evaluation of coherent structures in a stable boundary layer observed during a series of flights with an uncrewed aircraft system (DataHawk2) on 19 October 2016 at Oliktok Point, Alaska as part of the ERASMUS (Evaluation of Routine Atmospheric Sounding Measurements using Unmanned Systems) field campaign. During a sequence of five flights over a nine-hour period, 57 profiles of atmospheric properties (0–400 m a.g.l) were collected. Turbulence was identified using derived Richardson Number, temperature structure function parameter, and turbulence kinetic energy dissipation. Throughout all flights on this strongly stable day, intermittent turbulence was observed. These turbulent layers showed well-mixed potential temperature profiles embedded within otherwise stable potential temperature profiles; often resulting in a characteristic staircase pattern. Turbulent layers ranged from 1 to 30 m deep, with most individual layers being 1–2 m deep. Vertical propagation velocities of layers in the lower atmosphere were on the order of a few cm s −1 , typical of non-convective environments. In different regions of the profile, turbulence was driven by a different balance of buoyancy and shear forces, with turbulence in the near surface environment driven by strong shear forces overcoming strong resistance to buoyancy, while turbulence in elevated layers characterized by weaker shear forces overcoming weaker resistance to buoyancy. We discuss the potential of such datasets for improving subgrid parameterizations of small-scale turbulence embedded within stable boundary layers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Solidification modes during additive manufacturing of steel revealed by high-speed X-ray diffraction

Solidification during fusion-based additive manufacturing (AM) is characterized by high solidification velocities and large thermal gradients, two factors that control the solidification mode of metals and alloys. Using two synchrotron-based, in situ setups, we perform high-speed X-ray diffraction measurements to investigate the impact of the solidification velocities and thermal gradients on the solidification mode of a hot-work tool steel over a wide range of thermal conditions of relevance to AM of metals. The solidification mode of primary δ-ferrite is observed at a cooling rate of 2.12x10 4 K/s, and at a higher cooling rate of 1.5x10 6 K/s, δ-ferrite is sup-pressed, and primary austenite is observed. The experimental thermal conditions are evaluated and linked to a Kurz-Giovanola-Trivedi (KGT) based solidification model. The modelling results show that the predictions from the multicomponent KGT model agree with the experimental observations. This work highlights the role of in situ XRD measurements for a fundamental understanding of the microstructure evolution during AM and for vali-dation of computational thermodynamics and kinetics models, facilitating parameter and alloy development for AM processes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomic Force Microscopy-Based Force Spectroscopy and Multiparametric Imaging of Biomolecular and Cellular Systems

During the last three decades, a series of key technological improvements turned atomic force microscopy (AFM) into a nanoscopic laboratory to directly observe and chemically characterize molecular and cell biological systems under physiological conditions. In this work, we review key technological improvements that have established AFM as an analytical tool to observe and quantify native biological systems from the micro- to the nanoscale. Native biological systems include living tissues, cells, and cellular components such as single or complexed proteins, nucleic acids, lipids, or sugars. We showcase the procedures to customize nanoscopic chemical laboratories by functionalizing AFM tips and outline the advantages and limitations in applying different AFM modes to chemically image, sense, and manipulate biosystems at (sub)nanometer spatial and millisecond temporal resolution. We further discuss theoretical approaches to extract the kinetic and thermodynamic parameters of specific biomolecular interactions detected by AFM for single bonds and extend the discussion to multiple bonds. Finally, we highlight the potential of combining AFM with optical microscopy and spectroscopy to address the full complexity of biological systems and to tackle fundamental challenges in life sciences.

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

Geochemical Modeling of Celestite (SrSO 4 ) Precipitation and Reactive Transport in Shales

Celestite (SrSO 4 ) precipitation is a prevalent example of secondary sulfate mineral scaling issues in hydraulic fracturing systems, particularly in basins where large concentrations of naturally occurring strontium are present. Herein, we present a validated and flexible geochemical model capable of predicting celestite formation under such unconventional environments. Simulations were built using CrunchFlow and guided by experimental data derived from batch reactors. These data allowed the constraint of key kinetic and thermodynamic parameters for celestite precipitation under relevant synthetic hydraulic fracturing fluid conditions. Effects of ionic strength, saturation index, and the presence of additives were considered in the combined experimental and modeling construction. This geochemical model was then expanded into a more complex system where interactions between hydraulic fracturing fluids and shale rocks were allowed to occur subject to diffusive transport. We find that the carbonate content of a given shale and the presence of persulfate breaker in the system strongly impact the location and extent of celestite formation. The results of this study provide a novel multicomponent reactive transport model that may be used to guide future experimental design in the pursuit of celestite and other sulfate mineral scale mitigation under extreme conditions typical of hydraulic fracturing in shale formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Influence of Nd(II) and Nd(III) Ions on the Speciation and Kinetic Dynamics of Radiolytic Transients in Molten LiCl-KCl Eutectic Salt Mixtures

In the pyrochemical reprocessing of used nuclear fuel, understanding the effects of ionizing radiation on the chemical behavior of fission products in molten chloride eutectic salt is crucial. Here, this study investigates the impact of radiation on the chemical and kinetic speciation of Nd ions in molten LiCl-KCl eutectic, which are challenging to separate due to their multivalent oxidation states and comproportionation reactions. Electron pulse radiolysis techniques were used to determine transient absorption spectra, chemical kinetics, and activation parameters for the reaction of Nd ions with the primary products of molten LiCl-KCl eutectic radiolysis: the solvated electron (e S - ) and dichlorine radical anion (Cl 2 •- ). Initially, Nd(III) reacted rapidly with e S - , forming Nd(II) with a second-order rate coefficient of k = (4.54 ± 0.07) × 10 10 M -1 s -1 at 673 K. The resulting Nd(II) ions then reacted slowly with Cl 2 •- , regenerating Nd(III) with an estimated rate of k = (1.72 ± 0.04) × 10 10 M -1 s -1 at 673 K. Additionally, the comproportionation reaction of Nd(III) and Nd(s) to form Nd(II) was monitored chemically and electrochemically. This work suggests that the comproportionation reaction of Nd(III) with Nd(s) to form Nd(II) is unstable in molten LiCl-KCl eutectic.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

Comparisons of bpy and phen Ligand Backbones in Cr-Mediated (Co-)Electrocatalytic CO 2 Reduction

Due to the rise in atmospheric carbon dioxide (CO 2 ) concentrations, there is a need for the development of new strategies to enhance the selectivity and activity of the electrocatalytic conversion of CO 2 to value-added products. The incorporation of redox mediators (RMs) as co-catalysts to enhance the transfer of redox equivalents during catalysis has been gaining more attention in recent years across a variety of small molecule transformations. We have shown that using Cr-centered complexes with sulfone-based RMs leads to an enhancement of CO 2 reduction electrocatalysis under protic conditions via an inner-sphere mechanism. In these co-catalytic systems, an oxygen atom of the reduced RM binds to the Cr center to form a key intermediate stabilized by pancake bonding between the reduced aromatic components of the catalyst ligand backbone and the RM. This interaction facilitates the transfer of an electron and accesses a more kinetically favorable reaction pathway. Here, we show that expanding the aromatic character of the ligand backbone of the catalyst as well as the RM can cause a greater enhancement of co-electrocatalytic activity. These results suggest that further activity improvements can be achieved by focusing on the kinetic and thermodynamic parameters which control association between the catalyst and RM.

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

Temperature and Solvent Effects on H 2 Splitting and Hydricity: Ramifications on CO 2 Hydrogenation by a Rhenium Pincer Catalyst

The catalytic hydrogenation of carbon dioxide holds immense promise for applications in sustainable fuel synthesis and hydrogen storage. Mechanistic studies that connect thermodynamic parameters with the kinetics of catalysis can provide new understanding and guide predictive design of improved catalysts. Reported here are thermochemical and kinetic analyses of a new pincer-ligated rhenium complex ( tBu POCOP)Re(CO) 2 ( tBu POCOP = 2,6-bis(di-tert-butylphosphinito)phenyl) that catalyzes CO 2 hydrogenation to formate with faster rates at lower temperature. Because the catalyst follows the prototypical “outer sphere” hydrogenation mechanism, comprehensive studies of temperature and solvent effects on the H 2 splitting and hydride transfer steps are expected to be relevant to many other catalysts. Strikingly large entropy associated with cleavage of H2 results in a strong temperature dependence on the concentration of [( tBu POCOP)Re(CO) 2 H] – present during catalysis, which is further impacted by changing the solvent from toluene to tetrahydrofuran to acetonitrile. New methods for determining the hydricity of metal hydrides and formate at temperatures other than 298 K were developed, providing insight into how temperature can influence the favorability of hydride transfer during catalysis. These thermochemical insights guided the selection of conditions for CO 2 hydrogenation to formate with high activity (up to 364 h –1 at 1 atm or 3330 h –1 at 20 atm of 1:1 H 2 CO 2 ). In cases where hydride transfer is the highest individual kinetic barrier, entropic contributions to outer sphere H 2 splitting lead to a unique temperature dependence: catalytic activity increases as temperature decreases in tetrahydrofuran (200-fold increase upon cooling from 50 to 0 °C) and toluene (4-fold increase upon cooling from 100 to 50 °C). Furthermore, the ramifications on catalyst structure-function relationships are discussed, including comparisons between “outer sphere” mechanisms and metal–ligand cooperation mechanisms.

14 SOLAR ENERGY↗
DOE OSTI2022

A model of spatio-temporal regulation within biomaterials using DNA reaction–diffusion waveguides

In multi-cellular organisms, cells and tissues coordinate biochemical signal propagation across length scales spanning micrometres to metres. Designing synthetic materials with similar capacities for coordinated signal propagation could allow these systems to adaptively regulate themselves across space and over time. Here, we combine ideas from cell signalling and electronic circuitry to propose a biochemical waveguide that transmits information in the form of a concentration of a DNA species on a directed path. The waveguide could be seamlessly integrated into a soft material because there is virtually no difference between the chemical or physical properties of the waveguide and the material it is embedded within. We propose the design of DNA strand displacement reactions to construct the system and, using reaction–diffusion models, identify kinetic and diffusive parameters that enable super-diffusive transport of DNA species via autocatalysis. Finally, to support experimental waveguide implementation, we propose a sink reaction and spatially inhomogeneous DNA concentrations that could mitigate the spurious amplification of an autocatalyst within the waveguide, allowing for controlled waveguide triggering. Chemical waveguides could facilitate the design of synthetic biomaterials with distributed sensing machinery integrated throughout their structure and enable coordinated self-regulating programmes triggered by changing environmental conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A Parallel Kinetic Model for Surface and Bulk Charge Storage in ε -MnO 2 Pseudocapacitors

Pseudocapacitive materials such as manganese dioxide (MnO 2 ) are attractive for energy storage applications due to their ability to combine the fast kinetics of capacitors with the higher energy density of battery-type systems. However, the electrochemical behavior of MnO 2 remains difficult to interpret mechanistically, in part because existing models often fail to distinguish between surface-based redox processes and bulk intercalation mechanisms. In this work, we develop a physics-based model that represents MnO 2 pseudocapacitance as a linear combination of two independent, parallel electrochemical processes: (i) the surface or near-surface redox storage and (ii) lithium ion intercalation into the bulk material. These two processes are treated with distinct kinetic and thermodynamic parameters and are assumed to proceed independently. The total measured current is assumed to be the sum of these two partial currents. We validate the model using rate-dependent cyclic voltammetry experiments, demonstrating that it captures key trends and provides physically interpretable parameters reflecting the relative contributions of surface and bulk processes. By enabling a clear separation between these mechanisms, the model offers a useful framework for analyzing pseudocapacitive materials and can guide the rational design of high-performance energy storage electrodes.

Energy - Storage↗