Calculation of tensorial flow properties on pore level: Exploring the influence of boundary conditions on the permeability of three-dimensional stochastic reconstructions

Three-dimensional construction of hyperuniform, nonhyperuniform, and antihyperuniform disordered heterogeneous materials and their transport properties via spectral density functions

Rigorous theories connecting physical properties of a heterogeneous material to its microstructure offer a promising avenue to guide the computational material design and optimization. The spectral density function χ[over ̃]_{_{V}}(k), which can be obtained experimentally from scattering data, enables accurate determination of various transport and wave propagation characteristics, including the time-dependent diffusion spreadability S(t) and effective dynamic dielectric constant ε_{e} for electromagnetic wave propagation. Moreover, χ[over ̃]_{_{V}}(k) determines rigorous upper bounds on the fluid permeability K. Given the importance of χ[over ̃]_{_{V}}(k), we present here an efficient Fourier-space based computational framework to construct three-dimensional (3D) statistically isotropic two-phase heterogeneous materials corresponding to targeted spectral density functions. In particular, we employ a variety of analytical functional forms for χ[over ̃]_{_{V}}(k) that satisfy all known necessary conditions to construct disordered stealthy hyperuniform, standard hyperuniform, nonhyperuniform, and antihyperuniform two-phase heterogeneous material systems at varying phase volume fractions. We show that by tuning the correlations in the system across length scales via the targeted functions, one can generate a rich spectrum of distinct structures within each of the above classes of materials. Importantly, we present the first realization of antihyperuniform two-phase heterogeneous materials in 3D, which are characterized by autocovariance function χ_{_{V}}(r) with a power-law tail, resulting in microstructures that contain clusters of dramatically different sizes and morphologies. We also determine the diffusion spreadability S(t) and estimate the fluid permeability K associated with all of the constructed materials directly from the corresponding spectral densities. Although it is well established that the long-time asymptotic scaling behavior of S(t) only depends on the functional form of χ[over ̃]_{_{V}}(k), with the stealthy hyperuniform and antihyperuniform media, respectively, achieving the most and least efficient transport, we find that varying the length-scale parameter within each class of χ[over ̃]_{_{V}}(k) functions can also lead to orders of magnitude variation of S(t) at intermediate and long time scales. Moreover, we find that increasing the solid volume fraction ϕ_{1} and correlation length a in the constructed media generally leads to a decrease in the dimensionless fluid permeability K/a^{2}, while the antihyperuniform media possess the largest K/a^{2} among the four classes of materials with the same ϕ_{1} and a. These results indicate the feasibility of employing parameterized spectral densities for designing composites with targeted transport properties.

Evaluation of three-point correlation functions from structural images on CPU and GPU architectures: Accounting for anisotropy effects

Structures, or spatial arrangements of matter and energy, including some fields (e.g., velocity or pressure) are ubiquitous in research applications and frequently require description for subsequent analysis, or stochastic reconstruction from limited data. The classical descriptors are two-point correlation functions (CFs), but the computation of three-point statistics is known to be advantageous in some cases as they can probe non-Gaussian signatures, not captured by their two-point counterparts. Moreover, n-point CFs with n≥3 are believed to possess larger information content and provide more information about studied structures. In this paper, we have developed algorithms and code to compute S_{3},C_{3},F_{sss},F_{ssv}, and F_{svv} with a right-angle and arbitrary triangle pattern. The former was believed to be faster to compute, but with the help of precomputed regular positions we achieved the same speed for arbitrary pattern. In this work we also implement and demonstrate computations of directional three-point CFs-for this purpose right-triangular pattern seems to be superior due to explicit orientation and high coverage. Moreover, we assess the errors in CFs' evaluation due to image or pattern rotations and show that they have minor effect on accuracy of computations. The execution times of our algorithms for the same number of samples are orders of magnitude lower than in existing published counterparts. We show that the volume of data produced gets unwieldy very easily, especially if computations are performed in frequency domain. For these reasons until information content of different sets of correlation functions with different "n-pointness" is known, advantages of CFs with n>3 are not clear. Nonetheless, developed algorithms and code are universal enough to be easily extendable to any n with increasing computational and random access memory (RAM) burden. All results are available as part of open-source package correlationfunctions.jl [V. Postnicov et al., Comput. Phys. Commun. 299, 109134 (2024)10.1016/j.cpc.2024.109134.]. As described in this paper, three-point CFs computations can be immediately applied in a great number of research applications, for example: (1) flow and transport velocity fields analysis or any data with non-Gaussian signatures, (2) deep learning for structural and physical properties, and (3) structure taxonomy and categorization. In all these and numerous other potential cases the ability to compute directional three-point functions may be crucial. Notably, the organization of the code functions allows computation of cross correlation, i.e., one can compute three-point CFs for multiphase images (while binary structures were used in this paper for simplicity of explanations).

Measuring structural nonstationarity: The use of imaging information to quantify homogeneity and inhomogeneity

Heterogeneity is the concept we encounter in numerous research areas and everyday life. While "not mixing apples and oranges" is easy to grasp, a more quantitative approach to such segregation is not always readily available. Consider the problem from a different angle: To what extent does one have to make apples more orange and oranges more "apple-shaped" to put them into the same basket (according to their appearance alone)? This question highlights the central problem of the blurred interface between heterogeneous and homogeneous, which also depends on the metrics used for its identification. This work uncovers the physics of structural stationarity quantification, based on correlation functions (CFs) and clustering based on CFs different between image subregions. By applying the methodology to a wide variety of synthetic and real images of binary porous media, we confirmed computationally that only periodically unit-celled structures and images produced by stationary processes with resolutions close to infinity are strictly stationary. Natural structures without recurring unit cells are only weakly stationary. We established a physically meaningful definition for these stationarity types and their distinction from nonstationarity. In addition, the importance of information content of the chosen metrics is highlighted and discussed. We believe the methodology as proposed in this contribution will find its way into numerous research areas dealing with materials, structures, and measurements and modeling based on structural imaging information.

Computational design of anisotropic stealthy hyperuniform composites with engineered directional scattering properties

Disordered hyperuniform materials are an emerging class of exotic amorphous states of matter that endow them with singular physical properties, including large isotropic photonic band gaps, superior resistance to fracture, and nearly optimal electrical and thermal transport properties, to name but a few. Here we generalize the Fourier-space-based numerical construction procedure for designing and generating digital realizations of isotropic disordered hyperuniform two-phase heterogeneous materials (i.e., composites) developed by Chen and Torquato [Acta Mater. 142, 152 (2018)1359-645410.1016/j.actamat.2017.09.053] to anisotropic microstructures with targeted spectral densities. Our generalized construction procedure explicitly incorporates the vector-dependent spectral density function χ[over ̃]_{_{V}}(k) of arbitrary form that is realizable. We demonstrate the utility of the procedure by generating a wide spectrum of anisotropic stealthy hyperuniform microstructures with χ[over ̃]_{_{V}}(k)=0 for k∈Ω, i.e., complete suppression of scattering in an "exclusion" region Ω around the origin in Fourier space. We show how different exclusion-region shapes with various discrete symmetries, including circular-disk, elliptical-disk, square, rectangular, butterfly-shaped, and lemniscate-shaped regions of varying size, affect the resulting statistically anisotropic microstructures as a function of the phase volume fraction. The latter two cases of Ω lead to directionally hyperuniform composites, which are stealthy hyperuniform only along certain directions and are nonhyperuniform along others. We find that while the circular-disk exclusion regions give rise to isotropic hyperuniform composite microstructures, the directional hyperuniform behaviors imposed by the shape asymmetry (or anisotropy) of certain exclusion regions give rise to distinct anisotropic structures and degree of uniformity in the distribution of the phases on intermediate and large length scales along different directions. Moreover, while the anisotropic exclusion regions impose strong constraints on the global symmetry of the resulting media, they can still possess structures at a local level that are nearly isotropic. Both the isotropic and anisotropic hyperuniform microstructures associated with the elliptical-disk, square, and rectangular Ω possess phase-inversion symmetry over certain range of volume fractions and a percolation threshold ϕ_{c}≈0.5. On the other hand, the directionally hyperuniform microstructures associated with the butterfly-shaped and lemniscate-shaped Ω do not possess phase-inversion symmetry and percolate along certain directions at much lower volume fractions. We also apply our general procedure to construct stealthy nonhyperuniform systems. Our construction algorithm enables one to control the statistical anisotropy of composite microstructures via the shape, size, and symmetries of Ω, which is crucial to engineering directional optical, transport, and mechanical properties of two-phase composite media.

Robust surface-correlation-function evaluation from experimental discrete digital images

Correlation functions (CFs) are universal structural descriptors; surface-surface F_{ss} and surface-void F_{sv} CFs are a subset containing additional information about the interface between the phases. The description of the interface between pores and solids in porous media is of particular importance and recently Ma and Torquato [Phys. Rev. E 98, 013307 (2018)2470-004510.1103/PhysRevE.98.013307] proposed an elegant way to compute these functions for a wide variety of cases. However, their "continuous" approach is not always applicable to digital experimental 2D and 3D images of porous media as obtained using x-ray tomography or scanning electron microscopy due to nonsingularities in chemical composition or local solid material's density and partial volume effects. In this paper we propose to use edge-detecting filters to compute surface CFs in the "digital" fashion directly in the images. Computed this way, surface correlation functions are the same as analytically known for Poisson disks in case the resolution of the image is adequate. Based on the multiscale image analysis we developed a C_{0.5} criterion that can predict if the imaging resolution is enough to make an accurate evaluation of the surface CFs. We also showed that in cases when the input image contains all major features, but do not pass the C_{0.5} criterion, it is possible with the help of image magnification to sample CFs almost similar to those obtained for high-resolution image of the same structure with high C_{0.5}. The computational framework as developed here is open source and available within the CorrelationFunctions.jl package developed by our group. Our "digital" approach was applied to a wide variety of real porous media images of different quality. We discuss critical aspects of surface correlation functions computations as related to different applications. The developed methodology allows applying surface CFs to describe the structure of porous materials based on their experimental images and enhance stochastic reconstructions or super-resolution procedures, or serve as an efficient metrics in machine learning applications due to computationally effective GPU implementation.

Pore-network extraction using discrete Morse theory: Preserving the topology of the pore space

Pore-scale modeling based on the 3D structural information of porous materials has enormous potential in assessing physical properties beyond the capabilities of laboratory methods. Such capabilities are pricey in terms of computational expenses, and this limits the applicability of the direct simulations to a small volume and requires high-performance computational resources, especially for multiphase flow simulations. The only pore-scale technique capable of dealing with large representative volumes of porous samples is pore-network (PNM) based modeling. The problem of the PNM approach is that 3D pore geometry first needs to be simplified into a graph of pores and throats that conserve topological and geometrical properties of the original 3D image. While significant progress has been achieved in terms of geometry representation, no methodology provides full conservation of the topological features of the pore structure. In this paper we present a pore-network extraction algorithm for binary 3D images based on discrete Morse theory and persistent homology that by design targets topology preservation. In addition to methodological developments, we also clarify the relationship between topological characteristics of constructed Morse chain complex and pore-network elements. We show that the Euler numbers calculated for PNMs based on our methodology coincide with those obtained using the direct topological analysis. The characteristics of the extracted pore network are calculated for several 3D porous binary images and compared with the results of maximum inscribed balls-based and watershed-based approaches as well as a hybrid approach to support our methodology.

Improved recurrent generative model for reconstructing large-size porous media from two-dimensional images

Modeling the three-dimensional (3D) structure from a given 2D image is of great importance for analyzing and studying the physical properties of porous media. As an intractable inverse problem, many methods have been developed to address this fundamental problems over the past decades. Among many methods, the deep learning-(DL) based methods show great advantages in terms of accuracy, diversity, and efficiency. Usually, the 3D reconstruction from the 2D slice with a larger field-of-view is more conducive to simulate and analyze the physical properties of porous media accurately. However, due to the limitation of reconstruction ability, the reconstruction size of most widely used generative adversarial network-based model is constrained to 64^{3} or 128^{3}. Recently, a 3D porous media recurrent neural network based method (namely, 3D-PMRNN) (namely 3D-PMRNN) has been proposed to improve the reconstruction ability, and thus the reconstruction size is expanded to 256^{3}. Nevertheless, in order to train these models, the existed DL-based methods need to down-sample the original computed tomography (CT) image first so that the convolutional kernel can capture the morphological features of training images. Thus, the detailed information of the original CT image will be lost. Besides, the 3D reconstruction from a optical thin section is not available because of the large size of the cutting slice. In this paper, we proposed an improved recurrent generative model to further enhance the reconstruction ability (512^{3}). Benefiting from the RNN-based architecture, the proposed model requires only one 3D training sample at least and generates the 3D structures layer by layer. There are three more improvements: First, a hybrid receptive field for the kernel of convolutional neural network is adopted. Second, an attention-based module is merged into the proposed model. Finally, a useful section loss is proposed to enhance the continuity along the Z direction. Three experiments are carried out to verify the effectiveness of the proposed model. Experimental results indicate the good reconstruction ability of proposed model in terms of accuracy, diversity, and generalization. And the effectiveness of section loss is also proved from the perspective of visual inspection and statistical comparison.

Understanding degeneracy of two-point correlation functions via Debye random media

It is well known that the degeneracy of two-phase microstructures with the same volume fraction and two-point correlation function S_{2}(r) is generally infinite. To elucidate the degeneracy problem explicitly, we examine Debye random media, which are entirely defined by a purely exponentially decaying two-point correlation function S_{2}(r). In this work, we consider three different classes of Debye random media. First, we generate the "most probable" class using the Yeong-Torquato construction algorithm [Yeong and Torquato, Phys. Rev. E 57, 495 (1998)1063-651X10.1103/PhysRevE.57.495]. A second class of Debye random media is obtained by demonstrating that the corresponding two-point correlation functions are effectively realized in the first three space dimensions by certain models of overlapping, polydisperse spheres. A third class is obtained by using the Yeong-Torquato algorithm to construct Debye random media that are constrained to have an unusual prescribed pore-size probability density function. We structurally discriminate these three classes of Debye random media from one another by ascertaining their other statistical descriptors, including the pore-size, surface correlation, chord-length probability density, and lineal-path functions. We also compare and contrast the percolation thresholds as well as the diffusion and fluid transport properties of these degenerate Debye random media. We find that these three classes of Debye random media are generally distinguished by the aforementioned descriptors, and their microstructures are also visually distinct from one another. Our work further confirms the well-known fact that scattering information is insufficient to determine the effective physical properties of two-phase media. Additionally, our findings demonstrate the importance of the other two-point descriptors considered here in the design of materials with a spectrum of physical properties.

Three-dimensional reconstruction of porous media using super-dimension-based adjacent block-matching algorithm

As porous media play an essential role in a variety of industrial applications, it is essential to understand their physical properties. Nowadays, the super-dimensional (SD) reconstruction algorithm is used to stochastically reconstruct a three-dimensional (3D) structure of porous media from a given two-dimensional image. This algorithm exhibits superiority in accuracy compared with classical algorithms because it learns information from the real 3D structure. However, owing to the short development time of the SD algorithm, it also has some limitations, such as inexact porosity characterization, long run time, blocking artifacts, and suboptimal accuracy that may be improved. To mitigate these limitations, this study presents the design of a special template that contains two parts of data (i.e., adjacent blocks and a central block); the proposed method matches adjacent blocks during reconstruction and assigns the matched central block to the area to be reconstructed. Furthermore, we design two important mechanisms during reconstruction: one for block matching and the other for porosity control. To verify the effectiveness of the proposed method compared with an existing SD method, both methods were tested on silica particle material and three homogeneous sandstones with different porosities; meanwhile, we compared the proposed method with a multipoint statistics method and a simulated annealing method. The reconstructed results were then compared with the target both visually and quantitatively. The experimental results indicate that the proposed method can overcome the aforementioned limitations and further improve the accuracy of existing methods. This method achieved 4-6 speedup factor compared with the traditional SD method.

Application of percolation, critical-path, and effective-medium theories for calculation of two-phase relative permeability

There has been active development of numerical pore-network simulation of two-phase immiscible flows in porous media in recent years. These models allow for generation of capillary pressure and relative permeability curves. However, percolation models provide an efficient alternative, with reduced reliance on numerical techniques. Implementation of effective medium or critical path theory along with the percolation model allows for evaluation of the relative permeability curves. Both approximations failed to match the irreducible water saturation for water relative permeability. While the effective medium approximation poorly matches the pore network simulator, the critical path approximation is shown to match the result of the oil relative permeability. Despite the difference in end points, there is qualitative agreement between critical path approximation and the pore network simulator. Moreover, observed differences are not necessarily a drawback due to important boundary effects as discussed in the paper. Our results indicate that percolation-theory based predictions have the potential to become an efficient tool for upscaling by computing two-phase flow properties for numerous porosity subdomains.

Generation and structural characterization of Debye random media

In their seminal paper on scattering by an inhomogeneous solid, Debye and coworkers proposed a simple exponentially decaying function for the two-point correlation function of an idealized class of two-phase random media. Such Debye random media, which have been shown to be realizable, are singularly distinct from all other models of two-phase media in that they are entirely defined by their one- and two-point correlation functions. To our knowledge, there has been no determination of other microstructural descriptors of Debye random media. In this paper, we generate Debye random media in two dimensions using an accelerated Yeong-Torquato construction algorithm. We then ascertain microstructural descriptors of the constructed media, including their surface correlation functions, pore-size distributions, lineal-path function, and chord-length probability density function. Accurate semianalytic and empirical formulas for these descriptors are devised. We compare our results for Debye random media to those of other popular models (overlapping disks and equilibrium hard disks) and find that the former model possesses a wider spectrum of hole sizes, including a substantial fraction of large holes. Our algorithm can be applied to generate other models defined by their two-point correlation functions, and their other microstructural descriptors can be determined and analyzed by the procedures laid out here.

Probing information content of hierarchical n-point polytope functions for quantifying and reconstructing disordered systems

Disordered systems are ubiquitous in physical, biological, and material sciences. Examples include liquid and glassy states of condensed matter, colloids, granular materials, porous media, composites, alloys, packings of cells in avian retina, and tumor spheroids, to name but a few. A comprehensive understanding of such disordered systems requires, as the first step, systematic quantification, modeling, and representation of the underlying complex configurations and microstructure, which is generally very challenging to achieve. Recently, we introduced a set of hierarchical statistical microstructural descriptors, i.e., the "n-point polytope functions" P_{n}, which are derived from the standard n-point correlation functions S_{n}, and successively included higher-order n-point statistics of the morphological features of interest in a concise, explainable, and expressive manner. Here we investigate the information content of the P_{n} functions via optimization-based realization rendering. This is achieved by successively incorporating higher-order P_{n} functions up to n=8 and quantitatively assessing the accuracy of the reconstructed systems via unconstrained statistical morphological descriptors (e.g., the lineal-path function). We examine a wide spectrum of representative random systems with distinct geometrical and topological features. We find that, generally, successively incorporating higher-order P_{n} functions and, thus, the higher-order morphological information encoded in these descriptors leads to superior accuracy of the reconstructions. However, incorporating more P_{n} functions into the reconstruction also significantly increases the complexity and roughness of the associated energy landscape for the underlying stochastic optimization, making it difficult to convergence numerically.

Three-dimensional multiscale fusion for porous media on microtomography images of different resolutions

Accurately acquiring the three-dimensional (3D) image of a porous medium is an imperative issue for the prediction of multiple physical properties. Considering the inherent nature of the multiscale pores contained in porous media such as tight sandstones, to completely characterize the pore structure, one needs to scan the microstructure at different resolutions. Specifically, low-resolution (LR) images cover a larger field of view (FOV) of the sample, but are lacking small-scale features, whereas high-resolution (HR) images contain ample information, but sometimes only cover a limited FOV. To address this issue, we propose a method for fusing the spatial information from a two-dimensional (2D) HR image into a 3D LR image, and finally reconstructing an integrated 3D structure with added fine-scale features. In the fusion process, the large-scale structure depicted by the 3D LR image is fixed as background and the 2D image is utilized as training image to reconstruct a small-scale structure based on the background. To assess the performance of our method, we test it on a sandstone scanned with low and high resolutions. Statistical properties between the reconstructed image and the target are quantitatively compared. The comparison indicates that the proposed method enables an accurate fusion of the LR and HR images because the small-scale information is precisely reproduced within the large one.