Reconstruction of porous media from extremely limited information using conditional generative adversarial networks

Oral delivery of posaconazole-loaded phospholipid-based nanoformulation: Preparation and optimization using design of experiments, machine learning, and TOPSIS

Phospholipid-based nanosystems show promising potentials for oral administration of hydrophobic drugs. The study introduced a novel approach to optimize posaconazole-loaded phospholipid-based nanoformulation using the design of experiments, machine learning, and Technique for Order of Preference by Similarity to the Ideal Solution. These approaches were used to investigate the impact of various variables on the encapsulation efficiency (EE), particle size, and polydispersity index (PDI). The optimized formulation, with %EE of ∼ 74 %, demonstrated a particle size and PDI of 107.7 nm and 0.174, respectively. The oral pharmacokinetic profiles of the posaconazole suspension, empty nanoformulation + drug suspension, and drug-loaded nanoformulation were evaluated. The nanoformulation significantly increased maximum plasma concentration and the area under the drug plasma concentration-time curve (∼3.9- and 6.2-fold, respectively) and could be administered without regard to meals. MTT and histopathological examinations were carried out to evaluate the safety of the nanoformulation and results exhibited no significant toxicity. Lymphatic transport was found to be the main mechanism of oral delivery. Caco-2 cell studies demonstrated that the mechanism of delivery was not based on an increase in cellular uptake. Our study represents a promising strategy for the development of phospholipid-based nanoformulations as efficient and safe oral delivery systems.

Three-dimensional reconstruction of granular porous media based on deep generative models

Reconstruction of microstructure in granular porous media, which can be viewed as granular assemblies, is crucial for studying their characteristics and physical properties in various fields concerned with the behavior of such media, including petroleum geology and computational materials science. In spite of the fact that many existing studies have investigated grain reconstruction, most of them treat grains as simplified individuals for discrete reconstruction, which cannot replicate the complex geometrical shapes and natural interactions between grains. In this work, a hybrid generative model based on a deep-learning algorithm is proposed for high-quality three-dimensional (3D) microstructure reconstruction of granular porous media from a single two-dimensional (2D) slice image. The method extracts 2D prior information from the given image and generates the grain set as a whole. Both a self-attention module and effective pattern loss are introduced in a bid to enhance the reconstruction ability of the model. Samples with grains of varied geometrical shapes are utilized for the validation of our method, and experimental results demonstrate that our proposed approach can accurately reproduce the complex morphology and spatial distribution of grains without any artificiality. Furthermore, once the model training is complete, rapid end-to-end generation of diverse 3D realizations from a single 2D image can be achieved.

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

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.

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.

Super-resolution of real-world rock microcomputed tomography images using cycle-consistent generative adversarial networks

Digital rock imaging plays an important role in studying the microstructure and macroscopic properties of rocks, where microcomputed tomography (MCT) is widely used. Due to the inherent limitations of MCT, a balance should be made between the field of view (FOV) and resolution of rock MCT images-a large FOV at low resolution (LR) or a small FOV at high resolution (HR). However, large FOV and HR are both expected for reliable analysis results in practice. Super-resolution (SR) is an effective solution to break through the mutual restriction between the FOV and resolution of rock MCT images, for it can reconstruct an HR image from a LR observation. Most of the existing SR methods cannot produce satisfactory HR results on real-world rock MCT images. One of the main reasons for this is that paired images are usually needed to learn the relationship between LR and HR rock images. However, it is challenging to collect such a dataset in a real scenario. Meanwhile, the simulated datasets may be unable to accurately reflect the model in actual applications. To address these problems, we propose a cycle-consistent generative adversarial network (CycleGAN)-based SR approach for real-world rock MCT images, namely, SRCycleGAN. In the off-line training phase, a set of unpaired rock MCT images is used to train the proposed SRCycleGAN, which can model the mapping between rock MCT images at different resolutions. In the on-line testing phase, the resolution of the LR input is enhanced via the learned mapping by SRCycleGAN. Experimental results show that the proposed SRCycleGAN can greatly improve the quality of simulated and real-world rock MCT images. The HR images reconstructed by SRCycleGAN show good agreement with the targets in terms of both the visual quality and the statistical parameters, including the porosity, the local porosity distribution, the two-point correlation function, the lineal-path function, the two-point cluster function, the chord-length distribution function, and the pore size distribution. Large FOV and HR rock MCT images can be obtained with the help of SRCycleGAN. Hence, this work makes it possible to generate HR rock MCT images that exceed the limitations of imaging systems on FOV and resolution.

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

While it is well known that permeability is a tensorial property, it is usually reported as a scalar property or only diagonal values are reported. However, experimental evaluation of tensorial flow properties is problematic. Pore-scale modeling using three-dimensional (3D) images of porous media with subsequent upscaling to a continuum scale (homogenization) is a valuable alternative. In this study, we explore the influence of different types of boundary conditions on the external walls of the representative modeling domain along the applied pressure gradient on the magnitude and orientation of the computed permeability tensor. To implement periodic flow boundary conditions, we utilized stochastic reconstruction methodology to create statistically similar (to real porous media structures) geometrically periodic 3D structures. Stochastic reconstructions are similar to encapsulation of the porous media into statistically similar geometrically periodic one with the same permeability tensor. Seven main boundary conditions (BC) were implemented: closed walls, periodic flow, slip on the walls, linear pressure, translation, symmetry, and immersion. The different combinations of BCs amounted to a total number of 15 BC variations. All these BCs significantly influenced the resulting tensorial permeabilities, including both magnitude and orientation. Periodic boundary conditions produced the most physical flow patterns, while other classical BCs either suppressed crucial transversal flows or resulted in unphysical currents. Our results are crucial to performing flow properties upscaling and will be relevant to computing not only single-phase but also multiphase flow properties. Moreover, other calculation of physical properties such as some mechanical, transport, or heat conduction properties may benefit from the technique described in this study.