Numerical Modeling of Peridynamic Richards’ Equation with Piecewise Smooth Initial Conditions Using Spectral Methods

Abstract

In this paper, we introduce peridynamic theory and its application to Richards’ equation with a piecewise smooth initial condition. Peridynamic theory is a non-local continuum theory that models the deformation and failure of materials. Richards’ equation describes the unsaturated flow of water through porous media, and it plays an essential role in many applications, such as groundwater management, soil science, and environmental engineering. We develop a peridynamic formulation of Richards’ equation that includes the effect of peridynamic forces and a piecewise smooth initial condition, further introducing a non-standard symmetric influence function to describe such peridynamic interactions, which turns out to provide beneficial effects from a numerical point of view. Moreover, we implement a numerical scheme based on Chebyshev polynomials and symmetric Gauss–Lobatto nodes, providing a powerful spectral method able to capture singularities and critical issues of Richards’ equation with piecewise smooth initial conditions. We also present numerical simulations that illustrate the performance of the proposed approach. In particular, we perform a computational investigation into the spatial order of convergence, showing that, despite the discontinuity in the initial condition, the order of convergence is retained.