High order fast sweeping methods for static Hamilton-Jacobi equations

We construct high order fast sweeping numerical methods for computing viscosity solutions of static Hamilton-Jacobi equations on rectangular grids. These methods combine high order weighted essentially non-oscillatory (WENO) approximations to derivatives, monotone numerical Hamiltonians and Gauss-Seidel iterations with alternating-direction sweepings. Based on well-developed first order sweeping methods, we design a novel approach to incorporate high order approximations to derivatives into numerical Hamiltonians such that the resulting numerical schemes are formally high order accurate and inherit the fast convergence from the alternating sweeping strategy. Extensive numerical examples verify efficiency, convergence and high order accuracy of the new methods.