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.

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

Fabio V. Difonzo
;
Francesco Di Lena
2023-01-01

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.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11586/431940
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus ND
  • ???jsp.display-item.citation.isi??? ND
social impact