Search NASA⌕ Search

SEARCH · Search NASA

Results for “random microstructure”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Metallographic Characterization of Wrought Depleted Uranium

Metallographic characterization was performed on wrought depleted uranium (DU) samples taken from the longitudinal and transverse orientations from specific locations on two specimens. Characterization of the samples included general microstructure, inclusion analysis, grain size analysis, and microhardness testing. Comparisons of the characterization results were made to determine any differences based on specimen, sample orientation, or sample location. In addition, the characterization results for the wrought DU samples were also compared with data obtained from the metallographic characterization of cast DU samples previously characterized. No differences were observed in microstructure, inclusion size, morphology, and distribution, or grain size in regard to specimen, location, or orientation for the wrought depleted uranium samples. However, a small difference was observed in average hardness with regard to orientation at the same locations within the same specimen. The longitudinal samples were slightly harder than the transverse samples from the same location of the same specimen. This was true for both wrought DU specimens. Comparing the wrought DU sample data with the previously characterized cast DU sample data, distinct differences in microstructure, inclusion size, morphology and distribution, grain size, and microhardness were observed. As expected, the microstructure of the wrought DU samples consisted of small recrystallized grains which were uniform, randomly oriented, and equiaxed with minimal twinning observed in only a few grains. In contrast, the cast DU microstructure consisted of large irregularly shaped grains with extensive twinning observed in most grains. Inclusions in the wrought DU samples were elongated, broken and cracked and light and dark phases were observed in some inclusions. The mean inclusion area percentage for the wrought DU samples ranged from 0.08% to 0.34% and the average density from all wrought DU samples was 1.62E+04/cm 2 . Inclusions in the cast DU samples were equiaxed and intact with light and dark phases observed in some inclusions. The mean inclusion area percentage for the cast DU samples ranged from 0.93% to 1.00% and the average density from all wrought DU samples was 2.83E+04/cm 2 . The average mean grain area from all wrought DU samples was 141 μm 2 while the average mean grain area from all cast DU samples was 1.7 mm2. The average Knoop microhardness from all wrought DU samples was 215 HK and the average Knoop microhardness from all cast DU samples was 264 HK.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Microstructure evolution during low-strain thermo-mechanical processing and its repercussion on intergranular corrosion in alloy 600H

Highlights: • GBE microstructure in alloy 600H is realized through the sole activation of SIBM. • Higher TRD size and lower fractal dimension are characteristics of GBE microstructure. • Partially/fully recrystallized specimen shows higher DOS than as-received condition. • GBE specimen shows remarkable resistance against percolation and grain dropping. • Extent of percolation is governed by residual strain and random HAGBs connectivity. - Abstract: The current investigation deals with the evolution of microstructures as a consequence of varying processing parameters viz. extent of deformation, annealing temperature and time during grain boundary engineering (GBE) type thermo-mechanical processing in alloy 600H. Furthermore, the influence of the various product microstructures (viz. partially recrystallized, fully recrystallized and GBE) on intergranular corrosion is evaluated and compared with the as-received condition. Detailed characterization of the as-processed specimens have revealed that lower extent of deformation (7.5%) followed by a longer annealing treatment (i.e. 60 min) at a higher temperature (1373 K) could lead to a GBE microstructure with a higher proportion of Σ3{sup n} boundaries (~72%) due to the predominant activation of multiple twinning. In contrast, the specimens associated with moderately higher extent of deformation (10% and 15%) have revealed the occurrence of static recrystallization. The GBE specimen has demonstrated an excellent resistance against percolation essentially due to its discontinuous random high angle boundaries (HAGBs) network. In contrast, the nucleation of the finer grains in the recrystallized specimens caused the regeneration of random HAGBs network (as substantiated by the fractal analysis) that leads to higher percolation and severe grain dropping. In spite of the similar grain boundary character distribution, the severity of corrosion is relatively lesser in the as-received condition than the completely recrystallized specimen owing to its coarse-grained microstructure that commendably retarded the diffusion of chromium from the grain interior to the boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Text Mining for Process–Structure–Properties Relationships in Metals

With the advent of large language models (LLMs), the vast unstructured text within millions of academic papers is increasingly accessible for materials discovery—although significant challenges remain. While LLMs offer promising few- and zero-shot learning capabilities, particularly valuable in the materials domain where expert annotations are scarce, general-purpose LLMs often fail to address key materials-specific queries without further adaptation. To bridge this gap, fine-tuning LLMs on human-labeled data is essential for effective structured knowledge extraction (Liu in The Importance of Human-Labeled Data in the Era of LLMs, 2023). Here, in this study, we introduce a novel annotation schema designed to extract generic process–structure–properties relationships from scientific literature. We demonstrate the utility of this approach using a dataset of 128 abstracts, with annotations drawn from two distinct domains: high-temperature materials (Domain I) and uncertainty quantification in simulating materials microstructure (Domain II). Initially, we developed a conditional random field (CRF) model based on MatBERT—a domain-specific BERT variant—and evaluated its performance on Domain I. Subsequently, we compared this model with a fine-tuned LLM (GPT-4o from OpenAI) under identical conditions. Our results indicate that fine-tuning LLMs can significantly improve entity extraction performance over the BERT-CRF baseline on Domain I. However, when additional examples from Domain II were incorporated, the performance of the BERT-CRF model became comparable to that of the GPT-4o model. These findings underscore the potential of our schema for structured knowledge extraction and highlight the complementary strengths of both modeling approaches.

Materials science↗
DOE OSTI2025

COBRA:COMPUTED-TOMOGRAPHY BASED RANDOM-FIELD APPROXIMATION

SF-25-115 COBRA (COmputed-tomography Based Random-field Approximation) is a Python application for generating statistically equivalent random fields from CT-scan imagery. It leverages Karhunen–Loève expansions to model microstructural variability, enabling users to: Preprocess CT scans (filtering and Gaussian transformation); Fit covariance kernels fromempirical data; Solve eigenproblems to obtain KL modes; Sample random fields onsistent with fitted statistics; Postprocess samples back into the physical domain.

Hu, Tianchen↗
DOE OSTI2023

Regulating Gas Transport in Molecularly Engineered Polymer Membranes (Final Technical Report)

Energy-efficient separation processes are essential for a wide range of applications ranging from clean fuels (e.g., hydrogen purification) and petroleum refining (e.g., natural gas processing) to water purification and carbon capture. Membrane-mediated separations have shown tremendous promise in providing high productivity and high separation efficiency at significantly lower energy consumption, e.g., up to 90% less energy cost than traditional thermally driven processes such as distillation. Polymeric membranes–the dominant separation membrane materials–have yet to reach their full potential due to their limitations in long-term durability (e.g., productivity loss over the period of their lifetime due to physical aging) or insufficient stability under harsh conditions (e.g., high temperature, chemically complex feeds). This research seeks to establish a new paradigm in polymer membrane material design by harnessing crosslinked model networks with well-defined yet finely tailorable microstructure to facilitate fast and selective gas transport and simultaneously enhance membrane stability. Unlike traditional randomly crosslinked polymers, which suffer from structural inconsistencies and consequently suboptimal gas separation performance, crosslinked model network membranes prepared via a precisely controlled end-linking process enables the creation of previously unattainable microstructure tunability, which, in turn, results in versatile crosslinked membranes with high separation performance that meet the needs of various challenging gas separations. Using model network framework as a fundamental tool by applying this concept in diverse polymer categories, this work has led to the development of various innovative crosslinked membrane structures such as unimodal, bimodal and clustered model networks. These advanced crosslinked polymer membranes not only demonstrate exceptional gas separation performance that significantly outperform existing randomly crosslinked membranes, but also possess excellent long-term durability and robust stability under complex operating conditions. From a fundamental perspective, results from this research provide critical mechanistic insights into gas separation in crosslinked polymer membranes, addressing key knowledge gaps and opening new avenues for membrane design to meet various separation needs. The new membrane materials produced from this research enable the use of polymeric membranes for high temperature gas separations, offering substantial energy and cost savings by eliminating the need for repeated cooling-heating cycles in industrial processes.

02 PETROLEUM↗
DOE OSTI2024

A Finite Difference informed Random Walk solver for simulating radiation defect evolution in polycrystalline structures with strongly inhomogeneous diffusivity

Diffusivity of species and defects on grain boundaries is usually several orders of magnitude larger than that inside grains. Such strongly inhomogeneous diffusivity requires prohibitively high computational demands for modeling microstructural evolution. Here, this paper presents a highly-efficient numerical solver, combining the Finite Difference method and Random Walk model, designed for accurately modeling strongly inhomogeneous diffusion within polycrystalline structures. The proposed solver, termed Finite Difference informed Random Walk (FDiRW), integrates a customized Finite Difference (cFD) scheme tailored for fast diffusion along thin grain boundaries represented by a single-layer of nodes. Numerical experiments demonstrate that the FDiRW solver achieves an impressive efficiency gain of 1560x compared to traditional Finite Difference methods while maintaining accuracy, making it feasible for personal computer machines to handle diffusional systems with strongly inhomogeneous diffusivity across static polycrystalline microstructures. The model has been successfully applied to simulate radiation defect evolution, showcasing its scalability to engineering scales in both length and time dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design of novel refractory equiatomic multi-principal elemental alloys based on Mo-Nb-Ti system for Gen IV reactor applications

Excellent irradiation damage resistance demonstrated by multi-principal elemental alloys (MPEAs) has sparked significant interest among researchers, prompting exploration into their vast compositional space, to validate their suitability for nuclear applications. A combined approach of thermodynamic and empirical parameters calculations alongside CALPHAD (CALculation of PHAse Diagrams) for phase formation predictions enable high-throughput material selection for sophisticated applications like nuclear, overcoming laborious and time-consuming experiments. Key thermodynamic and empirical parameters for eight novel equiatomic MPEAs, based on seven low thermal neutron cross section refractory elements, for predicting phase formation were calculated, and equilibrium and non-equilibrium simulations in CALPHAD were employed to comprehensively model the systems. Pseudo binary phase diagram simulations showed that Zr, V or equiatomic CrV additions to the base MoNbTi alloy (MoNbTi-Zr, MoNbTi-V and MoNbTi-CrV alloys) favor the formation of isomorphous body-centered cubic (BCC) phase at high temperatures, while Cr, Al, equiatomic ZrV, or equiatomic CrAl additions (MoNbTi-Cr, MoNbTi-Al, MoNbTi-ZrV or MoNbTi-CrAl alloys) limit the solubility of them. Equilibrium CALPHAD simulations at 750 oC were consistent with XRD results on MoNbTi, MoNbTiZr and MoNbTiCr alloys, and partially for others. Notably, elemental segregation observed in the backscattered electron (BSE) scanning electron microscopy (SEM) images of the alloys was accurately simulated through non-equilibrium Scheil solidification calculations in CALPHAD, further verified by experiments. The precipitation of TiCr2 Laves phase in Cr containing MoNbTiCr and MoNbTiCrAl was accurately predicted while discrepancies were noted in MoNbTiCrV. The equilibrium simulations also provided insights into phase compositions at specific temperatures offering a pathway for tailoring the desired microstructure and properties of these systems. Empirical parameters calculations successfully predicted random solid solution in the base MoNbTi alloy, and with an exception in MoNbTiV and MoNbTiAl, predicted intermetallic precipitation in the rest, especially, Laves phase precipitation in Cr containing alloys.

36 - MATERIALS SCIENCE↗
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 OSTI2025

Blocky Selective Postpolymerization C–H Functionalization of Polyolefins

C–H functionalization of commodity polyolefins affords functional materials derived from a high-volume, low-cost resource. However, current postpolymerization modification strategies result in randomly distributed functionalization along the length of the polymer backbone, which has a negative impact on the crystallinity of the resultant polymers, and thus the thermomechanical properties. Here, we demonstrate an amidyl radical mediated C–H functionalization of polyolefins to access blocky microstructures, which exhibit a higher crystalline fraction, larger crystallite size, and improved mechanical properties compared to their randomly functionalized analogues. Taking inspiration from the site-selective C–H functionalization of small molecules, we leverage the steric protection provided by crystallites and target polymer functionalization to amorphous domains in a semicrystalline polyolefin gel. The beneficial outcomes of blocky functionalization are independent of the identity of the pendant functional group that is installed through functionalization. Here, the decoupling of functional group incorporation and crystallinity highlights the promise in accessing nonrandom microstructures through selective functionalization to circumvent traditional tradeoffs in postpolymerization modification, with potential impact in advanced materials and upcycling plastic waste.

Neidhart, Eliza K. [The University of North Caroli↗
DOE OSTI2023

Predictive machine learning approaches for the microstructural behavior of multiphase zirconium alloys

Abstract Zirconium alloys are widely used in harsh environments characterized by high temperatures, corrosivity, and radiation exposure. These alloys, which have a hexagonal closed packed (h.c.p.) structure thermo-mechanically degrade, when exposed to severe operating environments due to hydride formation. These hydrides have a different crystalline structure, than the matrix, which results in a multiphase alloy. To accurately model these materials at the relevant physical scale, it is necessary to fully characterize them based on a microstructural fingerprint, which is defined here as a combination of features that include hydride geometry, parent and hydride texture and crystalline structure of these multiphase alloys. Hence, this investigation will develop a reduced order modeling approach, where this microstructural fingerprint is used to predict critical fracture stress levels that are physically consistent with microstructural deformation and fracture modes. Machine Learning (ML) methodologies based on Gaussian Process Regression, random forests, and multilayer perceptrons (MLP) were used to predict material fracture critical stress states. MLPs, or neural networks, had the highest accuracy on held-out test sets across three predetermined strain levels of interest. Hydride orientation, grain orientation or texture, and hydride volume fraction had the greatest effect on critical fracture stress levels and had partial dependencies that were highly significant, and in comparison hydride length and hydride spacing have less effects on fracture stresses. Furthermore, these models were also used accurately predicted material response to nominal applied strains as a function of the microstructural fingerprint.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Machine Learning Approach to Quantitative Analysis of Enamel Microstructure from Scanning Electron Microscopy Images

Dental enamel, the outermost tissue of mammalian teeth, must withstand a lifetime of wear and cyclic contact. To meet this demand, enamel possesses a combination of high hardness and resistance to fracture, properties that are typically mutually exclusive. The impressive damage tolerance has been attributed largely to decussation of the enamel rods, the principal unit of its microstructure. As such, enamel is inspiring the design of next‐generation structural materials. However, quantitative descriptions of the decussated enamel rod microstructure remain limited due to challenges encountered in applying computed tomography and in acquiring quality images appropriate for traditional digital processing methods. Here, a machine learning segmentation method is applied to images of the enamel obtained using scanning electron microscopy to support quantitative analysis of the microstructure. A pretrained convolutional neural network is used to expand the input training image dataset to allow the training of a random forest classifier, which ultimately segments the image with a very small training set ( n = 3 images). A validation of this segmentation method is presented, in addition to its application to calculate relevant microstructural parameters for images of tooth enamel from selected mammalian species. The methodology applied here is equally applicable to other hard tissues.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Effects of TiB{sub 2} Particle and Short Fiber Sizes on the Microstructure and Properties of TiB{sub 2}-Reinforced Composite Coatings

In this study, particle and short fiber-reinforced titanium matrix composite coatings are prepared via laser in situ technique using (0.5 and 50 μm) TiB{sub 2} and Ti powder as cladding materials. The microstructure and properties of the composite coatings are studied, and the changing mechanism of the microstructure is discussed. The results reveal that particle agglomeration is prone to appear with using fine TiB{sub 2} particles. Decomposition of the particles preferentially occurs with using coarse TiB{sub 2} particles. The cracks and pores on the surface of the coating are formed at a lower laser energy density. With the increase in the laser energy density, cracking on the surface of the coating diminishes, but the coating exhibits depression behavior. The depression extent of the coating using fine TiB{sub 2} particle as the reinforcement is much less than that of the coating using coarse TiB{sub 2} particle. Moreover, the size of the aggregate and the tendency of cracking can be reduced with the increase in Ti addition. Meanwhile, short TiB fiber bundles are formed by the diffusion mechanism of rod aggregate, and randomly oriented TiB short fibers are formed mainly by the dissolution–precipitation mechanism of fine TiB{sub 2} particles. Moreover, the growth of short TiB fibers can be in an alternating manner between B27 and Bf structures. The micro-hardness and wear resistance of the coatings are evidently higher than that of the titanium alloy substrate. The wear resistance of the large size TiB{sub 2} coating is higher than that of the small size TiB{sub 2} coating under the condition of low load.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hierarchical reconstruction of 3D well-connected porous media from 2D exemplars using statistics-informed neural network

The relationships between porous microstructures and transport properties are of fundamental importance in various scientific and engineering applications. Due to the intricacy, stochasticity and heterogeneity of porous media, reliable characterization and modeling of transport properties often require a complete dataset of internal microstructure samples. However, it is often an unbearable cost to acquire sufficient 3D digital microstructures by purely using microscopic imaging systems. Herein this paper presents a machine learning-based technique to hierarchically reconstruct 3D well-connected porous microstructures from one isotropic or several anisotropic low-cost 2D exemplar(s). To compactly characterize the large-scale microstructural features, a Gaussian image pyramid is built for each 2D exemplar. Local morphology patterns are collected from the Gaussian image pyramids, and then they serve as the training data to embed the 2D morphological statistics into feed-forward neural networks at multiple length levels. By using a specially-developed morphology integration scheme, the 3D morphological statistics at different levels can be inferred from the statistics-informed neural networks. Gibbs sampling is adopted to hierarchically reconstruct 3D microstructures by using multi-level 3D morphological statistics, where the large-scale, regional and local morphological patterns are statistically generated and successively added to the same 3D random field. The proposed method is tested on a series of porous media with distinct morphologies, and the statistical equivalence between the reconstructed and the real microstructures is systematically evaluated by comparing morphological descriptors and transport properties. The results demonstrate that the proposed 2D-to-3D microstructure reconstruction method is a universal and efficient approach to generating morphologically and physically realistic samples of porous media.

42 ENGINEERING↗
DOE OSTI2018

Annealing-induced microstructural evolution and mechanical anisotropy improvement of the Mg-Gd-Y-Zr alloy processed by hot ring rolling

Highlights: • Growth of SRX grains in GBR is random, while that of TBR is anisotropic. • The new texture component in annealed MR weakens the anisotropy obviously. • T6 treatment leads to an excellent comprehensive mechanical property. - Abstract: Hot ring rolling (HRR) is widely used to fabricate large-size seamless ring-shape components. However, this process produces a severe mechanical anisotropy, which is not acceptable in practical industrial application and should be avoided. This paper reports the annealing-induced microstructural evolution and mechanical anisotropy improvement of HRRed Mg-8.12Gd-1.82Y-0.35Zr alloy ring. Results show that annealing treatment could improve the microstructural homogeneity of the ring owing to the occurrence of static recrystallization (SRX). SRX grains are preferentially formed at twin boundaries and initial grain boundaries due to high dislocation density, high level of stored energy, and stress concentration. Twin boundaries, to some extent, can suppress the migration of SRX grain boundaries. In addition, occurrence of SRX weakens the texture intensity and basal pole inclination, leading to a decrease in the mechanical anisotropy of the ring. Considering the weakened anisotropy, comparable strength, and improved ductility, annealing + aging (T6), compared with direct aging (T5), is a more suitable heat treatment schedule for the HRRed Mg-8.12Gd-1.82Y-0.35Zr alloy ring.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Installing Controlled Stereo-Defects Yields Semicrystalline and Biodegradable Poly(3-Hydroxybutyrate) with High Toughness and Optical Clarity

Stereo-defects present in stereo-regular polymers often diminish thermal and mechanical properties, and hence suppressing or eliminating them is a major aspirational goal for achieving polymers with optimal or enhanced properties. Here, we accomplish the opposite by introducing controlled stereo-defects to semicrystalline biodegradable poly(3-hydroxybutyrate) (P3HB), which offers an attractive biodegradable alternative to semicrystalline isotactic polypropylene but is brittle and opaque. We enhance the specific properties and mechanical performance of P3HB by drastically toughening it and also rendering it with the desired optical clarity while maintaining its biodegradability and crystallinity. This toughening strategy of stereo-microstructural engineering without changing the chemical compositions also departs from the conventional approach of toughening P3HB through copolymerization that increases chemical complexity, suppresses crystal-lization in the resulting copolymers, and is thus undesirable in the context of polymer recycling and performance. More specifically, syndio-rich P3HB (sr-P3HB), readily synthesized from the eight-membered meso-dimethyl diolide, has a unique set of stereo-microstructures comprising enriched syndiotactic [rr] and no isotactic [mm] triads but abundant stereo-defects randomly distributed along the chain. This sr-P3HB material is characterized by high toughness (UT = 96 MJ/m 3 ) as a result of its high elongation at break (>400%) and tensile strength (34 MPa), crystallinity (T m = 114 degrees C), optical clarity (due to its submicron spherulites), and good barrier properties, while it still biodegrades in freshwater and soil.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗