Convergence Analysis of a Spectral Numerical Method for a Peridynamic Formulation of Richards' Equation

We study the implementation of a Chebyshev spectral method with forward Euler integrator to investigate a peridynamic nonlocal formulation of Richards' equation. We prove the convergence of the fully-discretization of the model showing the existence and uniqueness of a solution to the weak formulation of the method by using the compactness properties of the approximated solution and exploiting the stability of the numerical scheme. We further support our results through numerical simulations, using initial conditions with different order of smoothness, showing reliability and robustness of the theoretical findings presented in the paper.