Properties-preserving high order numerical methods for a kinetic eikonal equation

For the BGK (Bhatnagar-Gross-Krook) equation in the large scale hyperbolic limit, the density of particles can be transformed as the Hopf-Cole transformation, where the phase function converges uniformly to the viscosity solution of an effective Hamilton-Jacobi equation, referred to as the kinetic eikonal equation. In this work, we present efficient high order finite difference methods for numerically solving the kinetic eikonal equation. The methods are based on monotone schemes such as the Godunov scheme. High order weighted essentially non-oscillatory techniques and Runge-Kutta procedures are used to obtain high order accuracy in both space and time. The effective Hamiltonian is determined implicitly by a nonlinear equation given as integrals with respect to the velocity variable. Newton's method is applied to solve the nonlinear equation, where integrals with respect to the velocity variable are evaluated either by a Gauss quadrature formula or as expansions with respect to moments of the Maxwellian. The methods are designed such that several key properties such as the positivity of the viscosity solution and the positivity of the effective Hamiltonian are preserved. Numerical experiments are presented to demonstrate the effectiveness of the methods.