Title

The Radial Basis Functions Method for Improved Numerical Approximations of Geological Processes in Heterogeneous Systems

Document Type

Article

Publication Date

2019

Publication Title

Mathematical Geosciences

ISSN

1874-8953

DOI

https://dx.doi.org/10.1007/s11004-019-09820-w

Abstract

A robust, high order modeling approach is introduced, based on the finite difference-based radial basis functions method, for solving the groundwater flow equation in the presence of an active well, in the case of a confined aquifer in a complex geological environment. The two important novelties of this work are the analytical handling of the wells’ singularities and the ability to do this accurately and efficiently in a heterogeneous medium. It is argued that the most commonly used methods for this type of problem have severe weaknesses in both the treatment of the singularities associated with the well, and in representing heterogeneities which commonly occur in geological processes. The method presented here is first applied to the groundwater flow problem in a homogeneous medium for which the analytical solution is known, to show its high order algebraic convergence. The method is then compared against the United States geological survey’s MODFLOW software on a quasi-realistic benchmark test case in a heterogeneous medium. It is shown that much fewer nodes are needed by the proposed method to yield a similar accuracy.

Comments

This is a post-print version of an article published in Mathematical Geosciences.

This document is currently not available here.

Share

COinS