A level set approach for computing discontinuous solutions of Hamilton-Jacobi equations

We introduce two types of finite difference methods to compute the L-solution and the proper viscosity solution recently proposed by the second author for semi-discontinuous solutions to a class of Hamilton-Jacobi equations. By regarding the graph of the solution as the zero level curve of a continuous function in one dimension higher, we can treat the corresponding level set equation using the viscosity theory introduced by Crandall and Lions. However, we need to pay special attention both analytically and numerically to prevent the zero level curve from overturning so that it can be interpreted as the graph of a function. We demonstrate our Lax-Friedrichs type numerical methods for computing the L-solution using its original level set formulation. In addition, we couple our numerical methods with a singular diffusive term which is essential to computing solutions to a more general class of HJ equations that includes conservation laws. With this singular viscosity, our numerical methods do not require the divergence structure of equations and do apply to more general equations developing shocks other than conservation laws. These numerical methods are generalized to higher order accuracy using weighted ENO local Lax-Friedrichs methods as developed recently by Jiang and Peng. We verify that our numerical solutions approximate the proper viscosity solutions obtained by the second author in a recent Hokkaido University preprint. Finally, since the solution of scalar conservation law equations can be constructed using existing numerical techniques, we use it to verify that our numerical solution approximates the entropy solution.

[1]  Olivier Guès.  Ondes oscillantes simples quasilinéaires , 1989 .

[2]  Thea Iberall,et al.  Dextrous robot hands , 1990 .

[3]  S. Osher,et al.  Regular Article: A PDE-Based Fast Local Level Set Method , 1999 .

[4]  E. Barron,et al.  Semicontinuous Viscosity Solutions For Hamilton–Jacobi Equations With Convex Hamiltonians , 1990 .

[5]  J. Łuczka,et al.  Application of statistical mechanics to stochastic transport , 1999 .

[6]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[7]  Danping Peng,et al.  Weighted ENO Schemes for Hamilton-Jacobi Equations , 1999, SIAM J. Sci. Comput..

[8]  Lawrence C. Evans.  A Geometric Interpretation of the Heat Equation with Multivalued Initial Data , 1996 .

[9]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[10]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[11]  M. Bardi,et al.  Optimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations , 1997 .

[12]  H. Ishii.  Perron’s method for Hamilton-Jacobi equations , 1987 .

[13]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[14]  Yoshikazu Giga,et al.  Very singular diffusion equations , 1999 .

[15]  S. Osher.  A level set formulation for the solution of the Dirichlet problem for Hamilton-Jacobi equations , 1993 .

[16]  H. Ishii,et al.  Comparison principle and convexity preserving properties for singular degenerate parabolic equations on unbounded domains , 1991 .

[17]  P. Lions,et al.  Two approximations of solutions of Hamilton-Jacobi equations , 1984 .

[18]  James A. Sethian,et al.  Fast-marching level-set methods for three-dimensional photolithography development , 1996, Advanced Lithography.

[19]  G. Barles.  Solutions de viscosité des équations de Hamilton-Jacobi , 1994 .

[20]  Phillip Colella,et al.  Two new methods for simulating photolithography development in 3D , 1996, Advanced Lithography.

[21]  S. Osher,et al.  A PDE-Based Fast Local Level Set Method 1 , 1998 .

[22]  H. Ishii,et al.  Existence and uniqueness of solutions of Hamilton-Jacobi equations , 1986 .

[23]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[24]  Yoshikazu Giga,et al.  Equations with Singular Diffusivity , 1998 .

[25]  J. Tsitsiklis,et al.  Efficient algorithms for globally optimal trajectories , 1994, Proceedings of 1994 33rd IEEE Conference on Decision and Control.

[26]  Rahul Pandharipande,et al.  EQUIVARIANT CHOW RINGS OF O (K), SO (2K + 1), AND SO (4) , 1998 .

[27]  Yoshikazu Giga.  Shocks and very strong vertical diffusion , 2000 .

[28]  Yoshikazu Giga,et al.  Crystalline and level set flow - convergence of a crystalline algorithm for a general anisotropic cu , 2000 .

[29]  Chi-Tien Lin,et al.  High-Resolution Nonoscillatory Central Schemes for Hamilton-Jacobi Equations , 1999, SIAM J. Sci. Comput..

[30]  Y. Tsai.  Rapid and accurate computation of the distance function using grids , 2002 .

[31]  Constantin Carathéodory,et al.  Calculus of variations and partial differential equations of the first order , 1965 .

[32]  H. Ishii.  Hamilton-Jacobi Equations with Discontinuous Hamiltonians on Arbitrary Open Sets , 1985 .

[33]  P. Souganidis.  Approximation schemes for viscosity solutions of Hamilton-Jacobi equations , 1985 .

[34]  S. Osher,et al.  High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations , 1990 .

[35]  Yoshikazu Giga,et al.  A LEVEL SET APPROACH TO SEMICONTINUOUS VISCOSITY SOLUTIONS FOR CAUCHY PROBLEMS , 1999 .

[36]  G. Barles.  Discontinuous viscosity solutions of first-order Hamilton-Jacobi equations: a guided visit , 1993 .

[37]  P. Lax.  Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves , 1987 .