Search NASA⌕ Search

SEARCH · Search NASA

Results for “microstructure scale modeling”

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 73 records · Page 4

Multiscale and Machine Learning Modeling for Process-informed Microstructure Prediction in Additively Manufactured Materials using MALAMUTE

The Advanced Materials and Manufacturing Technologies (AMMT) program under the Department of Energy Office of Nuclear Energy aims to develop and qualify additively manufactured materials for nuclear applications. One key challenge to this is the microstructural variability observed in the additively manufactured products and their impact on the properties and performance of the material in extreme environments. AMMT is using a combination of high-throughput experimental and modeling techniques to accelerate qualification. Conventionally, in-situ and ex-situ characterizations and testing are performed to correlate different aspects of the additive manufacturing process to the final product and its performance. However, adopting a trial-and-error approach to experimentally evaluate the vast range of process parameters required to capture microstructural variability is cost-prohibitive. Modeling and simulation provide a comparatively inexpensive way to understand and correlate the microstructural evolution to the processing conditions. The modeling and simulation work-packages within the AMMT program aims to use physics-based and machine learning models to develop a digital twin for additive manufacturing that can correlate the process conditions to the final product and establish a process-structure-property-performance (PSPP) correlation. The melting and subsequent solidification that occurs during the additive process is a complex phenomenon that requires multiscale multiphysics analysis. This work package focuses on understanding the role of process variabilities on the unique microstructural characteristics of additively manufactured materials. Microstructural features at the subgrain level, such as compositional micro-heterogeneity and dislocation cells, are of particular interest here since they can influence the creep properties and radiation performance. Idaho National Laboratory’s Multiphysics Object-Oriented Simulation Environment (MOOSE), specifically the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) software, provides an ideal platform for developing the multiphysics multiscale model to explore the intricacies of the microstructural evolution during the AM processes within a single framework. Furthermore, given that such full-fidelity simulations can be computationally intensive, reduced order models are necessary to explore the PSPP space for additively manufactured materials in an efficient, reliable, and cost-effective way. This work focuses on capturing the microstructural variabilities at the subgrain level that are often missing in the part-scale models. In fiscal year 2025, we significantly advanced upon our work in the last fiscal year, in terms of the predictive capabilities of the physics-based and ML models, by adding the capabilities to capture subgrain-level micro-segregation during solidification using phase-field model and to predict the time-dependent dynamics of the AM process through the MOGPAR model. The alloy solidification model in MOOSE incorporates the thermodynamic properties and free energy relevant to 316 stainless steel. The model demonstrates the Cr and Ni segregation that occurs during solidification, including that the rate of solidification. The microstructural evolution model is connected to the process conditions via the surrogate model developed in this work. This enables predictions of the final microstructure in conjunctions with the manufacturing process. This work supports AMMT's rapid qualification goals by laying the foundation for an efficient and cost-effective model establishing the PSPP correlation for AM. The generated microstructures and predicted micro-segregation can be used by other work packages under AMMT to evaluate the properties and environmental response of the material at the mesoscale. Thus, this work helps to identify the key microstructural features at the subgrain level that are significant in property and performance predictions of additively manufactured components. This work will also provide inputs to the large-scale process variability models to reevaluate and validate assumptions and simplifications made in the part-scale models. Furthermore, through active learning this work can help identify the data need from both modeling and experimental sides for development of a robust digital twin for additive manufacturing and accelerate the AMMT's qualification efforts.

36 - MATERIALS SCIENCE↗

SEM Milestone Report for MARMOT Validation Project

MARMOT is a meso-scale fuel code under development by NEAMS designed to predict the microstructure evolution of fuel materials to determine updated thermal and mechanical property models. These updated models can be used by engineering scale codes such as BISON. To support its development and validation the MARMOT Validation Project (MVP) was initiated and consists of extensive microstructure and thermal characterization of U-10Zr fuel material before and at discreet levels of irradiation. This report summarizes the SEM and EDS characterization of the lowest burnup samples (1.5% FIMA) and how they compare to pre-irradiation SEM and EDS.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗

Calibrating uncertain parameters in melt pool simulations of additive manufacturing

Melt pool scale numerical modeling of additive manufacturing (AM) processes can provide predictive capabilities and theoretical insight into the process-property-structure-performance relationships for AM parts. Despite capabilities of numerical models to solve complex multi-physics problems, it is often important to consider a tradeoff between detailed physics and computational cost. Therefore, sources of uncertainty in both experimental conditions and the parameters needed for modeling require models to be validated against empirical evidence. Here, a method is proposed to calibrate uncertain parameters used in continuum-scale melt pool models for powder bed fusion (PBF) AM. Both a simplified heat transfer model and a heat transfer and fluid flow model were investigated. A surrogate model and Markov chain-based optimization algorithm calibrated melt pool geometry for models within experimental variation of the target melt pool width and depth from the NIST AM-Bench 2018-02 dataset. The melt pool temperature distributions, solidification parameters, and simulated multi-layer solidification microstructures were compared between the two models. Similar results from both models indicate that calibrated, lower fidelity numerical models may be used in place of higher fidelity models to generate melt pool solidification data. Finally, these calibrated models therefore enable lower computational cost melt pool simulations without a noticeable decrease in simulation accuracy for grain-scale microstructure simulations.

36 MATERIALS SCIENCE↗

Towards inverse microstructure-centered materials design using generative phase-field modeling and deep variational autoencoders

The field of Integrated Computational Materials Engineering (ICME) combines a broad range of methods to study materials’ responses over a spectrum of length scales. A relatively unexplored aspect of microstructure-sensitive materials design is uncertainty propagation and quantification (UP/UQ) of materials’ microstructure, as well as establishing process-structure–property (PSP) relationships for inverse material design. In this study, an efficient UP technique built on the idea of changing probability measures and a deep generative unsupervised representative machine learning method for microstructure-based design of thermal conductivity of materials is proposed. Probability measures are used to represent microstructure space, and Wasserstein metrics are used to test the efficiency of the UP method. By using deep Variational AutoEncoder (VAE), we identify the correlations between the material/process parameters and the thermal conductivity of heterogeneous dual-phase microstructures. Through high-throughput screening, UP, and the deep-generative VAE method, PSP relationships that are too complex can be revealed by exploiting the materials’ design space with an emphasis on microstructures. As a last point, we demonstrate generative machine learning serves as a useful tool for inverse microstructure-centered materials design, and we demonstrate this by examining the inverse design of thermal conductivity in nano-structured materials. Here, the results reveal the effects of morphology, volume fraction, characteristic length scale, and the individual thermal diffusivity of phases on the thermal conductivity of dual-phase alloys. Our findings emphasize the advantages of high-throughput phase-field modeling and generative deep learning for linking PSP and inverse microstructure-centered materials design.

36 MATERIALS SCIENCE↗

A Gaussian process autoregressive model capturing microstructure evolution paths in a Ni–Mo–Nb alloy

Additive manufacturing is increasingly being employed to produce components of complex geometries in structural alloys because of the expected energy savings associated with the near-net-shape capability and the ability to build in novel internal features that are not possible with many conventional manufacturing approaches. However, because of the extreme thermal conditions encountered, the non-equilibrium microstructures produced during powder bed-based additive manufacturing processes must be subjected to custom post-heat treatment processes to recover the target mechanical properties. Phase-field models and simulation techniques have matured to a state where the microstructure evolution paths, and the morphologies of the resulting precipitate phases can be predicted reasonably accurately, considering alloy-specific thermodynamic and kinetic aspects of the nucleation and growth processes. However, phase-field simulations are computationally intensive, which precludes the ability to apply the simulations directly to the length scale of the entire component. Therefore, it is highly desirable to develop low-computational-cost surrogate models that effectively capture the physics at the microstructural length scale, while facilitating the design of optimized processing conditions resulting in location-specific targeted microstructures at the component scale. The work presented here demonstrates the application of the materials knowledge system framework to develop a surrogate model that effectively captures the microstructural path during annealing of a Ni–Mo–Nb alloy containing different Mo and Nb compositions known to segregate during solidification under additive manufacturing conditions. Specifically, the surrogate model built in this work is based on a Gaussian process autoregressive model informed by statistical representation of simulated microstructures using two-point correlations and dimensionality reduction through principal component analysis. In conclusion, this surrogate model is shown to capture the bifurcation of the microstructural path during precipitation, which yields a microstructure dominated by the $\gamma^{\prime\prime}$ phase at high Nb concentrations and the $\delta$ phase at low Nb concentrations.

36 MATERIALS SCIENCE↗

A generalizable machine learning-assisted fast Fourier transform algorithm to simulate the large strain phenomena in polycrystalline materials

Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.

36 MATERIALS SCIENCE↗

The Challenges of Modeling Defect Behavior and Plasticity across Spatial and Temporal Scales: A Case Study of Metal Bilayer Impact

Atomistic molecular dynamics (MD) and a microstructural dislocation density-based crystalline plasticity (DCP) framework were used together across time scales varying from picoseconds to nanoseconds and length scales spanning from angstroms to micrometers to model a buried copper–nickel interface subjected to high strain rates. The nucleation and evolution of defects, such as dislocations and stacking faults, as well as large inelastic strain accumulations and wave-induced stress reflections were physically represented in both approaches. Both methods showed similar qualitative behavior, such as defects originating along the impactor edges, a dominance of Shockley partial dislocations, and non-continuous dislocation distributions across the buried interface. The favorable comparison between methods justifies assumptions used in both, to model phenomena, such as the nucleation and interactions of single defects and partials with reflected tensile waves, based on MD predictions, which are consistent with the evolution of perfect and partial dislocation densities as predicted by DCP. This substantiates how the nanoscale as modeled by MD is representative of microstructural behavior as modeled by DCP.

36 MATERIALS SCIENCE↗

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗

Fused Deposition Modeling Additive Manufacturing of Carbonized Structures via Waste-Enhanced Filaments

This report details the development, characterization, and use of coal-enhanced composite materials in additive manufacturing applications. High coal loading formulations—containing up to 70 wt.% coal—were successfully extruded and processed using commercially available 3D printers. Extensive experimental testing was conducted to assess the mechanical, thermal, and microstructural properties of the composites. In parallel, multi-scale computational modeling was employed to elucidate atomistic interactions and evaluate the effects of printing-induced defects on structural performance. Large-scale printability trials demonstrated the feasibility of fabricating complex components for tooling and construction applications, including wind turbine blade molds and modular wall sections. Techno-economic analyses demonstrated the cost-effectiveness and scalability of coal-enhanced composites for large-scale additive manufacturing applications such as wind turbine blade tooling.

01 COAL, LIGNITE, AND PEAT↗

Micro-structural features and material properties impact on adhesive metal joints via computational modeling and machine learning

The quality of structural bonding in practical applications depends on various factors arising from materials, pre-processing conditions, and manufacturing. Understanding how these factors influence bonding performance and determining their relative importance are of significant interest. Thus, this study evaluates the effects of microstructural features and material properties on the structural strength of adhesively-bonded metal joints at the submillimeter scale, utilizing a combination of Finite Element Modeling (FEM) and Machine Learning (ML) with Gradient Boosting Regression (GBR). The microstructural features include adhesive thickness, internal voids within the adhesive, adherend-adhesive interfacial voids, void size and volume fraction, and surface roughness. The material properties include the constitutive behavior of the adhesive, as well as the adherend-adhesive interfacial strength and fracture energy. The changes in structural strength and morphologies of the bonded metal structures with respect to different microstructural features and material properties were clarified by FEM. By further leveraging ML-GBR, the sequence of importance of these factors affecting bonding performance across various scenarios was summarized. This work provides valuable insights into the development of improved structural bonding for adhesive joints in industries such as automotive , aerospace, and beyond.

36 MATERIALS SCIENCE↗

Microstructural-Level Fuel Performance Modeling of U-Mo Monolithic Fuel

As the physics that governs the microstructural evolution of nuclear fuel span various time and spatial scales, to fully understand the fuel behavior inevitably involves atomic to mesoscale resolution that can be difficult to determine experimentally. Microstructural-level modeling and simulations can be used to develop physics-based materials models that can provide physical understanding to inform fabrication process control, as well as a valuable feedback mechanism between post-irradiation examination (PIE) results and fabrication parameters. In accordance with the program schedule, the primary goals of the microstructure modeling effort are to: 1. Address critical microstructural questions and provide practical guidance to the fabricator via the fuel product specification; 2. Provide mechanistic inputs for the existing fuel performance code to improve its descriptive and predictive capability at the macroscopic scale. In fiscal year (FY)-21, the work scope consisted of six main facets: (1) the effect of carbides on fuel performance; (2) gas diffusivity in different phases; (3) integration of microstructural fuel performance modeling; (4) property degradation; (5) irradiation creep; and (6) historical analysis of microstructure data. Brief summaries of each are included below.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications

Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogeneities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications

Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogeneities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications

Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogeneities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications

Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogenities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications

Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogenities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A data-driven framework for permeability prediction of natural porous rocks via microstructural characterization and pore-scale simulation

Understanding the microstructure–property relationships of porous media is of great practical significance, based on which macroscopic physical properties can be directly derived from measurable microstructural informatics. However, establishing reliable microstructure–property mappings in an explicit manner is difficult, due to the intricacy, stochasticity, and heterogeneity of porous microstructures. In this paper, a data-driven computational framework is presented to investigate the inherent microstructure–permeability linkage for natural porous rocks, where multiple techniques are integrated together, including microscopy imaging, stochastic reconstruction, microstructural characterization, pore-scale simulation, feature selection, and data-driven modeling. A large number of 3D digital rocks with a wide porosity range are acquired from microscopy imaging and stochastic reconstruction techniques. A broad variety of morphological descriptors are used to quantitatively characterize pore microstructures from different perspectives, and they compose the raw feature pool for feature selection. Here high-fidelity lattice Boltzmann simulations are conducted to resolve fluid flow passing through porous media, from which reliable permeability references are obtained. The optimal feature set that best represents permeability is identified through a performance-oriented feature selection process, upon which a cost-effective surrogate model is rapidly fitted to approximate the microstructure-permeability mapping via data-driven modeling. This surrogate model exhibits great advantages over empirical/analytical formulas in terms of prediction accuracy and generalization capacity, which can predict reliable permeability values spanning four orders of magnitude. Besides, feature selection also greatly enhances the interpretability of the data-driven prediction model, from which new insights into the mechanism of how microstructural characteristics determine intrinsic permeability are obtained.

58 GEOSCIENCES↗