An ADER discontinuous Galerkin method on moving meshes for Liouville's equation of geometrical optics

• Published in: Journal of computational physics 2023-09, Vol.488, p.112208, Article 112208
• Main Authors: van Gestel, Robert A.M., Anthonissen, Martijn J.H., ten Thije Boonkkamp, Jan H.M., IJzerman, Wilbert L.
• Format: Article
• Language: English
• Subjects: ADER; Arbitrary-Lagrangian Eulerian; Discontinuous Galerkin; Geometrical optics; Liouville's equation
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2023.112208

Liouville's equation describes light propagation through an optical system. It governs the evolution of an energy distribution on phase space. This energy distribution is discontinuous across optical interfaces. Curved optical interfaces manifest themselves as moving boundaries on phase space. In this paper, an ADER discontinuous Galerkin (DG) method on a moving mesh is applied to solve Liouville's equation. In the ADER approach a temporal Taylor series is computed by replacing temporal derivatives with spatial derivatives using the Cauchy-Kovalewski procedure. The result is a fully discrete explicit scheme of arbitrary high order of accuracy. A moving mesh is not sufficient to be able to solve Liouville's equation numerically for the optical systems considered in this article. To that end, we combine the scheme with a new method we refer to as sub-cell interface method. When dealing with optical interfaces in phase space, non-local boundary conditions arise. These are incorporated in the DG method in an energy-preserving manner. Numerical experiments validate energy-preservation up to machine precision and show the high order of accuracy. Furthermore, the DG method is compared to quasi-Monte Carlo ray tracing for two examples showing that the DG method yields better accuracy in the same amount of computational time. •Explicit Taylor expansion on moving mesh for discontinuous Galerkin.•Non-local boundary conditions on phase space.•High order energy-preserving method for geometrical optics.•Outperforms quasi-Monte Carlo ray tracing in a few examples.