Search NASA⌕ Search

Engineering topics

Simunovic, Srdjan

Publications and source records attributed to Simunovic, Srdjan.

At least 19 records

Integrated top-down process and voxel-based microstructure modeling for Ti-6Al-4V in laser wire direct energy deposition process

Laser-wire metal additive manufacturing (AM) is one of the ideal direct energy deposition (DED) processes for creating large-scale parts with a medium level of complexity. However, the DED process involves complex thermal signatures and wide length scales making the fabrication of realistic AM components and part qualification often reliant on experimental trial-and-error optimization. While experimental measurements over the full volume of a part are valuable and necessary, measuring the entire area of a part is significantly laborious and practically infeasible, particularly for large parts in terms of cost and rapid qualification. Therefore, in this work, we developed an effective thermal and microstructure modeling framework based on the Johnson–Mehl-Avrami-Kolmogorov (JMAK) and Koistinen & Marburger (KM) models through a top-down approach that considers plate distortion-affected thermal profiles. A voxel-by-voxel simulation method is used to predict individual phase fractions of Ti-6Al-4 V. The predicted results were validated through detailed metallurgical measurements. A combined voxel-by-voxel approach with a sparse data reconstruction technique produced a near-perfect reconstruction of the original data. This approach anticipates a significant reduction in data points and computation time and resources. Lastly, we conclude with potential extensions of this work to other modeling efforts.

36 MATERIALS SCIENCE↗

Analysis and Mitigation of Cascading Failures Using a Stochastic Interaction Graph with Eigen-analysis

In studies on complex network systems using graph theory, eigen-analysis is typically performed on an undirected graph model of the network. However, when analyzing cascading failures in a power system, the interactions among failures suggest the need for a directed graph beyond the topology of the power system to model directions of failure propagation. To accurately quantify failure interactions for effective mitigation strategies, this paper proposes a stochastic interaction graph model and associated eigen-analysis. Different types of modes on failure propagations are defined and characterized by the eigenvalues of a stochastic interaction matrix, whose absolute values are unity, zero, or in between. Finding and interpreting these modes helps identify the probable patterns of failure propagation, either local or widespread, and the participating components based on eigenvectors. Then, by lowering the failure probabilities of critical components highly participating in a mode of widespread failures, cascading can be mitigated. Here, the validity of the proposed stochastic interaction graph model, eigen-analysis and the resulting mitigation strategies is demonstrated using simulated cascading failure data on an NPCC 140-bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Comparing Interaction Graphs on Cascading Outages under Different Loading Conditions

Interaction graphs on cascading outages of power systems provide valuable insights into how cascading outages evolve and propagate, and which components and links are critical to the propagation of cascading outages, enabling further development of mitigation strategies to support decision-making. However, the sensitivity of the interaction graph’s topology to the system’s loading condition has not been studied sufficiently. This paper compared interaction graphs under various loading conditions on the Northeastern Power Coordinating Council 140-bus system, and discovers the strong relationships between the graph topology, the cascade size distribution, and the load condition. Accordingly, three representative interaction graphs are constructed and illustrated.

Guo, Zhenping↗

Analysis and Mitigation of Cascading Outages Using an Interaction Graph Addressing Transient Stability

Cascading outages of power systems pose great threats to system security and reliability, potentially leading to large-scale blackouts. For analysis and mitigation of cascading outages, this paper proposes a transient stability-incorporated interaction graph. This graph statistically quantifies the interactions among line outages and instabilities of generators, which can model propagation paths and patterns of cascading outages. Compared with an interaction graph that only models line outages, this new interaction graph provides important insights on how transient instability occurs along with cascading outages. It also offers more effective strategies for mitigating outage propagation. The proposed interaction graph can be constructed from datasets of historical or simulated cascading events. It is demonstrated on an NPCC 140-bus system with mitigation strategies.

Guo, Zhenping↗

Material Model Parameters Optimization in Liquid Mercury Target Dynamics Simulation With Machine Learning Surrogates

A pulsed spallation target is subjected to very short (∼0.7μs) but intense loads (23.3 kJ) from repeated proton pulses, which knock away neutrons from the mercury atoms’ nuclei for a wide range application in physics, engineering, medicine, petroleum exploration, biology, chemistry, etc. The effect of this pulsed loading on the stainless-steel target module which contains the flowing mercury target material is difficult to predict not only due to its short but intense explosive-like physical reaction, but also the nonlinear material behavior of the liquid mercury in the structure. Injecting small helium bubbles in the mercury has been an efficient method of mitigating the pressure wave at high power level stage. However, prediction of the resultant loading on the target is more difficult when helium gas is intentionally injected into the mercury. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. A parameter sensitivity study was firstly employed to understand their impact on the simulation strains. The investigated parameters included E, μ, γ, σ, n, VFgas, and gas cumulative volume curve control parameters a and b. Verification and validation results from sparse polynomial expansions (SPE) method and directional Gaussian smoothing (DGS) optimization show that the surrogate model had training error of ∼7% and validation error of ∼15%, indicating that machine learning methods and surrogate models can help optimize the uncertain parameters in the complex 2-phase material model. This approach is expected to fill the knowledge gap between unknown liquid-gas mixture material model and measured vessel strain responses.

Lin, Lianshan↗

A Review on Simulation Models of Cascading Failures in Power Systems

Among various power system disturbances, cascading failures are considered the most serious and extreme threats to grid operations, potentially leading to significant stability issues or even widespread power blackouts. Simulating power systems' behaviors during cascading failures is of great importance to comprehend how failures originate and propagate, as well as to develop effective preventive and mitigative control strategies. The intricate mechanism of cascading failures, characterized by multi-timescale dynamics, presents exceptional challenges for their simulations. This paper provides a comprehensive review of simulation models for cascading failures, providing a systematic categorization and a comparison of these models. The challenges and potential research directions for the future are also discussed.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Accelerating Battery Simulations by Using High Performance Computing and Opportunities with Machine Learning

The US Department of Energy (DOE) estimates that the battery pack cost of $60 per kilowatt-hour while increasing the driving range to over 300 miles and vehicle charging under 15 min or less would enable mass penetration of electric vehicles in the USA by 2030. These projections are based on the currently available high-density cell chemistry combined with a system level design and optimized electrodes. Electrodes for current state-of-the-art lithium-ion cell technology are fabricated from electrode slurry comprising of active material, polymeric binder, and conductive diluent such as carbon black that are coated on metal current collectors such as copper and aluminum. Given these challenging requirements for development of electrical energy storage devices for future transportation needs, a predictive simulation capability which can accelerate design by considering performance and safety implications of different geometry, materials, and chemistry choices is required. In this chapter, we discuss our state-of-the-art three-dimensional modeling framework, providing examples on battery performance and safety simulations. We also present approaches to use machine learning in the context of large datasets that are becoming available.

Allu, Srikanth↗

Lithium-Ion Batteries with Safer Current Collectors

The US Department of Energy’s Oak Ridge National Laboratory and Soteria Battery Innovation Group collaborated to develop a metallized polymer film as a current collector for lithium-ion batteries and design slit patterns for electrodes to improve battery safety. The metallized polymer film serves as a fuse that will break under high temperature resulting from a short circuit. Consequently, the short circuit will be broken and the heat generated from the short circuit will be alleviated to avoid thermal runaway. The electrodes with slit patterns are expected to be broken into small segments upon mechanical impact, which can isolate the damaged electrodes from the rest. Thus, the heat generation associated with the mechanical impact will be reduced, along with the likelihood of thermal runaway. This effort aimed to improve lithium-ion battery safety by replacing traditional metal foils with metallized polymer films as current collector and introducing slit patterns to battery electrodes. The metallized polymer films were polyethylene terephthalate with a thin metal (aluminum for cathodes and copper for anodes) coating on both sides. Slit patterns with various geometries and dimensions were designed via simulation. Temperature distribution on cells with electrodes with and without slit patterns was also simulated via nemerical modeling. Electrodes on traditional metal foils and metallized polymer films were coated via a pilot-scale slot-die coater at the US Department of Energy Battery Manufacturing Facility at Oak Ridge National Laboratory. The electrodes were calendered to 35% porosity and punched out with and without slit patterns. Pouch cells were assembed with the electrodes inside a dry room for electrochemical and safety testing. Electrodes coated on metallized polymer films exhibited comparable electrochemical performance to their counterparts coated on traditional metal foils. They also showed improved safety in nail penetration and indentation tests.

25 ENERGY STORAGE↗

Influence of geometry on columnar to equiaxed transition during electron beam powder bed fusion of IN718

Correlation between spot-melt scan parameters (linear spot-density aka areal energy density), build geometry, and solidification microstructure evolution (columnar vs equiaxed) in a powder bed fusion technology is investigated. It is shown that to maintain the equiaxed solidification microstructure evolution in electron beam powder bed additive manufacturing (AM), the areal energy density per layer needs to be scaled with respect to the 2D cross-sectional area of the layer being melted. Samples with two different cross-sectional areas (40 × 40 mm and 20 × 20 mm) have been fabricated with varying areal energy densities. For a given square cross-section (20 × 20 mm), increasing the areal energy density (4.8 MJ/sq.m to 14.7 MJ/sq.m) transitioned the solidification microstructure from columnar to equiaxed. The observed microstructure data (Electron Back Scattered Diffraction - EBSD) is quantified by calculating the principal component (PC) score using a spatial statistics methodology. The sample with equiaxed grains is found to have a low PC score while the sample with columnar grain had a high PC score. A semi-analytical model is used to simulate the heat transfer and the local solidification conditions as a function of processing parameters (linear spot-density). The result from the heat transfer model is correlated with previously quantified microstructure data. Space-Time analysis of the melt pattern is done and correlated with the observed microstructure. In addition, from the findings, appropriate parameters have been used to additively manufacture a turbine blade with site-specific or hybrid solidification microstructure (traditional fabrication possible via a patented method of localized cold working and heat treatment).

36 MATERIALS SCIENCE↗

Implementation of Disruptive Designs for Gas Turbine Components Utilizing Additive Manufacturing

Electron beam melting (EBM) is an additive manufacturing process used to print metallic components with complex geometries for aerospace applications. A preheating characteristic enables this technology to avoid large accumulation of residual stress during printing. However, cracking mechanisms particularly in Ni-based superalloys are associated with not only residual stress but also with microstructure (i.e., solidification morphology, transient evolution of γ'). Siemens is committed to AM technology for gas turbine applications. The ability to novel predictive designs of turbine blades would reduce process development time and allow acquisition of key microstructure and performance data early in a design campaign. Siemens initiated a project to establish the design rules for the fabrication of non-printable Nibased superalloys. This phase 1 program demonstrates the feasibility for virtual fabrication of the actual blade geometries through EBM process.

42 ENGINEERING↗

Powder spreading, densification, and part deformation in binder jetting additive manufacturing

Binder jetting additive manufacturing (AM) can print complex structures in economical and scalable manner. Binder jetting AM comprises of deposition and weak binding of particles, known as green part, at room temperature and subsequent binder removal and sintering densification at high temperatures. However, during the densification (i.e., sintering), the part significantly deforms due to volume shrinkage. The deformation during sintering is difficult to predict, which prevents the widespread application of this technology. In this research, powder spreading simulation using discrete element method (DEM) was performed first to capture local variations in powder bed configuration. Second, finite element method (FEM) with a phenomenological continuum constitutive model was used to predict part shrinkage during the sintering process. DEM simulation showed variations in packing density, particle segregation, formation of uneven powder bed surface, and shift in particle size distribution (PSD). The sintering simulation modeled part deformation with a reasonable accuracy of 3% for solid-state sintering and intermediate liquid phase sintering. A demonstration case with non-uniform initial packing density showed that inhomogeneous green part density and PSD should be accounted for prediction of part deformation in binder jetting AM.

36 MATERIALS SCIENCE↗

Thermochemically-informed mass transport model for interdiffusion of U and Zr in irradiated U-Pu-Zr fuel with fission products

A new formulation for redistribution of constituents in metallic U-Pu-Zr nuclear fuel is presented that can incorporate the contributions from fission products. The formulation is based on the thermodynamic driving forces derived from the generalized chemical potential that includes effects of composition and temperature. As a result, the redistribution model can readily account for the composition changes due to the generation of fission products while using only a limited set of transport coefficients. The thermodynamic model for the metallic fuel and thermochemistry solver Thermochimica were coupled with the nuclear fuel performance code BISON to implement the redistribution model. The simulations reproduce an experimentally observed Zr-depletion zone in the mid-radius region of an irradiated fuel slug. The generation of fission products during burnup is shown to have a stabilizing effect on the fuel chemistry, slowing the rate of U-Zr interdiffusion and reducing the size of the Zr-depleted zone. Here, the overall objective of the work is to advance predictive capabilities of fission product behaviour in U-Zr metallic fuel in the context of fuel performance and safety.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗