SEARCH · Engineering Papers

Results for “materials knowledge systems”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Three-beam hybrid fs/ps coherent anti-Stokes Raman scattering of rotation–vibration non-equilibrium

In this work, we demonstrate time-resolved rotational and vibrational temperature measurements to probe plasma non-equilibrium through a simple three-beam hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (fs/ps CARS) system. A single pump/Stokes pair is employed to generate both the pure-rotational and ro-vibrational Raman coherence. A novel phase-matching scheme, to our knowledge, is employed to spatially overlap the two CARS signals, simplifying the optical design and allowing signal selection by tuning a single spectrometer grating. Measurements were performed near the electrodes in a N 2 DC glow discharge. Both the rotational temperature and the vibrational populations up to $v$ = 8 were calculated from the separately measured single-shot pure-rotational and ro-vibrational CARS spectra. Strong rotation–vibration non-equilibrium was observed at both electrodes, with the cathode showing higher vibrational and rotational temperatures. In addition, non-Boltzmann behaviors were observed at both electrodes. This simplified approach enables measurements of rotational and vibrational temperatures with high spatial resolution in non-equilibrium flows.

42 ENGINEERING↗

DOE OSTI2024

Universality of Collective Density Fluctuations in Liquids at and away from Equilibrium – An Integrated Neutron Scattering, Theoretical, and Computational Study (Final Technical Report)

Liquids, ubiquitous on earth, are prototypical disordered condensed matter. As the mediate phase between non-structured gases and ordered solids, liquids p ossess intriguing complexities that are absent on either extreme counter phase, especially when they are driven away from equilibrium. Despite decades of studies, the physics of liquids still keeps surprising us and challenges our understanding of condensed matter, emergent phenomena, and complex systems. Among the many interesting properties, the nature of collective dynamics in liquids remains elusive because of both the presence of strong interactions and the absence of translational invariance. This project studied the universal features of collective density fluctuations (or collective modes), which include both phonon-like excitations and collective relaxations, in the Q-dependent generalized hydrodynamic regime in liquids and liquid-like matter both at and away from equilibrium. We performed systematic studies of the fast collective dynamics (including both phonon-like excitations and collective relaxations) of four classes of liquids characterized by different fragilities and interactions using synergistically integrated coherent neutron scattering experiments (Inelastic Neutron Scattering (INS), Quasi-Elastic Neutron Scattering (QENS), and Neutron Spin Echo (NSE)), a ViscoElastic Hydrodynamics theory we have been developing, andMolecular Dynamics (MD) simulations and our analysis package LiquidLib. The transformative knowledge of the collective dynamics in liquids is important to understand and ultimately control the transport of mass, energy, and charge in liquids, both at and away from thermal equilibrium, which may will likely impact a wide range of disciplines.

36 MATERIALS SCIENCE↗

DOE OSTI2025

ORNL Design Engineering Library: Validated Thermoelectric Digital Twins

Radioisotope thermoelectric generators (RTGs) serve a crucial role in supplying thermal and electrical energy for reliable and long-duration power in remote and extreme environments. The thermal energy is provided by the decay of radioisotopes and is converted into electrical energy using the steady-state thermoelectric process by exploiting the Seebeck effect. Maintaining a thermal gradient through the material is necessary to drive current production. Modern predictive multi-physics tools are able to effectively describe the complex and strongly coupled phenomena needed for accurate model-based system engineering efforts needed to develop higher performing RTGs. Experimental validation of model simulations is essential to confirm the accuracy, reliability, and applicability of predictive multiphysics tools. Validation fosters confidence among users, helps meet regulatory requirements, and contributes to the ongoing improvement of simulation techniques, ultimately leading to better, safer, and more efficient designs and processes. Validated models can then be used as digital twins and can be interrogated to understand and quantify the performance gaps between theoretical and physical systems. These insights can be used to build higher performing systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Comparative Studies of the Structural and Transport Properties of Molten Salt FLiNaK Using the Machine-Learned Neural Network and Reparametrized Classical Forcefields

Despite surging interest in molten salt reactors and thermal storage systems, knowledge of the physicochemical properties of molten salts are still inadequate due to demanding experiments that require high temperature, impurity control, and corrosion mitigation. Therefore, the ability to predict these properties for molten salts from first-principles computations is urgently needed. Herein, we developed and compared a machine-learned neural network force field (NNFF) and a reparametrized rigid ion model (RIM) for a prototypical molten salt LiF–NaF–KF (FLiNaK). We found that NNFF was able to reproduce both the structural and transport properties of the molten salt with first-principles accuracy and classical-MD computational efficiency. Furthermore, the correlation between the local atomic structures and the dynamics was identified by comparing with RIMs, suggesting the significance of polarization of anions implicitly embedded in the NNFF. Furthermore, this work demonstrated a computational framework that can facilitate the screening of molten salts with different chemical compositions, impurities, and additives, and at different thermodynamic conditions suitable for the next-generation nuclear reactors and thermal energy storage facilities.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Connection between the physicochemical characteristics of amorphous carbon thin films and their electrochemical properties

Connecting a material's surface chemistry with its electrocatalytic performance is one of the major questions in analytical electrochemistry. This is especially important in many sensor applications where analytes from complex media need to be measured. Unfortunately, today this connection is still largely missing except perhaps for the most simple ideal model systems. Here we present an approach that can be used to obtain insights about this missing connection and apply it to the case of carbon nanomaterials. In this paper we show that by combining advanced computational techniques augmented by machine learning methods with x-ray absorption spectroscopy (XAS) and electrochemical measurements, it is possible to obtain a deeper understanding of the correlation between local surface chemistry and electrochemical performance. As a test case we show how by computationally assessing the growth of amorphous carbon (a-C) thin films at the atomic level, we can create computational structural motifs that may in turn be used to deconvolute the XAS data from the real samples resulting in local chemical information. Then, by carrying out electrochemical measurements on the same samples from which x-ray spectra were measured and that were further characterized computationally, it is possible to gain insight into the interplay between the local surface chemistry and electrochemical performance. To demonstrate this methodology, we proceed as follows: after assessing the basic electrochemical properties of a-C films, we investigate the effect of short HNO 3 treatment on the sensitivity of these electrodes towards an inner sphere redox probe dopamine to gain knowledge about the influence of altered surface chemistry to observed electrochemical performance. Overall, these results pave the way towards a more general assessment of electrocatalysis in different systems and provide the first steps towards data driven tailoring of electrode surfaces to gain optimal performance in a given application.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Corrosion of Containment Alloys in Molten Salt Reactors and the Prospect of Online Monitoring

The aim of this review is to communicate some essential knowledge of the underlying mechanism of the corrosion of structural containment alloys during molten salt reactor operation in the context of prospective online monitoring in future MSR installations. The formation of metal halide species and the progression of their concentration in the molten salt do reflect containment corrosion, tracing the depletion of alloying metals at the alloy salt interface will assure safe conditions during reactor operation. Even though the progress of alloying metal halides concentrations in the molten salt do strongly understate actual corrosion rates, their prospective 1 st order kinetics followed by near-linearly increase is attributed to homogeneous matrix corrosion. The service life of the structural containment alloy is derived from homogeneous matrix corrosion and near-surface void formation but less so from intergranular cracking (IGC) and pitting corrosion. Online monitoring of corrosion species is of particular interest for molten chloride systems since besides the expected formation of chromium chloride species CrCl 2 and CrCl 3 , other metal chloride species such as FeCl 2 , FeCl 3 , MoCl 2 , MnCl 2 and NiCl 2 will form, depending on the selected structural alloy. The metal chloride concentrations should follow, after an incubation period of about 10,000 hours, a linear projection with a positive slope and a steady increase of <1 ppm per day. During the incubation period metal concentration show 1 st order kinetics and increasing linearly with time. Ideally, a linear increase reflects homogeneous matrix corrosion, while a sharp increase in the metal chloride concentration could set a warning flag for potential material failure within the projected service life, e.g. as result of intergranular cracking or pitting corrosion. Continuous monitoring of metal chloride concentrations can therefore provide direct information about the mechanism of the ongoing corrosion scenario and offer valuable information for a timely warning of prospective material failure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Thermodynamic re-modelling of the Cu–Nb–Sn system: Integrating the nausite phase

Currently available Cu–Nb–Sn phase diagrams lack the recently discovered nausite phase (Cu,Nb)Sn 2 , which is an important intermediate in the course of thermal processing of superconducting Nb 3 Sn wires. Processing decisively determines the resulting microstructure of Nb 3 Sn and, thus, its superconducting properties. Lack of suitable and complete phase diagrams, however, obstructs rational design of such thermal processing procedures. To close this gap and to obtain valid knowledge of homogeneity and stability range of nausite, various Cu–Nb–Sn samples, which are heat-treated between 300 °C and 500 °C, are investigated. By means of energy-dispersive X-ray spectroscopy (EDX), a temperature-dependent homogeneity range of nausite is observed, which covers average mole fractions of Cu between 0.09 and 0.15. This is correlated with a change in the mean atomic volume and can be seen in the lattice parameters determined by X-ray diffraction (XRD). Additionally performed first-principles calculations on different CuSn 2 and NbSn 2 model structures confirm this trend. Furthermore, the peritectic decomposition of nausite to NbSn 2 and liquid at 586 °C is determined by means of in situ XRD and differential scanning calorimetry (DSC). By using the CALPHAD (CALculation of PHase Diagrams) approach, all these findings are used to extend a previous thermodynamic description of the Cu–Nb–Sn system by including the nausite as an additional phase. Finally, with this noteworthy integration, the updated modelling of the Cu–Nb–Sn system can be used for optimizing the multistage heat-treatment steps during processing superconducting Nb 3 Sn wires.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Synthesis of Nano-Oxide Precipitates by Implantation of Ti, Y and O Ions in Fe-10%Cr: Towards an Understanding of Precipitation in Oxide Dispersion-Strengthened (ODS) Steels

The properties of oxide dispersion-strengthened steels are highly dependent on the nature and size distribution of their constituting nano-oxide precipitates. A fine control of the processes of synthesis would enable the optimization of pertinent properties for use in various energy systems. This control, however, requires knowledge of the precise mechanisms of nucleation and growth of the nanoprecipitates, which are still a matter of debate. In the present study, nano-oxide precipitates were produced via the implantation of Y, Ti, and O ions in two different sequential orders in an Fe-10%Cr matrix that was subsequently thermally annealed. The results show that the oxides that precipitate are not necessarily favoured thermodynamically, but rather result from complex kinetics aspects related to the interaction between the implanted elements and induced defects. When Y is implanted first, the formation of nanoprecipitates with characteristics similar to those in conventionally produced ODS steels, especially with a core/shell structure, is evidenced. In contrast, when implantation starts with Ti, the precipitation of yttria during subsequent high-temperature annealing is totally suppressed, and corundum Cr 2 O 3 precipitates instead. Moreover, the systematic involvement of {110} matrix planes in orientation relationships with the precipitates, independently of the precipitate nature, suggests matrix restriction effects on the early stages of precipitation.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A Novel Dew Point Meter: Application to the Measurement of the Sulfuric Acid Dew Point for Combustion Flue Gas

Accurate knowledge of acid dew point is essential for industrial and applied combustion. Sulfur in the fuel or raw materials is converted to sulfur dioxide (SO2) during combustion, and a portion of the SO2 is oxidized to sulfur trioxide (SO3). The SO3 will react to form H2SO4 vapor when in the presence of water vapor. Even with just trace levels of H2SO4 vapor in the gas phase (1-10 ppm), the dew point can reach 100°C and higher. To avoid acid condensation and the resulting corrosion on heat recovery equipment, plant engineers must ensure that surface temperatures are above the acid dew point, but this decreases the efficiency of the thermal energy recovery. Thus, there is a trade-off between minimizing equipment corrosion and maximizing thermal energy recovery efficiency, and the acid dew point is a key parameter for this optimization. In this work, we analyze and discuss the theory of the sulfuric acid condensation process and develop a novel dew point meter based on this analysis. The meter consists of a novel optical instrument that is designed to monitor the slightest appearance of condensation on a hydrophobic window surface, and the surface temperature of the window is slowly decreased until condensation is first observed, yielding an accurate measurement of the dew point. The basis of the instrument is that a collimated beam from a diode laser generates scattered light when encountering surface condensate, and a sophisticated array detector is used to sensitively monitor the light scattering. The measurement procedures are established to rapidly find the acid dew point, while minimizing error. Further, we propose a calibration system based on a liquid bubbler, which can generate a stable gas flow with known sulfuric acid dew point, to test the dew point meter. Test results show that the dew point meter can accurately measure the acid dew point over a wide range.

Cheng, Mao↗

DOE OSTI2021

Advanced Syngas Cleanup Using Radically Engineered Modular Systems

The overall objective of this project was to address key knowledge gaps and develop modular sorbent-based warm syngas cleanup designs that will enable 1- to 5-MW REMS-based plants to be cost-competitive with large state-of-the-art commercial plants utilizing all of our abundant domestic coal reserves. We proposed to develop two potential new desulfurization process designs specifically targeted for REMS plants for combined heat and power or polygeneration from coal gasification. The first design is a modification of our commercially offered dual transport reactor design based Warm-gas Desulfurization Process (WDP) where the transport-reactor regenerator is replaced with a fluidized-bed regenerator. This design leverages the learnings obtained from our pre-commercial demonstration testing of the WDP technology for over 3,000 hours and conclusions from the lessons learned workshop. The second design is a fixed-bed process that leverages the successful composition and chemistry of our WDP sorbent to develop an effective fixed-bed sorbent formulation and fixed-bed process design. These two potential desulfurization process designs have a strong technical likelihood of success and offer the best opportunity to effectively employ design standardization, mass production and advanced manufacturing techniques to enable cost competitiveness for syngas cleanup at smaller scale.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2020

Library of Advanced Materials for Engineering (LAM) 4.56

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAME advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAME library in application, this effort seeks to document and verify the various models in the LAME library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Tight-Binding-Based Method for Simulations of Complex Materials with Accompanying Database

The NRL-Tight-Binding (NRL-TB) method has been applied successfully to many monatomic, binary, and a few ternary systems. The method consists of a fit to Density Functional Theory (DFT) results and has transferability beyond the fitted DFT database including performing tight-binding molecular dynamics simulations. We are working on a new website with a detailed guide of how to use the method in performing the fit to DFT results (Fitting Code) and obtaining non-orthogonal TB parameters expressed in a polynomial form as environment-dependent Slater-Koster parameters (on-site terms), and neighbor distance-dependent Hamiltonian and Overlap terms. These parameters are then used in a Static Code to verify or make predictions for large systems often beyond the computational capabilities of standard DFT codes. Our starting point was a website developed at NRL some years ago which NRL no longer supports. We have now been augmenting that website with new examples of expanding the original capabilities of the method such as incorporating in the codes f-electron states. Our website (esd.cos.gmu.edu/tb) contains two separate branches: One named NRL-TB with the full documentation of the method including source codes and examples. The other branch, named ESD-dB, is a library of electronic structure calculations for all the monatomic elements (a subset of the elements also found in a recent book by the author of this proposal. See reference 22) and an approximately 700 hundred compounds in specific crystal structures. We discuss below the progress made so far in improving and expanding the capabilities of our TB method and the plans for further extensions of our approach.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Mitigation documentation for the removal of hazardous materials from the engine installation vehicle and manned control car area 25, Nevada national security site, NYE county, Nevada

The purpose of this letter report is to document the mitigation of adverse effects of a proposed undertaking that would remove hazardous materials from the Engine Installation Vehicle (EIV) in Area 25 of the Nevada National Security Site (NNNSS) in compliance with Section 106 of the National Historic Preservation Act (NHPA) and the terms of the 2024 Programmatic Agreement among the U.S. Department of Energy and the Nevada State Historic Preservation Officer and the Advisory Council on Historic Preservation Concerning the Protection of Historic Properties on the Nevada National Security Site, Nye County, Nevada, hereafter referred to as the NNSS PA. The EIV (State Historic Preservation Office [SHPO] Resource No. S3057) has been determined individually eligible for listing in the National Register of Historic Places (NRHP) under Criteria A and C and as a contributing element to the Nuclear Rocket Development Station (NRDS) Historic District, which is eligible for the NRHP under all four Significance Criteria (Reno et al. 2023; Reed 2024).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Corrosion Testing of Refractory inContact with Molten Glasses Designed for Waste Vitrification - VSL Touchpoint Matrix Glasses

It is known that the predictive life of the refractory ceramic liner of nuclear waste glass melters is conservative, as demonstrated by performance of these materials such as in the Defense Waste Processing Facility (DWPF). The motivation for this task is to maximize the useful life of the melters that will be operated at the Waste Treatment and Immobilization Plant (WTP), which will in turn minimize procurement and disposal costs and melter outage times, as well as to identify maximum loadings in the waste glass of those species that corrode melter components. This task was initiated jointly with Pacific Northwest National Laboratory (PNNL) with the objective to develop a methodology and model to enable more accurate prediction of refractory service life under prototypic conditions from laboratory-scale material corrosion tests. Refractory corrosion is generally reported as physical material loss, measured in units of distance (e.g., inch) or as physical material loss rate, measured in units of distance per time (e.g., inch/day). In post-operational melters, the refractory corrosion is measured directly, sometimes reported as corrosion depth. Crucible tests are used in the laboratory to accelerate the refractory corrosion to facilitate a meaningful measurement in a commensurate amount of time. Crucible tests are particularly useful in understanding refractory corrosion across a large glass composition space, where operational testing would be prohibitive. Some of the critical parameters that are known to influence refractory corrosion by molten glass in a crucible test are temperature, system redox, molten salt phases, glass chemistry, and test duration. The majority of data collected for Monofrax® K-3 (hereafter referred to as K-3) corrosion is from crucible tests, but a small amount comes directly from scaled and production melters. Crucible test data has been collected under varying conditions, whereas data collected from operational melters is relatively fewer and represents conditions specific to the melter campaign. The result is that the published data can be grouped and analyzed in multiple ways, not all of which are readily comparable. The Standard Test Method for Isothermal Corrosion Resistance of Refractories to Molten Glass (ASTM C621) outlines the general guidelines used across industry. That method describes a sealed, static test in which the surface area of the refractory coupon and the volume of glass are fixed. A significant portion of the crucible data pertaining to nuclear waste glasses has been collected in a modified configuration; the most notable differences being the surface area of the refractory coupon to volume of the glass and use of a method for bubbling the melt. To our knowledge, the influence of those parameters on the refractory corrosion has not been quantified. In this work, it was determined that static tests and bubbled tests would be performed. Savannah River National Laboratory (SRNL) was tasked with setting up and performing static testing while PNNL was tasked with setting up and performing bubbled testing. Initial activities were performed to establish laboratory methods that reproduce data comparable to existing data sets of K-3 refractory corrosion by low activity waste (LAW) and high-level waste (HLW) glass compositions. Later activities were focused on refining the test parameters to establish a standard test practice to be used between Laboratories and collecting additional data to be used in the enhanced waste glass model development. This document serves primarily to convey the refractory loss measurement results from corrosion testing of K-3 refractory with waste glass compositions developed for use in the WTP melters. The data will be used in the enhanced property/composition models being developed for waste glass vitrification and melter operations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Prediction of the Al-rich part of the Al-Cu phase diagram using cluster expansion and statistical mechanics

In this work, the thermodynamic properties of α-Al and other phases (GP zones, θ'', θ' and θ) in the Al-rich part of the Al-Cu system have been obtained by means of the cluster expansion formalism in combination with statistical mechanics. This information was used to build the Al-rich part of the Al-Cu phase-diagram taking into account vibrational entropic contributions for θ', as those of the other phases were negligible. The simulation predictions of the phase boundaries between α-Al and either θ'', θ' or θ phases as a function of temperature are in good agreement with experimental data and extend the phase boundaries to a wider temperature range. The DFT calculations reveal the presence of a number of metastable Guinier-Preston-zone type configurations that may coexist with α-Al and θ'' at low temperatures. They also demonstrate that θ' is the stable phase below 550K but it is replaced by θ above this temperature due to the vibrational entropic contribution to the Gibbs energy of θ'. This work shows how the combination of cluster expansion and statistical mechanics can be used to expand our knowledge of the phase diagram of metallic alloys and to provide Gibbs free energies of different phases that can be used as input in mesoscale simulations of precipitation.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Library of Advanced Materials for Engineering (LAMÉ) 5.2

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗