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

Fabio V. Difonzo,Francesco Di Lena
DOI: https://doi.org/10.3390/sym15050960
2023-04-24
Symmetry
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.
multidisciplinary sciences
What problem does this paper attempt to address?