A domain decomposition approach to pore-network modeling of porous media flow

We propose a domain-decomposition pore-network method (DD-PNM) for modeling single-phase Stokes flow in porous media. The method combines the accuracy of finite-element discretizations on body-fitted meshes within pore subdomains with a sparse global coupling enforced through interface unknowns. Local Dirichlet-to-Neumann operators are precomputed from finite-element solutions for each pore subdomain, enabling a global Schur-complement system defined solely on internal interfaces. Rigorous mathematical analysis establishes solvability and discrete mass conservation of the global system. Moreover, we constructively recover classical pore-network models by fitting half-throat conductivities to local Dirichlet-to-Neumann maps, providing a principled bridge between mesh-based and network-based frameworks. Numerical results are presented to demonstrate the validity and effectiveness of the overall methodology.