Engineering Papers⌕ Search

Engineering topics

Devanathan, Ram

Publications and source records attributed to Devanathan, Ram.

Dynamic data-driven multiscale modeling for predicting the degradation of a 316L stainless steel nuclear cladding material

Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗

Monte Carlo Simulations of 347H Stainless Steel Aging for the Synthetic Generation of Microstructures Under Creep Conditions

Here, a Monte Carlo simulation method capable of replicating the kinetics of M 23 C 6 precipitation in 347H stainless steels was developed for the purpose of producing synthetic microstructures that approximate its microstructural evolution under aging periods of up to 10,000 hours at temperatures between 600 °C and 750 °C. To accomplish this, experimental data from the literature was used to parameterize simulations and replicate the nucleation and growth kinetics of M 23 C 6 particles within 347H and similar austenitic stainless steel alloys. These simulations were found to have considerable fidelity to previous efforts to study the precipitation of M 23 C 6 in other 300 series stainless steel alloys. Synthetic 347H microstructures were then generated that accounted the effects of aging temperature, duration, dislocation density, and the presence of boron within the microstructure. These simulations predict several key trends, those being that (1) the size of M 23 C 6 precipitates decreased with aging temperature and (2) the growth rate of M 23 C 6 particles decreased with aging temperature. Further, while (3) the addition of dislocation density due to creep conditions resulted in increasing intragranular nucleation of M 23 C 6 precipitates with increasing dislocation density and (4) B additions within the microstructure led to modest increases in precipitate size above 700 °C, which indicates that more complex physics are necessary to account for the presence of B.

36 MATERIALS SCIENCE↗

Finite Element Analysis and Machine Learning Guided Design of Carbon Fiber Organosheet-Based Battery Enclosures for Crashworthiness

Carbon fiber composite can be a potential candidate for replacing metal-based battery enclosures of current electric vehicles (E.V.s) owing to its better strength-to-weight ratio and corrosion resistance. However, the strength of carbon fiber-based structures depends on several parameters that should be carefully chosen. Here, in this work, we implemented high throughput finite element analysis (FEA) based thermoforming simulation to virtually manufacture the battery enclosure using different design and processing parameters. Subsequently, we performed virtual crash simulations to mimic a side pole crash to evaluate the crashworthiness of the battery enclosures. This high throughput crash simulation dataset was utilized to build predictive models to understand the crashworthiness of an unknown set. Our machine learning (ML) models showed excellent performance (R 2 > 0.97) in predicting the crashworthiness metrics, i.e., crush load efficiency, absorbed energy, intrusion, and maximum deceleration during a crash. We believe that this FEA-ML work framework will be helpful in down select process parameters for carbon fiber-based component design and can be transferrable to other manufacturing technologies.

36 MATERIALS SCIENCE↗

Chemical composition based machine learning model to predict defect formation in additive manufacturing

With a goal of exploiting additive manufacturing to improve the manufacturing of existing reactor materials, we developed a chemical composition-based machine learning model to predict the printability of any given alloy in laser powder bed fusion (L-PBF) using experimental data from peer-reviewed literature. We defined printability as the ability to avoid defects like cracking, balling, porosity, and lack of fusion, that are caused by thermal stresses (during solidification or liquation), molten pool disintegration into disconnected small beads or lack of heat input respectively. Our models predict the tendency of balling defect formation and porosity percentage for a given composition, under a given set of processing conditions. To predict the likelihood of balling defect, three models: a random forest classifier, a gradient boost regressor and a neural network were trained on a dataset containing both traditional alloys and high entropy alloys. The neural network model showed the highest accuracy of 92.3 % in predicting the balling defect formation. A random forest regressor, gradient boost regressor and neural network were trained and tested on a dataset of various alloys to predict porosity. The random forest regressor showed the best predictions with an R 2 score of 0.97. The models also revealed the relative importance of the input descriptors on defect-formation tendency. Of particular significance was the identification of carbon as an important element in determining the occurrence of balling and percent porosity in alloys like steel, as well as being moderately important to the percentage porosity in other alloys as well as steel. Manganese was also identified as a key descriptor for the percentage of porosity in steel and other alloys. Manganese’s low thermal conductivity and consistent presence in the dataset is the likely cause for its contribution. Carbon’s role is attributable to its relatively high specific heat and high melting temperature. In conclusion, our model serves as a swift, chemistry-based tool to design experiments and find modified compositions better suited for additive manufacturing.

36 MATERIALS SCIENCE↗

Cluster dynamics simulations of tritium and helium diffusion in lithium ceramics

Tritium (T) and He diffusion in LiAlO 2 and LiAl 5 O 8 phases influences the performance of tritium producing burnable absorber rods (TPBARs) by affecting the gas release, swelling and thermal conductivity of Li-bearing ceramic pellets. Frenkel pair defects and clusters created by irradiation can attract T and He interstitials and form clusters of the type He i x Li, He i x Al, He i x O, T i x Li, T i x Al, and T i x O, 1 ≤ x ≤ 4 in a Li, Al or O vacancy site (notation denotes x He or T atoms in a 1 Li, 1 Al or 1 O vacant site). The concentration and mobility of each of these clusters collectively contribute to the diffusion of the He and T gases in LiAlO 2 and LiAl 5 O 8 . Here, in this work, free energy cluster dynamics simulations implemented in the Centipede code, are used to obtain the concentration and diffusivities of these clusters which are then used to calculate the total diffusivity of T and He gases in LiAlO 2 and LiAl 5 O 8 . The results show that diffusivity of T is at least one order of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 whereas He diffusion is 2–13 orders of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 . There is a higher concentration of highly diffusive species (T interstitials and T i 03 Li for the case of tritium and He i 01 Li, He i 02 Li, and He i 03 Li for the case of He) in LiAlO 2 than in LiAl 5 O 8 which increase the total diffusion of T and He in LiAlO 2 .

36 MATERIALS SCIENCE↗

Molecular dynamics simulations of displacement cascades in LiAlO2 and LiAl5O8 ceramics

Abstract Molecular dynamics was employed to investigate the radiation damage due to collision cascades in LiAlO 2 and LiAl 5 O 8 , the latter being a secondary phase formed in the former during irradiation. Atomic displacement cascades were simulated by initiating primary knock-on atoms (PKA) with energy values = 5, 10 and 15 keV and the damage was quantified by the number of Frenkel pairs formed for each species: Li, Al and O. The primary challenges of modeling an ionic system with and without a core–shell model for oxygen atoms were addressed and new findings on the radiation resistance of these ceramics are presented. The working of a variable timestep function and the kinetics in the background of the simulations have been elaborated to highlight the novelty of the simulation approach. More importantly, the key results indicated that LiAlO 2 experiences much more radiation damage than LiAl 5 O 8 , where the number of Li Frenkel pairs in LiAlO 2 was 3–5 times higher than in LiAl 5 O 8 while the number of Frenkel pairs for Al and O in LiAlO 2 are ~ 2 times higher than in LiAl 5 O 8 . The primary reason is high displacement threshold energies (E d ) in LiAl 5 O 8 for Li cations. The greater E d for Li imparts higher resistance to damage during the collision cascade and thus inhibits amorphization in LiAl 5 O 8 . The presented results suggest that LiAl 5 O 8 is likely to maintain structural integrity better than LiAlO 2 in the irradiation conditions studied in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Automated Grain Boundary (GB) Segmentation and Microstructural Analysis in 347H Stainless Steel Using Deep Learning and Multimodal Microscopy

Austenitic 347H stainless steel offers superior mechanical properties and corrosion resistance required for extreme operating conditions such as high temperature. The change in microstructure due to composition and process variations is expected to impact material properties. Identifying microstructural features such as grain boundaries thus becomes an important task in the process-microstructure-properties loop. Applying convolutional neural network (CNN)-based deep learning models is a powerful technique to detect features from material micrographs in an automated manner. In contrast to microstructural classification, supervised CNN models for segmentation tasks require pixel-wise annotation labels. However, manual labeling of the images for the segmentation task poses a major bottleneck for generating training data and labels in a reliable and reproducible way within a reasonable timeframe. Microstructural characterization especially needs to be expedited for faster material discovery by changing alloy compositions. Here, in this study, we attempt to overcome such limitations by utilizing multimodal microscopy to generate labels directly instead of manual labeling. We combine scanning electron microscopy images of 347H stainless steel as training data and electron backscatter diffraction micrographs as pixel-wise labels for grain boundary detection as a semantic segmentation task. The viability of our method is evaluated by considering a set of deep CNN architectures. We demonstrate that despite producing instrumentation drift during data collection between two modes of microscopy, this method performs comparably to similar segmentation tasks that used manual labeling. Additionally, we find that naïve pixel-wise segmentation results in small gaps and missing boundaries in the predicted grain boundary map. By incorporating topological information during model training, the connectivity of the grain boundary network and segmentation performance is improved. Finally, our approach is validated by accurate computation on downstream tasks of predicting the underlying grain morphology distributions which are the ultimate quantities of interest for microstructural characterization.

36 MATERIALS SCIENCE↗

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗

Insights into radiation resistance of titanium alloys from displacement cascade simulations

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.

36 MATERIALS SCIENCE↗

Molecular Dynamics Simulation of Hygroscopic Aging Effects in Epoxy Polymer

The automobile industry is incorporating more lightweight content in car designs to boost fuel-economy. New structural adhesives are needed to mitigate the corrosion and thermal expansion issues associated with joining dissimilar lightweight materials, but adhesive developers lack a fundamental understanding of the chemistry that occurs in the adhesive as the joint ages. In this study, we developed structural adhesive molecular models and applied classical molecular dynamics simulations and density functional theory calculations to gain molecular insights into the influence of water molecules on the properties of epoxy-based adhesives (DGEBA + Jeffamine (JD230)). The simulations were complemented by experimental synthesis and characterization. Our work underscores the impact of water molecules on the local structure of the epoxy network as well as resulting mechanical properties. Water molecules were mainly coordinated with hydroxyls, primary amines and secondary amines, but also weakly coordinated with ether linkages, which were found most probable to be labile. Simulated stress–strain data indicates that increasing the water content deteriorates the mechanical properties. The Young’s modulus decreased by ~ 30% when the water content increased to 3 wt%. We conclude, this integration of molecular-level chemical insights with mechanical property simulations of the hydrated epoxy system and experimental validation holds the promise to advance lightweight joint technologies.

36 MATERIALS SCIENCE↗

Rapid discovery of high hardness multi-principal-element alloys using a generative adversarial network model

Multi-principal element alloys (MPEAs) continue to gain research prominence due to their promising high-temperature microstructural and mechanical properties. Recently, machine learning (ML) and materials informatics have been used extensively for screening MPEAs, however, most of these efforts were focused on constructing classification and regression models for predicting phase stability and mechanical properties of known compositions. These approaches may accelerate the screening process but optimizing new compositions with desirable properties within a practical time frame from an infinitely large design space of MPEA systems remains a grand challenge. To tackle this composition optimization challenge, a generative adversarial network coupled with a neural-network ML model was utilized to design MPEAs by filtering compositions that have high hardness. Even in a high-dimensional space with 18 elements as descriptors, the ML model was able to generate optimized compositions from which one composition was found to have 10% higher hardness (941 HV) than the maximum in the training data (857 HV). Density-functional theory was used to provide thermodynamic and electronic insights to higher hardness of the new MPEA found. The present work can optimize compositions from a wide design space of 18 elements (including W, Ta and Nb) that presents an opportunity to synthesize new compositions for applications ranging from corrosion-resistant alloys to nuclear materials. Here the findings suggest that generative ML can greatly accelerate materials discovery by identifying novel compositions, which can serve as a data-informed tool to guide experiments.

36 MATERIALS SCIENCE↗

Molecular dynamics study of primary damage in the near-surface region in nickel

In this work, we carried out a large-scale molecular dynamics (MD) simulation to elucidate the effect of the free surface on the defect production of displacement cascades in pure nickel. These MD simulations were performed both in the bulk and near-surface regions with primary knock-on atom (PKA) energies of E PKA = 1, 5, or 10 keV and at temperatures T = 300, 425, or 525 K. Additionally, for every (T, E PKA ), near-surface cascade simulations were performed as a function of depth. In the near-surface simulations, assuming isotropic neutron fluence, the PKAs were initiated in random directions, including toward and parallel to the free surface. For both the near-surface and bulk cascades, the effect of E PKA and T on the defect production is similar. In both cases, the defect production increases with E PKA , but the temperature has minimal effect. However, the production and clustering of vacancies are higher for near-surface cascades, and they decrease with increasing depth. In contrast, the production and clustering of self-interstitial atoms are lower and increase with depth. Eventually, the production and clustering of both defect types approach bulk cascade-like behavior, and the depth at which this occurs increases with E PKA , but is independent of temperature.

36 MATERIALS SCIENCE↗

Advanced Material Studies for High Intensity Proton Production Targets and Windows

The titanium alloy Ti6Al4V is widely used in accelerator facilities as beam windows, which are exposed to high cycle compressive stress waves from intense pulsed proton beams. Such materials interacting with the beam are subject to various forms of radiation damage, which can adversely affect their endurance limit. However, no fatigue data is currently available for high energy proton irradiated titanium alloy. Due to limitations in proton irradiation facilities, only miniature flat samples can be used for irradiation. To address this issue, we have developed a custom-made bend fatigue tester at Fermilab specifically for testing proton irradiated titanium alloy. In this presentation, we report on the fatigue test results obtained from this custom-fatigue tester using a non-ASTM standard specimen design. We plan to validate these sparse fatigue data with ASTM standard samples using standard fatigue samples. Furthermore, we have modified another commercial bend fatigue tester to accommodate miniature samples, and discuss some inherent deficiencies of the commercial fatigue tester to test miniature samples. To overcome this issue, a new fixture design is presented, which enables satisfactory fatigue testing on miniature samples over long periods. Finally, we present an upgrade to the custom-fatigue tester, featuring this new fixture design.

43 PARTICLE ACCELERATORS↗

Advanced Material Studies for High Intensity Proton Production Targets and Windows

The titanium alloy Ti6Al4V is widely used in accelerator facilities as beam windows, which are exposed to high cycle compressive stress waves from intense pulsed proton beams. Such materials interacting with the beam are subject to various forms of radiation damage, which can adversely affect their endurance limit. However, no fatigue data is currently available for high energy proton irradiated titanium alloy. Due to limitations in proton irradiation facilities, only miniature flat samples can be used for irradiation. To address this issue, we have developed a custom-made bend fatigue tester at Fermilab specifically for testing proton irradiated titanium alloy. In this presentation, we report on the fatigue test results obtained from this custom-fatigue tester using a non-ASTM standard specimen design. We plan to validate these sparse fatigue data with ASTM standard samples using standard fatigue samples. Furthermore, we have modified another commercial bend fatigue tester to accommodate miniature samples, and discuss some inherent deficiencies of the commercial fatigue tester to test miniature samples. To overcome this issue, a new fixture design is presented, which enables satisfactory fatigue testing on miniature samples over long periods. Finally, we present an upgrade to the custom-fatigue tester, featuring this new fixture design.

43 PARTICLE ACCELERATORS↗

Data Archive and Portal (DAP) Platform for Solid Phase Processing Technologies

The scale and speed of data generated by modern scientific experiments have constantly challenged the research community to store, curate, manage and optimally use it to drive scientific discoveries. In this work, we have developed a data archive and portal (DAP) platform including analytics capabilities to collect, curate, and manage data and metadata stream for solid phase processing (SPP) techniques. We successfully hosted around ~347K files of data related to processing parameters, microscopic images, and spectroscopic data related to solid phase processing. The DAP platform for SPP will establish an enduring capability to support machine learning and grow collaboration at the intersection of materials science and data science.

36 MATERIALS SCIENCE↗