An Efficient Data Structure and Accurate Scheme to Solve Front Propagation Problems

  • Authors:
  • O. Bokanowski;E. Cristiani;H. Zidani

  • Affiliations:
  • Laboratoire Jacques-Louis Lions (UMR 7598), Université Paris 6, Paris, France 75005 and UFR Mathématiques, Université Paris 7, Paris, France 75013;CEMSAC, Università di Salerno, Salerno, Italy and IAC-CNR, Rome, Italy;Projet Commands, Ensta - Inria Saclay, Paris Cx 15, France 75379

  • Venue:
  • Journal of Scientific Computing
  • Year:
  • 2010

Quantified Score

Hi-index 0.00

Visualization

Abstract

In this paper, we are interested in some front propagation problems coming from control problems in d-dimensional spaces, with d驴2. As opposed to the usual level set method, we localize the front as a discontinuity of a characteristic function. The evolution of the front is computed by solving an Hamilton-Jacobi-Bellman equation with discontinuous data, discretized by means of the antidissipative Ultra Bee scheme.We develop an efficient dynamic storage technique suitable for handling front evolutions in large dimension. Then we propose a fast algorithm, showing its relevance on several challenging tests in dimension d=2,3,4. We also compare our method with the techniques usually used in level set methods. Our approach leads to a computational cost as well as a memory allocation scaling as O(N nb ) in most situations, where N nb is the number of grid nodes around the front. Moreover, we show on several examples the accuracy of our approach when compared with level set methods.