Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grain”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

PRESOLAR GRAIN ABUNDANCE VARIATION IN THE MILLER RANGE 090019 CO3.1 CHONDRITE

Presolar grains condensed in the outflows of evolved red giant stars and the ejecta of supernovae (SNe) and novae. These grains have greatly anomalous isotopic compositions compared to solar system material, reflecting their stellar origins [1]. They have been identified in primitive meteorites, interplanetary dust particles (IDPs), Antarctic micrometeorites, and comet Wild 2 samples returned by NASA’s Stardust mission. Presolar silicates are one of the most abundant presolar phases and their concentrations extend up to 1.5% in primitive IDPs believed to derive from comets [2]. These grains are highly susceptible to alteration and destruction by secondary processing in the interstellar medium, nebula, and asteroid or comet parent body. Presolarsilicate abundance variations between primitive meteorites and chemical and mineralogical studies provide indications of the extent of secondary hydrothermal alteration [3, 4]. The abundance of presolar SiC grains is generally consistent among chondrites, but lower abundances in some meteorites have been attributed to thermal alteration [5, 6]. Presolar grain abundance variations attributed to localized alteration have also been reported within different regions of a chondrite [e.g., 7].Carbonaceous chondrites from the CO and CR groups have the highest abundances of presolar silicates among meteorites, attesting to their primitive nature. The CO3 chondrite Dominion Range (DOM) 08006 has the highest presolar O-rich grain abundance of ~260 ppm [8, 9]. MillerRange (MIL) 090019 is classified as a CO3.1 chondrite and has affinities to Acfer 094, DOM 08004/6 and Allan Hills (ALH) 77307. These chondrites have high presolar silicate abundances and contain high abundances of various types of refractory inclusions. We previously conducted detailed studies of CAIs in MIL 090019 [10, 11]. Here we evaluate its presolar grain inventory to assess the degree of parent body alteration and compare to other chondrites

A. N. Nguyen↗
NASA NTRS2020

Grain Formation around the AGB Star L2 Puppis Based On ALMA Observations

While models of grain formation in the outflows of carbon-rich stars have been relatively successful, models of outflows from oxygen-rich, asymptotic giant branch stars have been less fortunate. Under current modeling, it is difficult to produce sufficient amounts of silicate grains with high enough opacity to form a dust-driven wind from these stars. To investigate the cause of this difference, this work is a comparison between typical outflow model results and a model using input from Atacama Large Millimeter/submillimeter Array observations of L2Puppis. The temperatures from these observations are much lower than would typically be used in the standard outflow model. In addition, the observed gas densities are much higher than predicted from typical outflow models. Both of these differences make the formation of silicate grains much more favorable than predicted in current outflow models. We then explore the effects of other possible nonideal conditions including the efficiency of cluster growth prior to nucleation, the efficiency of grain growth following nucleation and the variation of grain coupling to stellar radiation during grain growth. Finally, we calculate the potential enhancement in grain production based on possible increased refractory abundances resulting from the vaporization of millimeter-scale and larger particles left over from the presence of a former planetary system.

Joseph A Nuth↗
DOE OSTI2026

High temperature stabilization of ultrafine grain tungsten alloys through synergistic compositional complexities

Thermally-stabilized fine-grained microstructures in tungsten provide a pathway to harnessing enhanced properties such as a reduced ductile-to-brittle transition temperature and improved strength while mitigating the adverse effects of grain growth and recrystallization. Here, in this study, we employ a material design strategy that relies on grain boundary segregation in the nanocrystalline state driven by alloy thermodynamics balanced with impurity scavenging through in situ formation of kinetically-stabilizing metal carbides. A W-Ti-Cr alloy is designed through lattice Monte Carlo methods and subsequently synthesized in a single-phase nanocrystalline state, which upon annealing, evolves into an ultrafine grained microstructure containing chromium grain boundary segregation collectively with a dispersed TiX (X=C,O) phase through the reaction of titanium with carbon and oxygen impurities. In situ synchrotron X-ray diffraction experiments demonstrate that increased Ti promotes stabilization across a larger temperature range but with diminishing returns above 10 at.% Ti. Long-term stability was confirmed through the retention of the ultrafine grained microstructure for a total of 8 days (192 h) at 1300 °C without grain/carbide growth and/or recrystallization. Our results demonstrate that, through strategic tailoring of composition and microstructure, one can harness the benefits of both thermodynamic and kinetic stabilization mechanisms, opening pathways for future alloy formulations that expand the window of stability.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore

Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.

08 HYDROGEN↗
DOE OSTI2025

Experimental characterization and atomistic simulation of grain boundary segregation in Mg-Y alloys

As a rare earth solute element in Mg alloys, Y has the beneficial effects of increasing both the strength and the ductility as well as weakening the crystallographic texture. To achieve a more fundamental understanding on how Y addition affects the microstructural evolution and mechanical properties, the Y segregation behavior at grain boundaries was investigated in Mg-1wt.%Y and Mg-7wt.%Y alloys at different conditions. The segregation intensity and its dependence on the grain boundary misorientation angle were experimentally characterized and computationally predicted. Strong segregation at grain boundaries was observed in both low and high Y-containing alloys. Y segregation was found to remain in alloy Mg-7Y after high-temperature annealing heat treatment at 540 °C. No direct correlation between the Y segregation intensity and the grain boundary misorientation angle could be established based on either the experimental characterization or the atomistic simulation with a spectral model. We thus conclude that grain boundary segregation of Y is independent of grain boundary misorientation angle.

Grain boundary↗
DOE OSTI2025

Meshfree simulation and prediction of recrystallized grain size in friction stir processed 316L stainless steel

Friction stir processing (FSP) is a promising solid-phase microstructural modification technique that can repair and enhance damaged stainless steel surfaces exposed to harsh environments. The quality of the repaired material is closely correlated to the recrystallized grain size in the stir zone (SZ), which is influenced by the thermomechanical conditions dictated by FSP process parameters. Thus, establishing a reliable relationship between these parameters and recrystallized grain size in the SZ is crucial for optimizing repair quality. However, existing experimental approaches often rely on indirect temperatures measured far from the SZ, along with rough strain rate estimations, which are imprecise and time-consuming. Meanwhile, existing mesh-based modeling methods usually face numerical challenges when dealing with the large material deformations inherent in FSP. Here, to address these issues, this study introduces a meshfree process model for FSP based on the smoothed particle hydrodynamics (SPH) method, aimed at predicting process conditions under different parameters. The model is validated using experimental data from 11 combinations of tool traverse and rotation speeds on 316 L stainless steel. Correlations between process parameters, material flow, temperature, strain, strain rate, and recrystallized grain size are revealed through SPH simulations and electron backscatter diffraction (EBSD) imaging. The results show that in situ SZ temperatures range from 1071 to 1322°C, which exceed the tool temperature by over 300°C. Furthermore, SZ temperature, strain rate, and grain size increase monotonically with higher tool temperature and faster traverse speed. A relationship is then established between the model-predicted Zener-Hollomon parameter and the recrystallized grain size based on EBSD data, expressed as ln(d) = -0.364 ln(Z) + 14.673. Finally, this relationship exhibits satisfactory accuracy with errors of less than 26.9% in predicting grain sizes at various SZ locations, which offers valuable insights for optimizing FSP repair processes for 316 L stainless steel.

316L stainless steel↗
DOE OSTI2026

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗
DOE OSTI2024

Nanoscale Imaging and Measurements of Grain Boundary Thermal Resistance in Ceramics with Scanning Thermal Wave Microscopy

Material thermal conductivity is a key factor in various applications, from thermal management to energy harvesting. With microstructure engineering being a widely used method for customizing material properties, including thermal properties, understanding and controlling the role of extended phonon-scattering defects, like grain boundaries, is crucial for efficient material design. However, systematic studies are still lacking primarily due to limited tools. In this study, we demonstrate an approach for measuring grain boundary thermal resistance by probing the propagation of thermal waves across grain boundaries with a temperature-sensitive scanning probe. The method, implemented with a spatial resolution of about 100 nm on finely grained Nb-substituted SrTiO 3 ceramics, achieves a detectability of about 2 × 10 –8 K m 2 W –1 , suitable for chalcogenide-based thermoelectrics. The measurements indicated that the thermal resistance of the majority of grain boundaries in the STiO 3 ceramics is below this value. While there are challenges in improving sensitivity, considering spatial resolution and the amount of material involved in the detection, the sensitivity of the scanning probe method is comparable to that of optical thermoreflectance techniques, and the method opens up an avenue to characterize thermal resistance at the level of single grain boundaries and domain walls in a spectrum of microstructured materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Physics-informed machine learning analysis for nanoscale grain mapping by synchrotron Laue microdiffraction

Understanding the grain morphology, orientation distribution and crystal structure of nanocrystals is essential for optimizing the mechanical and physical properties of functional materials. Synchrotron X-ray Laue microdiffraction is a powerful technique for characterizing crystal structures and orientation mapping using focused X-rays. However, when the grain sizes are smaller than the beam size, mixed peaks in the Laue pattern from neighboring grains limit the resolution of grain morphology mapping. We propose a physics-informed machine learning (PIML) approach that combines a convolutional neural network feature extractor with a physics-informed filtering algorithm to overcome the spatial resolution limits of X-rays, achieving nanoscale resolution for grain mapping. Our PIML method successfully resolves the grain size, orientation distribution and morphology of Au nanocrystals through synchrotron microdiffraction scans, showing good agreement with electron backscatter diffraction results. This PIML-assisted synchrotron microdiffraction analysis can be generalized to other diffraction-based probes, enabling the characterization of nanosized structures with micrometre-sized probes.

X-ray crystallography↗
NASA NTRS2021

An Investigation of the Grain Size and Adhesion Strength of High-Speed Impact Ice

Ice accumulation on aircraft from supercooled water droplets (impact icing) can increase drag and decrease lift, leading to a decrease in fuel efficiency and a potential for unsafe situations. Anti-icing materials and coatings are of high interest to prevent aircraft icing. Because no anti-icing material suitable for in-flight situations has been realized, icing models aim to predict the accumulation and shape of the ice to understand its effects. In both cases, understanding the physical mechanism behind ice adhesion and shedding is necessary. A recent literature review has shown that the apparent ice adhesion strength to any specific material varies over several orders of magnitude across studies and test methods. One of the reasons for this high variance is likely to be the varying microstructure of the ice used for different studies. Although the microstructure of impact ice has been studied in the past, those conditions were not relevant to that of in-flight aircraft icing. The microstructure of the ice formed will determine both its bulk and surface properties, yet it has been overlooked in the ice adhesion community. Here, we use optical microscopy to study the non-uniform ice crystal grain size in impact ice accreted in NASA Glenn Research Center’s Icing Research Tunnel at various icing conditions. Similar to other environmental conditions studied in the literature, we show an increase in grain size at larger distances from the growth interface; however, the data shows substantially larger grains than what was expected from the prior literature. Ice accrued at lower temperatures has a smaller grain size at all distances from the growth interface. No change in grain size was noticed when varying only the wind speed. A correlational analysis shows that smaller grain sizes at the growth interface corresponded to an increase in the apparent ice adhesion strength when the adhesive failure was interfacial, but not when the failure was cohesive. These results may explain the overall trend in the literature of increasing apparent adhesion strength with decreasing temperature. Incorporating microstructure trends for various environmental conditions in icing models may increase the certainty in predicting ice adhesion strength and shedding, and will provide a better basis for comparison between adhesive tests on ice produced in different facilities and at different storage times.

• Ice Fabric↗
DOE OSTI2026

Elucidating the impact of stress states on grain boundary passivation in Ni-5Cr using the Rhines pack method

The oxidation behavior of a Ni–5Cr (at.%) alloy was evaluated at 420 ° C using the Rhines pack method and simultaneous tensile and compressive stress states via a miniature four-point bending fixture. At this moderate temperature, grain boundaries dominate mass transport and the resulting oxidation response. Oxidation produced approximately 1μ⁢m-wide protective Cr 2 O 3 films capping some grain boundaries and penetrative, intergranular Cr-rich oxides at other grain boundaries. Externally applied tensile and compressive stress during oxidation increased the prevalence of Cr 2 O 3 cap formation compared to no applied stress, with tensile stress resulting in more Cr 2 O 3 caps than compressive stress. The observed oxide cap morphology was similar across all test conditions. Regions under compressive stress showed an order of magnitude greater Cr depletion depth along the grain boundary. Chromium nitrides (CrN), likely from N contamination of the Rhines pack cell, were observed both at the oxide–metal interface and intergranularly within the alloy. Collectively these results demonstrate that applied stress promotes localized protective oxide cap formation over grain boundaries, with compressive stress additionally promoting deeper Cr depletion. Furthermore the experimental approach helped separate out the effects of stress on local oxide formation and grain boundary passivation.

4 point bend↗
DOE OSTI2026

Dopant adsorption as a function of bulk concentration at a near 42º (100) twist grain boundary in SrTiO3

The enrichment of grain boundaries with dopant atoms is of critical importance for the macroscopic physical properties of materials. In thermodynamic equilibrium the Gibbs adsorption isotherm relates grain boundary excess of dopant atoms, their chemical potential in the adjacent bulk, and the respective interface energy. This study has used bicrystals with a near 42º (100) twist grain boundary in SrTiO3 to demonstrate that different kinetic pathways of Fe dopant atom additions converge towards comparable grain boundary configurations. Despite differences in bulk chemical potentials remarkably similar grain boundary excess quantities were observed. The experimental results indicate the general experimental feasibility to establish quantitative relationships between variations of grain boundary energy, interfacial excess, and overall dopant concentration. Improved experimental counting statistics are needed to distinguish solute interface excess as a function of bulk chemical potential.

Hahn, William [University of California, Davis]↗
DOE OSTI2026

Revealing the complex chemistry of grain boundaries in K-doped BaFe 2 As 2 with atom probe tomography

Iron-based superconductors have attractive properties for high-field applications, but there is a lack of understanding of the effect of grain boundary chemistry on the in-field performance. The near atomic-scale resolution, ppm sensitivity and 3D analysis offered by atom probe tomography make it a powerful tool to investigate the nanoscale structure and chemistry of these defects in fine-grained K-doped BaFe 2 As 2 samples. A computational method to systematically extract and compare the Gibbsian interfacial excess of chemical species across grain boundaries has been explored in this work. The robustness of the method has been tested by evaluating the effects of selected variables on simulated APT datasets. The accuracy and precision of the calculated Gibbsian interfacial excess were found to be stable over a range of analysis conditions: varying grain boundary widths and detection efficiencies, spatial precisions below 1.5 nm, and bin widths between 1.2 and 1.6 nm. For the K-doped BaFe 2 As 2 samples studied, segregation of As, Ba, K and impurities of O, Na, and Sb were found at grain boundaries. The Gibbsian excess values were found to vary widely between different boundaries, showing the complexity of the grain boundary chemistry in this material. Possible links between the observed critical current density (Jc) of these samples and their nano- and micro-structure have also been investigated and discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Oxygen Distribution and Segregation at Grain Boundaries in Nb and Ta-Encapsulated Nb Thin Films for Superconducting Qubits

We report on atomic-scale analyses of oxygen distribution and segregation at grain boundaries (GBs) of Nb and Ta-encapsulated Nb (Ta/Nb) thin films for superconducting qubits using atom probe tomography (APT) and transmission electron microscopy (TEM). We observe oxygen segregation at grain boundaries (GBs) relative to the oxygen concentration within the grains for both Nb and Ta-capped Nb thin films and find that a higher oxygen concentration in the interior of Nb grains leads to greater oxygen segregation levels at GBs. This finding reveals that the formation of a local equilibrium of oxygen concentration between GBs and grain interiors of Nb is the primary driving force of the oxygen segregation behaviors in Nb and Ta-capped Nb. The enrichment factors (CGB/Cgrain) for oxygen segregation at GBs in Nb and Ta-capped Nb range from 2.4 ± 0.3 to 2.7 ± 0.4. The current results also highlight that controlling oxygen impurities in Nb during film deposition and fabrication processing is important to concomitantly reducing the level of oxygen segregation at GBs in Nb. Finally, we find that increases in the oxygen concentration in both Nb grains and GBs correlate with a suppression in the critical temperature for superconductivity (Tc). Together, our comparative chemical and charge transport property analyses provide atomic-scale insights into a potential mechanism, contributing to the decoherence in superconducting qubits.

Lee, Jaeyel [Fermilab] (ORCID:0000000185429612)↗
DOE OSTI2024

Phase-field modeling of stored-energy-driven grain growth with intra-granular variation in dislocation density

Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.

Materials Science↗
DOE OSTI2024

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗
DOE OSTI2025

Grain boundary energy control in zinc aluminate nanoceramics

This study investigates the grain boundary energy dependence on segregated dopants in nanocrystalline zinc aluminate ceramics. Atomistic simulations of Σ3 and Σ9 grain boundaries showed that trivalent ions of varying ionic radii [Sc 3+ (74.5 pm), In 3+ (80.0 pm), Y 3+ (90.0 pm), and Nd 3+ (98.3 pm)] have a tendency to segregate to both interfaces, with Y 3+ presenting the highest segregation potentials. The connection between segregation and the reduction of interfacial energies was explored by measuring the grain boundary energy on nanoceramics fabricated via high‐pressure spark plasma sintering (HP‐SPS) using differential scanning calorimetry (DSC). The results revealed that Y 3+ doping at 0.5 mol% reduces the grain boundary energy in zinc aluminate nanoceramics from 1.1–1.3 J/m 2 to 0.6–0.8 J/m 2 ; the range correlates with the observed size dependence of the excess energy, with higher values observed for the smaller grain sizes (∼17 nm). The noted decrease in interfacial energies for doped samples suggests it is indeed possible to alter the stability of zinc aluminate grain boundaries via dopant segregation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling the Impact of Grain Size on Corrosion Behavior of Ni-Based Alloys in Molten Chloride Salt via Cellular Automata

Molten chloride salts hold significant promise as both thermal transfer and storage media for next-generation concentrated solar power (CSP) systems. However, molten chlorides pose a considerable corrosion risk to structural materials, particularly Ni-based alloys. One approach to enhancing corrosion resistance is through the optimization of grain structure; however, it remains uncertain whether increasing or decreasing grain size enhances corrosion resistance. A cellular automata (CA) program was developed to evaluate the interplay between grain size and corrosion in Ni-based alloy. Our CA program tracks alloy composition, surface roughness, and thickness loss via a graphical user interface, displaying corrosion and diffusion status, and multiple user input cards for tuning the simulation. CA simulations of Inconel 625 indicate enhanced corrosion resistance with increased grain size, with passivating oxides offering limited protection. Additionally, the temporal evolution of alloy surface roughness demonstrates notable fluctuations, with abrupt increases attributed to corrosion along vertical grain boundaries and sudden decreases to grain detachment from the protective film.

Materials Science↗