Popis: |
We describe a numerical implementation of a differential model for the simulation of self-organized criticality (SOC) phenomena arising from recent papers by Barbu (Annu Rev Control 34:52–61, 2010; Math Methods Appl Sci 36:1726–1733, 2013). In that singular nonlinear diffusion problem an initial supercritical state evolves in a finite time towards a given critical solution, progressively from the boundary towards the internal regions. The key elements are the Heaviside function which plays the role of a switch for the dynamics, and the initial boundary contact with the critical state. |