An AI-assisted Framework for Image-Based Multiphase Flow Simulations
Abstract. Image-based pore-scale models are critical for studying pore-scale multiphase displacement physics and connecting pore-scale processes to Darcy-scale constitutive relationships. However, systematic image-based modeling studies remain challenging due to the fragmentation of image acquisition and processing approaches, pore-scale simulators, and post-processing tools. To address this challenge, we develop a human-in-the-loop, AI-assisted workflow that integrates complex image-based modeling tasks into a state-aware framework for multiphase flow simulations. External or local AI tools (i.e., large language models) are used to translate natural-language instructions into validated model actions and parameter updates, while deterministic codes are implemented to execute image generation/import and processing, quasi-static and dynamic two-phase displacement simulations, postprocessing. We demonstrate the framework through deterministic two-phase displacement simulations, direct comparison with reported primary-drainage measurements, and ensemble analyses across porous-medium realizations. Our simulations show that: (1) Sample-scale constitutive relationships (e.g., capillary-pressure–saturation and fluid–fluid-interfacial-area–saturation curves) from dynamic simulations deviate from their quasi-static counterparts and do not collapse onto universal curves, revealing flow-regime-dependent nonequilibrium behavior. (2) Comparison with reported primary-drainage measurements further shows that model–experiment agreement is strongly affected by pore-body and pore-throat geometric representations. (3) Quasi-static capillary-pressure and interfacial-area relationships are relatively insensitive to microstructural variability, whereas relative permeability is substantially more variable because of its dependence on connected flow pathways. (4) For dynamic displacements, constitutive relationships remain comparatively consistent within each flow regime but distinctly separated among regimes across overlapping saturation ranges. These results highlight the importance of developing regime-aware Darcy-scale constitutive relationships and demonstrate a reproducible framework for model benchmarking, model–data evaluation, and uncertainty-aware pore-scale analysis.
The manuscript egusphere-2026-6229 reports on the development of an AI-driven workflow to automate a dynamic pore network model simulation for multiphase flow in porous media on the basis of a segmented 3D image of a porous structure by using AI/large language models.
The paper is well-written and plausible, numerical results validated against experimental data, and publication is supported.
The image processing itself is not part of the workflow. However, image segmentation and subsequent network extraction is a primary source of uncertainty in the whole workflow. Therefore, recently there are publications on grey-level LBM simulations which avoid this issue to some extent. It would be good if the authors could comment on that. Furthermore, for pore network extraction it depends a bit on the type of porous media, but for rocks the watershed method produced somewhat less reliable results than the maximum ball method by Dong & Blunt. The underlying reason is that the watershed method uses maxima in the Euclidean distance map as seed points but most pores in e.g. sandstone rock have already multiple maxima within a single pore, which leads to topological issues. Anyway, this is a detail, but it has quite some impact on the result, therefore, would be good to mention that.
For the dynamic PNM it would be good to provide some more literature references, e.g. the early work by
There are probably more, but these are the classical papers on dynamic PNM.
Fig. 3 – maybe it would be nice to map the results on the classical Lenormand diagram (mobility ratio vs. capillary number)
This is more of a suggestion: It might be nice to also show the Euler characteristic – because together with saturation, interfacial area and mean curvature we would have the complete Minkowski functionals, which are the capillary state variables (Phys. Rev. Fluids 3(8), 084306, 2018.). It would be in general very interesting to follow trajectories of displacement processes in Minkowski space, see. E.g. work by the R. T. Johns group at PennState (various papers, e.g. Transport in Porous Media 153 (8), 103, 2026), Parente et al. (International Journal of Hydrogen Energy 229, 154730, 2026.), Wildenschild group (Physics of Fluids 38 (2),2026).