Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational materials design”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A modular methodology for time-domain stochastic seismic wave propagation

Presented here is a modular methodology for time-domain stochastic seismic wave propagation analysis. Presented methodology is designed to analyse uncertain seismic motions as an input, propagating through uncertain material. Traditional approach for uncertain wave propagation relies on models that include deep bedrock, local soil site, and their random process and random field information. Such models can become quite large and computationally intractable. The modular approach proposed herein features two step approach that allows separate consideration of the deep bedrock and local site along with corresponding random field information. In this work, the first step considers an auxiliary stochastic motions problem in the bedrock. Stochastic local site response can then be simulated in a reduced domain within certain depth from the surface. Application of uncertain seismic motions at depth, for local uncertain site response is done using stochastic effective forces developed through the Domain Reduction Method. By using Hermite polynomial chaos expansion to represent the non-Gaussian random field of material parameters and non-stationary random process of seismic motion, the proposed modular methodology is formulated using intrusive stochastic Galerkin approach, as seen in the Stochastic Elastic–Plastic Finite Element Method (SEPFEM). Developed modular methodology is illustrated using a 1-D stochastic seismic wave propagation analysis with three cases, and simulation results are also verified with results from conventional approach.

58 GEOSCIENCES↗
DOE OSTI2021

UiO-66 Metal–Organic Framework as an Anode for a Potassium-Ion Battery: Quantum Mechanical Analysis

The natural abundance of potassium in the earth’s crust is 1000 times higher than that of lithium, so energy technologies built on potassium are more sustainable. Potassium-ion batteries have attracted considerable attention because of their relatively low cost and high operating potential, but questions remain about the best anode material for such batteries. Here, in this paper, we report first-principles computations based on density functional theory to investigate the performance of the UiO-66 metal–organic framework as an anode material for potassium-ion batteries; the goal is to provide a fundamental understanding of metal–organic framework (MOF)-based electrodes to guide the design and development of high-performance potassium-ion batteries. Our study includes the stability and electronic properties of potassiated structures and the mechanisms of potassium intercalation and diffusion in the framework. The results indicate that UiO-66 has a maximum specific capacity of 644 mAh/g as the anode of a potassium-ion battery. During potassiation, we observe charge transfer from potassium to carbon or oxygen of UiO-66 near the intercalated K. During K diffusion, the K migrates along the UiO-66 framework with a maximal migration energy barrier of 0.377 eV in the optimal pathway, which is much lower than the barriers for Li and Na diffusion in UiO-66. The diffusion coefficient of K in the anode is several orders of magnitude larger than those of Li and Na. This favors potassium ions over lithium ions or sodium ions when UiO-66 is the anode.

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

Tritium adsorption and absorption on (100) and (001) surfaces of pure and tin defective zirconium

Zirconium alloys such as zircaloy-4 are used as tritium (T) getter materials in tritium-producing burnable absorber rods (TPBARs) due to their ability to capture T, thereby forming metal hydrides. Developing an understanding of T adsorption onto zircaloy prior to diffusion into the subsurface is relevant for rational tritium getter and TPBAR design, to improve material properties for nuclear applications. Herein, density functional theory calculations revealed the preferred binding sites for T adsorption on Zr(001) and Zr(100). The energy barriers of T transfer, along the surface and from the surface to the subsurface were computed. The adsorption properties of Zr(001) were found to be superior to those of Zr(100). Surface tin impurities were found to strongly repel T. The presence of subsurface and surface tin resulted in higher absorption energy barriers for both the forward and reverse processes. Based on the calculated energy barriers, a surface to surface T diffusion coefficient of 9.53 × 10 -10 m 2 s -1 is expected for pristine Zr(001). Finally, a surface to subsurface T diffusion coefficient on the order of 10 -13 m 2 s -1 is predicted in pristine Zr, decreasing to 10 -19 m 2 s -1 for the transfer with a subsurface tin impurity.

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

Computer simulation of carbonization and graphitization of coal

Abstract This study describes computer simulations of carbonization and graphite formation, including the effects of hydrogen, nitrogen, oxygen, and sulfur. We introduce a novel technique to simulate carbonization, ‘Simulation of Thermal Emission of Atoms and Molecules (STEAM)’, designed to elucidate volatile outgassing and density variations in the intermediate material during carbonization. The investigation analyzes the functional groups that endure through high-temperature carbonization and examines the graphitization processes in carbon-rich materials containing non-carbon impurity elements. The physical, vibrational, and electronic attributes of impure amorphous graphite are analyzed, and the impact of nitrogen on electronic conduction is investigated.

Ugwumadu, C. (ORCID:0000000199207594)↗
DOE OSTI2020

Design Optimization of A generic Fissile Solution for Mo99 production using Electron Beam-based Neutron Generator using MCNP+CFD

Molybdenum-99 ( 99 Mo) is a medical isotope used in 80% of all medical imaging procedures today. However, 99 Mo is not directly administered to patients for imaging; its decay product, technetium-99m ( 99 mTc), is a metastable isotope with a half-life of ~6 hours. Such a short half-life does not allow for production, separation, and application to a patient before decaying further. Therefore, the solution is to produce the precursor to 99 mTc, 99 Mo. 99 Mo has a half-life of 66 hours, long enough to produce the material and send it to the desired location before it fully decays. Under National Nuclear Security Administration’s (NNSA) highly enriched uranium (HEU) minimization mission, the Material Management and Minimization (M3) office leads the Molybdenum-99 (Mo-99) program. A major goal of the Mo-99 program is to develop methods of producing 99 Mo using low enriched uranium (LEU) opposed to the traditionally used HEU domestically in the U.S. Internationally, LEU is the new standard for use in research reactors and isotope production facility. In this study, we introduce a generic fissile solution system for 99 Mo production facility meeting the design requirement of 1) LEU usage as fissile material and 2) subcritical fissile system coupled with a horizontal electron beam (E-beam) accelerator. The focus of this research is to optimize a fissile solution configuration for maximum 99 Mo production yield using MCNP model, and ultimately implement the optimized design with computational fluid dynamics (CFD) model to understand thermal hydraulic behavior and solution convection characteristic. The report consists of two parts: 1) isotope yield calculation and criticality analysis with a wide range of parametric test matrix, and 2) thermal hydraulic analysis for solution mixing and cooling assessment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Dissimilar joining of ZEK100 and AA6022 for automotive application

Mg and Al joining is problematic due to different high temperature flow characteristics that often results in weld defects and brittle intermetallic compounds (IMC). We present process development on lap joining of 1.27 mm 6022 (top sheet) and 1.5mm ZEK100 (bottom sheet) for an automotive door application. We conducted process development in linear welding using power control mode (modulating torque and spindle speed) with multiple tool designs, and welding parameters. Process variables were correlated to joint micro-structure and mechanical properties. A 2D computational model was created to understand the failure mechanism and contribution of mechanical interlocking towards joint strength. For a thin (< 2 mm) automotive sheet with softer material as a bottom sheet, adequate control of process variables is critical to maintain weld stability and avoid excessive IMC layer at the interface. A combination of mechanical anchoring and metallurgical reaction is likely responsible for bonding.

Das, Hrishikesh↗
DOE OSTI2023

2.3.4.501 - Synthesis and Analysis of Performance-Advantaged Bioproducts

This project focuses on the synthesis and analysis of performance-advantaged bioproducts (PABPs). We have established collaborations with other BETO-funded projects and academic and industrial collaborators to source new molecules that have promising manufacturing pathways and that could be serve as performance-advantaged biochemicals or biopolymers. We conduct synthesis and characterization of biochemicals and biopolymers alongside techno-economic analysis and life cycle assessment to estimate their cost and environmental impacts relative to incumbent materials. As part of the project, Profs. Linda Broadbelt and Brent Shanks are developing computational pathway prediction tools to identify optimal production pathways for bio-based compounds via biological and chemo-catalytic transformations. When coupled to the polyID tool from the Inverse Design project, these tools will ultimately enable a narrowing of design space for PABPs. From FY21-FY23, we described a framework for benchmarking PAPBs, estimated the energy and GHG emissions for commodity organic chemicals, developed performance-advantaged nylons and polyesters from beta-ketoadipic acid, and produced lignin-based plasticizers. We have shown that aromatic amines can be used in performance thermosets and that polyhydroxyalkanoates with crosslinked side chains can exhibit rubber-like properties, along with repair and degradability. We are working actively with industry partners on scale-up and validation of multiple PABPs.

BIOMASS FUELS↗
DOE OSTI2021

StressNet - Deep learning to predict stress with fracture propagation in brittle materials

Abstract Catastrophic failure in brittle materials is often due to the rapid growth and coalescence of cracks aided by high internal stresses. Hence, accurate prediction of maximum internal stress is critical to predicting time to failure and improving the fracture resistance and reliability of materials. Existing high-fidelity methods, such as the Finite-Discrete Element Model (FDEM), are limited by their high computational cost. Therefore, to reduce computational cost while preserving accuracy, a deep learning model, StressNet, is proposed to predict the entire sequence of maximum internal stress based on fracture propagation and the initial stress data. More specifically, the Temporal Independent Convolutional Neural Network (TI-CNN) is designed to capture the spatial features of fractures like fracture path and spall regions, and the Bidirectional Long Short-term Memory (Bi-LSTM) Network is adapted to capture the temporal features. By fusing these features, the evolution in time of the maximum internal stress can be accurately predicted. Moreover, an adaptive loss function is designed by dynamically integrating the Mean Squared Error (MSE) and the Mean Absolute Percentage Error (MAPE), to reflect the fluctuations in maximum internal stress. After training, the proposed model is able to compute accurate multi-step predictions of maximum internal stress in approximately 20 seconds, as compared to the FDEM run time of 4 h, with an average MAPE of 2% relative to test data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multi deep learning-based stochastic microstructure reconstruction and high-fidelity micromechanics simulation of time-dependent ceramic matrix composite response

A multi deep learning-based framework is developed for efficient, automated microstructure reconstruction and generation of stochastic representative volume elements (SRVEs) with periodic boundary conditions (PBCs) for accurate modeling of ceramic matrix composite (CMC) response. The methodology comprises a convolutional neural network coupled with regression layers to act as a vanilla regression network for semantic segmentation of the microstructure, allowing accurate characterization of the phases and their distributions at the microscale. Scanning electron microscope and confocal microscope are used to obtain C/SiNC and SiC/SiNC CMCs micrographs for vanilla regression testing. Microstructure variability in terms of fiber volume fraction and porosity are quantified through the output regression layer, ensuring accurate representation of material variability in SRVE construction. Generative adversarial network (GAN) and its variants are designed to produce high-fidelity SRVE, spanning CMCs microstructure variability space. A circular padding algorithm is developed to generate SRVEs with PBCs during training of GANs. The accuracy of the generated SRVEs is established through micromechanics simulations, where an efficient formulation of the high-fidelity generalized methods of cells (HFGMC) approach is used to compute the effective mechanical properties. Furthermore, an iterative algorithm is implemented in the HFGMC solver to simulate time-dependent deformation of SiC/SiNC subjected to creep loading conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Methane Adsorption on Heteroatom-Modified Maquettes of Porous Carbon Surfaces

Experimental and theoretical studies disagree on the energetics of methane adsorption on carbon materials. However, this information is critical for the rational design and optimization of the structure and composition of adsorbents for natural gas storage. The delicate nature of dispersion interactions, polarization of both the adsorbent and the adsorbate, interplay between H-bonding and tetrel bonding, and induced dipole/Coulomb interactions inherent to methane physisorption require computational treatment at the highest possible level of theory. In this study, we employed the smallest reasonable computational model, a maquette of porous carbon surfaces with a central site for substitution and methane binding. The most accurate predictions of methane adsorption energetics were achieved by electron-correlated molecular orbital theory CCSD(T) and hybrid density functional theory MN15 calculations employing a saturated, all-electron basis set. Here, the characteristic geometry of methane adsorption on a carbon surface (“lander approach”) arises due to bonding interactions of the adsorbent π-system with the proximal H–C bonds of methane, in addition to tetrel bonding between the antibonding orbital of the distal C–H bond and the central atom of the maquette (C, B, or N). The polarization of the electron density, structural deformations, and the comprehensive energetic analysis clearly indicate a ∼3 kJ mol –1 preference for methane binding on the N-substituted maquette. The B-substituted maquette showed a comparable or lower binding energy than the unsubstituted, pure C model, depending on the level of theory employed. The calculated thermodynamic results indicate a strategy for incorporating electron-enriched substitutions (e.g., N) into carbon materials as a way to increase methane storage capacity over electron-deficient (e.g., B) modifications. The thermochemical analysis was revised for establishing a conceptual agreement between the experimental isosteric heat of adsorption and the binding enthalpies from statistical thermodynamics principles.

03 NATURAL GAS↗
DOE OSTI2023

Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials

With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Thermophysical Property Measurements on Salt Mixture

Thermal conductivity, specific heat, viscosity, and density are critical properties for understanding thermal hydraulics of heat transfer fluids (HTFs) like molten salts. These properties are important for design and safe operation of any HTF loop. These characterizations are required to benchmark computational chemical models that are used to predict properties, but existent data are limited. These measurements are extremely difficult to make due to the operating temperature range and corrosive properties of the salt.Density measurements are needed to predict volumetric behavior. Thermal expansion mismatches can cause mechanical failure of structural materials and important components that are in contact with the molten salt. All of these properties are related to heat transfer in HTFs. The heat generated in a reactor must be high enough to derive power; however it should not be so hot that it cannot be safely removed via the HTF at an appropriate rate. It is important to understand the thermal hydraulics required to operate at the desired power generation level and to avoid situations in which thermal shifts can damage components,cause failure of structural materials, and/or lead to accidents. Understanding these thermal properties can aid in component and process design and can contribute to optimized operational and safety parameters.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

An in-situ view cell system for investigating swelling behavior of elastomers upon high-pressure hydrogen exposure

The transition to hydrogen as a clean and efficient energy carrier is impeded by challenges in the compatibility of hydrogen with materials used within hydrogen infrastructure. Elastomers, crucial in sealing components, often exhibit premature failures in high-pressure hydrogen environments due to excessive swelling. This study employs an innovative in-situ view cell system to assess the swelling behavior of hydrogenated nitrile butadiene rubber (HNBR) under various hydrogen conditions. The system, designed to withstand pressures up to 96.5 MPa, incorporates Digital Image Correlation (DIC) for strain measurements and volume estimation. Results reveal non-linear volume increases during depressurization, challenging conventional assumptions. Furthermore, investigations into peak hydrogen pressures and pressure-holding scenarios during decompression highlight complex swelling trends. The introduction of a novel computer vision (CV) method enhances precision in volume estimation, overcoming DIC limitations. The study provides insights into mitigating elastomer swelling, crucial for developing robust materials to support future hydrogen-driven energy systems.

Elastomer↗
DOE OSTI2024

Chapter 4: Physically informed deep learning networks for simulating microstructure evolution of 3D polycrystals

As discussed in the previous chapter, high energy diffraction microscopy (HEDM) is used to study the micromechanical evolution of a material during in situ loading. HEDM experiments have been used to verify crystal plasticity (CP) simulations [119, 91, 90, 120], for experimental planning, material design, and to further analyze experimental results. However, Fast Fourier transform-based CP (CP-FFT) or finite element-based CP (CP-FE) methods are often too slow to be used in real-time during an experiment. CP-FFT is faster than CP-FE simulations due to the absence of meshing, but can still take hours to simulate the response of a single volume depending on the size and number of strain steps [127]. Reducing computation time would create a larger exploration space in planning and design, and enable faster analysis of experimental results and real-time feedback during an experiment. This research expands upon previous works to develop a workflow for predicting the full-field evolution of a 3D polycrystal. The workflow is simplified from previous works to predict only orientation and elastic strain tensors (from which stress tensors are calculated). The network is physically informed through loss functions and network architecture for a more robust model. The orientation predictions are informed about the cubic crystal symmetry of the material by incorporating disorientation and misorientation information into the network architecture and loss. The Von Mises stress is used to enforce the correct stress-strain trends in the strain tensor predictions. Additional total strain steps from the elastic and elastoplastic region are included to better capture the stress-strain evolution at smaller total strain steps. Material and hardening parameters are additional inputs into the networks to further inform the network and to study the network’s ability to predict different materials other than those used for training.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Iron-based binary ferromagnets for transverse thermoelectric conversion

Thermoelectric generation using the anomalous Nernst effect (ANE) has great potential for application in energy harvesting technology because the transverse geometry of the Nernst effect should enable efficient, large-area and flexible coverage of a heat source. For such applications to be viable, substantial improvements will be necessary not only for their performance but also for the associated material costs, safety and stability. In terms of the electronic structure, the anomalous Nernst effect (ANE) originates from the Berry curvature of the conduction electrons near the Fermi energy. To design a large Berry curvature, several approaches have been considered using nodal points and lines in momentum space. Here we perform a high-throughput computational search and find that 25 percent doping of aluminium and gallium in alpha iron, a naturally abundant and low-cost element, dramatically enhances the ANE by a factor of more than ten, reaching about 4 and 6 microvolts per kelvin at room temperature, respectively, close to the highest value reported so far. The comparison between experiment and theory indicates that the Fermi energy tuning to the nodal web—a flat band structure made of interconnected nodal lines—is the key for the strong enhancement in the transverse thermoelectric coefficient, reaching a value of about 5 amperes per kelvin per metre with a logarithmic temperature dependence. We have also succeeded in fabricating thin films that exhibit a large ANE at zero field, which could be suitable for designing low-cost, flexible microelectronic thermoelectric generators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Organic Photovoltaic Efficiency Predictor: Data-Driven Models for Non-Fullerene Acceptor Organic Solar Cells

In the design of organic solar cells, there has been a need for materials with high power conversion efficiencies. Scharber’s model is commonly used to predict efficiency; however, it exhibits poor performance with new non-fullerene acceptor (NFA) devices, since it was designed for fullerene-based devices. In this work, an empirical model is proposed that can be a more accurate alternative for NFA organic solar cells. Additionally, many screening studies use computationally expensive methods. A model based on using semiempirical simplified time-dependent density functional theory (sTD-DFT) as an alternative method can accelerate the calculations and yield a similar accuracy. The models presented in this paper, termed organic photovoltaic efficiency predictor (OPEP) models, have shown significantly lower errors than previous models, with OPEP/B3LYP yielding errors of 1.53% and OPEP/sTD- DFT of 1.55%. As a result, the proposed computational models can be used for the fast and accurate screening of new high-efficiency NFAs/donor pairs.

14 SOLAR ENERGY↗
DOE OSTI2024

Rational Design of Novel Biomimetic Sequence-Defined Polymers for Mineralization Applications

Silica biomineralization is a naturally occurring process, wherein organisms use proteins and other biological structures to direct the formation of complex, hierarchical nanostructures. Discovery and characterization of such proteins and their underlying mechanisms spurred significant efforts to identify routes for biomimetic mineralization that reproduce the exquisite shapes and size selectivities found in nature. A common strategy has been the use of short peptide sequences with chemistry mimicking those found in natural systems, such as the use of the silaffin-derived R5 peptide. While progress has been made using this approach, there are many limitations that have prevented breakthroughs in biomimicry. To advance our ability to use charged macromolecules for silica formation, we propose to use sequence-defined synthetic polymers known as peptoids, or N-substituted polyglycines, which present significant capability for the precise tuning of sequence and structure beyond what can often be achieved with peptides alone. This study presents a computationally predicted design of these polymers that leads to the controlled formation of silica nanomaterials. We investigate surface adsorption and the mineralization process through analysis of binding mechanisms and energetics of the R5 system. Next, we synthesized two R5-inspired peptoids and validated our prediction in the design of mineralization polymers through characterization using surface plasmon resonance and electron microscopy. Here, this computationally guided study holds great promise for designing new sequences with unprecedented control of the placement of chemical functional groups, thus allowing for further unraveling of silicification mechanisms and the eventual design of sequence-defined synthetic polymers leading to the predictive synthesis of nanostructured functional materials.

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

Bimodal Visualization of Industrial X-Ray and Neutron Computed Tomography Data

Advanced manufacturing creates increasingly complex objects with material compositions that are often difficult to characterize by a single modality. Our collaborating domain scientists are going beyond traditional methods by employing both X-ray and neutron computed tomography to obtain complementary representations expected to better resolve material boundaries. However, the use of two modalities creates its own challenges for visualization, requiring either complex adjustments of bimodal transfer functions or the need for multiple views. Together with experts in nondestructive evaluation, we designed a novel interactive bimodal visualization approach to create a combined view of the co-registered X-ray and neutron acquisitions of industrial objects. Using an automatic topological segmentation of the bivariate histogram of X-ray and neutron values as a starting point, the system provides a simple yet effective interface to easily create, explore, and adjust a bimodal visualization. Here, we propose a widget with simple brushing interactions that enables the user to quickly correct the segmented histogram results. Our semiautomated system enables domain experts to intuitively explore large bimodal datasets without the need for either advanced segmentation algorithms or knowledge of visualization techniques. We demonstrate our approach using synthetic examples, industrial phantom objects created to stress bimodal scanning techniques, and real-world objects, and we discuss expert feedback.

image segmentation↗