A Meshless Solution Technique for the Solution of 3D Unsaturated Zone Problems, Based on Local Hermitian Interpolation with Radial Basis Functions

Abstract

Recent developments in meshless numerical methods have led to algorithms that can be used to solve arbitrarily large problems without the support of a connected mesh, and without the computational cost and numerical ill-conditioning issues usually associated with such solution techniques. This work applies the Local Hermitian Interpolation (LHI) method, based on local interpolation with Radial Basis Functions (RBFs), to the solution of 3D unsaturated porous media problems. The proposed implementation is capable of handling real soil properties, provided either as an analytical function or as a series of pointwise measurements. The technique is implemented with implicit and explicit timestepping, and is validated against two transient Richards’ equation models, of which one has a known analytical solution. In addition, a real-world infiltration problem based on a saturated–unsaturated formulation is modelled, using a realistic variation of soil properties with water-pressure.