For the computation of nonlinear solutions of Hamilton?Jacobi scalar equations in two space dimensions, we develop high order accurate numerical schemes that can be applied to complicated geometries. Previously, the recently developed essentially nonoscillatory (ENO) technology has been applied in simple domains like squares or rectangles using dimension-by-dimension algorithms. On arbitrary two dimensional closed or multiply connected domains, first order monotone methods were used. In this paper, we propose two different techniques to construct high order accurate methods using the ENO philosophy. Namely, any arbitrary domain is triangulated by finite elements into which two dimensional ENO polynomials are constructed. These polynomials are then differentiated to compute a high order accurate numerical solution. These new techniques are shown to be very useful in the computation of numerical solutions of various applications without significantly increasing CPU running times as compared to dimension-by-dimension algorithms. Furthermore, these methods are stable and no spurious oscillations are detected near singular points or curves.
[1]
E. Rouy,et al.
A viscosity solutions approach to shape-from-shading
,
1992
.
[2]
S. Osher,et al.
High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations
,
1990
.
[3]
S. Osher,et al.
Uniformly high order accurate essentially non-oscillatory schemes, 111
,
1987
.
[4]
P. Lions.
Generalized Solutions of Hamilton-Jacobi Equations
,
1982
.
[5]
S. Osher,et al.
Efficient implementation of essentially non-oscillatory shock-capturing schemes,II
,
1989
.
[6]
S. Osher,et al.
Uniformly High-Order Accurate Nonoscillatory Schemes. I
,
1987
.
[7]
J. Sethian,et al.
Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations
,
1988
.
[8]
S. Osher,et al.
Triangle based adaptive stencils for the solution of hyperbolic conservation laws
,
1992
.