Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random alloys”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Three-dimensional atomic structure and local chemical order of medium- and high-entropy nanoalloys

Medium- and high-entropy alloys (M/HEAs) mix several principal elements with near-equiatomic composition and represent a model-shift strategy for designing previously unknown materials in metallurgy, catalysis and other fields. One of the core hypotheses of M/HEAs is lattice distortion, which has been investigated by different numerical and experimental techniques. However, determining the three-dimensional (3D) lattice distortion in M/HEAs remains a challenge. Moreover, the presumed random elemental mixing in M/HEAs has been questioned by X-ray and neutron studies, atomistic simulations, energy dispersive spectroscopy and electron diffraction, which suggest the existence of local chemical order in M/HEAs. However, direct experimental observation of the 3D local chemical order has been difficult because energy dispersive spectroscopy integrates the composition of atomic columns along the zone axes and diffuse electron reflections may originate from planar defects instead of local chemical order. Here, in this work, we determine the 3D atomic positions of M/HEA nanoparticles using atomic electron tomography and quantitatively characterize the local lattice distortion, strain tensor, twin boundaries, dislocation cores and chemical short-range order (CSRO). We find that the high-entropy alloys have larger local lattice distortion and more heterogeneous strain than the medium-entropy alloys and that strain is correlated to CSRO. We also observe CSRO-mediated twinning in the medium-entropy alloys, that is, twinning occurs in energetically unfavoured CSRO regions but not in energetically favoured CSRO ones, which represents, to our knowledge, the first experimental observation of correlating local chemical order with structural defects in any material. We expect that this work will not only expand our fundamental understanding of this important class of materials but also provide the foundation for tailoring M/HEA properties through engineering lattice distortion and local chemical order.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First-principles study on the mechanical and thermodynamic properties of MoNbTaTiW

Refractory high-entropy alloys (RHEAs) are emerging as new materials for high temperature structural applications because of their stable mechanical and thermal properties at temperatures higher than 2273 K. In this study, the mechanical properties of MoNbTaTiW REDEA are examined by applying calculations based on first-principles density functional theory (DFT) and using a large unit cell with 100 randomized atoms. The phase calculation of MoNbTaTiW with CALPHAD method shows the existence of a stable body-centered cubic structure at a high temperature and a hexagonal closely packed phase at a low temperature. The predicted phase, shear modulus, Young’s modulus, Poisson’s ratio, and hardness values are consistent with available experimental results. The linear thermal expansion coefficient, vibrational entropy, and vibrational heat capacity of MoNbTaTiW RHEA are investigated in accordance with Debye-Grüneisen theory. These results may provide a basis for future research related to the application of RHEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Texture Memory in Hexagonal Metals and Its Mechanism

Texture memory is a phenomenon in which retention of initial textures occurs after a complete cycle of forward and backward transformations, and it occurs in various phase-transforming materials including cubic and hexagonal metals such as steels and Ti and Zr alloys. Texture memory is known to be caused by the phenomena called variant selection, in which some of the allowed child orientations in an orientation relationship between the parent and child phases are preferentially selected. Without such variant selection, the phase transformations would randomize preferred orientations. In this article, the methods of prediction of texture memory and mechanisms of variant selections in hexagonal metals are explored. The prediction method using harmonic expansion of orientation distribution functions with the variant selection in which the Burgers orientation relationship, {110}β//{0001} α-hex <1$\bar{1}$1>β//$\langle2\bar{11}0\rangle$ α-hex , is held with two or more adjacent parent grains at the same time, called “double Burgers orientation relation (DBOR)”, is introduced. This method is shown to be a powerful tool by which to analyze texture memory and ultimately provide predictive capabilities for texture changes during phase transformations. Variation in nucleation and growth rates on special boundaries and an extensive growth of selected variants are also described. Analysis of textures of commercially pure Ti observed in situ by pulsed neutron diffraction reveals that the texture memory in CP-Ti is indeed quite well predicted by consideration of the mechanism of DBOR. The analysis also suggests that the nucleation and growth rates on the special boundary of 90° rotation about $\langle2\bar{11}0\rangle$ α-hex should be about three times larger than those of the other special boundaries, and the selected variants should grow extensively into not only one parent grain but also other grains in α-hex(hexagonal)→β(bcc) transformation. The model calculations of texture development during two consecutive cycles of α-hex→β→α-hex transformation in CP-Ti and Zr are also shown.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Leveraging Artificial Intelligence to Predict Novel Eutectic Alloys

The goal of this project was to train an artificial neural network (ANN) to predict the fractional composition and melting point of eutectic alloys using fundamental atomic properties as inputs. The fundamental properties considered include atomic number, atomic weight, atomic radius, valence electron concentration, electronegativity, and electron affinity. The project involved several phases, starting with data preparation, where phase diagram data was harvested from the ASM International database. Approximately 1300 binary eutectics were collected and cleaned to ensure relevance and accuracy. A regression model was selected for training, utilizing a rectified linear unit as the activation function. Various model configurations were evaluated for predictive accuracy, with validation techniques employed to ensure robustness. The model demonstrated predictive capabilities above random guessing and was able to achieve up to 11% accuracy under certain conditions. An ablative test identified atomic radius and valence electron concentration as critical inputs for model performance. Incorporating the melting point of atomic constituents improved accuracy significantly, although ultimately the model’s predictive capability still fell short of the 80% target. This report details the methodology, results, and implications of the research, contributing to the understanding of employing artificial intelligence to predict the phase transition behavior of eutectic alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Energy and mechanical properties predictions in Fe-Ni binary system by ab initio calculations

Here, the fcc random solid solutions in Fe-Ni binary system were investigated using ab initio calculations. Density functional theory (DFT) was applied along with supercells generated by the Special Quasi-Random Structure (SQS) method to predict the enthalpy of mixing and magnetic and elastic properties. The results were compared with relevant experimental and computational data produced by other groups. The effect of SQS supercell size and relaxation procedure on the predicted properties of the alloys were also studied. Our calculations show that 16-atoms SQS model provides adequate approach to the mixing energy of the alloy. However, calculation of mechanical properties requires larger supercell. Excellent agreement between experimental and calculated elastic properties was found in fcc region of Fe-Ni system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Application of Machine Learning and Data Augmentation Algorithms in the Discovery of Metal Hydrides for Hydrogen Storage

The development of efficient and sustainable hydrogen storage materials is a key challenge for realizing hydrogen as a clean and flexible energy carrier. Among various options, metal hydrides offer high volumetric storage density and operational safety, yet their application is limited by thermodynamic, kinetic, and compositional constraints. In this work, we investigate the potential of machine learning (ML) to predict key thermodynamic properties—equilibrium plateau pressure, enthalpy, and entropy of hydride formation—based solely on alloy composition using Magpie-generated descriptors. We significantly expand an existing experimental dataset from ~400 to 806 entries and assess the impact of dataset size and data augmentation, using the PADRE algorithm, on model performance. Models including Support Vector Machines and Gradient Boosted Random Forests were trained and optimized via grid search and cross-validation. Results show a marked improvement in predictive accuracy with increased dataset size, while data augmentation benefits are limited to smaller datasets and do not improve accuracy in underrepresented pressure regimes. Furthermore, clustering and cross-validation analyses highlight the limited generalizability of models across different material classes, though high accuracy is achieved when training and testing within a single hydride family (e.g., AB2). The study demonstrates the viability and limitations of ML for accelerating hydride discovery, emphasizing the importance of dataset diversity and representation for robust property prediction.

augmentation↗
DOE OSTI2024

Insights into Prismatic Loop Formation in Irradiated Fe–Cr Alloys from Hypothesis-Driven Active Learning and Causal Analysis

Neutron and electron irradiation experimental studies conducted on body-centered cubic Fe and Fe–Cr alloys have established two prismatic dislocation loop populations, which have Burgers vectors of either a/2$\langle$111$\rangle$ or a$\langle$100$\rangle$. Here, the loop formation depends on factors such as dose (D), dose rate (D rt ), temperature (T), chromium content (Cr%), and other alloying elements. Hence, it is important to understand how irradiation-induced dislocation loops evolve conditional upon the loop characteristics, such as loop density (DD), average loop size d̅, and irradiation parameters (D, D rt , T, and irradiation type), which is still an active area of research. To understand these complex structure–property relationships, machine learning (ML) is employed in a three-step approach. This includes imputing missing data with a k-nearest neighbor, generating functionalized features, and assessing feature importance with random forest classification and regression. Physics-based features are incorporated in a hypothesis-driven active learning scheme to overcome data unavailability challenges. Insights obtained from ML models (i) to categorize dislocation loop types, show the highest correlation with d̅; (ii) Log(DD), obtained through mathematical formulations involving D, Cr%, d̅, and T (e.g., Log(DD) ~ D + exp(-Cr%) + 1/d̅ and log(DD) ~ D + exp(-Cr%) + 1/T). Hypothesis-driven active learning is able to predict Log(DD) in which the experimental date is not known. Causal models verify cause–effect relationships for dislocation loop classification and irradiation factors in FeCr alloys.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The Novel Hybrid Ab Initio Model of High-Performance Structural Alloys Design for Fossil Energy Power Plants

The current project developed a novel HT-CALPHAD/DFT approach, which can quickly design new high-performance structural alloys for the application of FE power plants. The PI will mainly take charge of high-throughput DFT simulations and computational thermodynamics of the selected multicomponent alloy systems for the FE power plant applications. At the end of the project, a novel hybrid model based on high throughput CALPHAD/DFT simulations and computational thermodynamics will be developed to provide guidance on how to identify multi-component new high-performance structural alloys with much less computational effort needed. It will address the extensive computation time needed for DFT on the new alloys design. In addition, it will also address the well-known headache of DFT, i.e. how to make the accurate prediction of the high-temperature equilibria. The novel hybrid model the PI proposed will not only be applied to the design of high-performance structural alloys in FE power plants but in many different applications, such as nuclear reactors. This hybrid modeling approach includes four sections: 1) HT-CALPHAD modeling of Al-Co-Cr-Ni-Fe system with FCC and BCC phase. In this section, 3561 non-equiatomic compositions were randomly generated in order to investigate the phase stability of single FCC and BCC phases. Meanwhile, we proposed a data screening procedure to screen out the good candidates within these compositions, considering the temperature range, average density, and melting temperature, etc. 2) Investigation of FCC-Cr lattice stability in Fe-Cr system. We systematically assessed the reliability of FCC-Cr lattice stability derived by DFT and CALPHAD approaches. Meanwhile, the Fe-Cr binary system was chosen as a case study to verify the Cr lattice stability obtained by both approaches. 3) High-throughput DFT modeling on elastic properties of Al-Co-Cr-Ni-Fe systems. We predicted and established the FCC quinary elastic constant database of the Al-Co-Cr-Fe-Ni systems at 0K by using special quasi-random structure (SQS) approach. The predictions will start with pure elements of Al-Co-Cr-Fe-Ni system and will be continued with binaries, ternaries, quaternaries, and finally the quinary compositions. 3) Modeling of temperature-dependent elastic properties in Al-Co-Cr-Ni-Fe systems. In this part, we predict the thermal expansion coefficient and elastic stiffness coefficient as a function of temperature by applying quasiharmonic approximation. With this approach, the elastic properties of HEAs at elevate temperature can be estimated.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Three-dimensional morphology of an ultrafine Al-Si eutectic produced via laser rapid solidification

Al-Si alloys processed by laser rapid solidification yield eutectic microstructures with ultrafine and interconnected fibers. Such fibrous structures have long been thought to bear resemblance to those formed in impurity-doped alloys upon conventional casting. Here, we show that any similarity is purely superficial. By harnessing high-throughput characterization and computer vision techniques, we perform a three-dimensional analysis of the branching behavior of the ultrafine eutectic and compare it against an impurity-modified eutectic as well as a random fractal (as a benchmark). Differences in the branching statistics point to different microstructural origins of the impurity- and quench-modified eutectic. Finally, our quantitative approach is not limited to the data presented here but can be used to extract abstract information from other volumetric datasets, without customization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 8. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 4. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 6. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 7. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 8. OUTCARL: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 10. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of hatch spacing and laser power on microstructure, texture, and thermomechanical properties of laser powder bed fusion (L-PBF) additively manufactured NiTi

This study systematically evaluates the effects of laser powder bed fusion additive manufacturing (L-PBF-AM) parameters (hatch spacing and laser power) on the thermomechanical behavior and microstructure of Ni 50.8 Ti 49.2 shape memory alloy. The samples were fabricated with hatch spacings from 40 to 240 µm and laser powers of 50 and 100 W at a constant scanning speed of 125 mm/s, resulting in parts with volumetric energy density levels from 55 to 666 J/mm 3 and two sets of linear energy densities of 0.4 and 0.8 J/mm. The results showed a reduced melt pool size and discontinuity of scan tracks with decreased laser power. Additionally, the porosity level was increased with larger hatch spacing and lower laser power. More notably, the transformation temperatures increased, and the critical stress, recoverable strain, and functional stability of samples improved with lower hatch spacing, where the recovery ratio of up to 90% was observed, regardless of the employed laser power. This study also discussed the relationship between the fabrication process and texture formation in the L-PBF-AM process. In conclusion, the advantage of L-PBF-AM was revealed in tailoring the microstructure from highly textured samples in [1 1 1] or [0 0 1] direction when hatch spacing lower than laser beam focused was employed, to the appearance of equiaxed solidification front with island grains and random orientations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Data-driven analysis and prediction of stable phases for high-entropy alloy design

High-entropy alloys (HEAs) represent a promising class of materials with exceptional structural and functional properties. However, their design and optimization pose challenges due to the large composition-phase space coupled with the complex and diverse nature of the phase formation dynamics. In this study, a data-driven approach that utilizes machine learning (ML) techniques to predict HEA phases and their composition-dependent phases is proposed. By employing a comprehensive dataset comprising 5692 experimental records encompassing 50 elements and 11 phase categories, we compare the performance of various ML models. Our analysis identifies the most influential features for accurate phase prediction. Furthermore, the class imbalance is addressed by employing data augmentation methods, raising the number of records to 1500 in each category, and ensuring a balanced representation of phase categories. The results show that XGBoost and Random Forest consistently outperform the other models, achieving 86% accuracy in predicting all phases. Additionally, this work provides an extensive analysis of HEA phase formers, showing the contributions of elements and features to the presence of specific phases. We also examine the impact of including different phases on ML model accuracy and feature significance. Notably, the findings underscore the need for ML model selection based on specific applications and desired predictions, as feature importance varies across models and phases. This study significantly advances the understanding of HEA phase formation, enabling targeted alloy design and fostering progress in the field of materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Identification of mechanisms driving heterogeneous void growth in ductile aluminum

Void growth plays a central role in ductile fracture, yet the specific mechanisms that control this remain obscure. Classical models, such as those proposed by Rice and Tracey in 1969, are able to capture average rates of void growth, but cannot capture the heterogeneity of individual void growth. Building on recent work, the present study employs laboratory-based diffraction contrast tomography and in-situ x-ray computed tomography to investigate the effect of grain structure and other microstructural factors on void growth in an Al-2219 alloy. Crystal plasticity finite element (CP-FE) modeling is used alongside experimental data to evaluate the contributions of local mechanical states, grain orientation, grain size, and neighboring microstructural features. No strong linear relationships are found with any of the considered descriptors and void growth rate. Potential complex nonlinear relationships are explored with the use of a random forest regression model, which identifies initial void volume, void aspect ratio, local normal stress state, local shear stress state, and local equivalent plastic strain (EQPS) as features that most improve void growth rate predictions. The combination of these analyses suggests that these features should be prioritized to improve models of void growth.

Diffraction contrast tomography (DCT)↗
DOE OSTI2022

Suppression of Superconductivity and Nematic Order in Fe 1–y Se 1–x S x (0 ≤ x ≤ 1; y ≤ 0.1) Crystals by Anion Height Disorder

Connections between crystal chemistry and critical temperature T c have been in the focus of superconductivity, one of the most widely studied phenomena in physics, chemistry and materials science alike. In most Fe-based superconductors materials chemistry and physics conspire so that T c correlates with the average anion height above the Fe plane, i.e. with the geometry of the FeAs 4 or FeCh 4 (Ch = Te, Se, S) tetrahedron. Here, by synthesizing Fe 1-y Se 1-x S x (0 ≤ x ≤ 1, y ≤ 0.1) we find that in alloyed crystals T c is not correlated with the anion height as most other Fe superconductors. Instead, changes in T c (x) and tetragonal-to-orthorombic (nematic) transition T s (x) on cooling are correlated with disorder in Fe vibrations in direction orthogonal to Fe planes, along the crystallographic c-axis. The disorder stems from the random nature of S substitution, causing deformed Fe(Se,S) 4 tetrahedra with different Fe-Se and Fe-S bond distances. Our results provide evidence of T c and T s suppression by disorder in anion height. The connection to local crystal chemistry may be exploited in computational prediction of new superconducting materials with FeSe/S building blocks.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Capturing the softening in T91 steel annealed with molten LBE through interfacial gradient plasticity

In order to study the effect of annealing with molten lead–bismuth eutectic alloy (mLBE) on the mechanical properties of T91 steel, micropillar compression tests were analyzed through the consideration of mechanical interface energy terms within gradient plasticity. The stress–strain curves of the micropillars showed significant stochastic effects, characterized by a different elastic and plastic behavior response for each pillar. Among the various pillars, some exhibited a similar trilinear behavior. Scanning electron microscopy (SEM) images of these specimens demonstrated localized slip deformation and slip planes after compression, while transmission electron microscopy (TEM) images showed the presence of a severe slip zone along the random grain boundaries (RGBs), indicating that the GBs played a dominant role in the deformation/slip. By employing interfacial gradient plasticity that can explicitly account for the presence of GBs, the trilinear response was captured, allowing for the determination of the mechanical interface parameter. As anticipated the pillars which underwent severe slip at the GBs had a lower value for the mechanical interface parameter. Finally, annealing with mLBE can, therefore, in some cases result in softening in the overall stress–strain, as it lowers the mechanical interface energy of GBs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A probabilistic creep model incorporating test condition, initial damage, and material property uncertainty

Uncertainty is prevalent in the creep resistance of alloys, where at elevated temperature and low pressure, rupture can range across logarithmic decades. In this study, a probabilistic continuum-damage-mechanics (CDM)-based model is derived to capture the uncertainty of creep resistance. To meet this objective, creep data for alloy 304 Stainless Steel is gathered. A constitutive model, “Sinh”, is calibrated deterministically to determine the statistical variability of the material properties. Three sources of uncertainty are injected into the model: test condition (stress and temperature), initial damage, and material properties. Probabilistic simulations are carried out by (a) calibrating probability distribution functions (pdfs) for each source of uncertainty (b) randomly sampling the pdfs using Monte Carlo methods and (c) executing simulations to replicate the uncertain creep behavior. A sensitivity analysis is performed to evaluate the relative effect of each source of uncertainty. In full probabilistic simulations, the cumulative uncertainty of creep behavior is evaluated. The probabilistic model accurately predicts the creep deformation and rupture of the available experiments. The probabilistic model is validated for interpolation but lacks extrapolation ability. Several future works are proposed to further improve the model.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A constitutive model incorporating grain refinement strengthening on metallic alloys

Surface nano-crystallization techniques have been recently developed as one of the most efficient ways to optimize materials’ structure, and therefore develop the local and global mechanical behavior as to increase strength without compromising ductility. Here, in this work, we present a constitutive model incorporating grain refinement hardening to simulate the nano-crystallization technique, specifically, surface mechanical attrition treatment. The computation is implemented using user-defined VUMAT subroutines. As an example of its application, a geometry model with full coverage of random impacts are employed. The results show that the model has rather precise predictability of grain size evolution during plastic deformation. The readily embedded with a computational code of material dynamics enables this novel model to be a promising tool to study the dynamic evolution of microstructures under plastic deformation.

36 MATERIALS SCIENCE↗