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.
Temporal analysis of products statistical analysis package (tapsap) is a package for reading, cleaning, analyzing and making inference of temporal analysis of products (TAP) information. This package is designed with the intent for collaboration within the TAP community for advancing TAP analysis.
The temporal analysis of products (TAP) reactor, a transient kinetic tool, provides users with information as the catalyst state evolves. However, the state of the art for TAP uncertainty quantification only considers the experimental noise present in the outlet flow signal. Additional sources of uncertainty, including, but not limited to, initial surface coverages, the catalyst zone location, the inert void fraction, and the gas pulse intensity and pulse delay, are not included. For this reason, a framework for quantifying all uncertainty sources present in TAP experiments is presented and applied to a carbon monoxide oxidation case study. Herein, two methods for quantifying these sources of uncertainty are introduced. The first utilizes initial state sensitivities to approximate the parameter variances, as well as to provide insights into the structural certainty of the model. The second generates parameter confidence distributions through an ensemble-based sampling algorithm. This initial state covariance matrix can ultimately be merged with the experimental noise covariance matrix, providing a unified description of the parameter uncertainties for a single TAP experiment.
Kinetics-informed neural networks improve fit quality for multi-pulse and noisy temporal analysis of products datasets.
Not Available
Surface science methodologies for understanding thermodynamic aspects of surface processes are at an advanced level. However, instrumentation and approaches for extracting kinetic parameters from elementary steps are far less accessible. In this work, we present an approach combining the use of a fast gas pulsing valve synchronized with data acquisition to enable surface transient kinetics studies using infrared reflection absorption spectroscopy. This methodology applies to the study of reversible processes and borrows concepts and ideas from molecular beam scattering and temporal analysis of products. Here, a temporal resolution of ~67 ms is achieved, and this is illustrated through the study of CO adsorption and desorption on a Pd(111) crystal in the presence and absence of background O 2 . The same approach can be extended to other surface spectroscopies, such as X-ray photoelectron spectroscopy, to obtain spectra with high temporal resolution and signal-to-noise ratio and enable future multimodal surface transient kinetic studies aiming at elucidating reaction mechanisms.
The water vapor is a relevant greenhouse gas in the Earth's climate system, and satellite products become one of the most effective way to characterize and monitor the columnar water vapor (CWV) content at global scale. Recently, a new product (MCD19) was released as part of MODIS (Moderate Resolution Imaging Spectroradiometer) Collection 6 (C6). This operational product from the Multi-Angle Implementation for Atmospheric Correction (MAIAC) algorithm includes a high 1-kilometer resolution CWV retrievals. This study presents the first global validation of MAIAC C6 CWV obtained from MODIS MCD19A2 product. This evaluation was performed using Aerosol Robotic Network (AERONET) observations at 265 sites (2000-2017). Overall, the results show a good agreement between MAIAC/AERONET CWV retrievals, with correlation coefficient higher than 0.95 and RMS (Root Mean Square) error lower than 0.250 centimeters. The binned error analysis revealed an underestimation (approximately 10 percent) of Aqua CWV retrievals with negative bias for CWV higher than 3.0 centimeters. In contrast, Terra CWV retrievals show a slope of regression close to unity and a low mean bias of 0.075 centimeters. While the accuracy is relatively similar between 1.0 and 5.0 centimeters for both sensor products, Terra dataset is more reliable for applications in humid tropical areas (less than 5.0 centimeters). The expected error was defined as plus or minus 15 percent, with less than 68 percent of retrievals falling within this envelope. However, the accuracy is regionally dependent, and lower error should be expected in some regions, such as South America and Oceania. Since MODIS instruments have exceeded their design lifetime, time series analysis was also presented for both sensor products. The temporal analysis revealed a systematic offset of global average between Terra and Aqua CWV records. We also found an upward trend (approximately 0.2 centimeters per decade) in Terra CWV retrievals, while Aqua CWV retrievals remain stable over time. The sensor degradation influences the ability to detect climate signals, and this study indicates the need for revisiting calibration of the MODIS bands 17-19, mainly for Terra instrument, to assure the quality of the MODIS water vapor product. Finally, this study presents a comprehensive validation analysis of MAIAC CWV over land, raising the understanding of its overall quality.
The Temporal Analysis of Products (TAP) pulse response technique provides valuable insights into catalytic function and reaction kinetics. However, complex fragmentation patterns in the TAP mass spectrometry signals can complicate precise quantification, particularly when analyzing transient gas flux data typical of TAP experiments. This work demonstrates a standard defragmentation method that deconvolves transient TAP signals while maintaining the temporal resolution of the experiment. First, the integrals of calibration gas fluxes are used to determine the fingerprint fragmentation pattern and construct a fragmentation matrix. This matrix is then used to defragment experimental flux data at each recorded time point via a non-negative least squares regression. The effectiveness of this method is demonstrated using virtual data and control experiments with a TAP reactor system. The defragmentation is then applied to the more complex propane dehydrogenation reaction on a chromia/alumina catalyst, which can contain up to ten significant gas species in the reactor outlet. Initial propane pulsing reveals an induction period during which propane is fully oxidized to CO2, followed by partial reduction to CO. Afterwards, there is a transition in chemistries towards coking and propylene production. Our example illustrates a practical method for the accurate determination of the time-dependent reactant/product concentrations and rates for a thorough analysis of the propane dehydrogenation kinetics. This approach can be broadly applied to any transient mass spectrometry experiment for a better understanding of catalyst-reaction dynamics.
This perspective reviews the current understanding of the Oxidative Coupling of Methane (OCM) reaction over the supported Mn-Na 2 WO 4 /SiO 2 catalyst, with a focus on recent insights gained from state-of-the-art in-situ and operando spectroscopic characterization and chemical probe experiments under controlled environments. The supported Mn-Na 2 WO 4 /SiO 2 catalyst exhibits dynamic structural changes during the OCM reaction, involving multiple reactive lattice and adsorbed oxygen species, each associated with different oxide phases. These oxygen species play distinct roles in various steps of the OCM mechanism. The catalytic active sites for activation of CH 4 are associated with isolated surface Na-WO x sites on the SiO 2 support and the role of surface MnO x sites on SiO 2 is to oxidatively dehydrogenate C 2 H 6 to C 2 H 4 . Furthermore, this paper provides a detailed discussion of these roles and also introduces new experimental data from Temporal Analysis of Products (TAP) studies to clarify the ongoing debate in the literature regarding the contributions of lattice versus adsorbed oxygen species in OCM reaction product formation. Additionally, recommendations are offered for optimizing the performance of supported Mn-Na 2 WO 4 /SiO 2 catalysts to enhance CH 4 activation and C 2 product selectivity.
Selective oxidation of propene to acrolein over industrial multicomponent bismuth molybdate (BMO) catalysts significantly depends on reaction conditions that include operating parameters and catalyst state. Here, this work investigates selective propene oxidation in the intrinsic kinetic paradigm of Temporal Analysis of Products (TAP) reactor by systematically varying catalyst redox state, temperature, and oxygen-to-propene feed ratio. A 93 % propene conversion with an acrolein yield of 80 % is achieved at elevated temperatures (450 °C) on oxidized catalysts under oxygen-rich conditions (O 2 :C 3 H 6 = 10). However, these conditions diminish acrolein-to-CO 2 selectivity due to enhanced total oxidation to CO 2 . In contrast, a reduced catalyst state, moderate temperature (350 °C), and lower oxygen feed (O 2 :C 3 H 6 = 1) nearly doubles the acrolein to-CO 2 selectivity, albeit at a lower acrolein yield (33 %). Transient kinetic studies together with a kinetic model that simplify the major product formation pathways in lumped non-elementary forms reveal that the availability of surface oxygen species plays a pivotal role in governing reaction pathways. Additionally, density functional theory (DFT) calculations on pure BMO catalysts inform the role of surface redox states on propene and oxygen activation barriers. Lattice oxygen at acrolein-selective sites drives both acrolein and CO 2 formation, while adsorbed oxygen at activation sites favors unselective CO 2 generation. This work establishes a critical relationship between transient product selectivity and surface oxygen availability, which is strongly influenced by catalyst redox state, feed ratio, and reaction temperature. These insights underscore the importance of dynamic reactor operation strategies and offer a foundation for designing next-generation propene oxidation processes with tunable acrolein selectivity.
Here, this paper introduces the use of statistical distributions based on transport differential equations for clear distinction of transport modes within transient kinetic experiments. More specifically, novel techniques are developed for the transient data obtained through the Temporal Analysis of Products (TAP) reactor and are applicable to experiments where pulse response into a gas flow is used. The methodology allows distinguishing between two domains of diffusion transport in heterogeneous catalytic systems, i.e., Knudsen and non-Knudsen diffusion, using statistical fingerprints, and finding the transition domain. Two distribution parameters were obtained that directly result in coefficients that correspond to the concentration and the rate of transport. Using a linear relationship between the rate and concentration coefficients, Knudsen diffusion is revealed when the rate of transport is constant and non-Knudsen diffusion is confirmed when the rate of transport coefficient is a function of the concentration coefficient. As a result, accurate transport information can be extracted from experimental data even in the presence of comprising instrument drift or noise particularly when analyzing higher pressure pulse responses with complex transport. This enables more direct investigation of experiments influenced by gas-phase reactions.
Oxidative coupling of methane (OCM) is a promising industrial process to upgrade natural gas to high value chemicals. In this study, Temporal Analysis of Products (TAP) and steady-state experiments were conducted to distinguish how the composition of surface and gas phase oxygen influence mechanistic details of the selective conversion of CH 4 to C 2 H 4 over the Mn-Na 2 WO 4 /SiO 2 catalyst. The results from TAP studies indicate that methane activation on this catalyst proceeds predominantly via a short-lived, transient surface oxygen species and there is a competition for this species to form either CO or methyl radicals on the surface. This active species has a total lifetime of 3 s and is identified to have a dioxygen (e.g. O 2 2- or O 2 -) form. We show that the concentration of the transient surface oxygen species significantly impacts the OCM performance. Oxygen attributed to the catalyst lattice (in a singular form e.g., O - ), is found to activate methane to a lesser degree, but exclusively forms CO 2 . Evidence for these surface pathways for methyl radical, CO and CO 2 formation identified by TAP are also validated through steady-state experiments. Finally, by distinguishing different catalyst oxygen species and their role in selective/nonselective pathways, important screening criteria have been identified for the advancement of superior catalyst formulations.
Oxidative coupling of methane (OCM) is an attractive direct route for upgrading methane to valuable chemicals. In this study, Temporal Analysis of Products (TAP) and steady state experiments are conducted to understand the role of individual oxide phases and their combinations in supported Mn-Na 2 WO 4 /SiO 2 catalysts for OCM. The results from TAP transient kinetic studies indicate that Mn plays an important role in promoting gas phase oxygen activation, while NaO x /SiO 2 and WO x /SiO 2 are relatively inert towards gas phase oxygen and methane activation. However, the supported catalyst combining Na and W in the form Na 2 WO 4 show enhanced gas phase oxygen activation exhibiting a much lower oxygen activation energy (148 kJ/mol) and enhanced activity toward methane activation as compared to the individual supported oxide catalysts. Addition of Mn to Na 2 WO 4 /SiO 2 further decreases the oxygen activation energy by 40 kJ/mol. Moreover, methane activation is also enhanced with CH 3 as the main intermediate but, with increasing Mn content, more CH 2 intermediates are observed. Different forms of oxygen (both dioxygen and atomic) are detected on the catalyst surface using isotopic pump/probe pulsing and their distribution is found to depend on the catalyst composition. An optimal Mn content in the Na 2 WO 4 /SiO 2 catalyst system is needed to enhance the amount of dioxide surface species (e.g., superoxide 16 O 2 - or peroxide 16 O 2 2- ) associated with the Na 2 WO 4 leading to high C 2 selectivity for OCM. When the Mn content is too high, the larger MnOx domains are shown to contribute to the formation of higher concertation of monoxide surface species that lead to nonselective OCM pathways. This insight from transient kinetic characterization using TAP combined with conventional steady state studies, provides a deeper understanding of the role of individual oxide phases and their combination on supported catalysts toward the formation of intermediate surface species and their impact on the OCM reaction mechanism. This knowledge is critical toward designing superior catalyst formulations for OCM.
We present a case study of catalytic propane dehydrogenation with our open-source multiphysics solver, Mesoflow. The solver was developed to simulate reactive flow coupled to heterogeneous catalytic reactions and deactivation in the context of complex, mesoscale geometry. The method leverages cartesian block-structured adaptive mesh refinement to capture realistic catalyst microstructural features acquired directly from X-ray computed tomography data. A kinetic model for propane dehydrogenation and catalyst deactivation was developed based on temporal analysis of products (TAP) reactor experiments. The TAP reactor experiments allow for precise characterization of intrinsic kinetic reaction steps which are implemented into Mesoflow simulations to model the spatial and temporal evolution of reactants, products, and catalyst active sites. The short-term and long-term deactivation behavior is studied by using XCT data collected from fresh and aged catalyst pellets, which exhibit different microstructural features. This study employs time-splitting algorithms to connect disparate reaction and flow timescales, enabling the simulations to achieve realistic deactivation timescales on the order of minutes while the flow time-scales for small particles (100 microns) are several milliseconds. We also introduce a flexible automated python script that writes the necessary files to construct a Mesoflow simulation from user-created chemical mechanisms. We will also introduce a few new features that are added to Mesoflow such as higher order schemes, implicit chemistry integrators and the ability to run on AMD and NVIDIA graphics-processing-units.
Syngas, a combination of carbon monoxide (CO) and hydrogen (H2), is a precursor to many chemicals and fuels, contributing to the billion-dollar global hydrocarbon industry. Chemical looping reforming (CLR) of the greenhouse gases methane (CH4) and carbon dioxide (CO2) allows for energy-efficient production of syngas. 5wt.% nickel (Ni) on ceria-zirconia (5wt%Ni@Ce0.6Zr0.4O2) mixed metal oxide catalyst was investigated here to explore pathways for enhanced syngas production on sustainable earth-abundant transition metal supported catalysts. The role of reduction-oxidation (redox) state of the catalyst, and carbon formation on the catalyst surface during chemical looping is explored to drive superior reaction kinetics, conversions, selectivity, and syngas ratios (H2/CO). The bulk and surface structure of the catalyst, along with carbon deposition features were characterized by electron microscopy, X-Ray diffraction, and ex situ Raman spectroscopy. The dynamic evolution of catalysts under CLR reaction conditions and the intrinsic reaction mechanisms were probed with in situ Raman spectroscopy, in situ Fourier transform spectroscopy, dynamic oxygen storage capacity (DOSC), and the Temporal Analysis of Products (TAP) reactor studies. Intrinsic kinetics of syngas production via CLR was correlated to the redox state of catalyst and the participation of nickel-catalyzed multiwalled carbon nanotube (CNT) growth, allowing enhanced CLR reaction performance.
This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.
Abstract The complex structure of the catalytic active phase, and surface‐gas reaction networks have hindered understanding of the oxidative coupling of methane (OCM) reaction mechanism by supported Na 2 WO 4 /SiO 2 catalysts. The present study demonstrates, with the aid of in situ Raman spectroscopy and chemical probe (H 2 ‐TPR, TAP and steady‐state kinetics) experiments, that the long speculated crystalline Na 2 WO 4 active phase is unstable and melts under OCM reaction conditions, partially transforming to thermally stable surface Na‐WO x sites. Kinetic analysis via temporal analysis of products (TAP) and steady‐state OCM reaction studies demonstrate that ( i ) surface Na‐WO x sites are responsible for selectively activating CH 4 to C 2 H x and over‐oxidizing CH y to CO and ( ii ) molten Na 2 WO 4 phase is mainly responsible for over‐oxidation of CH 4 to CO 2 and also assists in oxidative dehydrogenation of C 2 H 6 to C 2 H 4 . These new insights reveal the nature of catalytic active sites and resolve the OCM reaction mechanism over supported Na 2 WO 4 /SiO 2 catalysts.
Abstract The complex structure of the catalytic active phase, and surface‐gas reaction networks have hindered understanding of the oxidative coupling of methane (OCM) reaction mechanism by supported Na 2 WO 4 /SiO 2 catalysts. The present study demonstrates, with the aid of in situ Raman spectroscopy and chemical probe (H 2 ‐TPR, TAP and steady‐state kinetics) experiments, that the long speculated crystalline Na 2 WO 4 active phase is unstable and melts under OCM reaction conditions, partially transforming to thermally stable surface Na‐WO x sites. Kinetic analysis via temporal analysis of products (TAP) and steady‐state OCM reaction studies demonstrate that ( i ) surface Na‐WO x sites are responsible for selectively activating CH 4 to C 2 H x and over‐oxidizing CH y to CO and ( ii ) molten Na 2 WO 4 phase is mainly responsible for over‐oxidation of CH 4 to CO 2 and also assists in oxidative dehydrogenation of C 2 H 6 to C 2 H 4 . These new insights reveal the nature of catalytic active sites and resolve the OCM reaction mechanism over supported Na 2 WO 4 /SiO 2 catalysts.
The temporal analysis of products (TAP) reactor provides a vast amount of transient kinetic information that may be used to describe a variety of chemical features including residence time distributions, kinetic coefficients, number of active sites, reaction mechanism, etc. However, as with any measurement device, the TAP reactor signal is convoluted with noise and drift is common. In order to reduce the uncertainty of the kinetic measurement and any derived parameters or mechanisms, proper preprocessing must be performed prior to any advanced type of analysis. This preprocessing includes baseline correction, i.e., a shift in the voltage response, and calibration, i.e., a scaling of the flux response based on prior experiments. The traditional methodology of preprocessing requires significant user discretion and reliance on separate calibration experiments that may drift over time. Herein we use machine learning techniques combined with physical constraints to understand the noise and drift that is being generated within and between experiments for enhancement of the chemical kinetic signal. As such, the proposed methodology demonstrates clear benefits over the traditional preprocessing approach by eliminating the need for separate calibration experiments or heuristic input from the user.