High-order accurate discontinuous finite element solution of the 2D Euler equations
Journal of Computational Physics
The Runge-Kutta discontinuous Galerkin method for conservation laws V multidimensional systems
Journal of Computational Physics
Geometric optics in a phase-space-based level set and Eulerian framework
Journal of Computational Physics
High-order RKDG Methods for Computational Electromagnetics
Journal of Scientific Computing
Journal of Computational Physics
P2Q2Iso2D = 2D isoparametric FEM in Matlab
Journal of Computational and Applied Mathematics
A Level Set Framework for Capturing Multi-Valued Solutions of Nonlinear First-Order Equations
Journal of Scientific Computing
Hi-index | 31.45 |
We present an adaptive spectral/discontinuous Galerkin (DG) method on curved elements to simulate high-frequency wavefronts within a reduced phase-space formulation of geometrical optics. Following recent work, the approach is based on the use of level sets defined by functions satisfying the Liouville equations in reduced phase-space and, in particular, it relies on the smoothness of these functions to represent them by rapidly convergent spectral expansions in the phase variables. The resulting (hyperbolic) system of equations for the coefficients in these expansions are then amenable to a high-order accurate treatment via DG approximations. In the present work, we significantly expand on the applicability and efficiency of the approach by incorporating mechanisms that allow for its use in scattering simulations and for a reduced overall computational cost. With regards to the former we demonstrate that the incorporation of curved elements is necessary to attain any kind of accuracy in calculations that involve scattering off non-flat interfaces. With regards to efficiency, on the other hand, we also show that the level-set formulation allows for a space p-adaptive scheme that under-resolves the level-set functions away from the wavefront without incurring in a loss of accuracy in the approximation of its location. As we show, these improvements enable simulations that are beyond the capabilities of previous implementations of these numerical procedures.