Modeling groundwater flow with random hydraulic conductivity using radial basis function partition of unity method

Abstract

Simulating groundwater flows in heterogeneous aquifers is one of the most widely studied problems. The heterogeneity is modeled through random hydraulic conductivity fields log-normally distributed. In this paper, we aim to generate the realization of the log-normal hydraulic conductivity by summing up a finite number of random periodic modes with the Kraichnan algorithm. To address Neumann conditions, we propose a change of variables that results in Dirichlet boundary conditions along all boundaries of the computational domain. Then, we suggest using the radial basis function partition of unity method (RBFPUM) and RBFPUM with QR factorization (RBFPUM-QR) to analyze groundwater flow in the heterogeneous porous medium. We evaluate the performance of the proposed methods by increasing the number of random modes and the variance of the log-hydraulic conductivity fields using Gaussian and exponential correlation. In addition, we verify the effectiveness of the presented approaches in one and two dimensions using the method of manufactured solutions. Furthermore, we provide statistical inferences through Monte Carlo simulations. Numerical results demonstrate that RBFPUM-QR performs well in handling higher heterogeneity in both single realizations and Monte Carlo simulations.