Engineering Papers⌕ Search

Engineering topics

Bhan, Aditya

Publications and source records attributed to Bhan, Aditya.

18 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Reaction Pathways and Energy Consumption in NH 3 Decomposition for H 2 Production by Low Temperature, Atmospheric Pressure Plasma

Pathways for NH 3 decomposition to N 2 and N 2 H 4 by atmospheric pressure nonthermal plasma are analyzed using a combination of molecular beam mass spectrometry measurements and zero-dimensional kinetic modeling. Experimental measurements show that NH 3 conversion and selectivity towards N 2 formation scale monotonically with the specific energy input into the plasma with ~ 100% selectivity to N 2 formation achieved at specific energy inputs above 0.12 J cm −3 (3.1 eV (molecule NH 3 ) −1 ). The kinetic model recovers these trends, although it underpredicts N 2 selectivity at low specific energy input. These discrepancies can be explained by the underestimation of reaction rate coefficients for reactions that consume N 2 H x species in collisions with H radicals and/or radial nonuniformities in power deposition, gas temperature, and species concentrations that are not represented by the plug flow approximation used in the model. The kinetic model shows that N 2 formation proceeds through N 2 H x decomposition pathways rather than NH x decomposition pathways in low temperature, atmospheric pressure plasma. Higher selectivity toward N 2 production can be achieved by operating at higher NH 3 conversion and with a higher gas temperature. Furthermore, the high energy cost of NH 3 decomposition by atmospheric pressure nonthermal plasma found in this work (25–50 eV (molecule NH 3 converted) −1 ; 17–33 eV (molecule H 2 formed) −1 ) is a result of the energy requirement for electron-impact dissociation of NH 3 and the significant re-formation of NH 3 by three-body recombination reactions between NH 2 and H.

Nonthermal plasma↗
DOE OSTI2024

NO formation by N 2 /O 2 plasma catalysis: The impact of surface reactions, gas-phase reactions, and mass transport

Pathways and timescales relevant to facilitate plasma-assisted N 2 -O 2 reactions are assessed by measuring the consumption of plasma-derived N and the formation of NO in the gas phase and over Ag catalytic surfaces. These measurements are enabled by a setup that enables N 2 activation in an atmospheric pressure RF plasma jet, enables O 2 addition in the plasma afterglow, facilitates reactions over an Ag wire catalyst, and allows species density quantification by molecular beam mass spectrometry. Gas-phase reactions consume N but do not form NO with high selectivity. The presence of the non-porous Ag wire catalyst increases the rate of N conversion to NO, though mass transfer processes, not surface reactions, dictate the rate of N consumption. When O 2 concentrations and the ratio of the surface area of the catalyst to the void volume of the reactor are high (3–5 mol% O 2 , 10900 m –1 ), N conversion to NO reaches 100 % selectivity. When both N 2 and O 2 are fed through the plasma jet, gas-phase NO production increases 10×, although plasma and gas-phase processes do not exclusively produce NO. Above a threshold NO density, N cannot diffuse to the catalyst surface faster than it is consumed in the gas phase by reactions with NO. Furthermore, the use of heterogeneous catalysts to enhance plasma-driven N x O y formation and control N x O y product selectivity is limited to cases where diffusive transport of N from the gas phase to the catalyst surface is faster than consumption of N from gas-phase reactions with NO.

Engineering↗
DOE OSTI2023

Availability and reactivity of N 2 ( v ) for NH 3 synthesis by plasma catalysis

Production of vibrationally excited N 2 (N 2 (v)) in atmospheric pressure nonthermal plasma and loss of N 2 (v) by gas-phase reactions and reactions on catalytic surfaces are analyzed to examine the role of N 2 (v) in NH 3 formation by plasma catalysis. Vibrational state-to-state kinetic models complemented with molecular beam mass spectrometry (MBMS) measurements demonstrate that N 2 (v> 0) is produced with densities 100× greater than the density of N radicals by a radiofrequency atmospheric pressure plasma jet. The experimentally measured loss of N 2 (v) corresponds with a state-to-state kinetic model that describes loss of N 2 (v) by surface-mediated vibrational relaxation without consideration of reactions that convert N 2 (v) to NH 3 over the catalyst surface. Rate constants for vibrational relaxation of N 2 (v) on catalyst surfaces exceed upper bounds on proposed rate constants for NH 3 formation reactions from N 2 (v) over Fe when v < 9, Ni when v < 18, and Ag when v < 39, which indicates that only higher vibrational levels can possibly contribute to catalytic NH 3 formation faster than they undergo vibrational relaxation on the surface. Densities of N 2 (v> 8), vibrational levels that can possibly react over Fe to form NH 3 faster than they undergo vibrational relaxation, are less than or similar to N densities at the inlet of the catalyst bed and measured NH 3 formation for the investigated conditions in this work, while densities of N 2 (v> 17) and N 2 (v> 38) are orders of magnitude below the N density at the inlet of the catalyst bed and the measured NH 3 formation. The loss of N 2 (v) by vibrational relaxation on the surface limits the ability of N 2 (v) to contribute to catalytic NH 3 formation and explains why N 2 (v) does not produce NH 3 in quantities that are comparable to NH 3 formation from N even though N 2 (v > 0) is more abundantly produced by the plasma.

state-to-state kinetic modeling↗
DOE OSTI2023

Rates and reversibility of CO 2 hydrogenation on Cu-based catalysts

Kinetics of reaction pathways involved in the conversion of CO 2 to methanol and CO on Cu/ZnO/Al 2 O 3 are resolved using in situ chemical titration, steady-state kinetic measurements, and mathematical formalisms for reversibility to probe salient species governing methanol selectivity and yield during CO 2 hydrogenation. Across a range of H 2 :CO 2 = 1:1 to H 2 :CO 2 = 80.5:1, active site density determined from in situ chlorine uptake remained invariant; hence, observed trends in rates can be interpreted as only arising from reaction kinetics and not from changing active site density. Kinetic and thermodynamic contributions to rates are decoupled to evaluate forward and reverse rates of methanol synthesis and reverse water-gas shift (RWGS) reactions. These kinetic analyses show that the forward rates of methanol synthesis exhibit persistent first order dependence on hydrogen pressure and are inhibited by water more significantly than the forward rates of RWGS. In contrast, the reverse rates of methanol synthesis and RWGS are both inhibited by H 2 . Consequently, without any modifications to the Cu/ZnO/Al 2 O 3 catalyst formulation, methanol selectivity can be increased to > 80% by increasing inlet H 2 partial pressure and methanol yield can be enhanced by ~20% by adding water adsorbents even under conditions far from equilibrium. Here, the kinetic treatments presented herein demonstrate a dearth of H* species during catalysis, provide thermodynamic constraints precluding sequential RWGS and CO hydrogenation as the pathway for methanol synthesis, reveal P H2 and P H2O as salient in determining methanol selectivity and yield by impacting both the forward and reverse rates of CO 2 hydrogenation on Cu/ZnO/Al 2 O 3 , and explicate the fundamentals of novel sorption-enhanced methanol synthesis, which not only alleviates equilibrium constraints but also alters the intrinsic rate at which the system approaches equilibrium.

42 ENGINEERING↗
DOE OSTI2023

Validation of the Cossee–Arlman mechanism for propylene oligomerization on Ni/UiO-66

Steady state rate expressions can be derived to distinguish the Cossee–Arlman and metallacycle mechanisms postulated for propylene oligomerization on nickel-based catalysts based on product selectivities, where product selectivities for the former are a function of olefin pressure because sequential coordination and insertion steps lead to independent mechanistic pathways for different hexene isomers. In contrast, the metallacycle mechanism presents pressure-independent product selectivities due to successive coordination prior to the kinetically relevant steps in each mechanism. In this work, steady state propylene oligomerization rates and selectivities were measured in the absence of an activator on nickel functionalized UiO-66 metal organic framework (MOF), Ni/UiO-66, to validate the Cossee–Arlman mechanism for light olefin oligomerization. In situ NO titrations reveal that ~5% of nickel sites were active during the reaction, and thus, not all nickel sites are relevant for catalysis. Propylene dimerization was first order in propylene pressure from 5 to 500 kPa with an apparent activation energy of ~20 kJ mol -1 from 453 to 493 K. Calculated apparent activation energies with density functional theory (DFT) calculations on cluster models of Ni/UiO-66 are in agreement with experiment to corroborate the Cossee–Arlman mechanism. Selectivities of hexene products and the ratio of hexene product selectivities on Ni/UiO-66 are in accordance with selectivity expressions derived from the Cossee–Arlman mechanism. In conclusion, analysis of product selectivities can be used more extensively to demarcate the Cossee–Arlman and metallacycle mechanisms for olefin oligomerization on metal-based catalysts.

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

Comparing the reaction profiles of single iron catalytic sites in enzymes and in reticular frameworks for methane-to-methanol oxidation

The design of synthetic inorganic catalysts mimicking the first coordination spheres of enzymatic cofactors often results in lower yields and selectivity than their biological counterparts. In this study, we exploit Kohn-Sham density functional methods to compare the reaction profiles of four single iron-based catalysts for the direct oxidation of methane to methanol: two biomimetic models based on two enzymes (cytochrome P450 and taurine dioxygenase [TauD]) and two synthetic reticular frameworks (iron-BEA zeolite and tri-iron oxo-center-based metal-organic framework). Both the biomimetic and inorganic catalysts show almost zero selectivity toward methanol for methane conversions >1% at ambient temperature. Furthermore, this study highlights that iron’s first coordination shell can influence selectivity toward methanol but to a limited extent. In the absence of methanol protection strategies, high selectivity can be reached only by mimicking the reaction microenvironment of enzymes beyond the first coordination shell of iron.

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

Concepts Relevant for the Kinetic Analysis of Reversible Reaction Systems

The net rate of a reversible chemical reaction is the difference between unidirectional rates of traversal along forward and reverse reaction paths. In a multistep reaction sequence, the forward and reverse trajectories, in general, are not the microscopic reverse of one another; rather, each unidirectional route is comprised of distinct rate-controlling steps, intermediates, and transition states. Consequently, traditional descriptors of rate (e.g., reaction orders) do not reflect intrinsic kinetic information but instead conflate unidirectional contributions determined by (i) the microscopic occurrence of forward/reverse reactions (i.e., unidirectional kinetics) and (ii) the reversibility of reaction (i.e., nonequilibrium thermodynamics). This review aims to provide a comprehensive resource of analytical and conceptual tools which deconvolute the contributions of reaction kinetics and thermodynamics to disambiguate unidirectional reaction trajectories and precisely identify rate- and reversibility-controlling molecular species and steps in reversible reaction systems. The extrication of mechanistic and kinetic information from bidirectional reactions is accomplished through equation-based formalisms (e.g., De Donder relations) grounded in principles of thermodynamics and interpreted in the context of theories of chemical kinetics developed in the past 25 years. Here, the aggregate of mathematical formalisms detailed herein is general to thermochemical and electrochemical reactions and encapsulates a diverse body of scientific literature encompassing chemical physics, thermodynamics, chemical kinetics, catalysis, and kinetic modeling.

42 ENGINEERING↗
DOE OSTI2023

Species, Pathways, and Timescales for NH 3 Formation by Low-Temperature Atmospheric Pressure Plasma Catalysis

Species, pathways, and timescales for NH3 production by plasma catalysis over transition-metal wools are determined by measuring plasma-derived species densities [N, H, and N 2 (v)], quantitatively correlating consumption of these species with NH 3 formation, and measuring consumption of plasma-derived species at different residence times. These findings are enabled by a capillary flow through Ar/N 2 /H 2 plasma jet reactor setup that allows for the measurement of gas-phase species densities by molecular beam mass spectrometry. Surface-mediated reactions involving N radicals are responsible for NH 3 formation over Fe, Ni, and Ag surfaces. N reacts to form NH 3 with ~100% selectivity over Ni and Ag when H/N > 3 and % H 2 ≥ 0.5. The selectivity to ammonia drops as H and H 2 densities decrease for each catalyst. A comparison between amounts of NH 3 formed and N consumed with and without catalysts present shows that surface reactions enable higher and more selective conversion of N to NH 3 than gas-phase reactions alone. The conversion of N to NH 3 is negligible in the absence of H, demonstrating that H is required to produce NH 3 at these operating conditions. The consumption of N occurs on the same timescale as NH 3 formation, further confirming that reactions involving N contribute to NH 3 formation. Though vibrationally excited N 2 [N 2 (v)] is produced in quantities exceeding N by 100-fold, consumption of N 2 (v) on the catalytic surface does not contribute to NH 3 formation. Furthermore, these findings show that for low-temperature atmospheric pressure plasma catalysis, surface- mediated reactions among radical N and H species drive NH 3 formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Structure and Site Evolution of Framework Ni Species in MIL-127 MOFs for Propylene Oligomerization Catalysis

A mixed-valence oxotrimer metal–organic framework (MOF), Ni-MIL-127, with a fully coordinated nickel atom and two iron atoms in the inorganic node, generates a missing linker defect upon thermal treatment in helium (>473 K) to engender an open coordination site on nickel which catalyzes propylene oligomerization devoid of any cocatalysts or initiators. This catalyst is stable for ~20 h on stream at 500 kPa and 473 K, unprecedented for this chemistry. The number of missing linkers on synthesized and activated Ni-MIL-127 MOFs is quantified using temperature-programmed oxidation, 1 H nuclear magnetic resonance spectroscopy, and X-ray absorption spectroscopy to be ~0.7 missing linkers per nickel; thus, a majority of Ni species in the MOF framework catalyze propylene oligomerization. In situ NO titrations under reaction conditions enumerate ~62% of the nickel atoms as catalytically relevant to validate the defect density upon thermal treatment. Propylene oligomerization rates on Ni-MIL-127 measured at steady state have activation energies of 55–67 kJ mol –1 from 448 to 493 K and are first-order in propylene pressures from 5 to 550 kPa. Density functional theory calculations on cluster models of Ni-MIL-127 are employed to validate the plausibility of the missing linker defect and the Cossee–Arlman mechanism for propylene oligomerization through comparisons between apparent activation energies from steady-state kinetics and computation. Here this study illustrates how MOF precatalysts engender defective Ni species which exhibit reactivity and stability characteristics that are distinct and can be engineered to improve catalytic activity for olefin oligomerization.

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

Benchmarking Cu/BEA and HBEA catalysts for high-octane gasoline synthesis

We distinguish rates at which carbon deposition occurs during initiation, rates at which catalytic centers are lost during deactivation, and paraffin-to-olefin ratio during propagation as benchmarks that distinguish 5 wt% Cu/H-BEA and H-BEA (Si/Al = 13.5) catalysts during dimethyl ether (DME) homologation in the presence of hydrogen. Studies that systematically vary initial DME contact time (210, 94, and 45 mol H+, initial s (mol C ) -1 ), DME pressure (4 and 22 kPa), and H 2 pressure (1, 24, and 48 kPa) reveal that Cu enables lower carbon deposition rates (on a per proton basis) in the induction period, increases the effluent paraffin-to-olefin ratio during propagation, and decreases the instantaneous site-loss yields by a factor of ~ 1.5-2x (moles of active sites lost per mole of DME) during termination sequences thus affecting the degree of product saturation and catalyst stability during DME homologation. Furthermore, these results provide mechanistic insights revealing the critical role of Cu in facilitating DME homologation to high value, high-octane gasoline-range hydrocarbons with higher cumulative turnovers than proton form H-BEA.

09 BIOMASS FUELS↗
DOE OSTI2022

Acid Catalysis over Low-Silica Faujasite Zeolites

We report low-silica faujasite (FAU) zeolites (with Si/Al ratio of ca. 1.2–1.8) sustain framework integrity and porosity upon moderate ion exchange (0.01 M NH 4 NO 3 solution for 1 h at ambient temperature), which introduces two kinds of protons, distinctive in reactivity and coordination to the zeolite framework, within supercages (H SUP ). Moderate ion exchange limited within supercages transpires while maintaining full occupancy of Na + cations within associated sodalite cages; this in turn helps stabilize the framework of low-silica H-FAU zeolites. Protons located on site II (H 3630 ) and site III (H 3650 ) within supercages on low-silica FAU zeolites can be classified and enumerated by virtue of infrared spectroscopy, providing an opportunity to compare reactivities of these distinct protons for monomolecular protolytic reactions of propane. Protons on site II exhibit prominently higher reactivity for monomolecular propane dehydrogenation and cracking than protons on site III. Low-silica proton-form FAU zeolites (zeolite X) upon moderate ion exchange possess protons on site III that are unavailable on high-silica FAU zeolites (zeolite Y) and limit ion exchange within supercages, providing unprecedented high-temperature structural and chemical stability (>773 K) and enabling their application as solid-acid catalysts.

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

Rates and Reversibilities in Interconnected Reaction Networks

Mathematical relations prescribing unidirectional forward and reverse rates originally derived based on single-path reaction sequences do not apply to interconnected reaction networks. The presence of branches in reaction networks, leading to alternative stable products, decreases unidirectional rates in reference to those calculated by single-path functional forms, as shown by simulated isotopic exchange rates, but impacts the unidirectional forward and reverse rates equally such that the functional form of effective reversibility remains unchanged. Regardless of stoichiometric numbers and network connectivity, the application of the pseudo-steady-state hypothesis on reactive intermediates in conjunction with consideration of unidirectional rates toward the product of interest and all alternative stable products results in mathematical expressions that accurately reflect simulated isotopic exchange rates. Further analyses based on kinetic resistance, a property akin to electrical resistance, illustrate the manifestation of nodal resistances in addition to the single-path kinetic resistance for interconnected reaction networks. The generalized formalism for assessing rates and reversibilities in interconnected networks derived herein enables us to demonstrate that unidirectional rates cannot be assessed solely from effective reversibilities and net rates of generation of stable species in such networks. Furthermore, isotopic exchange rates, under the condition that each elementary step in the overall reaction sequence forms a unique reactive intermediate that is consumed solely by the subsequent step, can be utilized to determine unidirectional rates and can serve to validate postulated reaction pathways in highly interconnected reaction networks (e.g., CO x hydrogenation).

42 ENGINEERING↗
DOE OSTI2021

Best practices in catalysis: A perspective

Catalysis, from its roots in petrochemical refining and conversion, has emerged as a transdisciplinary field that now encompasses synthesis of materials and molecules that enable applications in energy conversion and storage, environmental remediation, medicine, plastics, and fertilizer production, among numerous others. A handful of disciplines can claim relevance and success over such an extended period of time and continue to claim a preeminent role in defining the state-of-the-art in science and technology. Syncretic and rapid advancements in formulation and spectroscopic characterization of materials and molecules useful as catalysts, high-level density functional and molecular orbital theory calculations, and a strong foundation in concepts of physical chemistry, thermodynamics, and chemical kinetics offer new and abundant opportunities at the present day for addressing the grand challenge of controlling chemical transformations using catalysis. Here, we examine what we have learned of concepts that underpin heterogeneous catalysis but more importantly, how we learned to archive our knowledge in context of a set of best practices and standards that have emerged in the course of our learnings—ones we seek to highlight herein. Our perspective emphasizes concepts in synthesis, characterization, kinetics, and theory, because these four elements combined enable description of molecular acts that happen on surfaces and how fast they occur. In revisiting best practices in heterogeneous catalysis in these sub-fields and in authorship and peer review we aspire to augment clarity, reproducibility, and rigor in the science and practice of catalysis.

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

Enhanced Reactivity of Accessible Protons in Sodalite Cages of Faujasite Zeolite

Abstract Faujasite (FAU) zeolites (with Si/Al ratio of ca. 1.7) undergo mild dealumination at moderate ion exchange conditions (0.01 to 0.6 M of NH 4 NO 3 solutions) resulting in protons circumscribed by sodalite cages becoming accessible for reaction without conspicuous changes to bulk crystallinity. The ratio of protons in sodalite cages (H SOD ) to supercages (H SUP ) can be systematically manipulated from 0 to ca. 1 by adjusting ammonium concentrations used in ion exchange. The fraction of accessible protons in the sodalite cages is assessed by virtue of infrared spectra for H‐D exchange of deuterated propane based on the band area ratio of OD 2620 /OD 2680 (OD SOD /OD SUP ). Protons in sodalite cages (H SOD ) show higher rate constants of propane dehydrogenation ( k D ) and cracking ( k C ) than protons in supercages (H SUP ) plausibly due to confinement effects being more prominent in smaller voids. Rate constants of dehydrogenation and cracking including k D / k C ratios are also augmented as the fraction of accessible protons in the sodalite cages is enhanced. These effects of accessibility and reactivity of protons in sodalite cages hitherto inconspicuous are revealed herein via methods that systematically increase accessibility of cations located in sodalite cages.

Li, Xinyu↗
DOE OSTI2021

Enhanced Reactivity of Accessible Protons in Sodalite Cages of Faujasite Zeolite

We report Faujasite (FAU) zeolites (with Si/Al ratio of ca. 1.7) undergo mild dealumination at moderate ion exchange conditions (0.01 to 0.6 M of NH 4 NO 3 solutions) resulting in protons circumscribed by sodalite cages becoming accessible for reaction without conspicuous changes to bulk crystallinity. The ratio of protons in sodalite cages (H SOD ) to supercages (H SUP ) can be systematically manipulated from 0 to ca. 1 by adjusting ammonium concentrations used in ion exchange. The fraction of accessible protons in the sodalite cages is assessed by virtue of infrared spectra for H-D exchange of deuterated propane based on the band area ratio of OD 2620 /OD 2680 (OD SOD /OD SUP ). Protons in sodalite cages (H- SOD ) show higher rate constants of propane dehydrogenation (k D ) and cracking (k C ) than protons in supercages (H SUP ) plausibly due to confinement effects being more prominent in smaller voids. Rate constants of dehydrogenation and cracking including k D /k C ratios are also augmented as the fraction of accessible protons in the sodalite cages is enhanced. These effects of accessibility and reactivity of protons in sodalite cages hitherto inconspicuous are revealed herein via methods that systematically increase accessibility of cations located in sodalite cages.

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

Catalytic site ensembles: A context to reexamine the Langmuir-Hinshelwood kinetic description

The Langmuir-Hinshelwood formalism describes catalytic reactions of Langmuirian surface species under the assumption that all adsorbates are randomly-distributed – enabling adjacency of surface-bound intermediates to be determined solely by coverages of single-site occupants. We demonstrate herein that this approximation is inappropriate even for simple catalytic reactions (e.g. A + A → A 2 ) and manifestly neglects islanding of slowly-consumed species and partitioning of highly-reactive species inherently engendered by ≥ two-site elementary steps (e.g. A*–A* → A 2(g) + *–*). Rigorous description of kinetically-consequential islanding/partitioning phenomena requires explicit description of the coverage and chemical dynamics of all multi-site ensembles. Higher-order, ensemble-specific rate terms identify the particular microscopic events relevant to each ensemble, and, in doing so, reveal that each elementary step (e.g. A (g) adsorption) describes not one event (e.g. A (g) + * → A*) , but a sum over all ensemble-specific paths (e.g. A (g) + *–* → A*–* and A (g) + A*–* → A*–A*). De-convoluting each elementary step into its constituent multi-site paths proffers kinetic detail otherwise inaccessible – enabling (i) identification of rate- and selectivity-determining site ensembles, (ii) calculation of rates and degrees of rate control of ensemble-specific elementary steps, (iii) incorporation of adsorbate surface diffusion, (iv) incorporation of lateral adsorbate interactions, and (v) quantitative description of catalysis of multi-site-occupying intermediates (e.g. *C n H m * species in hydrocarbon (de-)hydrogenation and C-C bond coupling/cleavage reactions) which we demonstrate here is inaccessible to the Langmuir-Hinshelwood formalism even if adsorbate surface diffusion is infinitely-fast.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗