Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state process”

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 541 records · Page 30

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Optimization strategies for produced water networks with integrated desalination facilities

Optimal management and desalination of produced water is a major challenge for U.S. oil and gas development. Integrating rigorous desalination models into multi-period produced water network optimization problems presents several hurdles, which need to be tackled using advanced optimization strategies. Here, in this work, a novel multi-period produced water network formulation with separate solid and liquid flows is introduced to avoid singularities at zero flows. Rigorous steady state desalination models based on mechanical vapor recompression are embedded at the desalination sites in the network model. An integrated optimization formulation is developed to co-optimize the design of desalination units along with the operation of the network. Furthermore, a more robust approach based on the trust region filter method is developed to efficiently integrate complex desalination models into the multi-period planning problem. Both optimization approaches are demonstrated on a produced water network from the PARETO library (Drouven et al., 2022) using thermal desalination units. Our results show that while the TRF and integrated approaches have comparable solve times, the TRF approach has better performance reliability in terms of solver convergence. Furthermore, the optimal solution obtained by embedding rigorous models into the network is significantly different than when desalination costs are approximated using simple cost models, which motivates further research in this field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Consideration of Nonideal Detonation Regimes Influenced by Wave Modes in a Water-Cooled Rotating Detonation Engine Using OH* Chemiluminescence

Although inherently unstable, existing research in rotating detonation combustion supports its application in notionally steady processes resulting in greater availability compared to conventional, constant pressure combustion. Further improvements rely on a more in-depth understanding of system losses and identifying conditions which optimize device performance. Within this study, the presence and proportion of ideal and nonideal combustion regimes are compared across a variety of process conditions and wave modes. Here, large-scale data analysis seeks to summarize proportional heat release associated with commensal, parasitic, and detonative combustion averaged across individual traces of OH* chemiluminescent data acquired at the detonation plane. Means of regime partitioning based on the anatomy of the time-resolved OH* signal are proposed to ensure consistent analysis throughout the current and future studies concerning combustion regimes. Of particular interest is the possible influence of wave on the nonideal combustion in relative proportion to the desired detonation. Results showed improved percent detonation with increasing significance for the following trends: decreasing equivalence ratio, increasing wave count, decreasing wave velocity, and increasing detonation time. Increased wave number, brought on by decreased equivalence ratios and wave velocities, is thought to decrease fill region surface area, and therefore, decrease nonideal contact burning. Additional analysis is performed to consider possible trend variation due to the presence of stable galloping waves, which were found to have minimal influence on relative percent detonation values. The outcome of this study suggests operational states, which correspond to increased wave quantities for increased proportions of reactants consumed by the targeted detonative combustion regime.

42 ENGINEERING↗

Multifunctional Catalysts for the Tandem Reactions of Oxygenates

Industrially-relevant catalytic reactions rarely consist of a simple sequence of elementary steps. Moreover, kinetic coupling of multiple reactions on a catalyst surface is highly desired for process intensification and improved energy efficiency for large scale chemical transformations. The shifting landscape of hydrocarbon chemical feedstocks in the US also motivates research on the selective conversion of more complex molecules. One desirable type of catalytic reaction is the reduction of carboxylic acids that are produced from biomass feedstocks to their corresponding alcohols. The proposed research explores the fundamental importance of hydrogen spillover on a multifunctional catalyst for carboxylic acid reduction with H 2 composed of metal particles coupled to metal oxide particles. Recent work has demonstrated the excellent performance of supported tungsten oxide clusters for carboxylic acid reduction, but only after they are promoted with a late transition metal such as palladium. Elucidating the active state of the catalyst and the associated reaction mechanism for acid reduction on that active state are the overall goals of the proposed project and successful completion will enable future design of efficient multifunctional catalysts. The critically important role of the metal promoter is hypothesized to be its ability to dissociate H2 and spillover atomic H to the support. Although hydrogen spillover is a well-recognized phenomenon in catalysis, its role in both catalyst activation and catalytic turnover are still unresolved. The study combined materials synthesis, characterization, reactivity testing, and molecular simulations, to explore the effect of hydrogen chemical potential on the formation of the active catalytic sites and on the steady state catalytic reduction of carboxylic acid. Varying the hydrogen chemical potential through modification of the gas conditions, support composition, and metal loading to modulated the structure and catalytic performance of the reducible metal oxide. Dual function catalysts containing supported Pd and WO x species co-located on a non-reducible carrier (silica) and a reducible carrier (titania) were synthesized and characterized by electron microscopy, temperature-programmed reduction, and chemisorption. Spectroscopic methods such as X-ray absorption and UV-vis were also used to evaluate the catalysts, which were used in the reduction of propionic acid to aldehyde and alcohol. Quantum chemical calculations, including ab initio phase diagrams provided molecular insights into the H spillover phenomenon.

09 BIOMASS FUELS↗

Marine ice sheet experiments with the Community Ice Sheet Model

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ~2 km, as opposed to ~0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ~2 km (arguably 4 km) when the connectivity is low or moderate and ~1 km (arguably 2 km) when the connectivity is strong.

54 ENVIRONMENTAL SCIENCES↗

Cryptic Methane-Cycling by Methanogens During Multi-Year Incubation of Estuarine Sediment

As marine sediments are buried, microbial communities transition from sulfate-reduction to methane-production after sulfate is depleted. When this biogenic methane diffuses into the overlying sulfate-rich sediments, it forms a sulfate-methane transition zone (SMTZ) because sulfate reducers deplete hydrogen concentrations and make hydrogenotrophic methanogenesis exergonic in the reverse direction, a process called the anaerobic oxidation of methane (AOM). Microbial participation in these processes is often inferred from geochemistry, genes, and gene expression changes with sediment depth, using sedimentation rates to convert depth to time. Less is known about how natural sediments transition through these geochemical states transition in real-time. We examined 16S rRNA gene amplicon libraries and metatranscriptomes in microcosms of anoxic sediment from the White Oak River estuary, NC, with three destructively sampled replicates with methane added (586-day incubations) and three re-sampled un-amended replicates (895-day incubations). Sulfate dropped to a low value (∼0.3 mM) on similar days for both experiments (312 and 320 days, respectively), followed by a peak in hydrogen, intermittent increases in methane-cycling archaea starting on days 375 and 362 (mostly Methanolinea spp. and Methanosaeta spp., and Methanococcoides sp. ANME-3), and a methane peak 1 month later. However, methane δ 13 C values only show net methanogenesis 6 months after methane-cycling archaea increase and 4 months after the methane peak, when sulfate is consistently below 0.1 mM and hydrogen increases to a stable 0.61 ± 0.13 nM (days 553–586, n = 9). Sulfate-reducing bacteria (mostly Desulfatiglans spp. and Desulfosarcina sp. SEEP-SRB1) increase in relative abundance only during this period of net methane production, suggesting syntrophy with methanogens in the absence of sulfate. The transition from sulfate reduction to methane production in marine sediments occurs through a prolonged period of methane-cycling by methanogens at low sulfate concentrations, and steady growth of sulfate reducers along with methanogens after sulfate is depleted.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling and analysis of synthetic liquid fuel production from CO 2 and nuclear energy using methanol-to-diesel process

Electrofuels (e-fuels) are synthetic fuels produced from carbon dioxide (CO 2 ) and electricity for blending with or replacing petroleum fuels. Nuclear energy is an attractive energy feedstock for e-fuel production because of its low environmental footprint and its ability to provide steady heat and power essential for e-fuels production. We modeled and evaluated the cost and environmental footprint of e-fuels production in the distillate range for three nuclear power scales, 100, 500, and 1000 MWe, through methanol and olefins intermediates leveraging commercial or high technology readiness level (TRL) processes. Compared to the commonly studied e-fuels from Fischer Tropsch process that has a distillate yield of <70% with the rest being low value naphtha, the proposed process via methanol intermediate increases the product selectivity with distillate yield of 96% and only 4% naphtha. The modeled process has a carbon conversion ratio of 98%, and a process energy efficiency of 56% relative to the total equivalent nuclear electricity input. The e-fuel plant economics and GHG emissions were estimated by considering CO 2 collected from ethanol plants adjacent to nuclear power plants. The estimated minimum fuel selling prices (MFSP) of e-fuel is in the range of $5.7-$9.1/gal depending on e-fuel plant scale, electricity cost, and CO 2 transportation distance. The corresponding e-fuels life cycle GHG emissions is estimated in the range of 5-6 gCO 2 e/MJ of liquid fuel using the R&D Greenhouse gases, Regulated Emissions, and Energy use in Technologies (R&D GREET) model.

10 SYNTHETIC FUELS↗

Exploring the impact of surface topography on Rayleigh-Bénard dry convection in the Pi cloud chamber using OpenFOAM: In cylindrical and rectangular geometries

The Pi convection-cloud chamber can generate steady-state turbulence in both rectangular and cylindrical shapes via Rayleigh-Bénard convection (RBC) by maintaining warm bottom and cold top surfaces. Although most experiments in the Pi chamber were conducted in cylindrical shapes, all previous Pi chamber simulations were conducted in a rectangular shape due to the limitations of those models to discretize a cylindrical domain when using the finite difference method therein. Here, we use OpenFOAM, an open-source finite-volume-based Computational Fluid Dynamics (CFD) software package, to conduct Large-Eddy Simulation (LES) of dry RBC in the Pi chamber at high Rayleigh numbers (10 8 to 10 9 ). Results show that large-scale circulation (LSC) direction varies in the chamber with a constant side wall temperature. Imposing a slight temperature imbalance at the side wall ranging from 0.1 to 0.7 degrees can lock the LSC, aligning better with Pi chamber observations, particularly at higher Rayleigh numbers. In addition, we examine the impact of surface topography on LSC and heat transfer in RBC systems within cylindrical and rectangular shapes under varying conditions. Results show that roughing top/bottom surfaces by adding bars of a few tens millimeters height can strengthen thermal plumes and enhance temperature fluctuations in the chamber. Furthermore, we observe that different bar height configurations lead to notable changes in LSC orientation and thermal stratification, highlighting the complex interactions between surface features and convection patterns. This finding highlights how surface topography and chamber geometry affect Rayleigh-Bénard convection, improving understanding of turbulent heat transfer and atmospheric boundary-layer processes. Direct Numerical Simulations (DNS) are also conducted to validate LES results. In conclusion, while LES effectively captures qualitative behaviors seen in DNS, it tends to underestimate velocity variances near walls, illustrating a trade-off between computational efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗

Models implemented in the methodological approach to design the initial STEP first wall contour

The official Spherical Tokamak for Energy Production mission aims to demonstrate the ability to generate net electricity from fusion with the STEP Prototype Power plant. One of the key technological and engineering challenges in fusion power plants is managing the loads on the first wall within acceptable limits. Therefore, the conceptual design development of the STEP Prototype Power plant needs to be based on load estimates derived using legitimate plasma physics assumptions through dynamic and flexible tools. The current design foresees the STEP main chamber first wall to withstand steady-state heat loads of up to ~1 MW/m 2 , excluding critical regions expected to receive higher heat loads such as the baffle regions approaching the divertors. These critical areas will require ad hoc assessments and will be designed with the presence of limiters. This article focuses on the models and methodology adopted for designing the 2-D poloidal contour of the STEP first wall, based on the anticipated charged particle and radiation heat loads during normal operation. Firstly, the models adopted for calculating the charged particle and radiation heat loads are introduced. The first model is validated through benchmarking against the particle tracing code SMARDDA, while the second model is verified by comparing it with data from the MAST-U experiment. Secondly, the model used to design the 2-D first wall contour according to the heat loads is explained. We acknowledge that this preliminary design stage assumes certain simplifications, notably an axisymmetric geometry, for computational efficiency and clarity in presentation. It is understood that subsequent design phases will address the complexities of real-world engineering, including non-axisymmetric effects, transient plasma scenarios, and the impact of disruptions on the first wall design. Finally, an automatic procedure based on these models is presented for defining the 2-D poloidal contour of the STEP first wall to minimize heat loads, taking into account the need to radiate most of the alpha-particle and auxiliary heating power. Here, by providing an overview of the models, methodology, and an automatic procedure, this paper contributes to the design process of the STEP first wall, addressing the engineering challenges associated with fusion power plant development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isoenergetic two-photon excitation enhances solvent-to-solute excited-state proton transfer

Two-photon excitation (TPE) is an attractive means for controlling chemistry in both space and time. Since isoenergetic one- and twophoton excitations (OPE and TPE) in non-centrosymmetric molecules are allowed to reach the same excited state, it is usually assumed that they produce similar excited-state reactivity. We compare the solvent-to-solute excited-state proton transfer of the super photobase FR0-SB following isoenergetic OPE and TPE. We find up to 62% increased reactivity following TPE compared to OPE. From steady-state spectroscopy, we rule out the involvement of different excited states and find that OPE and TPE spectra are identical in non-polar solvents but not in polar ones. We propose that differences in the matrix elements that contribute to the two-photon absorption cross sections lead to the observed enhanced isoenergetic reactivity, consistent with the predictions of our high-level coupledcluster-based computational protocol. We find that polar solvent configurations favor greater dipole moment change between ground and excited states, which enters the probability for TPE as the absolute value squared. This, in turn, causes a difference in the Franck– Condon region reached via TPE compared to OPE. We conclude that a new method has been found for controlling chemical reactivity via the matrix elements that affect two-photon cross sections, which may be of great utility for spatial and temporal precision chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Microemulsions in the driven Widom-Rowlinson lattice gas

An investigation of the two-dimensional Widom-Rowlinson lattice gas under an applied drive uncovered a remarkable nonequilibrium steady state in which uniform stripes (reminiscent of an equilibrium lamellar phase) form perpendicular to the drive direction [R. Dickman and R. K. P. Zia, Phys. Rev. E 97, 062126 (2018)]. Here we study this model at low particle densities in two and three dimensions, where we find a disordered phase with a characteristic length scale (a “microemulsion”) along the drive direction. Further we develop a continuum theory of this disordered phase to derive a coarse-grained field-theoretic action for the nonequilibrium dynamics. The action has the form of two coupled driven diffusive systems with different characteristic velocities, generated by an interplay between the particle repulsion and the drive. We then show how fluctuation corrections in the field theory may generate the characteristic features of the microemulsion phase, including a peak in the static structure factor corresponding to the characteristic length scale. This work lays the foundation for understanding the stripe phenomenon more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Single-particle chemical analysis reveals organic-rich detonation soot products

The detonation of high explosives produces a wide variety of particulate matter (PM) with distinct properties, not all of which are traditionally studied for chemical composition, formation processes, and forensic applications. We report particle-resolved measurements of Composition B detonation soot using soot particle aerosol mass spectrometry (SP-AMS), identifying carbonaceous species and metals not previously characterized on a single-particle basis. Results are combined with scanning electron microscopy with energy dispersive X-ray spectroscopy (SEM-EDS) to enhance source-dependent signatures. Black carbon, including graphitic carbon and detonation nanodiamonds, contributed 50.5-71.4% of PM 2.5 mass, while 22.5-43.4% was non-refractory organic carbon, a previously overlooked component that exhibited a complex and varying composition. Detonations were performed with and without PMMA confinement and under steady and overdriven conditions. Distinct particles enriched with polycyclic aromatic hydrocarbons (PAHs) were observed in experiments utilizing confinement, with quantities dependent on manufacturing method. SEM-EDS validated SP-AMS findings of metals internally mixed with carbonaceous species and extended the particle size range to 100 µm. This work makes detecting non-refractory organics using single-particle techniques more feasible for detonation forensics and understanding high-energy soot formation. While this analysis uses offline sample collection, SP-AMS could be deployed for in-situ measurements of detonation PM 2.5 transported in the atmosphere.

36 MATERIALS SCIENCE↗

Characterization of a Low-Temperature Plasma (LTP) Ambient Ionization Source Using Temporally Resolved Monochromatic Imaging Spectrometry

The low-temperature plasma (LTP) probe is a common plasma-based source used for ambient desorption–ionization mass spectrometry (MS). While the LTP probe has been characterized in detail with MS, relatively few studies have used optical spectroscopy. In this paper, two-dimensional (2D) imaging at selected wavelengths is used to visualize important species in the LTP plasma jet. First, 2D steady-state images of the LTP plume for N 2 + (391.2 nm), He I (706.5 nm), and N 2 (337.1 nm) emissions were recorded under selected plasma conditions. Second, time-resolved 2D emission maps of radiative species in the LTP plasma jet were recorded through the use of a 200 ns detection gate and varying gate delays with respect to the LTP trigger pulse. Emission from He I, N 2 + , and N 2 in the plasma jet region was found to show a transient behavior (often referred to as plasma bullets) lasting only a few microseconds. The N 2 + and He I maps were highly correlated in spatial and temporal structure. Further, emission from N 2 showed two maxima in time, one before and one after the maximum emission for N 2 + and He I, due to an initial electronic excitation wave and ion–electron recombination, respectively. Third, the interaction of the LTP probe with a sample substrate and an electrically grounded metallic needle was studied. Emission from a fluorophore on the sample substrate showed an initial photon-induced excitation from plasma-generated photons followed by electronic excitation by other plasma species. The presence of a grounded needle near the plasma jet significantly extended the plasma jet lifetime and also generated a long-lived corona discharge on the needle. The effect of LTP operating parameters on emission spectra was correlated with mass-spectral results including reagent-ion signals. Lastly, five movies provide a side-by-side comparison of the temporal behavior of emitting species and insights into the interactions of the emission clouds with a sample surface as well as an external needle. Temporally and spatially resolved imaging provided insights into important processes in the LTP plasma jet, which will help improve analyte ion sampling in LTP–MS.

Chan, George C. -Y.↗

Aldol Condensation and Esterification over Ti-Substituted *BEA Zeolite: Mechanisms and Effects of Pore Hydrophobicity

Aldol condensation and esterification reactions provide paths to upgrade ethanol and acetaldehyde to higher- value molecules useful as fuels or intermediates for the synthesis of polymers. Transition-metal-substituted BEA zeolites (M-BEA) catalyze these reactions; however, the mechanisms for these processes in M-BEA and the effects of incidental or purposefully included silanol groups are not reported. Here, we combine kinetic and spectroscopic measurements obtained during catalytic reactions of acetaldehyde (CH 3 CHO), ethanol (C 2 H 5 OH), and hydrogen (H 2 ) mixtures over a series of Ti-BEA catalysts that possess a known range of silanol group densities to examine the kinetic relevance of intervening steps and the impact of silanol groups on catalytic rates. Across all Ti-BEA, rates for aldol condensation and esterification increase with the pressure of CH 3 CHO; however, C 2 H 5 OH and H 2 O weakly inhibit the rates of these reactions. The substitution of CD 3 CDO for CH 3 CHO decreases aldol condensation rates slightly (~10%) but leads to greater esterification rates (2- to 5-fold). The kinetic isotope effects together with the measured dependence of rates on reactant pressures suggest that aldol condensation and esterification occur on unoccupied Ti sites and involve multiple kinetically relevant steps. CH 3 CHO deprotonates irreversibly, and the kinetically relevant nucleophilic attack of the enolate to CH 3 CHO* (i.e., adsorbed CH 3 CHO on Ti sites) leads to aldol products, while the nucleophilic attack of the enolate to C 2 H 5 OH* gives esters. Selectivities toward aldol condensation increase with the ratio of CH 3 CHO to C 2 H 5 OH pressure and with increases in the silanol density of the as-synthesized Ti-BEA. During catalysis, in situ infrared spectroscopy demonstrates that these silanol groups react with C 2 H 5 OH to form ethoxysilane groups (i.e., SiOC 2 H 5 ) that modify the polarity of the environment near Ti active sites. As initial silanol densities increase, steady-state turnover rates for aldol condensation and esterification increase by factors of 5 and 2, respectively. The changes in rates and selectivities among Ti-BEA catalysts likely reflect changes in excess free energies of transition states for enolization and nucleophilic attack of the enolate to adsorbed coreactants. The differences in excess stability report on the interactions among reactive intermediates at framework Ti atoms and the ethoxysilane and remaining silanol groups present. Furthermore, the in situ modification of these pore environments confers changes in the stability of reactive species in a manner that contradicts intuition when considering the initial state of the catalyst but can be reconciled after accounting for the formation of persistent alkoxy surface moieties in the pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling Fundamental Activity–Stability Relationships in Rutile Oxides

The oxygen evolution reaction (OER) is a key anodic half-cell reaction that accompanies several critical electrochemical reduction reactions of interest to a variety of applications. Despite steady advances in understanding and qualitatively predicting OER activity and selectivity trends, a comprehensive description or prediction of material aqueous (in)stability and degradation mechanisms remains elusive, even though these processes critically influence device lifetime and economic feasibility. In this work, we investigate the interplay, or lack thereof, between OER activity and material aqueous stability across rutile oxides, with a particular focus on iridium oxide (IrO 2 ). By applying a Born–Haber cycle, we calculate the thermodynamic driving force for metal dissolution as a function of the applied bias and electrolyte conditions. We apply interpretable machine learning techniques, including principal component analysis and symbolic regression, to analyze trends across rutile oxides and find that key thermodynamic descriptors for OER activity and surface stability are only very weakly correlated. Instead, the local atomic environment─especially electronic structure signatures for interactions between the active site and its neighbors─plays a more important role in predicting material stability. Leveraging these insights, we investigate the impact of doping IrO 2 with a range of transition metals and show that the stability of Ir active sites can be tuned largely independently of its predicted OER activity. These insights lay the foundation for material design to improve stability with respect to corrosion, with the ultimate aim to enhance long-term stability without sacrificing catalytic performance in the OER.

evolution reactions↗

The Monolithic Heat Pipe Microreactor Reference Plant Model

This work introduces a reference plant model for a generic monolithic heat-pipe-cooled microreactor. The model will serve as a springboard to develop future evaluation models in the licensing process of similar microreactor designs at the U.S. Nuclear Regulatory Commission. This model has been developed with the Comprehensive Reactor Analysis Bundle and its specifications are based on open literature publications for the eVinci TM design. BlueCRAB is the U.S. Nu- clear Regulatory Commission non-light-water reactor analysis system based on the Multiphysics Object-Oriented Simulation Environment framework, which can couple the Griffin, BISON, and Sockeye applications to resolve the various physics that are essential for the safety analysis of this type of reactor system. The core specifications includes tristructural isotropic fuel, graphite monolith, graphite reflectors, and drums composed of graphite and B 4 C. No moderator or burnable poison pins are used in the design. The fuel enrichment is reduced to control excess reactivity in the core. This core design is not optimized and only serves for testing purposes, since the primary objective of this work is to exercise the multiphysics coupling for this type of reactor system. A three dimensional (3D) core heterogeneous Griffin discrete ordinates (SN) transport model allows the precise calculation of the flux distribution and pin powers. Griffin transfers the power density distribution and obtains a temperature distribution to and from BISON. The BISON model com- putes the 3D core temperature distribution and is coupled to 876 Sockeye subapplications running a heat pipe model. This 3D conduction model is coupled to the various heat pipes via heat flux boundary conditions. The model includes a small gap between the heat pipe and the monolith. Convective heat transfer boundaries with either ambient temperature or condenser temperature as heat sinks are imposed at the model boundaries. The 2D Sockeye heat pipe model uses a vapor- only methodology, which provides the needed resolution for transient calculations and allows the determination of various heat pipe limits. This approach is superior to the superconductor model traditionally used in steady-state calculations. BlueCRAB computes steady-state power and temperature distributions that serve as the initial condition for a loss-of-heat-sink transient simulation. The steady-state results show significant peaking due to the position of the control drum, but this is a characteristic of the particular design used, which is not optimized at this stage. The transient results show the reactor power slowly stabilizing towards a 3% power level after the partial loss of secondary heat removal. Several recriticalities are observed due to cooling through the secondary system but the reactor is self-stabilizing and behaves as expected.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Evaluation of BISON metallic fuel performance modeling against experimental measurements within FIPD and IMIS databases

Simulations were conducted using the BISON fuel performance code on an automated process to read initial and operating conditions from two databases—the Fuels Irradiation and Physics Database (FIPD) and Integral Fast Reactor Materials Information System (IMIS) database. These databases contain metallic fuel data from the Experimental Breeder Reactor-II (EBR-II) and the Fast Flux Test Facility (FFTF). The work demonstrates use of an integrated framework to access EBR-II fuel pin data for evaluating fuel performance models contained within BISON to predict fuel performance of next-generation metallic fuel systems. Between IMIS and FIPD, there is enough information to conduct 1,977 unique EBR-II metallic fuel pin histories from 29 different experiments, and 338 pins from FFTF MFF-3 and MFF-5 with varying levels of details between the two databases. Each of these fuel performance histories includes a high-resolution power history, flux history, coolant channel flow rates, and coolant channel temperatures, and new model developments in BISON since the initial demonstration of this integrated framework. Fission gas release (FGR), cumulative damage fraction, fuel axial swelling, FCCI wastage thickness, cladding profilometry, and burnup were all simulated in BISON and compared to post-irradiation examination (PIE) results to evaluate BISON fuel performance modeling. Implementation of new fuel performance models into a generic BISON input file coupled with IMIS and FIPD yielded results with a better representation of physics than the initial evaluation of the integrated framework. Cladding profilometry, FGR, and fuel axial swelling were found to be in good agreement with PIE measurements for most of the pins simulated. The chosen mechanical contact solver was found to significantly impact the axial fuel swelling and cladding strain predictions when used in conjunction with the U-Pu-Zr hot-pressing model since it bound the fuel to prevent further swelling and increased hydrostatic stresses. This work suggests that fuel performance modeling in BISON under steady-state conditions represents the PIE data well and should be reassessed when new PIE data become available in IMIS and FIPD databases and when improved physical models to better capture fuel performance are added to BISON.

Paaren, Kyle M.↗