Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Enthalpy”

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

Implications of rootless geothermal models: Missing processes, parameter compensation, and imposter convection

Numerical models of geothermal systems commonly capture only the top of a reservoir. Deeper areas of the reservoir are simplified to a boundary at the base of the model domain. Commonly, the basal boundary is given either a heat source or a source of mass and enthalpy. Here we developed and present simple numerical experiments which demonstrate that these approaches do not produce the correct model behavior in comparison to a model that captures the entire convecting domain with a heat flux only. We describe a variety of incorrect types of model behavior that arise directly from the choice of boundary condition, independent of the specific parameterization of the model. Heat sources are sensitive to the thickness of the domain and parameters take unphysical values to compensate for the reduced height. The combined mass/heat boundary can produce temperatures that show similarities to circulating geothermal systems, but with incorrect fluid flow and a strong boundary layer focused at the base of the simulated clay cap. These errors likely cause parameters to adjust their values to compensate for the incorrect physics. We highlight these issues and show an example from a developed reservoir model. Initial calibrations to natural state temperature were unsuitable for history matching. The 3D model required parameter adjustments to achieve a more realistic production model. Appropriate mitigation measures should be considered to reduce parameter compensation and improve decision-support models.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Use of tritium-rich fuel to improve the yield of layered deuterium/tritium inertial fusion capsules

In deuterium–tritium (DT) ice layered implosions, nearly all hot spot mass at peak burn comes from the dense fuel. Accurate prediction of the fuel mass ablation, including the enthalpy associated with mass inflow into the hot spot from the dense fuel, is essential to understanding the energetics and ignition of the hot spot in layered implosions. A recently published boundary layer analysis (Daughton et al., 2023) indicates a faster mass ablation rate than in previous analyses of layered implosions. Inclusion of this effect provides a better match to simulations and leads to a new ignition threshold where the temperature of the dense fuel plays a critical role. This analysis motivates possible new directions for improved capsule performance. Here, the authors present evidence in support of one such approach: the use of tritium-rich ice to decrease 14 MeV neutron scattering and heating of the dense fuel, resulting in less mass ablation and more robust burn of the hot spot. It is found from numerical simulations that despite a less favorable D:T ratio in the ice, the use of a 40:60 D:T ratio leads to an increase in capsule yield of 17% percent compared with that of a 50:50 D:T ratio fuel for capsules resembling those of the recent N210808 ignition experiment on the NIF (Abu-Shawareb et al., 2022) and an increase of 74% compared with that of a 60:40 D:T ratio fuel capsule. Finally, these results are potentially important for modeling all layered implosions, since some degree of DT fractionization may arise naturally during the beta layering process. In addition, this physics is important for the feasibility of high-gain capsule designs that seek to minimize tritium usage, as in some inertial fusion energy concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental measurements and mathematical modeling of cold plate for aviation thermal management

Herein, this study, which has been motivated by the recent applications of the cold plate device in aviation thermal management, reports on physics-based mathematical models derived from the conservation laws of mass, momentum, and energy, and empiricism-based models. One of the objectives of the present work is to report on an elaborate and successful experimental work carried out on an additively manufactured device for the purpose of rigorously validating the numerical predictions. The excellent agreement between the numerical predictions and measured performance provides much needed confidence in the implementation, in the software package, of the offset-strip fin passage correlations, as well as in the software implementation of user-defined wavy fin correlations for aerospace heat exchangers and cold plates operating with ram air at Reynolds numbers below 8000. The contributions of this work can also be found in the development of a new and accurate thermal-hydraulic analysis procedure, referred to in this paper as plate-fin analogy. Results from this procedure are compared with those from thermal resistance network. The comparative study in this paper of the bulk and discrete enthalpy flux method is also new, as is the relative assessment of four off-set strip fin thermal-hydraulic models.

42 ENGINEERING↗
DOE OSTI2025

Nature-inspired lotus-shaped fins combined with hybrid nanoparticles and metal foam for high-performance latent heat thermal energy storage

Latent heat thermal energy storage (LHTES) systems play a critical role in renewable energy integration by providing high energy density and nearly isothermal operation during phase transitions. However, their performance is often limited by slow melting/charging rates, which motivates the search for enhanced heat transfer designs. This study investigates the melting behavior of RT-82 phase change material (PCM) using novel lotus-shaped fins combined with copper metal foam and conductive graphene nanoparticles and carbon nanotubes. A two-dimensional enthalpy-porosity model in ANSYS Fluent was developed to simulate the charging/melting process, capturing non-thermal equilibrium between the foam and PCM/nano-PCM. In this study, effects of fin geometry, nanoparticle concentration, and foam porosity on melting dynamics and cost-performance trade-offs were investigated. Results showed that natural convection accelerated melting by ~12% compared to conduction-only scenarios. Optimized lotus-shaped fins with higher fin density (T3F4 and T3F10) achieved up to 63% faster melting relative to sparse configurations. Graphene nanoparticles improved thermal conductivity, with a 6% volume fraction, by reducing melting time by ~6.9%, while their combination with 75% porosity foam achieved a maximum reduction in the melting time of ~51% compared to pure PCM. Cost-performance analysis identified T3F4 as the most balanced design, offering rapid thermal response without excessive material costs, while moderate-density designs like T3S6 provided economical alternatives with acceptable performance. These results highlight the performance enhancement that can be achieved by integrating bio-inspired fins, nanoparticles, and foams, into compact and efficient LHTES for solar heating, building thermal management, and industrial waste-heat recovery applications.

25 ENERGY STORAGE↗
DOE OSTI2026

Technoeconomic analysis of hydrogen storage using 1,4-butanediol (BDO)/γ-butyrolactone (GBL) as a stationary backup power system

Liquid organic hydrogen carriers (LOHCs) are compounds that store and release hydrogen in stable forms at high density. While one-way carriers such as methanol and ammonia have gained attention, their economic advantages are often realized from their use as an export product and direct use as a fuel. Liquid carrier materials that can instead be cycled for energy storage have promise for stationary power applications. In particular, the reversible LOHC system 1,4-butanediol (BDO, H 2 -rich) and gamma-butyrolactone (GBL, H 2 -lean) has a lower enthalpy of dehydrogenation compared to conventional cyclic hydrocarbons and can utilize non-precious metal copper-based catalysts. Here, in this study, BDO/GBL system capital and operating expenses are characterized for vapor phase versus liquid phase hydrogenation and dehydrogenation in a 10 MW backup power application corresponding to sizing of Tier 2 datacenters as well as other critical infrastructure such as hospitals. Costs are benchmarked against two incumbent technologies: a well-established methylcyclohexane/toluene carrier system and compressed gas storage. BDO/GBL storage costs are found to differ substantially between operating modes, with liquid phase hydrogenation coupled with liquid phase dehydrogenation leading to the lowest LCOS of $\$$4.58/kg H 2 in the absence of byproduct formation. In this bounding case, LCOS for the BDO/GBL system is lower than for MCH/TOL ($\$$6.97/kg H 2 ) and compressed gas ($\$$8.48/kg H 2 at 170 bara and $\$$12.05/kg H 2 at 350 bara). However, escalating costs of carrier replacement due to byproduct formation (ranging from an added $\$$6–11/kg H 2 ) illustrate the need for highly selective catalysts to ensure BDO/GBL carrier viability.

BDO/GBL↗
DOE OSTI2020

Development of a CPU/GPU portable software library for Lagrangian–Eulerian simulations of liquid sprays

The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Approximation of refrigerant thermophysical properties using neural networks to speed up transient thermofluid simulations

Accurate and efficient evaluations of refrigerant thermophysical properties and their partial derivatives are essential for transient simulations of thermofluid systems, where several computations need to be executed at each integration time step. Since the utilization of an Equation of State for retrieving properties based on a pair of independent inputs typically involves numerical iterations in solution procedures, when the input variables differ from the refrigerant state variables employed in dynamic models, a variety of approaches including lookup table interpolation and curve fitting have been developed to explicitly approximate these properties based on the state variables, and consequently eliminate internal iterations. This paper presents an alternative method that exploits derivative-informed neural networks to model refrigerant properties explicitly from inputs of pressure and enthalpy, while ensuring consistent partial derivatives generated by differentiating the neural networks. Computational speed and accuracy of the proposed approach are demonstrated via transient simulations of a discretized heat exchanger model in Modelica, and comparisons against other property evaluation routines. Simulation results indicate that the proposed approach can realize a significant speedup with negligible discrepancies in predicted transients. The method is implemented in an open-source Modelica library.

Ma, Jiacheng↗
DOE OSTI2022

Product specific thermal degradation kinetics of bisphenol F epoxy in inert and oxidative atmospheres using evolved gas analysis–mass spectrometry

Knowledge of the degradation kinetics for polymer materials is important for understanding thermal stability. In this study, evolved gas analysis–mass spectrometry and pyrolysis gas-chromatography-mass spectrometry were evaluated for the potential to deliver additional insight into thermal degradation kinetics of diglycidal ether of bisphenol F (DGEBF) epoxy thermoset under inert and oxidative atmospheres. Degradation products of selected precursor ions were evaluated for their uniqueness to the specific precursor using extracted ion thermographs. Unique mass peaks, solely attributed to a single reaction pathway of a specific product, were determined from extracted ion thermographs and used to determine both activation energy (E a ) and pre-exponential factors for the specific primary reaction pathways. These primary reaction pathways for DGEBF epoxy degradation were then evaluated in the context of transition state theory (TST) and related transition state enthalpies (ΔH ‡ ) and entropies (ΔS ‡ ) of activation to further elucidate the degradation process. It was determined under pyrolysis conditions, as suggested by the E a , the formation of bisphenol F monomer was the rate-limiting step toward the formation of xanthene and phenol. In contrast, under thermo-oxidative conditions, reactions involving oxygen containing species were identified as the rate-limiting step for all observed products based on the large negative ΔS ‡ calculated from TST. This work demonstrates a powerful combination of technique and theory that can provide new insight into the degradation of polymer materials.

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

Primary solidification of ternary compounds in Al-rich Al–Ce–Mn alloys

Primary solidification of ternary compounds Al 20 Mn 2 Ce and Al 10 Mn 2 Ce were analyzed through the coupling of the thermodynamic modeling and classic nucleation theory. Thermodynamic models of Al 20 Mn 2 Ce and Al 10 Mn 2 Ce were developed using the CALPHAD approach based on first-principles calculated enthalpy of formation and experimental data obtained from this work and the literature. The analysis suggested that despite the larger thermodynamic driving force for nucleation of Al 10 Mn 2 Ce, nucleation is preferred for the Al 20 Mn 2 Ce phase in the highly undercooled liquid due to its smaller interfacial energy. Therefore, manufacturing methods with rapid cooling rates will favor primary solidification of Al 20 Mn 2 Ce for Al-rich Al–Ce–Mn alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Probing thermodynamics of radiogenic helium and defects in $δ$-plutonium alloys and interactions with adsorbed environmental gases

Differential scanning calorimetry coupled with simultaneous evolved gas analysis (DSC-EGA) on aged δ-Pu samples shows that most radiogenic helium remains trapped within the Pu matrix at temperatures very close to, or slightly above, the melting temperature. Our results indicate that helium release from 50-year-old δ-Pu occurs as a burst just below the melting temperature (>0.994 T m ), with subsequent pressure oscillations as temperature increases. Subordinate quantities of H 2 were also released along with helium. Here, the helium emission tails off and ceases above ~720 – 750°C. In a δ-Pu alloy aged 6 years, the initial helium burst occurs slightly above melting (~1.015 – 1.042 T m ), with a discrete, larger helium spike occurring between 670 and 686°C. The proximity of helium release to the liquidus transition presented challenges in the deconvolution of overlapping process enthalpies, the liquidus endotherm, and the exotherm resulting from bubble collapse, annealing and gas expulsion. Helium’s strong affinity for vacancy binding in a 2He-vac configuration is predicted by Density Functional Theory (DFT) modeling. The measured stored energy associated with the He release events in a 50-year-old δ alloy is on the order of ~ 10–11 J/g, which is significantly higher than stored energies measured in the sub-solidus regimes (~2 J/g) that are related to the solid-state annealing of processing- and radiation-induced defects. This implies that aged δ Pu alloys have a remarkable resilience to accommodate the lattice strain produced by the internal pressure of the helium bubbles and provides further insight into the thermodynamic behavior of aged δ Pu.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Degree of rate control and De Donder relations – An interpretation based on transition state theory

We describe sensitivities of elementary steps as salient parameters in determining the rate determining character of elementary steps within a reaction network and develop a formalism wherein the overall composite reaction is described by an apparent rate-determining step that is a sensitivity-weighted average of the elementary steps that comprise the reaction. Reaction parameters—apparent reaction orders, apparent enthalpy and entropy of activation—of the composite reaction network are determined within the framework we develop by application of transition state theory to the apparent rate-determining step. From this formalism we develop methods for determination of surface coverages by measuring only reaction orders, novel methods for discrimination between proposed mechanisms, a proof that the kinetic degrees of rate control sum to unity, and a proof for the relationship between fractional coverages and thermodynamic degrees of rate control. Two complementary formalisms for identifying rate-limiting transition states are broadly employed in chemical kinetics—De Donder relations based on measuring the thermodynamic driving forces of elementary steps during reaction and degrees of rate control based on knowledge of the kinetics of elementary reaction steps. Here, the formalism developed herein unifies these two strategies to elucidate in the most general case a relationship between the thermodynamic driving forces of elementary steps captured by the reversibilities and the kinetic and thermodynamic degrees of rate control.

42 ENGINEERING↗
DOE OSTI2021

Effects of water on the kinetics of acetone hydrogenation over Pt and Ru catalysts

We employed an approach combining reaction kinetics measurements at steady state conditions, electronic structure calculations employing density functional theory, and microkinetic modeling for acetone hydrogenation to provide insights into the effects of water on metal catalyst surfaces for the hydrogenation of oxygenates over a wide range of reaction conditions. Elucidation of the repulsive interactions due to adsorbed water molecules at various reaction conditions provides a basis to formulate rate expressions for heterogeneous catalytic processes of biomass oxygenates. Reaction kinetics experiments were carried out at partial pressures of H 2 , acetone, water and helium in the range of 0.51-0.79, 0.02-0.13, 0.08-0.23, 0-0.28 atm, respectively. We show that the addition of water enhances the hydrogenation rate at 353 K and 1 atm on oxophilic metal catalysts such as Ru/C, whereas the same promotional effect of water is not observed for Pt-based catalysts. Microkinetic model predictions for the hydrogenation of acetone on Ru in the absence and presence of water, using enthalpies and entropies obtained from DFT calculations, were in agreement with the experimentally observed reaction orders and activation barriers. The model shows that a water-assisted hydroxypropyl path is expected to be the favored path on Ru with a rate-determining step of H-OH-mediated hydrogenation of C 3 H 6 OH (i.e., the hydroxypropyl intermediate formed by H 2 O-mediated initial hydrogenation of acetone) to produce isopropyl alcohol (IPA). Furthermore, hydrogen, acetone, hydroxypropyl intermediate and hydroxyl species were predicted to be abundant on the Ru surface with a high coverage of nearly 85%. Lastly, the combined studies of computational and experimental catalysis on hydrogenation reactions help to elucidate the mechanistic role of water on metal catalyzed reactions for producing chemical building blocks from biomass-derived oxygenates.

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

Experimental and computational kinetics study of the liquid-phase hydrogenation of C=C and C=O bonds

Solvent effects on adsorption equilibrium and reaction kinetics are evaluated for hydrogenation reactions catalyzed by Pd/alumina in a series of different solvents. Three reactants – cyclohexene, benzene and benzaldehyde – and three solvents – n-heptane (HEP), methylcyclohexane (MCHA), decalin (DL) – have been investigated. Kinetic analysis of hydrogenation of cyclohexene reveals that hydrogen adsorbs on different sites from those where cyclohexene and the solvents adsorb; however, the presence of hydrogen on these separate sites affects the heats of adsorption of the hydrocarbons. When the solvent is a weakly interacting linear alkane (HEP), the rate determining step of the reaction is the first hydrogenation of adsorbed cyclohexene. Furthermore, this conclusion is supported by DFT calculations that show a higher enthalpy barrier for the first hydrogenation than for the second, while statistical thermodynamics analysis validates the physical significance of the entropy of adsorption parameters derived from the kinetic fitting of experimental data. By contrast, with solvents such as MCHA and DL, which interact more strongly with the metal surface and compete for active sites with the reactant and the surface intermediate, the rate limiting step seems to shift to the second hydrogenation step.

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

H 2 O-assisted O 2 reduction by H 2 on Pt and PtAu bimetallic nanoparticles: Influences of composition and reactant coverages on kinetic regimes, rates, and selectivities

Hydrogen peroxide (H 2 O 2 ) can replace hazardous oxidants in industrial processes but is currently too expensive for many such applications. While direct synthesis of H 2 O 2 (H 2 + O 2 → H 2 O 2 ) may reduce costs in comparison to incumbent technology, current catalysts lack the requisite stability and selectivity. Here, we examine the direct synthesis of H 2 O 2 on bimetallic Pt 1 Au x (0 ≤ x ≤ 230) and Pt catalysts at steady-state in pure water and relate kinetic parameters for H 2 O 2 and H 2 O formation to possible active site structures informed by complementary characterization methods. X-ray photoelectron spectra show significant Pt surface enrichment compared to the bulk composition. Analysis of infrared spectra of mixed monolayers of 12 CO* and 13 CO* indicate that Pt and Au form substitutional surface alloys. The Pt 1 Au x nanoparticles with the greatest mole fractions of Au predominantly expose Pt monomers (i.e., isolated Pt atoms), yet Pt atoms exposed upon all these nanoparticles possess electronic structures distinct from bulk Pt. Despite these differences, rate measurements are consistent with product formation through proton-electron transfer pathways for all Pt 1 Au x catalysts. In situ XAS indicate that Pt remains metallic during H 2 O 2 synthesis. Under the most oxidizing conditions, selectivities toward H 2 O 2 increase strongly with the Au to Pt ratio from 2% for monometallic Pt to 85% for Pt 1 Au 170 . However, selectivities are similar among all catalysts within reducing conditions. Comparisons of apparent activation enthalpies for the formation of H 2 O 2 and H 2 O across these catalysts and the range of conditions suggest that Pt monomers within Au provide the greatest selectivities for H 2 O 2 formation, because these active sites present high barriers for O-O bond rupture. Further, selectivities decrease with increasing ratios of H 2 to O 2 pressures, because Pt atoms aggregate and form oligomers that readily dissociate dioxygen intermediates. The combined use of spectroscopy, kinetics, and concepts employed in reaching these conclusions take inspiration from the legacy of Prof. Michel Boudart, and specifically his elegant methods for interrogating bimetallic catalysts.

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

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

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

Intrinsic activation energies for ring contraction of allylic cations in zeolites

Cyclic carbocations are important intermediates in zeolite-catalyzed chemistries such as methanol-to-hydrocarbons conversion, naphthenes ring opening, or coke formation. While information about their thermodynamic stability exists, little is known about the kinetics of formation and transformation of cyclic carbocations in zeolites. To fill this knowledge gap, ring contraction of the 1,3,5,5-tetramethylcyclohexenyl cation (C 10 H 17 + ), a representative of 6-membered ring allylic cations, was investigated by in situ UV–vis and IR spectroscopy. Protonic forms of zeolites served as catalysts, at temperatures from 80 °C to 135 °C. Significant oligomerization and hydride transfer in BEA and FAU hampered kinetics analysis, whereas ring contraction dominated in the channels of MOR. The reactant cation, characterized by an electronic absorption at 314 nm and an allylic stretch at 1549 cm −1 , contracted to both a 1,3-alkyl-substituted cyclopentenyl cation (287 nm and 1506 cm −1 ) and a 1,2,3-alkyl-substituted cyclopentenyl cation (297 nm and 1489 cm −1 ). Collection of time-resolved IR spectra and fitting of the intensities with various kinetic models revealed a third transformation, which is expected from thermodynamics: the 1,3-alkyl-substituted cyclopentenyl cation isomerizes to the 1,2,3-alkyl-substituted cyclopentenyl cation. Series of IR spectra recorded at different temperatures delivered intrinsic activation enthalpies (entropies) in MOR of 67 ± 2 kJ mol −1 (−130 ± 6 J mol −1 K −1 ) and 90 ± 3 kJ mol −1 (−70 ± 8 J mol −1 K −1 ) for the contraction to 1,3- and 1,2,3-substituted species, and of 88 ± 5 kJ mol −1 (−84 ± 12 J mol −1 K −1 ) for the isomerization. The findings characterize one path – via contraction of larger rings – to different cyclopentenyl species in zeolites; and the associated, moderate activation energies suggest such transformations contribute to many complex hydrocarbon reaction networks.

Acid catalysis↗
DOE OSTI2021

Molecular-level understanding of gibbsite particle aggregation in water

Using molecular dynamics simulations, we investigate the molecular scale origin of crystal face selectivity when one gibbsite particle attaches to another in water. A comparison of the free energy per unit surface area of particle–particle attachment indicates that particle attachment through edge surfaces, where the edge surfaces are either (1 0 0) or (1 1 0) crystal faces, is more energetically favorable compared to attachment between two basal surfaces (i.e., (0 0 1) crystal faces) or between the basal surface of one particle and the edge surface of another. This result suggests that gibbsite crystals with low basal/edge surface area ratio will preferentially attach through edge surfaces, potentially helping the crystals grow laterally. However, for larger gibbsite particles (high basal/edge surface area ratio) the total free energy, not normalized by surface area, of particle attachment through the basal surfaces is lower (more negative) than attachment through the edge surfaces, indicating that larger gibbsite particles will preferentially aggregate through basal surface attachments. The short-range electrostatic interactions including the interparticle hydrogen bonds from surface –OH groups drive particle attachment, and the dominant contribution to the free energy minimum is enthalpic rather than entropic. However, the enthalpy of basal-edge attachment is significantly offset by the entropy leading to a higher free energy (less negative) compared to that of basal-basal attachment. Study of the free energy for a few imperfect attachments of two particles indicates a higher free energy (i.e., less negative, less stable), compared to a perfect attachment

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

Probing surface-adsorbate interactions through active particle dynamics

Adsorbate molecules present in a reaction mixture may bind to and block catalytic sites. Measurement of the surface coverage of these molecules via adsorption isotherms is critical for modeling and design of catalytic reactions on surfaces. However, it is challenging to measure isotherms in solution in a way that is directly relevant to catalytic activity under reaction conditions, particularly since adsorbates may bind with an enormous range of surface affinity parameters. Here we used the motion of self-propelled catalytic Janus particles, which employ the decomposition of hydrogen peroxide fuel as a propulsion mechanism, to determine the effective surface coverage of thioglycerol, furfural, and ethanol on a platinum surface as a function of concentration in aqueous solution by measuring the decrease in active motion due to the blocking of active sites. For strongly adsorbing thioglycerol, this effective coverage was compared and contrasted to the total adsorbed amount measured using inductively-coupled plasma analysis. Demonstrating the broad applicability of this approach, the surface affinity of the three adsorbates spanned more than four orders of magnitude. For each species, the adsorbate-mediated attenuation of active motion occurred over a wide concentration range and was well-described by a Langmuir isotherm. The strongly interacting thioglycerol had the highest affinity towards the surface (K a = 15.5 ± 4.3 mM –1 ) and fully deactivated the active particle motion at surface saturation. Furfural had an intermediate affinity (K a = 0.42 ± 0.07 mM –1 ) but did not fully block H 2 O 2 access to the surface at apparent saturation, consistent with a maximum fractional surface coverage of θ max = 0.67. Ethanol exhibited even lower affinity (K a = 0.0025 ± 2x10 -4 mM –1 ) and its coverage saturated at only θ max = 0.38. Analysis of isotherms at elevated temperatures enabled direct extraction of the enthalpies of adsorption. The degree of surface coverage at adsorbate saturation appeared to correlate with the relative energies of adsorption for the different adsorbate species and was consistent with adsorbate saturation of one of multiple active site populations towards H 2 O 2 decomposition. Furthermore, computational investigations into solvent effects on furfural adsorption showed good quantitative agreement with the experimental results. This work leverages unique properties of active particles to explore fundamental catalysis questions and demonstrates a novel paradigm for significant and experimentally accessible multidisciplinary research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗