Numerical Approximations for a Nonlocal Evolution Equation

In this paper we study numerical approximations of continuous solutions to the nonlocal $p$-Laplacian type diffusion equation, $u_t (t,x) = \int_{\Omega} J(x-y)|u(t,y) - u(t,x)|^{p-2}(u(t,y) - u(t,x)) \, dy$. First, we find that a semidiscretization in space of this problem gives rise to an ODE system whose solutions converge uniformly to the continuous one as the mesh size goes to zero. Moreover, the semidiscrete approximation shares some properties of the continuous problem: it preserves the total mass and the solution converges to the mean value of the initial condition as $t$ goes to infinity. Next, we also discretize the time variable and present a totally discrete method which also enjoys the above mentioned properties. In addition, we investigate the limit as $p$ goes to infinity in these approximations and obtain a discrete model for the evolution of a sandpile. Finally, we present some numerical experiments that illustrate our results.

[1]  Mikhail Feldman Growth of a sandpile around an obstacle , 1999 .

[2]  P. Fife,et al.  Spatial effects in discrete generation population models , 2005, Journal of mathematical biology.

[3]  J. D. Rossi,et al.  Boundary fluxes for non-local diffusion , 2006, math/0607057.

[4]  Noemi Wolanski,et al.  How to Approximate the Heat Equation with Neumann Boundary Conditions by Nonlocal Diffusion Problems , 2006, math/0607058.

[5]  I. Ekeland,et al.  Convex analysis and variational problems , 1976 .

[6]  Xuefeng Wang,et al.  Metastability and Stability of Patterns in a Convolution Model for Phase Transitions , 2002 .

[7]  Fuensanta Andreu,et al.  The limit as p → ∞ in a nonlocal p-Laplacian evolution equation: a nonlocal approximation of a model for sandpiles , 2009 .

[8]  Noemi Wolanski,et al.  Boundary fluxes for nonlocal diffusion , 2007 .

[9]  Jean-Michel Morel,et al.  Neighborhood filters and PDE’s , 2006, Numerische Mathematik.

[10]  Julio D. Rossi,et al.  A Nonlocal Diffusion Equation whose Solutions Develop a Free Boundary , 2005 .

[11]  P. Bates,et al.  Traveling Waves in a Convolution Model for Phase Transitions , 1997 .

[12]  Paul C. Fife,et al.  Some Nonclassical Trends in Parabolic and Parabolic-like Evolutions , 2003 .

[13]  A Stochastic Model for Growing Sandpiles and its Continuum Limit , 1998 .

[14]  Fuensanta Andreu,et al.  The Neumann problem for nonlocal nonlinear diffusion equations , 2008 .

[15]  Emmanuel Chasseigne,et al.  Asymptotic behavior for nonlocal diffusion equations , 2006 .

[16]  Xinfu Chen,et al.  Existence, uniqueness, and asymptotic stability of traveling waves in nonlocal evolution equations , 1997, Advances in Differential Equations.

[17]  Peter W. Bates,et al.  An Integrodifferential Model for Phase Transitions: Stationary Solutions in Higher Space Dimensions , 1999 .

[18]  Linghai Zhang,et al.  Existence, uniqueness and exponential stability of traveling wave solutions of some integral differential equations arising from neuronal networks , 2004 .

[19]  Mikhail Feldman,et al.  Fast/slow diffusion and collapsing sandpiles , 1997 .

[20]  Mario Milman,et al.  Monge Ampère equation : applications to geometry and optimization : NSF-CBMS Conference on the Monge Ampère Equation : Applications to Geometry and Optimization, July 9-13, 1997, Florida Atlantic University , 1999 .

[21]  J. M. Mazón,et al.  A nonlocal p-Laplacian evolution equation with Neumann boundary conditions , 2008 .

[22]  H. Brezis Opérateurs maximaux monotones et semi-groupes de contractions dans les espaces de Hilbert , 1973 .

[23]  H. Brezis,et al.  Convergence and approximation of semigroups of nonlinear operators in Banach spaces , 1972 .

[24]  U. Mosco Convergence of convex sets and of solutions of variational inequalities , 1969 .

[25]  Stanley Osher,et al.  Deblurring and Denoising of Images by Nonlocal Functionals , 2005, Multiscale Model. Simul..

[26]  J. L. Webb OPERATEURS MAXIMAUX MONOTONES ET SEMI‐GROUPES DE CONTRACTIONS DANS LES ESPACES DE HILBERT , 1974 .

[27]  Luis Silvestre,et al.  Holder estimates for solutions of integra-differential equations like the fractional laplace , 2006 .

[28]  Paul C. Fife,et al.  A convolution model for interfacial motion: the generation and propagation of internal layers in higher space dimensions , 1998, Advances in Differential Equations.