Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Scalable molecular dynamics on CPU and GPU architectures with NAMD

NAMD is a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

high-performance computing↗
DOE OSTI2023

Permutationally Invariant Networks for Enhanced Sampling (PINES): Discovery of Multimolecular and Solvent-Inclusive Collective Variables

The typically rugged nature of molecular free energy landscapes can frustrate efficient sampling of the thermodynamically relevant phase space due to the presence of high free energy barriers. Enhanced sampling techniques can improve phase space exploration by accelerating sampling along particular collective variables (CVs). A number of techniques exist for data-driven discovery of CVs parameterizing the important large scale motions of the system. A challenge to CV discovery is learning CVs invariant to symmetries of the molecular system, frequently rigid translation, rigid rotation, and permutational relabeling of identical particles. Of these, permutational invariance has proved a persistent challenge in frustrating the data-driven discovery of multi-molecular CVs in systems of self-assembling particles and solvent-inclusive CVs for solvated systems. In this work, we integrate Permutation Invariant Vector (PIV) featurizations with autoencoding neural networks to learn nonlinear CVs invariant to translation, rotation, and permutation, and perform interleaved rounds of CV discovery and enhanced sampling to iteratively expand sampling of configurational phase space and obtain converged CVs and free energy landscapes. Here, we demonstrate the Permutationally Invariant Network for Enhanced Sampling (PINES) approach in applications to the self-assembly of a 13-atom Argon cluster, association/dissociation of a NaCl ion pair in water, and hydrophobic collapse of a C 45 H 92 n-pentatetracontane polymer chain. We make the approach freely available as a new module within the PLUMED2 enhanced sampling libraries.

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

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Computational modeling predicts the stability of both Pd + and Pd2 + ion-exchanged into H-CHA

Passive NO x adsorbers (PNA) using Pd/zeolites have emerged as a promising solution for the reduction of cold-start emissions from vehicle exhaust. However, the nature of the active sites and the mechanisms underlying NO x adsorption in Pd/zeolites remain a subject of ongoing investigation. In this study, we employ quantum chemical simulations to investigate the structure of Pd species in cation-exchange sites at isolated Al and Al pairs in the 6-ring and 8-ring of the CHA framework before the introduction of NO x . Our calculations show that the speciation of Pd in these exchange sites strongly depends on the precise Al arrangement within the framework, as well as the operating conditions. Ionically dispersed Pd is found to be the most favorable species over a wide range of oxidizing and reducing conditions. Small oligomers of PdO and metallic Pd do not appear to be competitive at either isolated Al or Al pairs. Notably, our calculations show that ion exchange sites other than next–next–nearest neighbor Al pairs in the 6-ring will be preferentially occupied by Pd + instead of Pd 2+ . The stability of Pd + in the zeolite environment is an interesting contrast with its rareness in molecular Pd compounds. Nonetheless, a detailed analysis of the electronic structure shows that predicted Pd oxidation states are consistent with chemical intuition for all complexes investigated in this study. We also discuss the potential ambiguity in Pd characterization provided by typical experimental techniques such as XANES, EXAFS and UV-Vis, and highlight the need for additional EPR spectroscopy studies to further elucidate the initial Pd speciation in zeolites for PNA applications.

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

Classification of Cloud Particle Imagery and Thermodynamics (COCPIT): A New Databasing Tool for the Characterization of Cloud Particle Images Captured During DOE Field Campaigns

The Department of Energy for decades has explored the earth system and atmosphere through research and deployment of in-situ and remote sensing platforms during field campaigns. Among these datasets exists a vast supply of cloud particle images that provide visual insight into the complex microphysics in the clouds that span our globe. The millions of images collected over decades of deployments provides a unique opportunity to further our understanding of our atmosphere down to the crystal size. This work over the past 5 years has sought to organize these images into digestible datasets that can then be used by scientists to further our understanding of microphysics. A machine learning model was developed that categorizes over 1.5 million images across 11 weather events with over 90% accuracy according to particle type. The database was then extended to include dimensional characteristics of the particle as well as co-location of environmental properties, such as temperature and water content. Then, to initialize the connection between these data and our understanding of how crystals form and grow, weather research and forecasting simulations were run to generate the growth histories of the classified crystals. This research culminates with 2 databases per event: (1) a database of all classified crystals and their dimensional and environmental properties and (2) simulated growth histories of each crystal. Finally, a user interface was created to allow researchers to explore data statistics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A statistical model for the product energy distribution in reactions leading to prompt dissociation

Direct dynamics calculations have been performed for three reactions: C 3 H 8 + H → i-C 3 H 7 + H 2 , C 3 H 8 + H → n-C 3 H 7 + H 2 , and C 2 H 3 + O 2 → HCO + CH 2 O. The fraction of the population for the radical products that promptly dissociates is computed. Here the results for C 3 H 8 + H are qualitatively similar to previous results for C 3 H 8 + OH, but the new results exhibit a slightly higher branching fraction for prompt dissociation products, owing to the fact that a greater fraction of the internal energy in the transition state ends up in the radical. For C 2 H 3 + O 2 → HCO + CH 2 O, the fraction of HCO that promptly dissociates is in excess of 99%. Consequently, the main product for C 2 H 3 + O 2 at lower temperatures should be written as H + CO + CH 2 O and not HCO + CH 2 O. These results are then compared with four previous systems: CH 2 O + H → HCO + H 2 , CH 2 O + OH → HCO + H 2 O, C 3 H 8 + OH → i-C 3 H 7 + H 2 O, and C 3 H 8 + OH → n-C 3 H 7 + H 2 O. Based upon these seven system, several statistical models are presented. The goal of these statistical models is to predict the fraction of the transition state energy that ends up in the rovibrationally excited radical. On average, these statistical models provide an excellent prediction of product energy distribution. Consequently, these models can be used instead of costly trajectory simulations for predicting prompt radical dissociation for larger species.

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

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗
DOE OSTI2021

Adaptive variational quantum eigensolvers for highly excited states

Highly excited states of quantum many-body systems are central objects in the study of quantum dynamics and thermalization that challenge classical computational methods due to their volume-law entanglement content. In this work, we explore the potential of variational quantum algorithms to approximate such states. We propose an adaptive variational quantum eigensolver (VQE) for excited states (X) that self-generates a variational ansatz for arbitrary eigenstates of a many-body Hamiltonian H by attempting to minimize the energy variance with respect to H. We benchmark the method by applying it to an Ising spin chain with integrable and nonintegrable regimes, where we calculate various quantities of interest, including the total energy, magnetization density, and entanglement entropy. We also compare the performance of adaptive VQE-X to an adaptive variant of the folded-spectrum method. For both methods, we find a strong dependence of the algorithm's performance on the choice of operator pool used for the adaptive construction of the ansatz. In particular, an operator pool including long-range two-body gates accelerates the convergence of both algorithms in the nonintegrable regime. Here, we also study the scaling of the number of variational parameters with system size, finding that an exponentially large number of parameters may be necessary to approximate individual highly excited states. Nevertheless, we argue that these methods lay a foundation for the use of quantum algorithms to study finite-energy-density properties of many-body systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Identifying Green Solvent Mixtures for Bioproduct Separation Using Bayesian Experimental Design

Liquid–liquid extraction (LLE) is a widely used technique for the separation and purification of liquid-phase products with applications in various industries, including pharmaceuticals, petrochemicals, and renewable chemistry. A critical step in the design of an LLE process is the selection of appropriate solvents. This study presents a new methodology for identifying solvent mixtures for bioproduct separation using Bayesian experimental design (BED). Motivated by the need for environmentally friendly and effective separation methods, we address the challenge of selecting solvent systems that balance separation efficiency, selectivity, and environmental impact while also tackling the difficulty of separating multiple bioproducts using complex solvent systems. Our approach specifically seeks to predict product partition coefficients (log10 Kp values) as thermodynamic parameters underlying solvent selection. The iterative approach integrates Bayesian optimization with experimental measurements to guide solvent selection and leverages COSMO-RS simulations to enhance high-throughput experimentation. Using the design of solvent systems for the separation of lignin-derived aromatic products via centrifugal partition chromatography (CPC) as a case study, we show that within seven iterations/cycles of the methodology, we can identify new mixtures of green solvents that align with CPC design principles. Furthermore, these results demonstrate the efficacy of the BED framework in optimizing green solvent systems for complex separations, highlighting the potential of this method to advance the field of green chemistry and contribute to the development of sustainable industrial processes.

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

Evaluation of semi-empirical soot models for nonpremixed flames with increased stoichiometric mixture fraction and strain

Due to the size and complexity of industrial-scale systems and computing limitations, semi-empirical soot models are often employed in CFD simulations rather than computationally-expensive detailed models. Many of these models were developed for specific applications with unique characteristic timescales and/or validated only under fuel-air combustion conditions. Hence, their use in different contexts, such as oxy-combustion, could lead to inaccurate predictions. In this study, twelve semi-empirical models (1-step or 2-step) are evaluated on their ability to respond to changes in stoichiometric mixture fraction (Z st ) and strain in a series of ethylene counterflow flames which span the experimental sooting-to-non-sooting (yellow to blue) transition. A unique approach of plotting soot formation rate in normalized, local equivalence ratio space is introduced to aid in the analysis. Results show that no existing model is able to predict a blue flame when Z st is increased beyond the experimentally-measured sooting limit. Many models give the counter-intuitive result of increasing peak and/or integrated soot volume fraction (svf) as Z st increased, contrary to experimental observations. Two models are able to successfully predict a blue flame upon increase of strain to the sooting limit value. For two-step models, the strong dependence of growth rate on surface area results in high sensitivity of svf to changing flame boundary conditions, but with mixed accuracy. There is still a significant need for a robust semi-empirical model which can accurately predict soot fraction in systems with oxygen-enrichment or variable strain.

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

Influence of NOx chemistry on the prediction of natural gas end-gas autoignition in CFD engine simulations

Natural gas (NG) represents a promising low-cost/low-emission alternative to diesel fuel when used in high-efficiency internal combustion engines. Advanced combustion strategies utilizing high EGR rates and controlled end-gas autoignition can be implemented with NG to achieve diesel-like efficiencies; however, to support the design of these next-generation NG ICEs, computational tools, including single- and multi-dimensional simulation packages will need to account for the complex chemistry that can occur between the reactive species found in EGR (including NOx) and the fuel. Research has shown that NOx plays an important role in the promotion/inhibition of large hydrocarbon autoignition and when accounted for in CFD engine simulations, can significantly improve the prediction of end-gas autoignition for these fuels. However, reduced NOx-enabled NG mechanisms for use in CFD engine simulations are lacking, and as a result, the influence of NOx chemistry on NG engine operation remains unknown. Here, we analyze the effects of NOx chemistry on the prediction of NG/oxidizer/EGR autoignition and generate a reduced mechanism of a suitable size to be used in engine simulations. Results indicate that NG ignition is sensitive to NOx chemistry, where it was observed that the addition of EGR, which included NOx, promoted NG autoignition. The modified mechanism captured well all trends and closely matched experimentally measured ignition delay times for a wide range of EGR rates and NG compositions. Here, the importance of C2-C3 chemistry is noted, especially for wet NG compositions containing high fractions of ethane and propane. Finally, when utilized in CFD simulations of a Cooperative Fuels Research (CFR) engine, the new reduced mechanism was able to predict the knock onset crank angle (KOCA) to within one crank angle degree of experimental data, a significant improvement compared to previous simulations without NOx chemistry.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Steric Mapping, Ligand Dynamics, and Cycloisomerization Catalysis with Redox Robust Mn I/0/-I Dicarbenes

Manganese is perhaps the most electronically versatile element, yet the redox properties, reactivity, and catalytic applications of low-valent Mn 0 /Mn −I complexes remain underexplored due to the propensity for Mn 0 to dimerize, quenching highenergy metalloradicals. We report a series of redox-active monometallic Mn I , Mn 0 and Mn −I complexes containing a BH 2 - bridged dicarbene, characterized using a suite of experimental and cutting-edge computational (DFT) methods. Slow electron transfer kinetics at Mn I/0 are observed, with computations and electrochemical simulations in excellent agreement with experimental values. Despite the lack of steric bulk at the BH 2 -bridged Mn 0 , the t Bu groups at the dicarbene provide adequate steric protection to prevent dimerization, with percent buried volume (% V bur ) serving as a valuable steric ranking tool. We also show that a %V bur > 83% prevents dimerization for a diverse array of Mn 0 complexes from the literature. Ligand sterics of BPh 2 -and BH 2 -bridged complexes dictate reaction outcomes when Mn I and Mn −I are exposed to nucleophiles and electrophiles, respectively, while Mn0 facilitates the radical cycloisomerization catalysis of 6-iodo-1- hexene at room temperature. Furthermore, this work underscores the importance of ligand sterics in rationalizing reactivity patterns at Mn and provides valuable insights for designing chelating ligands that can selectively leverage Mn I/0/‑I states in redox-mediated catalytic reactions.

Ligands↗
DOE OSTI2021

Experimental and computational investigations of ethane and ethylene kinetics with copper oxide particles for Chemical Looping Combustion

In this work, reaction pathways for the oxidation of methane, ethane, and ethylene with CuO was obtained by ReaxFF Molecular Dynamics (MD) simulations between temperatures of 1000 K and 2000 K. Experiments in a fixed-bed flow reactor were preformed with methane, ethane, and ethylene at temperatures ranging from 500 K to 1000 K with time-dependent species measurements from an Electron-Ionization Molecular Beam Mass Spectrometer (MBMS), and species validation with Gas Chromatography (GC) for detection of complete and intermediate combustion products. The MBMS and GC allow for the detection of oxygenated species and larger species produced from radical reformation. The simulation and experiment agree on the production of such species as CH 3 CHO, CH 2 O, CO, and H 2 O, which allow for the creation of simple C1 and C2 reaction pathways, which can be used in kinetic models of C2 species and larger fuels such as biofuels, which inherently depend on C1 and C2 kinetics and reaction pathway. The simulation and experiment disagree on the formation of C 2 H 2 , CH 3 OH, and CO 2 with large amounts of C 2 H 2 being measured in the ethylene oxidation simulations and CH 3 OH being formed in methane oxidation simulations, while neither species were experimentally found. In the case of CO 2 large amounts of CO 2 are rapidly produced in experiments with C2 fuels at 800 K, while little-to-no CO 2 was observed in simulations. This is believed to be resulting from the extremely short timescale of the simulations, preventing total oxidation of the fuel. Here, the differences in products produced between simulation and experiment allow for the potential to modify the ReaxFF potential functions to more accurately model the experimental products of Cu–H–O–C reaction kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A New Nitrogen-Rich Energetic Material with So Many Tautomers

A new nitrogen-rich triazolo-triazole compound, (3-(6-methyl-1H-[1,2,4]triazolo[4,3-b][1,2,4]triazol-3-yl)-1H-1,2,4-triazol-5-amine), TTT1, has been prepared, and its acid–base and tautomeric behavior has been investigated. In the pH range of 0.3–12, TTT1 can accept up to two protons, forming a monocation and a dication, and can deliver one proton, forming a monoanion. The tautomeric behavior is particularly rich for the monocation, for which computational analysis predicts four different tautomers in a narrow energy range of 2 kcal/mol. Two of these tautomers (2H-7H-8H and 3H-7H-8H) have been isolated in salts of the monocation with suitable counterions (chloride, bromide, perchlorate). Surprisingly, the most stable predicted tautomer, 1H-3H-7H, has not been found in the four crystallized salts of the monocation. The energetic perchlorate salt of the monocation (3H-7H-8H tautomer) shows good thermal stability and good stability to impact, friction, and electric discharge. The packing of this compound shows the formation of H-bonded dimers with interactions between N8–H···N1. The crystal structure of this energetic salt was studied experimentally up to 2.8 GPa; no phase change or decomposition was observed.

Chemical structure↗
DOE OSTI2025

Revealing the Origin and Nature of the Buried Metal‐Substrate Interface Layer in Ta/Sapphire Superconducting Films

Abstract Despite constituting a smaller fraction of the qubit's electromagnetic mode, surfaces and interfaces can exert significant influence as sources of high‐loss tangents, which brings forward the need to reveal properties of these extended defects and identify routes to their control. Here, we examine the structure and composition of the metal‐substrate interfacial layer that exists in Ta/sapphire‐based superconducting films. Synchrotron‐based X‐ray reflectivity measurements of Ta films, commonly used in these qubits, reveal an unexplored interface layer at the metal‐substrate interface. Scanning transmission electron microscopy and core‐level electron energy loss spectroscopy identified an intermixing layer (≈0.65 ± 0.05 nm) at the metal‐substrate interface containing Al, O, and Ta atoms. Density functional theory modeling reveals that the structure and properties of the Ta/sapphire heterojunctions are determined by the oxygen content on the sapphire surface prior to Ta deposition for two atomic terminations of sapphire. Using a multimodal approach, we gained deeper insights into the interface layer between the metal and substrate, which suggests that the orientation of deposited Ta films depend on the surface termination of sapphire. The observed elemental intermixing at the metal‐substrate interface influences the thermodynamic stability and electronic behavior of the film, which may also affect qubit performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational Fermi Level Engineering and Doping-Type Conversion of Ga2O3 Via Three-Step Processing

Ga2O3 is being actively explored for power electronics, deep-ultraviolet optoelectronics, and other applications due to its ultra-wide bandgap and low projected fabrication cost of large-size and high-quality crystals. N-type doping of Ga2O3 can be achieved and tuned, but p-type doping faces fundamental obstacles due to deep character of acceptor levels and polaron transport of resulting holes. However, successful engineering of Ga2O3 based devices requires critical control of doping density, Fermi level position, and free carrier concentration, providing opportunities for predictive process simulation. We use first-principles defect theory and defect equilibrium calculations to simulate a 3-step growth-annealing-quench protocol for hydrogen assisted Mg doping in Ga2O3, taking into account the hydrogen-oxygen-water gas phase equilibrium. We predict type conversion to a net p-type regime following O-rich annealing after growth under reducing conditions in the presence of H2. This process is similar to the Mg acceptor activation by H removal in GaN. We show that there is an optimal temperature that maximizes the net acceptor density during the equilibrium annealing step for a given Mg doping level. Quenching of non-equilibrium annealed samples then results in a Fermi level EF below mid-gap down to about EV +1.5 eV, creating a significant number of uncompensated neutral MgGa0 acceptors. The resulting free hole concentration in quenched samples is very low (109 – 1012 cm-3) due to deep energy level of these Mg acceptors, but this type converted Ga2O3 material can create a significant built-in field in a p-n junction with an adjoining n-type material. Additionally, the electron concentration is greatly suppressed in such acceptor-doped Ga2O3, which could enable the use as a current blocking layer to fabricate normally-off (enhancement-mode) vertical Ga2O3 based metal-oxide-semiconductor field effect transistors (MOSFETs), and as a guard ring for edge termination in MOSFETs and Schottky barrier diodes (SBDs) with increased breakdown voltage.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗
DOE OSTI2020

Phase Equilibria and Thermochemistry of Advanced Fuels: Modeling Burnup Behavior (Final Report)

Achieving the goal of developing advanced fuel concepts that meet the DOE objectives of being robust, demonstrating high performance, and are more tolerant of accident conditions than current fuel systems will require a thorough understanding of the thermophysical and thermochemical properties of the constituent materials. Non-oxide fuel systems are being explored under the Advanced Fuels Program that hold significant promise for improved performance and accident tolerance, including the uranium silicide-based system considered in the current work. Prospective cladding materials currently considered that contribute to improved accident tolerance include silicon carbide composites and ferritic alloys (Fe-Cr-Al base compositions). Thus, the effort developed thermochemical models and values, supported with targeted experiments, to evaluate the ferritic alloy and silicon carbide composite cladding systems in contrast to current zirconium alloy cladding. The developed detailed understanding will serve to aid in relatively early screening of candidate systems to avoid wasted effort, guide development of new fuel forms, and to provide a basis for predicting and modeling fuel performance. Major deliverables for the project included: Thermochemical assessment and models of phases in the U-Si and U-Si-N systems; thermochemical evaluation supported by experimental measurements of fuel-cladding interactions of silicide fuel with baseline zirconium, ceramic composite, and ferritic alloy cladding; thermochemical assessment and models of phases supported by experimental measurements for silicide fuel with key fission products provided in a dataset and reported in refereed publications. Within the project a significantly refined U-Si phase diagram was developed and reported that now includes homogeneity ranges for key phases, such at the U 3 Si 2 proposed fuel phase, and settles issues with regard to uncertainty in the stability of some phases. Computational efforts together with key experiments has determined phase formation in interactions between U 3 Si 2 and Zircaloy-4 cladding material, a ferritic FeCrAlY alloy of interest as an advanced cladding material, and silicon carbide, also of interest as a fiber-reinforced composite cladding. As expected, very significant reactions occur between U3Si2 and Zircaloy-4, with much less interaction at higher temperatures for the ferritic alloy, and finally interactions with SiC only in the region of contact. A potentially major issue is the stability of U 3 Si 2 fuel that has undergone significant burnup. The result is the loss of the uranium metal, liberating silicon, and the formation of concomitant fission product elements that either dissolve in the U 3 Si 2 phase or form independent, and possibly silicide phases. A combination of experimental determinations of phase formation of U 3 Si 2 reacted with representative fission products yttrium, gadolinium, cerium, zirconium, and molybdenum and first principles calculations has helped understand the fuel chemistry. The behavior of the U 3 Si 2 phase and the partitioning of silicon to possible fission product phases with burnup was thus determined, with significant dissolution in U 3 Si 2 of cerium, gadolinium, zirconium, and plutonium predicted along with independent phase formation of a U-Mo-Si ternary phase, yttrium silicide, and elemental selenium. It can be concluded that at significant burnup there will be a very minor amount of the U 3 Si 2 fuel phase that will decompose to a lower silicide or a uranium alloy as silicon preferentially forms a secondary phase.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Using X-ray radiography to study oxygen flow in a proton exchange membrane electrolyzer operating under balanced pressure conditions

Of the various water electrolyzer technologies, the proton exchange membrane electrolyzer (PEMWE) is one of the best solutions for producing clean hydrogen without releasing CO 2 . In order to allow for widespread use of clean hydrogen, it is necessary to decrease its cost, which is intrinsically related to system operation. Current PEMWE plants operate in differential mode, directly pressurizing hydrogen and benefiting from thermodynamic compression, which increases overall system efficiency. However, high differential pressure above 30 bar can cause membrane stress, resulting in membrane creeping and failure. Pressurizing the water and operating at balanced pressure allows hydrogen to be produced at higher pressures while preserving the integrity of the membrane and porous layers. Nevertheless, the impact of pressurizing water on PEMWE performance must be better understood to maximize performance under balanced pressure conditions. Here, this study examined the impact of water pressure on two-phase flow. A high-pressure electrolyzer setup was developed to perform operando X-ray radiography and examine oxygen transport with high temporal resolution. The imaging segmentation process, developed to capture bubble properties in the channel, was applied to a specific experiment. The results clearly showed that as pressure increased up to 30 bars, the initial bubbly flow transitioned to slug flow, which led to channel saturation with oxygen. This work demonstrates that two-phase flow in an electrolyzer can be studied using X-ray radiography, which has the advantages of fast measurements and the ability to probe dense materials, such as those required for pressurized electrolyzers.

Balanced high pressure operation↗