In this Note, we discuss the numerical solution of a system of Eikonal equations with Dirichlet boundary conditions. Since the problem under consideration has infinitely many solutions, we look for those solutions which are nonnegative and of maximal (or nearly maximal)  -norm. The computational methodology combines penalty, biharmonic regularization, operator splitting, and finite element approximations. Its practical implementation requires essentially the solution of cubic equations in one variable and of discrete linear elliptic problems of the Poisson and Helmholtz type. As expected, when the spatial domain is a square whose sides are parallel to the coordinate axes, and when the Dirichlet data vanishes at the boundary, the computed solutions show a fractal behavior near the boundary, and particularly, close to the corners. To cite this article: B. Dacorogna et al., C. R. Acad. Sci. Paris, Ser. I 336 (2003).

Dans cette Note, on étudie la résolution numérique d'un système d'équations eiconales avec conditions aux limites du type Dirichlet. Dans la mesure, où le problème considéré a une infinité de solutions on recherche celles qui sont non-négatives et de norme   maximale (on presque maximale). La méthodologie numérique combine pénalité, régularisation biharmonique, décomposition d'opérateurs, et approximations par éléments finis. Son implémentation demande essentiellement la résolution d'équations à une variable du troisième degré et de problèmes linéaires elliptiques discrets pour le Laplacien et l'opérateur d'Helmholtz. Comme prévu, quand le domaine spatial est un carré de côtés parallèles aux axes de coordonnées les solutions calculées montrent un comportement fractal au voisinage de la frontière et plus particulièrement des coins. Pour citer cet article : B. Dacorogna et al., C. R. Acad. Sci. Paris, Ser. I 336 (2003).