Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials simulations”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Canister Centerline Cooling Experiments for DPF5-336 Reference Material Made with ERV3b Salt Simulant

This report provides experimental details and results on the elemental distributions (compositions), crystalline/amorphous phase distribution, and microstructures for DPF5-336 reference materials made with ERV3b salt that were subjected to the canister centerline cooling (CCC) heat treatment process. To fabricate these materials, ammonium dihydrogen phosphate was mixed with ERV3b salt simulant and a dechlorination procedure was run in the generation-2 dechlorinator at the Pacific Northwest National Laboratory up to 600°C in an alumina crucible. After the dechlorination process, the product was removed from the crucible, lightly crushed, loaded into a new alumina crucible, Fe 2 O 3 was added, and this was vitrified at 1100°C for 1 hour and quenched. The quenched material was mostly amorphous but showed some small diffraction peaks attributed to Li 3 Fe 2 (PO 4 ) 3 . This quenched material was ground to a fine particle size in a milling chamber, and this was subjected to the CCC cooling profile in duplicate within alumina and fused quartz crucibles. The CCC experiment is designed to simulate the slowest cooling profile (the vertical centerline) of a 0.61-m diameter canister where the slower cooling profile often results in the production of a variety of different crystalline phases upon formation of which result in a residual glass of unknown composition and properties. The information gleaned from these types of experiments will lead to improved composition-property predictions for materials with similar compositions. The results of the alumina and fused quartz experiments show different microstructures but similar diffraction patterns providing evidence of similar crystalline phases being present in each sample. Both samples had monazite [i.e., (Ce,Nd)PO 4 ] and mixtures of phases containing P, Fe, and alkali metal elements. This report completes the milestone M4FT-23PN030104042 with details provided in Appendix B.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Metal–ceramic composite structures for fabrication of high power density plasmonic devices

The recent decade brought many advances to plasmonics, but high power density plasmonic antennas designed to behave as heaters or operate in high temperature environments are still facing material stability challenges preventing their ultimate use. Gold has been the optimal choice among plasmonic materials but experiences morphology changes at temperature that result in device efficiency reduction and failure. Bulk titanium nitride has been explored as a solution but has deal-breaking tradeoffs in device quality factor. In this paper, we explore via proof-of-concept the use of a metal–ceramic composite structure to determine whether a bulk Au nanorod can provide strong plasmonic resonances while coated with an ultrathin conformal layer of titanium nitride or silica to provide morphological stability and sufficient plasmonic activity without excessive resonance quality degradation. We show SEM-level morphological stability for temperatures up to 500 °C with coatings below 4 nm. Computer modeling suggests the ultrathin titanium nitride has detrimental effects on the strong plasmonic resonances of a Au nanorod. Here, we then looked into other possible coatings for solutions to stabilize high power density plasmonic antennas including plasmonic oxides, metal adhesion layers, and silica, the latter appearing to be the best option while lowering the overall peak electric field intensity, the silica increases the electric field intensity at its boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sampling lattices in semi-grand canonical ensemble with autoregressive machine learning

Calculating thermodynamic potentials and observables efficiently and accurately is key for the application of statistical mechanics simulations to materials science. However, naive Monte Carlo approaches, on which such calculations are often dependent, struggle to scale to complex materials in many state-of-the-art disciplines such as the design of high entropy alloys or multi-component catalysts. To address this issue, we adapt sampling tools built upon machine learning-based generative modeling to the materials space by transforming them into the semi-grand canonical ensemble. Furthermore, we show that the resulting models are transferable across wide ranges of thermodynamic conditions and can be implemented with any internal energy model U, allowing integration into many existing materials workflows. We demonstrate the applicability of this approach to the simulation of benchmark systems (AgPd, CuAu) that exhibit diverse thermodynamic behavior in their phase diagrams. Finally, we discuss remaining challenges in model development and promising research directions for future improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

VERA-Grizzly Ex-Core Calculations: Watts Bar Unit 1 Cycles 1-2

The critical structures that comprise light-water reactor (LWR) nuclear power plants are subjected to operating environments that can challenge their integrity. Structures in close proximity to the reactor core, such as the reactor pressure vessel (RPV) and the biological shield wall, are subjected to high levels of radiation emanating from the core, as well as elevated temperatures. As the US fleet of operating LWRs ages, the effects of these operating environments on the integrity of these structures must be considered to ensure their continued safe operation. Extending the lifetime of commercial reactors and maintaining the aging reactor fleet require accurate prediction of the exposure of ex-core components to neutron and photon radiation. In particular, concrete degradation studies must be performed to evaluate the safety and long-term operation of reactors with lifetime extensions. The concrete reactor bioshield is important for providing radiological protection during operation and must last for the entire lifetime of the reactor. Recent interest in lifetime extensions furthers the need to accurately simulate concrete material degradation in the reactor bioshield. As a result of this need, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded this study to couple its tools, Virtual Environment for Reactor Applications (VERA) and Grizzly. VERA allows users to set up models to calculate time-dependent and fully coupled solutions (with thermal feedback) for ex-core quantities of interest such as vessel and coupon fluence and detector responses for multiple statepoints and cycles. Grizzly is a finite-element application based on the Multiphysics Object Oriented Simulation Environment (MOOSE) framework that is used to enable aging materials calculations. This report highlights the work performed to calculate the fluence in the vessel and concrete for Watts Bar Nuclear Plant Unit 1 (WBN1) Cycles 1 and 2. The fluences obtained from VERA were successfully transferred to Grizzly using a Python script. Four simulations were run with Grizzly: (1) the Mazars model with the initial Young’s modulus being the instantaneous modulus, (2) the Mazars model with the initial Young’s modulus being the delayed modulus, (3) the Mazars model with the initial Young’s modulus being the delayed modulus with the addition of the effects of micro-damage caused by irradiation, and (4) the Mazars model with the initial Young’s modulus being the instantaneous modulus, and with the addition of micro-damage and creep. Details regarding the methods used to obtain the fluence and the statistical errors associated with the VERA Monte Carlo Shift calculations are discussed in greater detail in this report. The results obtained from the four Grizzly models are also presented in this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
DOE OSTI2022

Extensive Attention Mechanisms in Graph Neural Networks for Materials Discovery

We present our research where attention mechanism is extensively applied to various aspects of graph neural net- works for predicting materials properties. As a result, surrogate models can not only replace costly simulations for materials screening but also formulate hypotheses and insights to guide further design exploration. We predict formation energy of the Materials Project and gas adsorption of crystalline adsorbents, and demonstrate the superior performance of our graph neural networks. Moreover, attention reveals important substructures that the machine learning models deem important for a material to achieve desired target properties. Our model is based solely on standard structural input files containing atomistic descriptions of the adsorbent material candidates. We construct novel methodological extensions to match the prediction accuracy of state-of-the-art models some of which were built with hundreds of features at much higher computational cost. We show that sophisticated neural networks can obviate the need for elaborate feature engineering. Our approach can be more broadly applied to optimize gas capture processes at industrial scale.

Cong, Guojing↗
DOE OSTI2023

A probabilistic-phase field model for the fracture of brittle materials

Abstract We develop a computational method to determine the failure probability of brittle materials under general mechanical loading conditions. The method is a combination of two parts: (1) numerical simulations of materials with multiple cracks using phase field theory, where the complete fracture process is viewed as ‘damage percolation’ along critical paths or clusters of cracks, rather than the traditional weak-link failure mechanism of Weibull, and (2) an extension of the Batdorf statistical theory of fracture to finite domains, where it is implemented within the finite element framework. The results of phase-field simulations at the ‘percolation threshold’ are used as failure data in the Batdorf theory to determine the overall probability of failure. The input to this approach is the size distribution of cracks in a pristine material. An example is shown, where alumina samples that were previously tested by Abe and coworkers (Abeet al2003J. Am. Ceram. Soc.861019–21) in four-point loading are compared to the results of our numerical simulations. The approach developed here has the advantage of being extendable to more complex thermomechanical loading.

Materials Science↗
DOE OSTI2025

Machine Learning-Driven Optimization of Building Enclosures for Moisture Durability and Thermal Performance

The design of moisture-durable building enclosures with low embodied carbon often involves an iterative process of selecting the materials for the specific exposure conditions to meet the performance requirements. While hygrothermal simulations are commonly used to evaluate moisture durability, they often require advanced expertise for proper implementation. Machine learning (ML) provides a promising alternative by streamlining the design process and minimizing the reliance on complex simulations. This study presents a machine learning-based approach for predicting moisture durability in residential wall assemblies. The ML model was trained to estimate the mold index and maximum moisture content of various layers under typical exposure conditions. The model achieved a high predictive accuracy, with a coefficient of determination (R²) exceeding 0.90 when compared to traditional hygrothermal simulations on materials that were not part of training the ML model. Building on these results, the ML model was developed into a practical tool for optimizing wall assembly designs. This tool allows users to automatically optimize material selections based on energy, moisture, and carbon performance criteria. By incorporating multi-objective optimization, the tool identifies configurations that minimize embodied carbon while maintaining moisture safety and code-compliant thermal performance. Additionally, it provides insights into how material choices influence assembly durability, energy efficiency, and carbon reduction. The tool will be implemented in the Building Science Advisor (BSA) to enhance its performance and provide more granularity on the results. This research highlights the potential for ML-driven tools to simplify the design of high-performance building enclosures, offering architects and engineers a faster, more efficient way to balance critical performance factors.

Salonvaara, Mikael [ORNL] (ORCID:0000000318991554)↗
DOE OSTI2026

Optimizing spin qubit coherence through materials codesign

The evolution of defect-based spin qubit systems is currently transitioning from fundamental studies and proof-of-concept demonstrations into applications in the burgeoning field of quantum technology. Within this context, new challenges emerge, in particular, the need to understand and engineer the fundamental materials that form the hardware building blocks critical for the scalability and wide-scale adoption of such technologies. While earlier discussions have often focused on qubits within idealized systems, major limitations on spin coherence and optical properties arise from effects imposed by the nonideality of the surrounding host matrix. Decoherence can stem from a variety of sources, including other qubits, nuclear spins, and parasitic point- and extended defects, which interact with the qubit via magnetic and electric fields, photons, phonons, and strain. In this article, we focus on the relevant sources and mechanisms through which decoherence occurs and provide potential mitigation strategies via the synergistic integration of first-principles simulations and materials synthesis and engineering. We aim to provide a tangible link between material properties and material functions thereby enabling materials-by-design.

ab initio simulations↗
DOE OSTI2021

Sensitivity Analysis of Tunable Equation of State Material Model In Pulsed Mercury Target Simulation

A pulsed neutron spallation target is subjected to very short but intense loads from repeated proton pulses. Approximately 60% of the energy from each proton pulse is deposited into the mercury target material and the stainless-steel target structure, leading to a high-pressure region in both the stationary target structure and the flowing mercury. The high-pressure region propagates and leads to fluid-structure interaction. The resultant loading on the target structure containing liquid mercury is difficult to predict, although various simulation approaches and material models for the mercury have been tried. To date, the best match of simulation to experimental data is obtained by using an equation of state (EOS) material model with a specified tensile cutoff pressure, which simulates the cavitation threshold. The inclusion of a threshold to represent cavitation is key to the successful predictions of stress waves triggered by the high-energy pulse striking the mercury and vessel. However, recent measurements of target structure strain show that significant discrepancies remain between the measured and simulated strain values in the EOS mercury model. These differences grow when noncondensable helium gas is intentionally injected into the flowing mercury to reduce the loading on the structure. An EOS-based proportional–integral–derivative (PID) mercury model has been proposed to reduce the gap between the measured and simulated vessel strain responses for targets with gas injection. The conceptual and numerical description and initial investigation of the PID model are presented in previous work. Further studies of this PID model — including the sensitivity of the structure’s strain response to model parameters (the tensile cutoff, PID parameters Kp, Ki, and Kd) — are reported in this article. Results show the strain response is more sensitive to changes in the tensile cutoff value than to changes in the model parameters Kp, Ki, and Kd. These results will aid in future work where the model parameters will be optimized to match simulation data to strain measurements.

Lin, Lianshan↗
DOE OSTI2020

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

DFT-based insight into finite-temperature properties of ferroelectric perovskites with lone-pair: the case of CsGeX 3 (X = Cl, Br, I)

Ferroelectrics remain in the focus of scientific attention for decades owing to their fundamental and practical appeal. Recently, ferroelectricity has been demonstrated in semiconducting halide perovskites (Zhang et al 2022 Sci. Adv. 8 eabj5881), offering both a rare combination of ferroelectricity and semiconductivity in the same material and a possible alternative to the prevailing perovskite oxide ferroelectrics. We propose a route to simulating such materials at finite temperatures capable of reproducing key experimental and first-principle data, such as Curie temperature, phase transition sequence, spontaneous polarization, and soft mode frequencies. The key methodological finding is the superior performance of hybrid exchange correlation functionals in parametrization of effective Hamiltonians for ferroelectrics with lone pair. The parametrization for effective Hamiltonians for CsGeX 3 (X = Cl, Br, I) is reported. The application of methodology to study polarization reversal in CsGeX 3 allows for the development of a ‘minimalistic’ model for polarization reversal in ferroelectrics that provides an insight into the mechanisms of polarization reversal and its key features, such as the relationship between the coercive field, temperature, and AC field frequency. Importantly, the model reveals the origin of the well-known and ever-puzzling overestimation of coercive fields in computations. Furthermore, we report a variety of finite-temperature properties of CsGeX 3 ferroelectrics, such as dielectric susceptibility, pyroelectric coefficients, and energy storage density, which reveal that these halide perovskites possess properties comparable to their oxide counterparts. Here, we believe that our work provides significant methodological advancements, deepens fundamental understanding of ferroelectrics, and reveals the potential of halide perovskite ferroelectrics.

effective Hamiltonian↗
DOE OSTI2025

Investigation of Electrocatalytic CO 2 Reduction on MXene Materials via First-Principles Simulations

Computational studies of CO 2 reduction to yield various products were carried out on the basal plane and edges of different MXene materials. The impact of vacancies upon Mo 2 TiC 2 T x , W 2 TiC 2 T x , and Ti 3 C 2 T x (T x = O and OH) was also examined for both edge and basal sites. Initial calibrations were carried out to generate surfaces with optimal oxide/hydroxide ratios and proper termination sites upon which various vacancy sites were explored to ensure an accurate model of the surfaces’ resting states. From this work, Mo 2 TiC 2 O was determined to exhibit the lowest theoretical overpotential for methane as determined by volcano plot analyses. At a large enough vacancy concentration, the CO 2 reduction reaction (CO 2 RR) is predicted to outcompete the hydrogen evolution reaction (HER) as the predominant reaction on the surface. In conclusion, when examining the edge of Mo 2 TiC 2 O, stronger CO 2 binding was exhibited to such an extent that the reaction was predicted to terminate after the generation of formate on the edge.

Lund, Colton [Argonne National Laboratory (ANL), A↗
DOE OSTI2021

Tusas: A fully implicit parallel approach for coupled phase-field equations

In this study, we develop a fully-coupled, fully-implicit approach for phase-field modeling of solidification in metals and alloys. Predictive simulation of solidification in pure metals and metal alloys remains a significant challenge in the field of materials science, as microstructure formation during the solidification process plays a critical role in the properties and performance of the solid material. Our simulation approach consists of a finite element spatial discretization of the fully-coupled nonlinear system of partial differential equations at the microscale, which is treated implicitly in time with a preconditioned Jacobian-free Newton-Krylov method. The approach is algorithmically scalable as well as efficient due to an effective preconditioning strategy based on algebraic multigrid and block factorization. We implement this approach in the open-source Tusas framework, which is a general, flexible tool developed in C++ for solving coupled systems of nonlinear partial differential equations. The performance of our approach is analyzed in terms of algorithmic scalability and efficiency, while the computational performance of Tusas is presented in terms of parallel scalability and efficiency on emerging heterogeneous architectures. We demonstrate that modern algorithms, discretizations, and computational science, and heterogeneous hardware provide a robust route for predictive phase-field simulation of microstructure evolution during additive manufacturing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Quantum Theory, Quantum Materials, Quantum Computing

The Sanibel Symposium series is renowned amongst materials theorists, quantum chemists, and condensed matter physicists as meetings driving progress on theory, mod eling, and simulation of materials and their molecular and nano-scale constituents. The Symposia are highly unusual (perhaps unique) in their priority emphasis on theory and computation, in having no parallel sessions, in cultivating well-attended Hot-Topic contributed oral sessions, and accessible poster sessions. These provide highly visible, influential platforms for cross-fertilization among specialist investigators, hence are strong contributors to the advance of quantum information sciences (QIS) research of strategic importance to the Office of Basic Energy Sciences (BES). As part of a five-year plan to highlight QIS challenges and opportunities and foster progress on them, each of the pre ceding three Sanibel Symposia had a thematic focus, Quantum Theory, Quantum Materi als, Quantum Computing, as a major program component. Emphasis was on quantum materials and their molecular constituents. The award for 2024 was for year four of that sustained thematic focus.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Realization of fermionic Laughlin state on a quantum processor

Strongly correlated topological phases of matter are central to modern condensed matter physics and quantum information technology but often challenging to probe and control in material systems. The experimental difficulty of accessing these phases has motivated the use of engineered quantum platforms for simulation and manipulation of exotic topological states. Among these, the Laughlin state stands as a cornerstone for topological matter, embodying fractionalization, anyonic excitations, and incompressibility. Although its bosonic analogs have been realized on programmable quantum simulators, a genuine fermionic Laughlin state has yet to be demonstrated on a quantum processor. Here, we realize the ν = 1/3 fermionic Laughlin state on IonQ’s trapped-ion quantum computer using an efficient and scalable Hamiltonian variational ansatz with 369 two-qubit gates on a 16-qubit circuit. Employing symmetry-verification error mitigation, we extract key observables that characterize the Laughlin state, including correlation hole, bulk-edge correspondence, and topological entanglement entropy, with strong agreement to exact diagonalization benchmarks. This work demonstrates an end-to-end workflow to simulate material-intrinsic topological orders and provides a starting point to explore its dynamics and excitations on digital quantum processors.

Shen, Lingnan [Univ. of Washington, Seattle, WA (U↗
DOE OSTI2023

Using Homogenized Models to Explore the Effect of Fracture Densities on Weathering

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to predict depth of weathering for a system where the fracture density decreases downward, the model will overestimate regolith depth. Learning how to develop accurately homogenized models could thus enable better conceptual models and predictions of weathering advance in natural systems.

58 GEOSCIENCES↗