In this paper, we introduce a new numerical procedure for simulations in geometrical optics that, based on the recent development of Eulerian phase-space formulations of the model, can deliver very accurate, uniformly resolved solutions which can be made to converge with arbitrarily high orders in general geometrical configurations. Following previous treatments, the scheme is based on the evolution of a wavefront in phase-space, defined as the intersection of level sets satisfying the relevant Liouville equation. In contrast with previous work, however, our numerical approximation is specifically designed: (i) to take full advantage of the smoothness of solutions; (ii) to facilitate the treatment of scattering obstacles, all while retaining high-order convergence characteristics. Indeed, to incorporate the full regularity of solutions that results from the unfolding of singularities, our method is based on their spectral representation; to enable a simple high-order treatment of scattering boundaries, on the other hand, we resort to a discontinuous Galerkin finite element method for the solution of the resulting system of equations. The procedure is complemented with the use of a recently derived strong stability preserving Runge-Kutta (SSP-RK) scheme for the time integration that, as we demonstrate, allows for overall approximations that are rapidly convergent.
Bibliographical noteFunding Information:
B.C. gratefully acknowledges support from NSF through Grant No. DMS-0107609. J.Q. gratefully acknowledges support from ONR through Grant No. N00014-02-1-0720. F.R. gratefully acknowledges support from NSF through Grant No. DMS-0311763, from AFOSR through Contract No. F49620-02-1-0052 and from the Army High Performance Computing Research Center (AHPCRC) under Army Research Laboratory cooperative agreement number DAAD19-01-2-0014.
Disclaimer. Effort sponsored by the Air Force Office of Scientific Research, Air Force Materials Command, USAF, under Grant No. F49620-02-1-0052, and by AHPCRC under the auspices of the Department of the Army, Army Research Laboratory cooperative agreement number DAAD19-01-2-0014. The US Government is authorized to reproduce and distribute reprints for governmental purposes notwithstanding any copyright notation thereon. The views and conclusions contained herein are those of the author and should not be interpreted as necessarily representing the official policies or endorsements, either expressed or implied, of the Air Force Office of Scientific Research, the Army Research Laboratory or the US Government.
- Discontinuous Galerkin
- Eikonal equation
- Geometrical optics
- Liouville equation
- Spectral methods
- Wave equation