A numerical study for growing sandpiles on flat tables with walls

We continue our study on the approximation of a system of partial differential equations recently proposed by Hadeler and Kuttler to model the dynamics of growing sandpiles on a flat bounded table. The novelty here is the introduction of (infinite) walls on the boundary of the domain and the corresponding modification of boundary conditions for the standing and for the rolling layers. An explicit finite difference scheme is introduced and new boundary conditions are analyzed. We show experiments in ID and 2D which characterize the steady-state solutions. © 2006 International Federation for Information Processing.