Search NASA⌕ Search

SEARCH · Search NASA

Results for “defects simulation”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Multiscale Failure Analysis of a 3D Woven Unit Cell Containing Defects

Multiscale failure simulations have been performed for a Three-dimensional woven composite unit cell considering five, or more, length scales spanning the woven composite mesoscale to the sub-microscale voids. The multiscale recursive micromechanics approach, which enables recursive integration of general micromechanics theories over an arbitrary number of length scales, has been employed within the NASA Multiscale Analysis Tool. The multiscale model uses both the generalized method of cells and Mori-Tanaka micromechanics theories, and considers failure in the constituent materials using a simple damage model. Baseline results, containing distributed voids, are compared to uniaxial experimental data for an AS4 carbon fiber/ RTM6 epoxy matrix 3D orthogonal woven composite with good agreement in terms of global stiffness and global failure stress. The simulations demonstrate that the 3D woven composite exhibits damage tolerance through sustaining increasing axial load far beyond the first initiation of damage. The multiscale model is used to examine the nonlinear response of the material to other loading conditions. Case studies, motivated by X-ray computed tomography data, are presented on the effects of manufacturing induced voids and cracks.

3D woven↗
DOE OSTI2024

pymatgen-analysis-defects: A Python package for analyzing point defects in crystalline materials

Point defects can often determine the properties of semiconductor and optoelectronic materials. Due to the large simulation cell and the higher-cost density functionals required for defect simulations, the computational cost of defect calculations is often orders of magnitude higher than that of bulk calculations. As such, managing and curating the results of the defect calculations generated by a single user has the potential to save a significant amount of computational resources. Moreover, eventually building a high-quality, persistent defects database will significantly reduce the computational cost of defect calculations for the entire community.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗
DOE OSTI2021

Computational Fermi Level Engineering and Doping-Type Conversion of Ga2O3 Via Three-Step Processing

Ga2O3 is being actively explored for power electronics, and other applications due to its ultra- wide bandgap and low projected fabrication cost of high-quality crystals. N-type doping of Ga2O3 can be achieved and tuned, but p-type doping faces fundamental obstacles. Successful engineering of Ga2O3 based devices requires critical control of doping density, Fermi level position, providing opportunities for predictive process simulation. We use first-principles defect equilibrium calculations to simulate a 3-step growth- annealing-quench protocol for hydrogen assisted Mg doping in Ga2O3, taking into account the hydrogen- oxygen-water gas phase equilibrium. We predict type conversion to a net p-type regime following O-rich annealing after growth under reducing conditions in the presence of H2. We show that there is an optimal temperature that maximizes the net acceptor density during the equilibrium annealing step for a given Mg doping level. Quenching of non-equilibrium annealed samples then results in a Fermi level EF below mid-gap down to about EV +1.5 eV, creating a significant number of uncompensated neutral MgGa0 acceptors. This type converted Ga2O3 can create significant built-in field in a p-n junction with an adjoining n-type material and could enable vertical MOSFETs.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗
DOE OSTI2021

Computational Fermi Level Engineering and Doping-Type Conversion of Ga2O3 Via Three-Step Processing

Ga2O3 is being actively explored for power electronics, deep-ultraviolet optoelectronics, and other applications due to its ultra-wide bandgap and low projected fabrication cost of large-size and high-quality crystals. N-type doping of Ga2O3 can be achieved and tuned, but p-type doping faces fundamental obstacles due to deep character of acceptor levels and polaron transport of resulting holes. However, successful engineering of Ga2O3 based devices requires critical control of doping density, Fermi level position, and free carrier concentration, providing opportunities for predictive process simulation. We use first-principles defect theory and defect equilibrium calculations to simulate a 3-step growth-annealing-quench protocol for hydrogen assisted Mg doping in Ga2O3, taking into account the hydrogen-oxygen-water gas phase equilibrium. We predict type conversion to a net p-type regime following O-rich annealing after growth under reducing conditions in the presence of H2. This process is similar to the Mg acceptor activation by H removal in GaN. We show that there is an optimal temperature that maximizes the net acceptor density during the equilibrium annealing step for a given Mg doping level. Quenching of non-equilibrium annealed samples then results in a Fermi level EF below mid-gap down to about EV +1.5 eV, creating a significant number of uncompensated neutral MgGa0 acceptors. The resulting free hole concentration in quenched samples is very low (109 – 1012 cm-3) due to deep energy level of these Mg acceptors, but this type converted Ga2O3 material can create a significant built-in field in a p-n junction with an adjoining n-type material. Additionally, the electron concentration is greatly suppressed in such acceptor-doped Ga2O3, which could enable the use as a current blocking layer to fabricate normally-off (enhancement-mode) vertical Ga2O3 based metal-oxide-semiconductor field effect transistors (MOSFETs), and as a guard ring for edge termination in MOSFETs and Schottky barrier diodes (SBDs) with increased breakdown voltage.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗
DOE OSTI2021

Porous Covalent Organic Polymers for Efficient Fluorocarbon‐Based Adsorption Cooling

Abstract Adsorption‐based cooling is an energy‐efficient renewable‐energy technology that can be driven using low‐grade industrial waste heat and/or solar heat. Here, we report the first exploration of fluorocarbon adsorption using porous covalent organic polymers (COPs) for this cooling application. High fluorocarbon R134a equilibrium capacities and unique overall linear‐shaped isotherms are revealed for the materials, namely COP‐2 and COP‐3. The key role of mesoporous defects on this unusual adsorption behavior was demonstrated by molecular simulations based on atomistic defect‐containing models built for both porous COPs. Analysis of simulated R134a adsorption isotherms for various defect‐containing atomistic models of the COPs shows a direct correlation between higher fluorocarbon adsorption capacities and increasing pore volumes induced by defects. Combined with their high porosities, excellent reversibility, fast kinetics, and large operating window, these defect‐containing porous COPs are promising for adsorption‐based cooling applications.

Zheng, Jian↗
DOE OSTI2020

Defect Energetics in Pseudo-Cubic Mixed Halide Lead Perovskites from First-Principles

Owing to the increasing popularity of lead-based hybrid perovskites for photovoltaic (PV) applications, it is crucial to understand their defect energetics and its influence on their optoelectronic properties. In this work, we simulate various point defects in pseudocubic structures of mixed iodide-bromide and bromide-chloride methylammonium lead perovskites with the general formula MAPbI 3-y Br y or MAPbBr 3-y Cl y (where y is between 0 and 3), and use first-principles based density functional theory computations to study their relative formation energies and charge transition levels. We identify vacancy defects and Pb on MA antisite defect as the lowest energy native defects in each perovskite. Here, we observe that while the low energy defects in all MAPbI 3-y Br y systems only create shallow transition levels, the Br or Cl vacancy defects in the Cl-containing pervoskites have low energy and form deep levels which become deeper for higher Cl content. We examine the structures and density of states of pure and defect-containing perovskite systems to obtain an understanding of the nature of defect levels. Further, we study extrinsic substitution by different elements at the Pb site in MAPbBr 3 , MAPbCl 3 , and the 50-50 mixed halide perovskite, MAPbBr 1.5 Cl 1.5 , and identify some transition metals that create lower energy defects than the dominant intrinsic defects and also create midgap charge transition levels.

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

Point defect evolution under irradiation: Finite size effects and spatio-temporal correlations

The evolution of point defect concentrations under irradiation is controlled by their diffusion properties, and by their formation and elimination mechanisms. The latter includes the mutual recombination of vacancies and interstitials, and the elimination of point defects at sinks. Two models are traditionally used to predict the defect concentration evolution, the standard rate theory (SRT) and the atomistic kinetic Monte Carlo (AKMC). We show in this work a large discrepancy in the defect concentrations between both methods when the average number of defects in the AKMC simulation box is close or lower than one. The reason is that AKMC naturally captures strong space and time correlations between vacancies and interstitials generated by the finite size of the periodic simulation box. These correlations strongly affect the recombination rate and the point defect concentrations. SRT fails to predict such correlations, and the corresponding solution deviates from the more accurate solution given by the AKMC simulation under similar conditions. These finite size correlation effects are strong when the elimination of point defects occurs by recombination only, but can still be significant in the presence of sinks. In order to account for the spacio-temporal correlation in a continuum framework, we introduce a Correlated Pair Theory (CPT). This theory fully takes into account the correlations between vacancy and interstitial pairs and predicts point defect concentrations in good agreement with AKMC simulations. Inversely, if the goal is for the AKMC to reproduce defect concentrations in bulk using small simulation boxes, the naturally occurring correlations need to be corrected. We show here that the CPT can be used to modify the elimination rates in the AKMC simulations, so as to yield point defect concentrations in agreement with SRT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Effective parameterization of phase-field models of fission gas bubble growth

Fission gas bubbles are one of the most important microstructural features of ceramic nuclear fuels. As gas bubbles grow and interconnect, they allow release of gases, with important consequences for fuel performance. Phase-field modeling has been increasingly used to simulate the evolution of fission gas bubble microstructural because of its capability to capture complex microstructural features. However, computational performance limitations have made it difficult to simulate all the defects present in fuels during operation. For this reason, phase-field models have often simulated only vacancies and used multiple approaches to include the effect of vacancy-interstitial recombination and sinks in a simplified way. Here, we compare some of the most prevalent approaches, including source-only and source/sink. The kinetics of bubble growth using these approaches are analyzed analytically, and simulations with these approaches are compared to a full vacancy-interstitial model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Simulation-based methodology to assess the lattice defects creation as energy storing process

In nuclear industry, damages caused to in-core structural materials by high-energy neutrons have been studied for years. Recently, the possibility of storing energy through the creation of defects in crystalline materials via heavy particle bombardment was investigated. To evaluate the potential of the Wigner effect as energy-storing process, a graphite-moderated Molten Salt Reactor was adopted as a test case. Neutronics simulations performed with OpenMC Monte Carlo code allowed estimating the amount of energy released during fission events, the fraction that is deposited in graphite, and the fraction that has the potential to cause radiation damages. Additionally, molecular dynamics simulations were then performed to study the evolution of cascades triggered by high-energy neutrons. As a major outcome, a more accurate estimate of the radiation damage with respect to the traditional approaches was obtained, and the fraction of the deposited neutron energy that can be stored through long-lived defects was quantified.

25 ENERGY STORAGE↗
DOE OSTI2022

Analyzing the effect of pressure on the properties of point defects in $\gamma\mathrm{U–Mo}$ through atomistic simulations

Uranium–molybdenum (U–Mo) alloys in monolithic fuel foil are the primary candidate for the conversion of high-performance research reactors in the USA. Monolithic fuel is utilized in a plate-type design with a zirconium diffusion barrier and aluminum cladding. These fuel types are unique in that they contain no plenum for the release of fission gases, which, in conjunction with the aluminum cladding, can lead to large stress states within the fuel. The nature of how fundamental processes of radiation damage, including the evolution of point defects, under such stresses occur is unknown. In this work, we present molecular dynamics simulations of the formation energy of point defects under applied stress. We report this work will allow for the implementation of stress-dependent microstructural evolution models of nuclear fuels, including those for both fission gas bubble growth and creep, which are critical to ensure the stable and predictable behavior of research reactor fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Development of Large Scale Extrusion Deposition for Structural Applications

Large Scale Extrusion Deposition (LSED) is an evolving additive manufacturing (AM) technology that research, such as that taking place at Oak Ridge National Laboratory (ORNL) and companies like Local Motors, are continuing to utilize and develop new applications. A major LSED application of interest has been molds and tooling as it allows for much shorter production time and lower cost. However, interest has been expanding into using LSED for more structural applications due to more frequent use of high performing polymer composites as feedstock. The use of LSED for structural applications is of particular interest to Local Motors as it is currently being used to create commercially viable, energy efficient electric vehicles. LSED offers a unique opportunity when compared to traditional manufacturing methods as it can significantly reduce the number of components necessary, while decreasing embodied energy and carbon emissions. Even with these benefits however, it is important that LSED is properly understood from a structural aspect as this field has not been as heavily researched as tooling. To ensure a high level of safety and repeatability is an essential responsibility of a manufacturer whose products are structural in nature. As such it is important to understand the materials and LSED process to make structural objects that the manufacturer can be confident in. The main goals of this project were to: develop and investigate materials that are of interest for structural LSED applications, further develop and understand the current machines used in LSED and develop simulations tools of LSED and the mechanical properties of the created structure. Material development focused on composite materials that have high mechanical properties and are stable in a variety of environments. The materials were tested to determine their as printed mechanical and thermal properties as these are necessary for simulations. Once the materials were investigated thoroughly, it allowed for simulations to be performed to compare to experimental data with simulations. Materials were also vetted to determine candidates for multi-material printing. The machine development focused around the areas of process monitoring, non-destructive evaluation, and quality control. Finally, the goal of the simulations was to develop a realistic model of printed structures, including in-process simulation and dynamic simulation. The routes to get to some of these goals and the depth in which they were investigated changed throughout the project due to personnel changes and the COVID-19 pandemic. This project resulted in many valuable results such as the development of an nondestructive evaluation (NDE) technique for interlayer defects, proof of simulation for warpage in simple parts, the development and utilization of a profilometer to monitor a print for defects or inconsistencies, thorough investigation of a material used commercially for structural LSED applications and valuable experimental data on the applicability and advantages of multi-material crush structures versus their neat counterparts by creation and testing of samples.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A WISER way for simulating the performance of gratings

Soft X-ray monochromators for synchrotron radiation sources have been continuously developed over the years, improving energy resolution and stability. Great effort has been made in improving the surface quality of the optics involved, reaching values permitting diffraction-limited images. Still, one problem has not been solved, nor fully understood, yet: groove placing errors. Nowadays, these are one of the major factors responsible for resolving the power reduction of diffraction-gratings-based X-ray monochromators. Despite decades of use of gratings, there is not yet a universally established method for predicting and simulating the effect of groove placing errors on monochromator performance. This is especially important in the new scenario of high-coherent X-ray sources, i.e. diffraction-limited storage rings and free-electron lasers. To address this problem, in this article an approach based on WISER (Wavefront propagatIon Simulation codE libRary) is presented. WISER is a physical optics simulation package, also available in the user-friendly Orange Synchrotron Radiation Suite – OASYS. Even though it was originally conceived to assess the focusing performance of X-ray mirrors in the presence of height defects, it perfectly simulates the performance of a periodic (or quasi-periodic) structure like a diffraction grating. In this article, the way to use WISER and its application to a specific case, e.g. the design of a monochromator for the upgrade of a beamline at the Advanced Light Source, are shown. A simple rule for estimating how well the grooves are placed on a grating, based on calculation of the Strehl ratio, is also presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Direct observation of the formation and stabilization of metallic nanoparticles on carbon supports

Direct formation of ultra-small nanoparticles on carbon supports by rapid high temperature synthesis method offers new opportunities for scalable nanomanufacturing and the synthesis of stable multi-elemental nanoparticles. However, the underlying mechanisms affecting the dispersion and stability of nanoparticles on the supports during high temperature processing remain enigmatic. In this work, we report the observation of metallic nanoparticles formation and stabilization on carbon supports through in situ Joule heating method. We find that the formation of metallic nanoparticles is associated with the simultaneous phase transition of amorphous carbon to a highly defective turbostratic graphite (T-graphite). Molecular dynamic (MD) simulations suggest that the defective T-graphite provide numerous nucleation sites for the nanoparticles to form. Furthermore, the nanoparticles partially intercalate and take root on edge planes, leading to high binding energy on support. This interaction between nanoparticles and T-graphite substrate strengthens the anchoring and provides excellent thermal stability to the nanoparticles. These findings provide mechanistic understanding of rapid high temperature synthesis of metal nanoparticles on carbon supports and the origin of their stability.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Experiments in monthly mean simulation of the atmosphere with a coarse-mesh general circulation model

The Hansen atmospheric model was used to compute five monthly forecasts (October 1976 through February 1977). The comparison is based on an energetics analysis, meridional and vertical profiles, error statistics, and prognostic and observed mean maps. The monthly mean model simulations suffer from several defects. There is, in general, no skill in the simulation of the monthly mean sea-level pressure field, and only marginal skill is indicated for the 850 mb temperatures and 500 mb heights. The coarse-mesh model appears to generate a less satisfactory monthly mean simulation than the finer mesh GISS model.

Lutz, R. J.↗
DOE OSTI2021

Porous Covalent Organic Polymers for efficient Fluorocarbon-based Adsorption Cooling

Adsorption-based cooling is an energy-efficient renewable-energy technology that can be driven using low-grade industrial waste heat and/or solar heat. Here, we report the first exploration of fluorocarbon adsorption using porous covalent organic polymers (COPs) for this cooling application. High fluorocarbon R134a equilibrium capacities and unique overall linear-shaped isotherms are revealed for the materials, namely COP-2 and COP-3. The key role of mesoporous defects on this unusual adsorption behavior was demonstrated by molecular simulations based on atomistic defect-containing models built for both porous COPs. Analysis of simulated R134a adsorption isotherms for various defect-containing atomistic models of the COPs shows a direct correlation between higher fluorocarbon adsorption capacities and increasing pore volumes induced by defects. Combined with their high porosities, excellent reversibility, fast kinetics, and large operating window, these defect-containing porous COPs are promising for adsorption-based cooling applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2021

Hot Isostatic Pressing for Reduction of Defects in Additively Manufactured Inconel-718 Preliminary Simulation Approach and Results

Selective laser melting (SLM) of Inconel-718 (IN718) is a promising advanced manufacturing technology; however, challenges exist for fracture-critical parts where pores induced during the manufacturing process reduce fatigue life of the material. Hot isostatic pressing (HIP) is a common post processing technique used to reduce the part porosity by exposing parts to high temperature and pressures. Pore closure depends on the pore shape, trapped gas, and HIP parameters, and the process can only be optimized by understanding the effect of these variables. An analysis capability is described to better understand the conditions necessary for pore closure, and the effect of several variables on porosity are explored and compared to a baseline case. The importance of valid material response models is identified as a significant factor in model accuracy; especially the creep material response model which, in the baseline case, predicts that creep strain accounts for 87.7% of the total strain at the conclusion of the HIP cycle. Lessons learned are described, in addition to future work including model validation.

Joshua M Fody↗