Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Materials 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 1,693 records · Page 94

Accelerating high-strain continuum-scale brittle fracture simulations with machine learning

Failure in brittle materials under dynamic loading conditions is a result of the propagation and coalescence of microcracks. Simulating this discrete crack evolution at the continuum level is computationally expensive or, in some cases, intractable, resulting in the need to make broad assumptions or neglect key physics. In this work, we have developed an approach using machine learning that overcomes the current inability to represent meso-scale physics at the macro-scale. Our approach leverages damage and stress data from a computationally expensive high-fidelity model that explicitly resolves microcrack behavior to build an inexpensive machine learning emulator. Once trained, the machine learning emulator is used to predict the evolution of crack length statistics, which then informs a continuum-scale constitutive model. This results in a significant speed-up of the workflow by four orders of magnitude. Both the machine learning emulator and the continuum-scale model are validated against the high-fidelity model and experimental data, respectively, showing excellent agreement. There are two key findings. The first is that we can reduce the dimensionality of the problem, establishing that the machine learning emulator only needs the length of the longest crack and one of the maximum stress components to capture the necessary physics. Another compelling finding is that the emulator can be trained in one experimental setting and transferred successfully to predict behavior in a different setting.

36 MATERIALS SCIENCE↗

Physical and mechanical metallurgy of zirconium alloys for nuclear applications: a multi-scale computational study

In the post-Fukushima world, the stability of materials under extreme conditions is an important issue for the safety of nuclear reactors. Because the nuclear industry is going to continue using advanced zirconium cladding materials in the foreseeable future, it become critical to gain fundamental understanding of the several interconnected problems. First, what are the thermodynamic and kinetic factors affecting the oxidation and hydrogen pick-up by these materials at normal, off-normal conditions, and in long-term storage? Secondly, what protective coatings (if any) could be used in order to gain extremely valuable time at off-normal conditions, e.g., when temperature exceeds the critical value of 2200°F? Thirdly, the kinetics of oxidation of such protective coating or braiding needs to be quantified. Lastly, even if some degree of success is achieved along this path, it is absolutely critical to have automated inspection algorithms allowing identifying defects of cladding as soon as possible. This work strives to explore these interconnected factors from the most advanced computational perspective, utilizing such modern techniques as first-principles atomistic simulations, computational thermodynamics of materials, diffusion modeling, and the morphological algorithms of image processing for defect identification. Consequently, it consists of the four parts dealing with these four problem areas preceded by the introduction and formulation of the studied problems. In the 1st part an effort was made to employ computational thermodynamics and ab initio calculations to shed light upon the different stages of oxidation of ziraloys (2 and 4), the role of microstructure optimization in increasing their thermal stability, and the process of hydrogen pick-up, both in normal working conditions and in long-term storage. The 2nd part deals with the need to understand the influence and respective roles of the two different plasticity mechanisms in Zr nuclear alloys: twinning (at low T) and crystallographic slip (higher T’s). For that goal, a description of the advanced plasticity model is outlined featuring the non-associated flow rule in hcp materials including Zr. The 3rd part describes the kinetic theory of oxidation of the several materials considered to be perspective coating materials for Zr alloys: SiC and ZrSiO 4 . In the 4th part novel and advanced projectional algorithms for defect identification in zircaloy coatings are described. In so doing, the author capitalized on some 12 years of his applied industrial research in this area. Our conclusions and recommendations are presented in the 5th part of this work, along with the list of used literature and the scripts for atomistic, thermodynamic, kinetic, and morphological computations.

36 MATERIALS SCIENCE↗

Materials technology for Stirling space power converters

This program was conducted in support of the NASA LeRC development of the Stirling power converter (SPC) for space power applications. The objectives of this contract were: (1) to perform a technology review and analyses to support the evaluation of materials issues for the SPC; (2) to evaluate liquid metal compatibility issues of the SPC; (3) to evaluate and define a transient liquid phase diffusion bonding (TLPDB) process for the SPC joints to the Udimet 720 heater head; and (4) to evaluate alternative (to the TLPDB) joining techniques. In the technology review, several aspects of the current Stirling design were examined including the power converter assembly process, materials joining, gas bearings, and heat exchangers. The supporting analyses included GLIMPS power converter simulation in support of the materials studies, and system level analysis in support of the technology review. The liquid metal compatibility study evaluated process parameters for use in the Stirling power converter. The alternative joining techniques study looked at the applicability of various joining techniques to the Stirling power converter requirements.

Baggenstoss, William↗

Real-Time Investigation of Solidification of Metal Matrix Composites

Casting of metal matrix composites can develop imperfections either as non- uniform distributions of the reinforcement phases or as outright defects such as porosity. The solidification process itself initiates these problems. To identify or rectify the problems, one must be able to detect and to study how they form. Until, recently this was only possible by experiments that employed transparent metal model organic materials with glass beads to simulate the reinforcing phases. Recent results obtained from a Space Shuttle experiment (using transparent materials) will be used to illustrate the fundamental physics that dictates the final distribution of agglomerates in a casting. We have further extended this real time investigation to aluminum alloys using X-ray microscopy. A variety of interface-particle interactions will be discussed and how they alter the final properties of the composite. A demonstration of how a solid-liquid interface is distorted by nearby voids or particles, particle pushing or engulfment by the interface, formations of wormholes, Aggregation of particles, and particle-induced segregation of alloying elements will be presented.

Kaukler, William↗

Machine learning for collective variable discovery and enhanced sampling in biomolecular simulation

Classical molecular dynamics simulates the time evolution of molecular systems through the phase space spanned by the positions and velocities of the constituent atoms. Molecular-level thermodynamic, kinetic, and structural data extracted from the resulting trajectories provide valuable information for the understanding, engineering, and design of biological and molecular materials. The cost of simulating many-body atomic systems makes simulations of large molecules prohibitively expensive, and the high-dimensionality of the resulting trajectories presents a challenge for analysis. Driven by advances in algorithms, hardware, and data availability, there has been a flare of interest in recent years in the applications of machine learning – especially deep learning – to molecular simulation. Furthermore, these techniques have demonstrated great power and flexibility in both extracting mechanistic understanding of the important nonlinear collective variables governing the dynamics of a molecular system, and in furnishing good low-dimensional system representations with which to perform enhanced sampling or develop long-timescale dynamical models. It is the purpose of this article to introduce the key machine learning approaches, describe how they are married with statistical mechanical theory into domain-specific tools, and detail applications of these approaches in understanding and accelerating biomolecular simulation.

74 ATOMIC AND MOLECULAR PHYSICS↗

In situ measurement of low-Z material coating thickness on high Z substrate for tokamaks

Rutherford backscattering of energetic particles can be used to determine the thickness of a coating of a low-Z material over a heavier substrate. Simulations indicate that 5 MeV alpha particles from an {sup 241}Am source can be used to measure the thickness of a Li coating on Mo tiles between 0.5 and 15 μm thick. Using a 0.1 mCi source, a thickness measurement can be accomplished in 2 h of counting. This technique could be used to measure any thin, low-Z material coating (up to 1 mg/cm{sup 2} thick) on a high-Z substrate, such as Be on W, B on Mo, or Li on Mo. By inserting a source and detector on a moveable probe, this technique could be used to provide an in situ measurement of the thickness of Li coating on NSTX-U Mo tiles. A test stand with an alpha source and an annular solid-state detector was used to investigate the measurable range of low-Z material thicknesses on Mo tiles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

BISON Fuel Performance Analysis of IFA-796 Rod 3 & 4 and Investigation of the Impact of Fuel Creep

In order to improve the accident tolerance of light water reactor (LWR) fuel, alternative cladding materials have been proposed to replace the currently used zirconium (Zr)-based alloys. Of these materials, there is a particular focus on iron-chromiumaluminum (FeCrAl) alloys because they exhibit slower oxidation kinetics in high-temperature steam than Zr-alloys. This should decrease the energy release due to oxidation and slow cladding consumption in the presence of high temperature steam. These alloys should also exhibit increased “coping time” in the event of an accident scenario by improving the mechanical performance at high temperatures, allowing greater flexibility to achieve core cooling. As a continuation of the development of these alloys, in-reactor irradiation testing of FeCrAl cladded fuel rods has started. In order to provide insight on the possible behavior of these fuel rods as they undergo irradiation in the Halden Boiling Water Reactor, engineering analysis has been performed using FeCrAl material models implemented into the BISON fuel performance code. This milestone report provides an update on the ongoing development of modeling capability to predict FeCrAl cladding fuel performance and to provide an early look at the possible behavior of planned in-reactor FeCrAl cladding experiments. In particular, this report consists of two separate analyses. The first analysis consists of fuel performance simulations of IFA-796 rod 4 and two segments of rod 3. These simulations utilize previously implemented material models for the C35M FeCrAl alloy and UO2 to provide a bounding behavior analysis corresponding to variation of the initial fuel cladding gap thickness within the fuel rod. The second analysis is an assessment of the fuel and cladding stress states after modification of the fuel creep model that is currently implemented in the BISON fuel performance code. Effects from modifying the fuel creep model were identified for the BISON simulations of the IFA-796 rod 4 experiment, but show that varying the creep model (within the range investigated here) only provide a minimal increase in the fuel radius and maximum cladding hoop stress. Continued investigation of fuel behavioral models will include benchmarking the modified fuel creep model against available experimental data, as well as an investigation of the role that fuel cracking will play in the compliance of the fuel. Correctly calculating stress evolution in the fuel is key to assessing fuel behavior up to gap closure and the subsequent deformation of the cladding due to PCMI. The inclusion of frictional contact should also be investigated to determine the axial elongation of the fuel rods for comparison with data from this experiment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Artificial neural network approach for multiphase segmentation of battery electrode nano-CT images

The segmentation of tomographic images of the battery electrode is a crucial processing step, which will have an additional impact on the results of material characterization and electrochemical simulation. However, manually labeling X-ray CT images (XCT) is time-consuming, and these XCT images are generally difficult to segment with histographical methods. We propose a deep learning approach with an asymmetrical depth encode-decoder convolutional neural network (CNN) for real-world battery material datasets. This network achieves high accuracy while requiring small amounts of labeled data and predicts a volume of billions voxel within few minutes. While applying supervised machine learning for segmenting real-world data, the ground truth is often absent. The results of segmentation are usually qualitatively justified by visual judgement. We try to unravel this fuzzy definition of segmentation quality by identifying the uncertainty due to the human bias diluted in the training data. Further CNN trainings using synthetic data show quantitative impact of such uncertainty on the determination of material’s properties. Nano-XCT datasets of various battery materials have been successfully segmented by training this neural network from scratch. We will also show that applying the transfer learning, which consists of reusing a well-trained network, can improve the accuracy of a similar dataset.

25 ENERGY STORAGE↗

Predicting Failure Progression and Failure Loads in Composite Open-Hole Tension Coupons

Failure types and failure loads in carbon-epoxy [45n/90n/-45n/0n]ms laminate coupons with central circular holes subjected to tensile load are simulated using progressive failure analysis (PFA) methodology. The progressive failure methodology is implemented using VUMAT subroutine within the ABAQUS(TradeMark)/Explicit nonlinear finite element code. The degradation model adopted in the present PFA methodology uses an instantaneous complete stress reduction (COSTR) approach to simulate damage at a material point when failure occurs. In-plane modeling parameters such as element size and shape are held constant in the finite element models, irrespective of laminate thickness and hole size, to predict failure loads and failure progression. Comparison to published test data indicates that this methodology accurately simulates brittle, pull-out and delamination failure types. The sensitivity of the failure progression and the failure load to analytical loading rates and solvers precision is demonstrated.

Arunkumar, Satyanarayana↗

Hypervelocity impact tests on Space Shuttle Orbiter thermal protection material

Hypervelocity impact tests were conducted to simulate the damage that meteoroids will produce in the Shuttle Orbiter leading edge structural subsystem material. The nature and extent of the damage is reported and the probability of encountering meteoroids with sufficient energy to produce such damage is discussed.

Humes, D. H.↗

Dielectric discharge characteristics in a two-electron simulation environment

The results of a series of experiments in which flux levels representative of the space electron environment were used are presented and compared to the results of high flux tests. The simulation approach was to partition the space electron spectrum into two parts, those electrons which do not penetrate a material and therefore contribute to charging and those which completely penetrate the material. The nonpenetrating electrons were simulated using 25 keV electrons and the penetrating electrons by 350 keV electrons. The materials included in this investigation were Kapton, optical solar reflectors (OSRs), and a ground test satellite surface potential monitor which contained Kapton, astroquart, OSRs and teflon.

Treadaway, M.↗

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

We report regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

36 MATERIALS SCIENCE↗

Computational Design and Experimental Validation of New Thermal Barrier Systems

This project (10/01/2010-9/30/2014), “Computational Design and Experimental Validation of New Thermal Barrier Systems”, originates from Louisiana State University (LSU) Mechanical Engineering Department and Southern University (SU) Department of Computer Science. This project will directly support the technical goals specified in DE-FOA-0000248, Topic Area 3: Turbine Materials, by addressing key technologies needed to enable the development of advanced turbines and turbine-based systems that will operate safely and efficiently using coal-derived synthesis gases. In this project, the focus is to develop and implement novel molecular dynamics method to improve the efficiency of simulation on novel TBC materials; perform high performance computing (HPC) on complex TBC structures to screen the most promising TBC compositions; perform material characterizations and oxidation/corrosion tests; and demonstrate our new thermal barrier coating (TBC) systems experimentally under integrated gasification combined cycle (IGCC) environments.

36 MATERIALS SCIENCE↗

Improved Axisymmetric and High Temperature Material Structural Modeling in MOOSE and NEML

This report describes improvements made to the solid mechanics formulation in the MOOSE open source finite element simulation environment and the open source Nuclear Engineering Material model Library (NEML) for mechanical constitutive models. The focus of these improvements is to improve the usability and performance of simulations involving one or both pieces of software. Specifically, this work completes a new system for solid mechanics simulations in the MOOSE ecosystem providing exact linearizations and optimal (quadratic) convergence, for a variety of coordinate systems and material types, including large deformation simulations. This work then provides users a framework to build highly efficient mechanical simulations of structures or materials or to couple in additional MOOSE physics modules to build complex, scalable multiphysics simulations.

36 MATERIALS SCIENCE↗

Thermochemical Ablation Analysis of the Orion Heatshield

The Orion Multi-Purpose Crew Vehicle will one day carry astronauts to the Moon and beyond, and Orion's heatshield is a critical component in ensuring their safe return to Earth. The Orion heatshield is the structural component responsible for absorbing the intense heating environment caused by re-entry to Earth's atmosphere. The heatshield is primarily composed of Avcoat, an ablative material that is consumed during the re-entry process. Ablation is primarily characterized by two processes: pyrolysis and recession. The decomposition of in-depth virgin material is known as pyrolysis. Recession occurs when the exposed surface of the heatshield reacts with the surrounding flow. The Orion heatshield design was changed from an individually filled Avcoat honeycomb to a molded block Avcoat design. The molded block Avcoat heatshield relies on an adhesive bond to keep it attached to the capsule. In some locations on the heatshield, the integrity of the adhesive bond cannot be verified. For these locations, a mechanical retention device was proposed. Avcoat ablation was modelled in CHAR and the in-depth virgin material temperatures were used in a Thermal Desktop model of the mechanical retention device. The retention device was analyzed and shown to cause a large increase in the maximum bondline temperature. In order to study the impact of individual ablation modelling parameters on the heatshield sizing process, a Monte Carlo simulation of the sizing process was proposed. The simulation will give the sensitivity of the ablation model to each of its input parameters. As part of the Monte Carlo simulation, statistical uncertainties on material properties were required for Avcoat. Several properties were difficult to acquire uncertainties for: the pyrolysis gas enthalpy, non-dimensional mass loss rate (B´c), and Arrhenius equation parameters. Variability in the elemental composition of Avcoat was used as the basis for determining the statistical uncertainty in pyrolysis gas enthalpy and B´c. A MATLAB program was developed to allow for faster, more accurate and automated computation of Arrhenius reaction parameters. These parameters are required for a material model to be used in the CHAR ablation analysis program. This MATLAB program, along with thermogravimetric analysis (TGA) data, was used to generate uncertainties on the Arrhenius parameters for Avcoat. In addition, the TGA fitting program was developed to provide Arrhenius parameters for the ablation model of the gap filler material, RTV silicone.

Sixel, William↗

Friction and wear of carbon-graphite materials for high-energy brakes

Caliper type brake simulation experiments were conducted on seven different carbon graphite materials formulations against a steel disk material and against a carbon graphite disk material. The effects of binder level, boron carbide (B4C) additions, SiC additions, graphite fiber additions, and graphite cloth reinforcement on friction and wear behavior were investigated. Reductions in binder level, additions of B4C, and additions of SiC each resulted in increased wear. The wear rate was not affected by the addition of graphite fibers. Transition to severe wear and high friction was observed in the case of graphite-cloth-reinforced carbon sliding against a disk of similar composition. The transition was related to the disruption of a continuous graphite shear film that must form on the sliding surfaces if low wear is to occur.

Bill, R. C.↗

Assessing the Impact of Point Defects on Molecular Diffusion in ZIF-8 Using Molecular Simulations

Because defects are ubiquitous in materials, they may play an important role in affecting the performance of materials in practical applications. In this paper, we investigate the influence of point defects on the diffusion of molecules including water, hydrocarbons, and acid gases in zeolitic imidazolate framework-8 (ZIF-8) using molecular simulations. To make these simulations possible, we introduce a force field that extends previous descriptions of pristine ZIF-8 to include experimentally relevant point defects. In general, the point defects we examined increase the local hopping rate for molecular diffusion, suggesting that low concentrations of these defects will not dominate long-range molecular diffusion in ZIF-8.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extending the density functional embedding theory to finite temperature and an efficient iterative method for solving for embedding potentials

A key element in the density functional embedding theory (DFET) is the embedding potential. We discuss two major issues related to the embedding potential: (1) its non-uniqueness and (2) the numerical difficulty for solving for it, especially for the spin-polarized systems. To resolve the first issue, we extend DFET to finite temperature: all quantities, such as the subsystem densities and the total system’s density, are calculated at a finite temperature. This is a physical extension since materials work at finite temperatures. We show that the embedding potential is strictly unique at T > 0. To resolve the second issue, we introduce an efficient iterative embedding potential solver. We discuss how to relax the magnetic moments in subsystems and how to equilibrate the chemical potentials across subsystems. The solver is robust and efficient for several non-trivial examples, in all of which good quality spin-polarized embedding potentials were obtained. We also demonstrate the solver on an extended periodic system: iron body-centered cubic (110) surface, which is related to the modeling of the heterogeneous catalysis involving iron, such as the Fischer-Tropsch and the Haber processes. This work would make it efficient and accurate to perform embedding simulations of some challenging material problems, such as the heterogeneous catalysis and the defects of complicated spin configurations in electronic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗