SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Nearly thirty years since its inception, the combined finite-discrete element method (FDEM) has made remarkable strides in becoming a mainstream analysis tool within the field of Computational Mechanics. FDEM was developed to effectively “bridge the gap” between two disparate Computational Mechanics approaches known as the finite and discrete element meth-ods. At Los Alamos National Laboratory (LANL) researchers developed the Hybrid Optimization Software Suite (HOSS) as a hybrid multi-physics platform, based on FDEM, for the simulation of solid material behavior complemented with the latest technological enhancements for full fluid–solid interaction. Furthermore, in HOSS, several newly developed FDEM algorithms have been implemented that yield more accurate material deformation formulations, inter-particle interaction solvers, and fracture and fragmentation solutions. Additionally, an explicit computational fluid dynamics solver and a novel fluid–solid interaction algorithms have been fully integrated (as opposed to coupled) into the HOSS’ solid mechanical solver, allowing for the study of an even wider range of problems. Advancements such as this are leading HOSS to become a tool of choice for multi-physics problems. Finally, HOSS has been successfully applied by a myriad of researchers for analysis in rock mechanics, oil and gas industries, engineering application (structural, mechanical and biomedical engineering), mining, blast loading, high velocity impact, as well as seismic and acoustic analysis. This paper intends to summarize the latest development and application efforts for HOSS.
A plasma-sprayed aluminum-silicon (AlSi)/polyester coating is applied in modern gas turbine engines as an abradable sealcoating to maintain tight clearances between the rotating blades and the static casing. While running the engine, the rotating blades “rub” with the abradable coating, which results in extreme strain rate (up to 10 6 ) dynamics and a high-temperature environment. Due to the difficulty of collecting direct measurements, predictive computational models are important for analyzing the deformation and failure of the abradable material, to help meet the design target of avoiding damage to the blade tip and maintaining high fuel efficiency. In this research, a microstructure-based finite element (FE) computational model was developed to capture the complex mechanical behavior of the AlSi/polyester microstructure. The model is based on a virtual representative-volume-element (RVE) of a metal-polymer microstructure, reconstructed from x-ray computed tomography. It models the plastic deformation of, and damage to, each AlSi and polyester constituents, as well as the failure at their interface. The model was calibrated and validated with uniaxial tension and compression experiments, conducted at two temperatures (298 K and 533 K) at an applied strain rate of 10 3 - 10 4 s -1 . The material exhibited strongly asymmetric tension-compression behavior and a sensitivity to temperature, which was well captured by the model. The model was further applied to investigate changes in mechanical behavior due to variations in constituents’ volume fractions, which provides guidance to the microstructural design of AlSi/polyester abradable materials. The model is expected to facilitate the development of improved abradable materials by bypassing the conventional trial-and-error approach and extensive testing requirements.
Predicting the mechanical response and failure of heterogeneous materials has proven difficult. Multiscale numerical methods based in higher order continuum theories attempt to bridge the gap between microscale and macroscale structural behavior. Micromorphic continuum theories have shown promise. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling workflows. These include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This work discusses micromorphic upscaling efforts of the University of Colorado Boulder PSAAP III multidisciplinary simulation center (MSC). Verification studies are presented that compare the accuracy of the upscaling workflow with analytical solutions for a trivial stress state and homogeneous material using direct numerical simulations (DNS) conducted in the Ratel finite element method (FEM) and GEOS material point method (MPM) codes. These verification studies consider upscaling using the Micromorphic Filter for a “single filter domain”. Finally, DNS of a heterogeneous composite material is upscaled using a “multiple filter domain” method.
The Fusion Nuclear Science Facility (FNSF) is a proposed tokamak reactor with the mission to investigate operation of a fusion reactor in a nuclear environment. The high neutron fluence component of the FNSF mission requires steady-state operation for extremely long pulses (t_{pulse} ∼ months) at full power. Plasma sustainment and current drive will be critical components of a successful FNSF. COMSOL Multiphysics software is used for combined radiofrequency (RF) and thermal simulations of the lower hybrid current drive (LHCD) antenna system. These simu- lations consider the resistive RF losses in the antenna including realistic surface roughness and a range of potential materials. The thermal analysis adds volumetric nuclear heating, plasma heat flux on leading edges, and electromagnetic radiation from the plasma to the RF heating calculated by COMSOL. Additional neutronics calculations have been performed to determine the impact of these antenna designs on activated waste disposal for the materials considered. The simulations show that it is technically feasible to implement a fully-active multi-junction (FAM) rather than a passive-active multi-junction (PAM) style of antenna if the septum between adjacent waveguides is sufficiently wide and the thermal conductivity of the structural material is sufficiently high. The FAM has the benefit of higher achievable power density with respect to the PAM, which results in a more compact antenna with potentially lower impact on neutron shielding and tritium breeding. These considerations point to tungsten rather than steel as the preferred structural material in constructing the antenna.
This review discusses atomistic modeling techniques used to simulate radiation damage in crystalline materials. Radiation damage due to energetic particles results in the formation of defects. The subsequent evolution of these defects over multiple length and time scales requiring numerous simulations techniques to model the gamut of behaviors. This work focuses attention on current and new methodologies at the atomistic scale regarding the mechanisms of defect formation at the primary damage state.
In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.
Photovoltaic (PV) modules and batteries can either be recycled or disposed of in landfills at end-of-life (EoL). This work focuses on disposal since the benefit of recycling PV modules and batteries is well established. This study characterizes the potential toxicity due to metals leaching from selected PV modules and batteries through both the Toxicity Characteristic Leaching Procedure (TCLP) and in-house batch leaching protocols to probe the impact on metallic ion mobility as a function of the different types of e-waste entering the waste stream, the magnitude of device damage when placed into the waste stream, and the simulated municipal solid waste (SMSW) composition. Our results showed that one PV module and three battery types in this study should be classified as hazardous waste within the U.S. However, for some of the other e-wastes, metals of concern including Cr, Cu, Hg, Ni, Pb, and Zn leached during the batch tests but not in the TCLP regulatory method. For most of the waste types, the amounts of metals that leached in the TCLP test and batch tests were much lower than the total extractable amounts, demonstrating the potential for additional unaccounted for amounts of metals to leach. Our results demonstrate that the TCLP regulatory method might fail at predicting potential leaching and at capturing the complexity of e-waste leaching in landfill conditions. It confirms that additional work is needed urgently to develop appropriate EoL procedures for MSW with PV and battery e-waste.
Continuum-scale material deformation models, such as crystal plasticity (CP), can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e. mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a CP flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization (BO). We demonstrate a semi-autonomous BO loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model CP flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.
Not Available
Radiation damage in beam window materials limits the use of high power proton beams in high energy physics research. The alloy Ti-6Al-4V is presently used as a beam window material but a prospective alternative, Ti-15V-3Cr-3Sn-3Al has been proposed. Since both these alloys contain dual phases at room temperature, we compare the radiation damage in the α and ß-phases of these two materials via primary knock-on atom (PKA) cascade simulations in the 10-40keV energy range. At PKA energies 30 and 40keV, the number of Frenkel pairs in the ballistic stage is higher in the ß-phase of Ti-6Al-4V than that in the ß-phase of Ti-15V-3Cr-3Sn-3Al almost by a factor of 2. The α-phase, of both these alloys, by far outperforms the ß-phases of the two alloys, both in terms of damage during the ballistic stage and in terms of the surviving defects. The average displacement threshold energy (E d ) in the α-phase of both alloys was found to be 66eV while that in the ß-phase of Ti-15-3 was 55 eV and in the ß-phase of Ti-6-4 was 46. Further, while the number of surviving defects is almost equal in both alloys, the vacancy and interstitial clustering mechanisms differ notably, which can impact the degree of radiation hardening and loss of ductility. Our simulations show larger vacancy and interstitial clusters form in Ti-6Al-4V as compared to that in Ti-15-3-3-3 alloy. These results indicate that Ti-15-3-3-3 alloys may be a promising candidate for next generation beam window material with a higher radiation tolerance than the existing Ti-6-4 alloy.
Not provided.
Throughout U.S. Department of Energy (DOE) complexes, safety engineers employ the five-factor formula to calculate the source term (ST) that includes parameters of airborne release fraction (ARF), respirable fraction (RF) and damage ratio (DR). Limited experimental data on fragmentation of solids, such as ceramic pellets (i.e., PuO2), and container breach due to mechanical insults (i.e., drop and forklift impact), can be supplemented by modeling and simulation using high fidelity computational tools to estimate these parameters. This paper presents the use of Sandia National Laboratories' SIERRA solid mechanics (SM) finite element code to investigate the behavior of the widely utilized waste container (such as 7A Drum) subject to a range of free fall impact and puncture scenarios. The resulting behavior of the container is assessed, and the estimates are presented for bounding DRs from calculated breach areas for the various accident conditions considered. This paper also describes a novel multiscale constitutive model recently implemented in SIERRA/SM that simulates the fracture of brittle materials such as PuO 2 and determines ARF during the fracture process. Furthermore, comparisons are made between model predictions and simple bench-top experiments.
Understanding fuel system performance during anticipated transients without scram (ATWSs) in boiling water reactors (BWRs) is necessary for refining current and future safety limits. High-fidelity material models and simulations are fundamental to rigorous assessment of zirconium-based cladding performance. However, experimental thermomechanical data during simulated ATWSs to validate these modes are limited. To provide relevant in-situ data, Zircaloy-4 cladding was subjected to cyclic heating in a steam environment to simulate an out-of-pile BWR ATWS. Digital image correlation was used to capture the cladding strain behavior in-situ for comparison against simulations using the BISON finite element code. Conventional high-temperature models were compared using multiple schemes to gain a better understanding of the applicability of three BISON models to BWR ATWS: (1) the default combination of creep models in BISON, (2) the high-temperature Erbacher model alone, and (3) the low-temperature Limback-Andersson model alone. The cases run with the Limback-Andersson model alone produced the lowest root mean square error (RMSE). The lowest RMSE for the Limback-Andersson model alone was 0.659%, and the highest RMSE reported was 4.22%. A data gap within the model in the temperature regime of interest was also identified, and to account for this gap, the current model in BISON is linearly interpolated between two separate datasets. In conclusion, this evaluation highlights the need to either develop a new model or to improve the existing model to capture transient creep effects resulting from a cyclic temperature transient.
A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.
Nuclear explosive debris simulants are an important material in training and validating aspects of post-detonation nuclear forensic processes. Realistic simulants should replicate several aspects of nuclear explosive debris such as the size, shape, color, density, and chemical and radiological properties. Silica particles produced via sol-gel synthesis have recently been found to successfully reproduce many of these parameters including the controllable incorporation of radionuclide content. However, to be useful as a benchmarking material for validation and verification of laboratory methodologies, radionuclide content from batch-to-batch must be reproducible. Here, in this work, we explore the variance in radionuclide distribution incorporated into sol-gel benchmarking materials with respect to sample subdivision. Results will help inform the sample sizes required to minimize variance between samples.
The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.
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.