Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compositional reasoning”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Compositional Reasoning for Hierarchical State Machines

Harel statecharts and its derivatives are popular graphical languages for specifying discrete control systems via hierarchical state machines. Separately, there has been a long line of work on specifying concurrent systems with process calculi which come equipped with an algebraic theory, the ability reason compositionally about various temporal properties, and strong type systems. While these two approaches to modeling systems are tantalizingly similar, the integrated reasoning principles that exist for process calculi have not been demonstrated in hierarchical state machines. A key issue is that operational theories for process calculi do not behave like control systems, and thus, there is virtually no tool support for modeling control systems with such languages. For a control system designer, bringing the integrated, more scalable reasoning from the process calculi to state-machine languages would enable the specification of more complex systems and a more modular systems development process. Our insight is that we can recover many important results from the process calculi in hierarchical state machines with local scope. We employ a structural operational semantics, which is ubiquitous in process and 𝜆-calculi but uncommon in hierarchical statemachine formalizations, to enable inductive reasoning about behavior. Taking inspiration from the structure of process calculi metatheories, we define a calculus of refinement and equivalence that we prove sound with respect to local notion of (bi)simulation. Furthermore, we prove that the calculus preserves the behavioral properties of reactivity, observational determinism, traces, and linear temporal properties. Our results are mechanized in the Rocq proof assistant.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Boron Coordination in Multicomponent Glasses: Analytical Models and Machine Learning With Uncertainty

Borosilicate glasses are extensively used in a variety of applications from kitchenware to nuclear waste immobilization due to the strong network formed by the Si-O-B bond that makes it resistant to chemical corrosion and gives it a low thermal expansion. Boron, however, exists in both trigonal BO3 and tetrahedral BO4 bonds in glass systems, which impacts the chemical durability and thermal resistance of the glass, amongst other properties. Boron coordination (N4), or the ratio of the amount of BO4 to BO3 within a glass, may aid in predicting these properties but is difficult to derive without experimental data due to the complexity of impacts from varied glass compositions and processing factors. For this reason, compositional models have been developed to predict boron coordination, but the models typically include a limited number of glass components. To help fill this gap in the models, in this work, a diverse multicomponent glass dataset of 809 glasses is compiled from a literature search, and then a number of analytical and machine learning (ML) models are trained on the dataset. Previously developed modified Bernstein and modified Du Stebbins analytical models were fitted to update parameters with the new dataset. Then, partially Bayesian neural networks, Gaussian process regressor, and heteroskedastic deterministic neural networks were evaluated. The ML models examined all have different strategies to overcome the potential for overfitting as a result of a limited training dataset, and return results that account for model uncertainty, which can be valuable for understanding model reliability. For the first time, cooling rate is introduced as an input parameter for ML models, showing consistent improvements in performance and solidifying the importance of including parameters outside of composition alone for N4 prediction. The machine learning models examined here show promise in accurate predictions of boron coordination in borosilicate glasses, all achieving R2 values of 0.91.

boron coordination↗
DOE OSTI2025

Evaluating High-Halide Waste Form Options for Salt-Based Nuclear Waste Simulants

In this study, waste forms were being explored for an electrochemical salt simulant, , referred to as ERV3, which is a high-LiCl/KCl salt containing simulated fission products (i.e., Sr, Cs, Nd) and Na to represent bond sodium from Experimental Breeder Reactor-II metallic fast reactor fuel. The goal was to find glassy systems that could be used to immobilize the salt in a single-step process. The envisionment of this process would be find a frit glass that could be added to the salt waste, heat treated, poured into waste canisters, and then stored for disposal. To perform this study, a literature review was conducted, the most promising seven systems were fabricated without the salt, and then mixed with salt simulant and heat treated under different processes. The criteria that were used to screen potential compositions included demonstrated alkali incorporation, could be melted at reasonably low temperature ( T ≤ 1000°C), and if the compositions had some demonstrated data for waste-form-related properties, that was a benefit. High marks were given for compositions that showed amorphous nature after slow cooling of samples containing ERV3.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

High-Speed Layup and Forming of Automotive Composite Components

This Project is focused on the design and manufacture of automotive components that meet functional and environmental requirements of an existing automotive application at a cost of ≤ $\$$11.00 per kilogram weight reduction. This project fosters the development of composite material technologies suitable for high volume automotive processes and run rates as well as industry workforce development with these technologies. Current automotive manufacturing involves utilizing steel or aluminum in sheet form which is rapidly stamped into components at rates up to 3600 per hour. The metallic sheets are available in many different thicknesses, strength levels, and manufacturing rates are reasonable independent of part size. While composite materials are available for use in automotive applications, the material cost, labor to manufacture and the processing of the waste far exceed the cost compared to metallic designs. Typical composite layer by layer layup procedures don’t meet the desired 60 second layup time that current automotive processes require and are also restricted by part size. Due to these factors, composites have not yet made advances into today’s high volume automotive applications. Industry partners DURA, BASF, Ford, and IACMI core innovation partner MSU collaborated to develop a manufacturing process technology that is capable of manufacturing composite blanks at high volume and independent of part size. IACMI core innovation partner Purdue provided FEA analysis and cost modelling. The objective of this project was to demonstrate a composite sheet layup and consolidation process that can be commercialized for high volume requirements, identify potential layup equipment suppliers, and develop a process of 60 second layup, forming, and trimming of a continuous fiber automotive component for the mainstream market.

42 ENGINEERING↗
DOE OSTI2025

Multimuons in cosmic-ray events as seen in ALICE at the LHC

ALICE is a large experiment at the CERN Large Hadron Collider. Located 52 meters underground, its detectors are suitable to measure muons produced by cosmic-ray interactions in the atmosphere. In this paper, the studies of the cosmic muons registered by ALICE during Run 2 (2015–2018) are described. The analysis is limited to multimuon events defined as events with more than four detected muons (N μ > 4) and in the zenith angle range 0° < θ < 50°. The results are compared with Monte Carlo simulations using three of the main hadronic interaction models describing the air shower development in the atmosphere: QGSJET-II-04, EPOS-LHC, and SIBYLL 2.3d. The interval of the primary cosmic-ray energy involved in the measured muon multiplicity distribution is about 4 × 10 15 < E prim < 6 × 10 16 eV. In this interval none of the three models is able to describe precisely the trend of the composition of cosmic rays as the energy increases. However, QGSJET-II-04 is found to be the only model capable of reproducing reasonably well the muon multiplicity distribution, assuming a heavy composition of the primary cosmic rays over the whole energy range, while SIBYLL 2.3d and EPOS-LHC underpredict the number of muons in a large interval of multiplicity by more than 20% and 30%, respectively. The rate of high muon multiplicity events (N μ > 100) obtained with QGSJET-II-04 and SIBYLL 2.3d is compatible with the data, while EPOS-LHC produces a significantly lower rate (55% of the measured rate). For both QGSJET-II-04 and SIBYLL 2.3d, the rate is close to the data when the composition is assumed to be dominated by heavy elements, an outcome compatible with the average energy E prim ∼ 10 17 eV of these events. This result places significant constraints on more exotic production mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Electrolyte and Cutoff Potential Effects on Cycle Life of Li4Ti5O12/LiNi0.9Mn0.1O2 Batteries for Behind-the-Meter Storage Applications

Behind-the-Meter Storage (BTMS) is a stationary battery energy storage system that is connected to the electrical distribution system on the customer's side of the utility's service meter. BTMS systems are used to store electrical energy from the grid as well as inconstant, renewable energy, such as local solar and wind generation. A successful BTMS system will allow the customer to pair their energy generation and storage to optimize electrical consumption from the grid, improving reliability and minimizing cost. For BTMS applications, batteries must be designed and optimized with different set of criteria from other leading segments of the Li-ion battery market, like transportation, due the system being stationary and proximal to the residential or commercial building it's benefitting. BTMS applications prioritize safety, cost (low/no-critical materials), reliability (20-year calendar life), and durability (10,000 cycle life), while having the ability to (minimally) compromise energy density and rate capability. Lithium titanate (Li4Ti5O12-, LTO) is a promising anode candidate for BTMS applications due to its high safety and capacity retention, while maintaining a reasonable 160 mAhg-1 reversable capacity and composition of relatively abundant materials. (1) Specifically, LTO has a high working voltage which helps to prevent Li dendrite formation, improving safety. Furthermore, LTO also has negligible lithiation-based volume change, leading to less mechanical pulverization, or loss of active material, upon cycling. For the cathode, materials with little or no Co are of high interest due to the high cost and low abundance of Co. LiMn2O4 (LMO) has been paired with LTO for BTMS applications in the past due to its safety, low cost (abundancy), and reasonably high operating voltage. (2-4) However, the low capacity of LMO limits energy density and specific energy. While not the highest priority for BTMS applications, increasing energy density will enable deployment in space constrained BTMS applications and decrease total cost. LiNi0.9Mn0.1O2 (LN-MO) is a recently developed material with promise due to its high operating voltage and relatively low price. (5) However, Ni-rich layered oxides, including LNMO, tend to struggle with capacity retention during high-voltage cycling due to mechanical pulverization, irreversible phase transitions, and unstable solid-electrolyte interphase. The study presented here focuses on building an understanding of how electrolyte solvent and varied cutoff potentials will impact the cycle life of LTO/LN-MO cells. Specifically, a comparison is provided between ethylene carbonate (EC), ethyl methyl carbonate (EMC), fluoroethylene carbonate (FEC), and Gen2 electrolyte solvents with 1M Lithium hexafluorophosphate (LiPF6) salt, cycling to two upper termination potentials, 2.6V and 2.7V. Electrochemical testing and diagnostics (e.g., differential capacity analysis, area specific impedance, constant voltage hold, and rate capability) and post-mortem characterization will be used to understand the aging behavior and failure mechanisms of the 8 cell combinations (four electrolytes and two voltage cutoffs). Cells with FEC electrolyte showed a lower initial capacity compared to cells with Gen2, EMC, and EC cycling at both voltages; however, the cells with FEC showed consistent trends in capacity retention with 2.6V and 2.7V termination potentials, while the cells with the other electrolytes showed much higher rates of capacity loss when cycling to the higher voltage. These results indicate that FEC may play a role in improving durability of high-voltage, Ni-rich electrode systems for use in high-cycle applications, such as BTMS.

electrolyte↗
DOE OSTI2025

NitroNet: Smart System to Quantify Nitrous Oxide Emissions

Agricultural croplands are the largest anthropogenic source of nitrous oxide (N 2 O), the third most important greenhouse gas. Emissions are characterized by “hot spots” and “hot moments”, meaning emissions are highly heterogeneous in space and time. Emissions are driven by nitrification and denitrification processes which depend upon complex factors such as soil characteristics (type, compaction, pH, moisture, topography), management practices (fertilizer type and application, tillage, crop type, field history, irrigation), biogeochemistry (organic carbon, microbial composition), and meteorology (precipitation, temperature, and wind). For these reasons, quantifying cropland emissions by either measurements or modeling is extremely challenging with large uncertainties.

09 BIOMASS FUELS↗
DOE OSTI2026

Two-stage formation-energy correction (NbZr, TaZr, VZr)

This bundle contains the scripts, the raw and corrected per-structure data, and the manuscript plots for the NbZr / TaZr / VZr BCC binary formation energies and the associated RMSDs. Why a two-stage correction is necessary: The "raw" formation energy of every relaxed VASP configuration is computed in the usual way, FE_raw(c) = E_alloy(c) - sum_i x_i * E_pure_i , where E_pure_i are the per-atom total energies of the pure-element reference structures (Nb, Ta, V, Zr in the same BCC supercell, with identical INCAR / KPOINTS / PAW choices). With perfectly consistent reference runs the raw FE should vanish at the two pure-element endpoints (x = 0 and x = 1) by construction. In practice this does not hold for two reasons that are present in our dataset: 1. Reference-energy inconsistency (composition-dependent bias). Even with identical input parameters, the pure-element runs (stored in `corrected_DFT_pure_element_runs/`) differ slightly from the values that would be implied by the alloy runs at near-pure compositions (a few meV/atom). This bias is approximately linear in concentration, because the residual error in E_pure_Nb (or E_pure_Ta / E_pure_V) propagates into FE_raw(c) as (1 - x) * dE_pure_1, and the corresponding error in E_pure_Zr propagates as x * dE_pure_2. Left uncorrected, this produces a non-physical "tilt" of FE_raw(x) and shifts the entire FE-vs-x cloud away from zero at the endpoints. 2. Endpoint anchoring against the audited true endpoints. The strict endpoint values (FE_x0_meVatom, FE_x1_meVatom in `corrected_fe_strict_endpoints_20260518/strict_endpoint_check_20260518.csv`) were re-derived from an independent cross-check of the pure-element runs. After stage 1 removes the linear bias, the near-pure compositions in the alloy dataset still extrapolate to values that differ slightly from these audited endpoints — because stage 1 is fit from a few near-end alloy bins, not from the audited pure-element references themselves. The README.txt file discusses how these issues are addressed by the two-stage correction, and describes folder layout, pipeline summary, and how to re-run.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulating Marine Ecosystem Dynamics and Biogeochemical Cycling With Multiple Plankton Functional Types

Current representations of marine ecosystems in Earth System Models are greatly simplified, neglecting key interactions between dynamic food webs, biogeochemistry, and climate change. We use the Marine Biogeochemistry Library code base within the Community Earth System Model 2.2.2 to create an expanded ecosystem model with eight phytoplankton groups and four zooplankton size classes (MARBL-8P4Z). Incorporating more specific plankton types and size classes has the potential to capture a wider range of possible behaviors of the ecosystem, its complex interactions with biogeochemistry, and its feedback to climate change. It also permits stronger observational constraints, including in situ group-specific biomass and various observational estimates of plankton community composition. MARBL-8P4Z broadly captures observed global-scale patterns in biomass and community composition for both phytoplankton and zooplankton, with a good performance in simulating broad biogeochemistry fields. The model shows comparable spatial patterns and magnitudes to the observed picophytoplankton biomass (Prochlorococcus, Synechococcus, picoeukaryotes), and captures the seasonal cycle of mesozooplankton biomass. Picophytoplankton groups and microzooplankton dominate biomass and production in oligotrophic, subtropical regions, while nano-phytoplankton, diatoms and the larger zooplankton groups prevail at higher latitudes and within upwelling zones. The model simulates reasonable energy transfer efficiency through the food web, with tight linkages between the phytoplankton community composition, zooplankton grazing, and carbon export, with the potential to link to fisheries models. Thus, MARBL-8P4Z has the potential to account for key climate-driven ecological shifts in the plankton that will modify ocean biogeochemistry in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Design of lightweight BCC multi-principal element alloys with enhanced hydrogen storage using a machine learning-driven genetic algorithm

Body-centered cubic (BCC) based multi-principal element alloy (MPEA) hydrides have demonstrated significant potential for compact and efficient hydrogen storage. In this work, we first leverage machine learning (ML) models to predict the hydrogen affinity, storage capacity and phase stability of BCC MPEAs, creating a unique hydrogen-to-metal (H/M) predictor for materials with unprecedented performance. We developed a metaheuristic optimizer high-throughput framework by interfacing ML models with a genetic algorithm for the accelerated search of {Mg, Al, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Nb, Mo} based lightweight BCC MPEAs with improved hydrogen storage characteristics. We report five new MPEAs with a predicted gravimetric hydrogen storage capacity of around 3.5 wt% or more, including Cr 0.09 Mg 0.73 Ti 0.18 (4.25 wt% H) and Cr 0.21 Nb 0.11 Ti 0.35 V 0.33 (3.5 wt% H). The electronic structure of the top-performing composition, Cr 0.09 Mg 0.73 Ti 0.18 , was analyzed using density functional theory (DFT) to understand the reasons for its improved hydrogen storage properties compared to TiFe (1.90 wt% H), LaNi 5 (1.37 wt% H) or BCC MPEAs like TiVNbCr (3.70 wt% H). Temperature-dependent molecular dynamics (MD) studies were further performed on optimized BCC MPEAs to qualitatively study hydrogen mobility and analyze the effect of different elemental composition on bulk hydrogen diffusion. Our findings demonstrate how a ML assisted genetic algorithm framework can be used for efficient search of stable, lightweight and cost-effective MPEAs while minimizing the need for expensive ab initio calculations.

DFT↗
DOE OSTI2025

Development of a two-zone zero-dimensional mixing model for actively fueled prechamber engines

In this numerical work, a zero-dimensional, two-zone mixing model was developed and validated to capture the mixture preparation process of an actively fueled prechamber engine. The model was developed for ethanol-gasoline fuel blends but is applicable to various fuels of interest relevant to prechamber applications. The model was exercised to demonstrate the effects of boundary conditions, fuel composition, and residual content on the mixture stoichiometry predictions. The 0D model showed reasonable agreement in prechamber equivalence ratio predictions with higher fidelity computational fluid dynamics modeling, but the accuracy of model was constrained to the governing assumptions of the injection process. The modeling assumptions surrounding the non-dimensional injection of mass into the prechamber were scrutinized to identify key areas of model improvement. The results indicate heat transfer effects play a crucial role in the evaporation of the fuel within the confined auxiliary chamber, especially when considering fuels with high latent heat of vaporization. Additionally, the presence of fuel films that collect on the prechamber surfaces complicated the modeled evaporative cooling loses as the modeling suggests the films evaporate primarily by means of conduction from the boundary. Furthermore, the developed model serves as a simplified framework to generate calibration maps of prechamber fueling or run real time as an open loop control strategy such that desired prechamber mixture proportions are effectively targeted.

Zero-dimensional↗
DOE OSTI2024

A FeCrAl-Al2O3 Composite Produced via Laser Powder Bed Fusion of a Mixed Powder for Porous Catalyst Scaffolds

This study proposes a novel approach for synthesizing and etching bicontinuous FeCrAl-Al2O3 composites as a means for replacing FeCrAl foams as catalyst scaffolds in bio-driven alcohol reactors for jet-fuel production. Conventional FeCrAl foams suffer from poor availability and consequent high costs. New additive manufacturing techniques provide an opportunity to produce tailored foams in reasonable times and at acceptable costs. This research aimed to generate a porous FeCrAl structure by etching a bicontinuous FeCrAl-Al2O3 composite produced by laser powder bed fusion of amalgamated FeCrAl and Al2O3 powders. The composite powder for laser powder bed fusion is created by ball-milling FeCrAl and Al2O3 powders. This research focuses on achieving a bi-continuous FeCrAl-Al2O3 structure, essential for the selective removal of the ceramic phase. The influence of laser processing parameters on the microstructure was examined across a range of laser powers (60-120 W) and scan speeds (100-400 mm/s), showing that higher powers and speeds produce finer metal struts. A bi-continuous microstructure was consistently obtained, marking a key achievement. The Al2O3 removal process involved a two-step etching method using hydrochloric and phosphoric acids, tested across various etching times. The alumina phase was reduced from 36 vol% to 17 vol% (corresponding to an increase in porosity from 24 vol% to 43 vol%), showing the potential for use as a porous catalyst scaffold. This research demonstrates the potential for using additive manufacturing to produce porous FeCrAl structures capable of replacing hard-to-source FeCrAl foams.

Son, Kwangtae↗
DOE OSTI2024

Hydrogen sensing in different hydrogen-carrying gases using composites of PdSnO 2 and halloysite nanotubes

Here, in this work, we developed composites of palladium-decorated tin dioxide (PdSnO 2 ) and halloysite nanotubes (HNTs) by adding different amounts of HNTs as additives into PdSnO 2 , and investigated how the different amounts of HNTs in the composites and the different hydrogen (H 2 ) carrying gases (air and helium (He)) could affect the H 2 sensing performance of such composites (PdSnO 2 -HNT). Through the sensing-performance characterization using H 2 carried by air, it was found that PdSnO 2 with an appropriate small amount of HNTs (i.e., 2 % of PdSnO 2 mass) can improve the sensing performance with respect to limit of detection (LOD) and response/recovery time. Further, the optimal PdSnO2-HNT composite as a sensor was tested to detect H 2 in He. The testing results indicated that the composite can detect H 2 in He, but its performance parameters (i.e., the profile of calibration curve, LOD, and response time) are different from those of such a sensor in detecting H 2 in air. Moreover, the composite still presented better sensing performance than PdSnO 2 without HNTs in detecting H 2 in He. Possible reasons for the effects of HNT and H 2 -carrying gas on the sensing performance of PdSnO 2 -HNT-based sensors were discussed. We believe that this study provides valuable insights into the functionality and the adaptability of PdSnO 2 -HNT-based H 2 sensors in diverse operational conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Model of Hydrogen Solubility in Palladium-Silver Alloys

Fusion fuel cycle designs depend on palladium and palladium-silver alloys. The selective permeability of these is used to separate hydrogen isotopes from other elements. Predicting the solubility of hydrogen isotopes in palladium-silver alloys is important for the design of key unit operations for tritium processing. This paper presents a model of hydrogen isotope solubility in palladium and palladium-silver alloys based on Guggenheim’s quasichemical lattice theory. The model parameters were determined as functions of isotope molecular weight and the silver content of the alloy. The model predicts reasonable values for the solubility of protium, deuterium, and tritium over a wide range of temperatures, pressures, and alloy compositions.

08 HYDROGEN↗
DOE OSTI2025

Development of fiber orientation in injection molding: Comparison of glass fiber, carbon fiber and their hybrid composites

Fiber orientation distribution (FOD) in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for glass fiber (GF), carbon fiber (CF), and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varied with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Theorems in Service of Sound Composition, Rapid Modeling and Scalable Analysis

This project extends the state of the art in formal verification modeling with modules and automatically checkable data-sharing patterns such that component modules can retain their assurance case when composed within a larger system. For users, smaller models make reasoning easier and help to ensure they accurately reflect text specifications. For automated methods, smaller models give exponential benefits for verification algorithm execution time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Analysis of the Surface Morphology and Chemical Composition of Zr-Nb3Sn Alloys with Different Zr Concentrations

In the field of superconducting science, Nb has been the material of choice for Radio Frequency (RF) Cavities for years. However, in recent years Nb3Sn has shown to be a promising alternative to Nb. Its combination of a significantly high superconducting transition temperature and a much higher value of upper critical magnetic field (Hc2), has been the reason the alloy has been so widely studied. This project aims to develop materials to further improve the high magnetic field performance of Nb3Sn by improving its upper critical magnetic field and critical current density (Jc). One way to improve the Jc and the Hc2 is to create more grain boundaries in the Nb3Sn alloy allowing them to trap more magnetic flux lines in the materials. To do this the grains must be made smaller which will naturally cause the formation of more grain boundaries. A common way to create smaller grains is by doping an impurity in the alloy. In this study, Nb3Sn was doped with Zr to form Zr-Nb3Sn alloys. This project investigates the surface properties and chemical composition of Zr-Nb3Sn alloys with varying zirconium concentrations (~0.5% and ~24% Zr). Using Scanning Electron Microscopy (SEM), Energy Dispersive Spectroscopy (EDS), and X-ray Photoemission Spectroscopy (XPS), allowed for the characterization of the differences in surface morphology and elemental composition between samples. SEM images were quantified using ImageJ giving insight into grain distribution. The results indicate that increased Zr concentration creates smaller grains in the Zr-Nb3Sn alloy. Furthermore, the ~ 24% Zr grew nanoscale particles on the surface that are possible oxides. XPS data revealed that as Zr concentration increased, the total Nb oxides (NbOx) thickness decreased. Conversely, as Zr concentration increased, the surface Sn oxides (SnOx) also increased slightly. Quantification via ImageJ of SEM data showed that the ~ 0.5% Zr sample had a grain size of ~ 5300 nm2. These results show that Zr inclusion leads to smaller grain formations, which is promising in the world of superconducting RF cavities for high magnetic field applications.

Sue, Micah↗