Local Discontinuous Galerkin Methods for the Khokhlov–Zabolotskaya–Kuznetzov Equation

Khokhlov–Zabolotskaya–Kuznetzov (KZK) equation is a model that describes the propagation of the ultrasound beams in the thermoviscous fluid. It contains a nonlocal diffraction term, an absorption term and a nonlinear term. Accurate numerical methods to simulate the KZK equation are important to its broad applications in medical ultrasound simulations. In this paper, we propose a local discontinuous Galerkin method to solve the KZK equation. We prove the L² stability of our scheme and conduct a series of numerical experiments including the focused circular short tone burst excitation and the propagation of unfocused sound beams, which show that our scheme leads to accurate solutions and performs better than the benchmark solutions in the literature.