ESAIM: M2AN 42 (2008) 535-563
DOI: 10.1051/m2an:2008016

Fully adaptive multiresolution schemes for strongly degenerate parabolic equations in one space dimension

Raimund Bürger¹, Ricardo Ruiz¹, Kai Schneider² and Mauricio Sepúlveda<sup>1

1</sup>  Departamento de Ingeniería Matemática, Universidad de Concepción, Casilla 160-C, Concepción, Chile. rburger@ing-mat.udec.cl; rruiz@ing-mat.udec.cl; mauricio@ing-mat.udec.cl
²  Centre de Mathématiques et d'Informatique, Université de Provence, 39 rue Joliot-Curie, 13453 Marseille Cedex 13, France. kschneid@cmi.univ-mrs.fr

Received January 31, 2007. Revised September 17, 2007. Published online May 27, 2008.

Abstract
We present a fully adaptive multiresolution scheme for spatially one-dimensional quasilinear strongly degenerate parabolic equations with zero-flux and periodic boundary conditions. The numerical scheme is based on a finite volume discretization using the Engquist-Osher numerical flux and explicit time stepping. An adaptive multiresolution scheme based on cell averages is then used to speed up the CPU time and the memory requirements of the underlying finite volume scheme, whose first-order version is known to converge to an entropy solution of the problem. A particular feature of the method is the storage of the multiresolution representation of the solution in a graded tree, whose leaves are the non-uniform finite volumes on which the numerical divergence is eventually evaluated. Moreover using the L¹ contraction of the discrete time evolution operator we derive the optimal choice of the threshold in the adaptive multiresolution method. Numerical examples illustrate the computational efficiency together with the convergence properties.

Mathematics Subject Classification. 35L65, 35R05, 65M06, 76T20.

Key words: Degenerate parabolic equation, adaptive multiresolution scheme, monotone scheme, upwind difference scheme, boundary conditions, entropy solution.


© EDP Sciences, SMAI 2008

What is OpenURL?