Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Integrated Computational Materials Engineering”

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 271 records · Page 15

Derivation of physical equations for high-speed laser welding using large language models

It is challenging to formulate complex physical phenomena that occur in a manufacturing process, particularly when the available data are limited, rendering conventional data-driven approaches ineffective. This study aims to predict humping onset in high-speed laser welding by introducing a novel framework, namely text-to-equations generative pre-trained transformer (T2EGPT). This method leverages the capabilities of large language models (LLMs), in combination with sparse experimental data and enriched literature data, to derive an interpretable and generalizable equation for predicting humping initiation. By capturing key correlations among physical parameters, T2EGPT generates a compact and dimensionless expression that accurately predicts hump formation. The equation reveals that humping arises from the interplay between inertia-driven backward melt flow and capillary-driven surface stabilization, where inertial forces drive molten metal backward and capillary forces resist surface deformation. Furthermore, compared to traditional data-driven models, T2EGPT demonstrates enhanced predictive accuracy and cross-material transferability. More broadly, this study highlights the potential of LLMs to integrate textual information with data-driven discovery, enabling the extraction of physical laws in data-scarce scientific domains.

36 MATERIALS SCIENCE↗

MyCrunchGPT: A LLM Assisted Framework for Scientific Machine Learning

Scientific machine learning (SciML) has advanced recently across many different areas in computational science and engineering. Here, the objective is to integrate data and physics seamlessly without the need of employing elaborate and computationally taxing data assimilation schemes. However, preprocessing, problem formulation, code generation, postprocessing, and analysis are still time- consuming and may prevent SciML from wide applicability in industrial applications and in digital twin frameworks. Here, we integrate the various stages of SciML under the umbrella of ChatGPT, to formulate MyCrunchGPT, which plays the role of a conductor orchestrating the entire workflow of SciML based on simple prompts by the user. Specifically, we present two examples that demonstrate the potential use of MyCrunchGPT in optimizing airfoils in aerodynamics, and in obtaining flow fields in various geometries in interactive mode, with emphasis on the validation stage. To demonstrate the flow of the MyCrunchGPT, and create an infrastructure that can facilitate a broader vision, we built a web app based guided user interface, that includes options for a comprehensive summary report. The overall objective is to extend MyCrunchGPT to handle diverse problems in computational mechanics, design, optimization and controls, and general scientific computing tasks involved in SciML, hence using it as a research assistant tool but also as an educational tool. While here the examples focus on fluid mechanics, future versions will target solid mechanics and materials science, geophysics, systems biology, and bioinformatics.

97 MATHEMATICS AND COMPUTING↗

Development of Numerical Model of Metal Foam with PCM for the Estimation of Effective Thermal Conductivity

Global warming due to climate change is a threat to humankind. Nuclear energy is one of the promising solutions to reduce fossil fuel usage. Nuclear energy can handle the base load, compensating for the volatility of renewable energy. If nuclear energy could achieve load following capability, the combination with renewable energy would be more suitable. Thermal energy storage (TES) is one of the options for enabling load following of nuclear reactors. The TES makes it possible to store surplus nuclear thermal energy and release it later as needed. In Idaho National Laboratory (INL), a new concept of latent heat TES integrated with high-temperature heat pipe has been proposed and is under development, which is called Heat pipe-Integrated Thermal Battery (HITB). HITB exchanges thermal energy between the reactor system and TES via heat pipe. The heat transferred to TES medium, made of phase change material (PCM), stores energy as sensible heat and/or latent heat. As PCM typically has poor thermal conductivity, however, various heat transfer enhancement techniques are required to achieve a rapid charging cycle. There are many techniques to enhance the heat transfer ability of TES medium such as disk, fin, and metal foam. Among them, metal foam is an appropriate option to enhance the heat transfer because it maximizes the heat transfer area through metal wicks. Metal foam is a lightweight metal structure that has a high porosity of over 0.9. The typical materials for metal foam are Aluminum, Copper, Nickel, and Silicon Carbide (SiC). Metal foam not only enhances heat transfer via conduction but also increases contact surface area. In the HITB design , the metal foam is being considered as one of the options to enhance the heat transfer of TES medium (PCM) [1]. To predict the enhanced thermal performance of TES, one should properly estimate the effective thermal conductivity of metal foam combined with PCM material or calculate heat transfer in distributed model. There are many experimental works that provides effective thermal conductivity of metal foam with various PCM [2,3]. Also, many theoretical models were developed based on the unit cell model of metal foam [4,5]. With a distributed model, on the other hand, detail heat transfer characteristics between metal foam and PCM material can be analyzed considering the geometry or buoyancy effect. However, due to the complex geometry of metal foam pores, the computational cost for three-dimensional modeling highly increases. Therefore, if metal foam structure can be modeled in simple and repetitive design, the computational cost would decrease Among the various metal foam models [2], lattice model is one of the simple and extendable design. The porosity and pores per inch (PPI) can be characterized by the size and spatial distance of lattice structure. If the three-dimensional metal foam model consists of lattice structure could properly estimate the heat transfer, which is characterized by effective thermal conductivity, it would be a good option to assess the thermal performance of metal foam with PCM. In this study, a three-dimensional numerical model was developed to simulate conductive heat transfer between metal foam and PCM. The three-dimensional lattice structure of square pillars was selected as a basic structure of the metal foam. The calculation result was characterized by the effective thermal conductivity of the whole domain. A sensitivity study was conducted for mesh size, domain size, and PPI to check whether the calculation result gives a converged result or not. Lastly, the effective thermal conductivity from the lattice model was compared with existing experimental data to validate the model result

25 ENERGY STORAGE↗

ALD with Alternative Coreactants: Which Work, Which Do Not, and Why

In this paper, we present a combined experimental and theoretical study that systematically evaluates a series of alcohols (i.e., primary, secondary, and tertiary) as coreactants for trimethylaluminum (TMA)-based atomic layer deposition (ALD) of Al 2 O 3 . We employed in situ quartz crystal microbalance techniques and ex situ X-ray photoelectron spectroscopy to probe growth across a range of temperatures from T = 120 to 285 °C. Dispersion-inclusive hybrid density functional theory was employed to identify potential reaction pathways and compute corresponding Gibbs free-energy barriers and rate coefficients. Experimentally on an Al 2 O 3 surface, sustained thin-film growth, where the rate is constant over many cycles, was only observed with tertiary alcohols (tert-butanol and 2-methyl-2-butanol) at elevated temperatures (e.g., T = 285 °C); for tert-butanol specifically, sustained growth also occurred at T = 240 °C, with no sustained growth at or below T = 210 °C. Growth using primary and secondary alcohols decayed to a negligible value after only 2–3 cycles on an Al 2 O 3 surface at all tested temperatures. Intentional trace additions of H 2 O to anhydrous alcohols restored sustained growth but at reduced deposition rates relative to pure H 2 O. Theoretical analysis supports a bimolecular alkoxy β-H elimination mechanism as the primary reaction pathway: initiated by the formation of a bound alkoxy species, the subsequent β-H elimination step is rate-determining, with a barrier that systematically decreases as the degree of the alcohol coreactant is increased from primary to secondary to tertiary. Computed rate coefficients (k) also indicate a strong temperature dependence (with k increasing by ∼4–5 orders of magnitude from T = 120 to 285 °C) and are consistent with the experimental observation that only tertiary alcohols lead to steady thin-film growth at elevated temperatures. This integrated approach establishes why only certain alcohols can sustain Al 2 O 3 ALD and delivers a predictive framework for identifying effective alcohol coreactants.

36 MATERIALS SCIENCE↗

Unraveling the Mechanics: Computational Modeling of Residual Stresses in U10Mo Fuel with Aluminum Cladding

LEU (Low Enriched Uranium) plate-type fuel elements consist of a high-density, low-enrichment U–Mo alloy-based fuel foil encapsulated in an aluminum cladding and is fabricated through the Hot Isostatic Pressing (HIP) technique. Understanding the mechanism of possible failure modes of nuclear fuel is critical to mitigate potential consequences. One of the major contributing factors in various failure modes is stress and it greatly affects the integrity under temperature, pressure, and irradiation. These stresses could originate from various sources, including manufacturing, and operation, and can lead to deformation, cracking, or even complete failure of fuel elements. This study focuses on the development of a computational model that accurately predicts the residual stresses generated in the U10Mo fuel during fabrication. It has been observed that cladding creep has a substantial impact on the residual stresses in the U10Mo fuel post the HIP fabrication process. Furthermore, during the HIP bonding process the fuel plate system is heated to a temperature of 560 oC and as a result, the aluminum cladding transitions from (Al 6061-T6 to -O). This presents a challenge in capturing the change in material properties accurately. Therefore, a temperature dependent creep model such as hyperbolic sine creep model is considered to estimate the creep properties of the Aluminum cladding. The proposed calibrated creep model accurately predicts the residual stresses in the U10Mo fuel foil and agrees well with the experiments.

42 ENGINEERING↗

A CALPHAD-informed approach to modeling constituent redistribution in Zr-based metallic fuels using BISON

Here, a CALPHAD-informed (Computer Coupling of Phase Diagrams and Thermochemistry) constituent redistribution model was developed for Zr-based metallic fuels and incorporated into the BISON fuel performance code. Three uncertain model parameters associated with β and γ phase kinetics were calibrated using integral test data from U-Zr fuel elements irradiated in Experimental Breeder Reactor II. The calibrated constituent redistribution model was shown to predict the behavior of U-Zr fuels with excellent accuracy. Model predictions for U-Pu-Zr fuels were physically reasonable but less accurate. Reduction of uncertainties in the ternary phase transition temperatures and collection of kinetic data for the ζ phase are expected to improve the model’s ternary predictions. Finally, the new model was coupled to existing thermomechanics models in BISON to simulate irradiation of an entire U-Zr fuel element, demonstrating its ability to accurately predict the behavior of U-Zr fuels with realistic geometries and mesh resolutions at the engineering scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Phase-field predictions of the influence of cooling rates during AM on the Evolution of Microstructures in Nickel-Based Single Crystal Superalloys

Additive manufacturing of single crystals made of Ni-based superalloys offers major cost savings for gas turbine engines with the inclusion of internal cooling channels. However, the lack of understanding of the effect of transient thermal conditions on solidification grain structure during additive manufacturing hinders the potential for process control to maintain the single crystal quality. The use of high-fidelity simulations through high performance computing to predict the evolution of the solidification microstructure will enhance the abilities to tailor the microstructures through process optimization. Phase field simulations are used to determine the effect local thermal conditions and defects on the stability of the solidification morphology, specifically with respect to the onset of columnar-to-equiaxed transition that results in the loss of the single crystal. The results are expected to be instrumental for developing future surrogate models to speed up the integration of design and manufacturing of turbine blades under the harsh in-service conditions.

36 MATERIALS SCIENCE↗

Initial Fracture Propagation Modeling of Graphite Components with Grizzly

Graphite has historically been extensively used in power reactor cores and will be used in multiple types of advanced reactors currently under development. These graphite structural components can experience significant stresses due to nonuniform volumetric strains induced by irradiation and thermal expansion, which can lead to fracture. Robust tools for predicting fracture initiation and propagation in graphite structural components in nuclear reactors are important for evaluating component integrity, developing design standards, and interpreting experimental results to characterize graphite performance. The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation program has been developing degradation models for other structural components in nuclear reactors within the Grizzly and BlackBear codes. This report documents an effort to develop initial capabilities for modeling graphite fracture within these codes, building on prior efforts to model fracture in other materials. Major elements of this effort include developing a new system for modeling fracture nucleation and growth in two dimensions using the extended finite element method and incorporating a damage and plasticity model. These capabilities are applied here to model a representative graphite component and a splitting disc experiment used to obtain tensile strength.

36 MATERIALS SCIENCE↗

Multiscale mechanical design of the lightweight, stiff, and damage-tolerant cuttlebone: A computational study

Cuttlebone, the endoskeleton of cuttlefish, offers an intriguing biological structural model for designing low-density cellular ceramics with high stiffness and damage tolerance. Cuttlebone is highly porous (porosity ~93%) and lightweight (density less than 20% of seawater), constructed mainly by brittle aragonite (95 wt%), but capable of sustaining hydrostatic water pressures over 20 atmospheres and exhibits energy absorption capability under compression comparable to many metallic foams (~4.4 kJ/kg). Here, in this work, we computationally investigate how such remarkable mechanical efficiency is enabled by the multiscale structure of cuttlebone. Using the common cuttlefish, Sepia Officinalis, as a model system, we first conducted high-resolution synchrotron micro-computed tomography (µ-CT) and quantified the cuttlebone's multiscale geometry, including the 3D asymmetric shape of individual walls, the wall assembly patterns, and the long-range structural gradient of walls across the entire cuttlebone (ca. 38 chambers). The acquired 3D structural information enables systematic finite-element simulations, which further reveal the multiscale mechanical design of cuttlebone: at the wall level, wall asymmetry provides optimized energy absorption while maintaining high structural stiffness; at the chamber level, variation of walls (number, pattern, and waviness amplitude) contributes to progressive damage; at the entire skeletal level, the gradient of chamber heights tailors the local mechanical anisotropy of the cuttlebone for reduced stress concentration. Our results provide integrated insights into understanding the cuttlebone's multiscale mechanical design and provide useful knowledge for the designs of lightweight cellular ceramics.

36 MATERIALS SCIENCE↗

RADAI: A Large-Scale Realistic Dataset for Radiation Detection Algorithm Development

Open, realistic datasets are essential for developing and benchmarking radiation detection algorithms, yet they remain scarce. The Radiological Anomaly Detection and Identification (RADAI) project was develop to create datasets that meet the training and testing needs for sophisticated radiation detection algorithms. The RADAI dataset is a large-scale synthetic resource that integrates high-fidelity Monte Carlo simulations with realistic urban scenarios to capture both background variability and source signatures. RADAI models construction-material NORM, people and vehicles, urban clutter, and dynamic environmental effects such as cosmic-ray and rain-induced transients, and they provide list-mode detector data with motion and response modeling suitable for algorithm training and evaluation. The RADAI project resulted in three publicly-released complementary datasets together with an online scoring portal for standardized performance assessment and an open software toolkit that supports data access, augmentation, model development, and evaluation. These resources enable reproducible comparisons across methods and promote rigorous studies at the scale required by contemporary machine learning. By grounding algorithm development in realistic, well-documented conditions, RADAI supports progress toward more robust detection, identification, and localization in complex urban environments.

Ghawaly, James M. [Division of Computer Science an↗

PyOECP: A flexible open-source software library for estimating and modeling the complex permittivity based on the open-ended coaxial probe (OECP) technique

Here, we present PyOECP, a Python-based flexible open-source software for estimating and modeling the complex permittivity obtained from the open-ended coaxial probe (OECP) technique. The transformation of the measured reflection coefficient to complex permittivity is performed based on three different methods. The software library contains the dielectric spectra of common reference liquids, which can be used to transform the reflection coefficient into the dielectric spectra. Several Python routines that are commonly employed (e.g., SciPy and NumPy) in the field of science and engineering are required only so that the users can alter the software structure depending on their needs. The modeling algorithm exploits the Markov Chain Monte Carlo method for the data regression. The discrete relaxation models can be built by a proper combination of well-known relaxation models. In addition to these models, electrode polarization, a typical measurement artifact for interpreting dielectric spectra, can be incorporated into the modeling algorithm. A continuous relaxation model, which solves the Fredholm integral equation of the first kind (a mathematically ill-posed problem), is also included. This open-source software enables users to freely adjust the physical parameters to obtain physical insight into their materials under test and will be consistently updated for more accurate measurement and interpretation of dielectric spectra in an automated manner. This work describes the theoretical and mathematical background of the software, lays out the workflow, and validates the software functionality based on both synthetic and empirical data included in the software.

97 MATHEMATICS AND COMPUTING↗

Non-Nuclear Test Facility Development to Support Verification and Validation Needs

This study focuses on developing a non-nuclear test facility to help meet the verification and validation (V&V) requirements of a water-cooled small modular reactor (SMR) system. The scaling analysis necessary for reactor system experimentation and assessment model V&V considers various reactor system accident scenarios such as main steam line break (MSLB), steam generator tube rupture (SGTR), and containment condensation and heat transfer. The goal is to identify key physics based phenomena of interest (POI), specific figures of merit (FOMs), and the state of knowledge (SOK) by employing phenomena identification and ranking table (PIRT) studies. Important physics phenomena such as heat transfer, fluid dynamics, structural integrity, and material selection must be addressed when scaling a thermal-hydraulics system. Alongside these considerations, modifications to the facility may be needed to align the scaled system with the available footprint (height and space) of the facility. These modifications may include infrastructure upgrades, changes to the piping layout, implementation of safety measures, and installation of new instrumentation and control systems. The test facility aims to provide a scaled-down version of prototypic facilities for comprehensive V&V activities by overcoming these challenges and making the necessary modifications. This approach minimizes scaling distortion, simulates prototypic plant conditions, and accurately models accident event progressions. Thus, the facility will be capable of adequately modeling a scaled SMR system to support the design, development and licensing needs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Melting Behavior of Hydrogen Direct Reduced Iron in Molten Steel and Slag: An Integrated Computational and Experimental Study

Direct reduced iron (DRI) and hot briquetted iron (HBI) are essential feedstocks for tramp element control in the electric arc furnace (EAF). Due to greenhouse gas (GHG) concerns related to CO2 emissions, hydrogen as a substitute for natural gas and a reductant in DRI production is being widely explored to reduce GHG emissions in ironmaking. This study examines the melting behavior of hydrogen DRI (H-DRI) pellets in the EAF containing low-carbon (0.1 wt.%) molten steel and molten slag. A computational heat transfer model was developed to predict the melting behavior of H-DRI pellets. To validate the model, a set of experimental laboratory simulations was conducted by immersing H-DRI in a molten steel bath and slag. The temperature history at the center of the pellet during melting and the shell thickness at different melting stages were utilized to validate the model. The simulation results agree with the experimental measurements of steel balls and H-DRI in different metallic molten steel and slag baths.

Materials Science↗

Fracture mechanics approach to TRISO fuel particle failure analysis

Weibull stress-based methods for failure probability assessment have been developed and analyzed to assess the integrity of tristructural isotropic (TRISO) fuel particles during fuel life cycles and accident operating conditions. While simple, these methods entail a number of drawbacks when stress concentrates near crack tips, including finite element mesh size dependency when the Weibull stress is averaged over the finite element domain. Fracture mechanics approaches involving the use of interaction integrals eliminate this lack of mesh convergence and produce consistent fracture predictions. In this work, we use an interaction integral approach to computing stress intensity factors for a crack in the inner pyrolytic carbon layer perpendicular to the silicon carbide layer, which is simplified representation of a failure mode in TRISO particles. Further, the interface between these two TRISO layers has been shown to become porous, which we simulate by considering a transition of mechanical properties over such porous length, i.e. the layers are modeled as a functionally graded material. Aspects such as porosity and thermal and irradiation eigenstrains are considered in computing the stress intensity factor from a fracture mechanics approach and compared with the known Weibull stress failure approach. The methodology introduced in this paper enables a more general fracture probability assessment in TRISO particles and eliminates the need to identify best suited parameters when using local or averaged stress-based failure criterion. Finally, the numerical sensitivity studies show how parameters such as the porous transition zone length, the material stiffness, and creep affect the probability of TRISO fuel particle failure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ultra-thick three-dimensional interpenetrating graphene electrode architectures for high volumetric density energy storage

For electrochemical energy storage, increasing the electrode thickness is an effective approach to achieving higher energy density from a given material. However, this often compromises ion transport, leading to diminished performance. Here, in this study, we present a novel platform for fabricating complex 3D interpenetrating electrode structures via photo-polymerization 3D printing, integrated with computational structural optimization for energy storage. The platform employs an acrylate resin system infused with graphene oxide (GO), enabling high-fidelity printing of optimized porous structures and facilitating efficient electron and ion transport in ultra-thick electrodes. The optimized 3D layouts substantially enhance energy and power densities compared to conventional configurations, ensuring superior material utilization and minimal ohmic losses. Supercapacitors fabricated using this approach achieved an exceptional energy density of 4.7 Wh L−1 at a power density of 1689.0 W L−1, surpassing traditional designs. This work underscores the transformative role of structural optimization in advancing electrochemical performance and establishes a versatile pathway for developing next-generation energy storage systems with exceptional efficiency and functionality.

Wang, Zhen [University of California, Berkeley, CA↗

Quantum Engineering With Hybrid Magnonic Systems and Materials (Invited Paper)

Quantum technology has made tremendous strides over the past two decades with remarkable advances in materials engineering, circuit design and dynamic operation. In particular, the integration of different quantum modules has benefited from hybrid quantum systems, which provide an important pathway for harnessing the different natural advantages of complementary quantum systems and for engineering new functionalities. This review focuses on the current frontiers with respect to utilizing magnetic excitations or magnons for novel quantum functionality. Magnons are the fundamental excitations of magnetically ordered solid-state materials and provide great tunability and exibility for interacting with various quantum modules for integration in diverse quantum systems. The concomitant rich variety of physics and material selections enable exploration of novel quantum phenomena in materials science and engineering. In addition, the relative ease of generating strong coupling and forming hybrid dynamic systems with other excitations makes hybrid magnonics a unique platform for quantum engineering. We start our discussion with circuit-based hybrid magnonic systems, which are coupled with microwave photons and acoustic phonons. Subsequently, we are focusing on the recent progress of magnon-magnon coupling within confined magnetic systems. Next we highlight new opportunities for understanding the interactions between magnons and nitrogen-vacancy centers for quantum sensing and implementing quantum interconnects. Lastly, we focus on the spin excitations and magnon spectra of novel quantum materials investigated with advanced optical characterization.

42 ENGINEERING↗

PVDeg: Enhancing Usability and AI-Driven Multi-Mechanism Degradation Modeling

PVDeg version 0.7.0, released in December 2025, introduced major enhancements to improve usability and performance. This update reorganized tutorials and tool notebooks to create a more intuitive experience, enabling users to easily follow and adapt workflows for their specific analyses. In addition to structural improvements, both the notebooks and core logic underwent significant optimization for efficiency, robustness, and style. These refinements were supported by new testing frameworks built on nbval and pytest, adherence to PEP8 standards, and extensive code refactoring, which collectively simplify onboarding for new developers. Looking ahead, version 0.8.0 will deliver advanced AI-driven capabilities. The primary focus is to further develop and automate the degradation workflow, designed to analyze PV module degradation across diverse locations and system configurations. By integrating large language models (LLMs) to scan literature and compile a comprehensive database of materials and degradation rates, this feature will enable modeling of multiple materials and mechanisms within a single, streamlined workflow. Users will be able to evaluate degradation impacts on different system architectures under varying environmental conditions, facilitating informed decisions on bill-of-materials optimization for specific deployment scenarios. These advancements position PVDeg as a powerful, user-friendly tool for accelerating PV reliability research and system design.

14 SOLAR ENERGY↗

Understanding Reliability Trade-Offs in 1T-nC and 2T-nC FeRAM Designs

Ferroelectric random access memory (FeRAM) is a promising candidate for energy-efficient nonvolatile memory, particularly for logic-in-memory and compute-in-memory (CIM) applications. Among the available cell architectures, One-Transistor–n-Capacitor (1T-nC) and two-transistor–n-capacitor (2T-nC) FeRAMs each offer distinct trade-offs in density, scalability, and reliability. In this work, we present a comparative study of these two architectures under both dimensional scaling ( XY/Z shrinkage) and vertical integration (increasing stacked capacitors per cell). Using technology computer-aided design (TCAD) and circuit-level simulations, we analyze how scaling impacts ferroelectric capacitance, parasitic coupling, and floating-node (FN) dynamics, which together dictate sense margin (SM) and read stability. A key mitigation strategy—floating unselected capacitors—is applied to both architectures, effectively decoupling the SM from the number of stacked capacitors and enabling tractable analysis across scaling regimes. Results show that 1T-nC suffers more from charge sharing with the bitline (BL), while 2T-nC benefits from transistor isolation and stronger low-voltage sensing at the cost of increased area. By systematically evaluating these behaviors across scaling directions, this work establishes the reliability trade-offs of 1T-nC and 2T-nC cells and provides design guidelines for high-density, vertically integrated FeRAM systems.

1T-nC↗