An adaptive spectral/DG method for a reduced phase-space based level set approach to geometrical optics on curved elements

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.

[1]  Håvar Gjøystdal,et al.  Estimation of multivalued arrivals in 3D models using wavefront construction , 1993 .

[2]  William E. Schiesser,et al.  Linear and nonlinear waves , 2009, Scholarpedia.

[3]  Gilles Lambaré,et al.  Can We Image Quantitatively Complex Models With Rays , 1998 .

[4]  Carsten Carstensen,et al.  P2Q2Iso2D = 2D isoparametric FEM in Matlab , 2006 .

[5]  Markus Voelter,et al.  State of the Art , 1997, Pediatric Research.

[6]  Fernando Reitich,et al.  An accurate spectral/discontinuous finite-element formulation of a phase-space-based level set approach to geometrical optics , 2005 .

[7]  S. Rebay,et al.  High-Order Accurate Discontinuous Finite Element Solution of the 2D Euler Equations , 1997 .

[8]  S. Osher,et al.  COMPUTATIONAL HIGH-FREQUENCY WAVE PROPAGATION USING THE LEVEL SET METHOD, WITH APPLICATIONS TO THE SEMI-CLASSICAL LIMIT OF SCHRÖDINGER EQUATIONS∗ , 2003 .

[9]  P. Lions,et al.  User’s guide to viscosity solutions of second order partial differential equations , 1992, math/9207212.

[10]  Hailiang Liu,et al.  A Level Set Framework for Capturing Multi-Valued Solutions of Nonlinear First-Order Equations , 2006, J. Sci. Comput..

[11]  B. Engquist,et al.  Computational high frequency wave propagation , 2003, Acta Numerica.

[12]  Håvar Gjøystdal,et al.  Part II: Tracing and interpolation1 , 1996 .

[13]  Stanley Osher,et al.  A level set method for the computation of multi-valued solutions to quasi-linear hyperbolic PDE's and Hamilton-Jacobi equations , 2003 .

[14]  S. Osher,et al.  A level set-based Eulerian approach for anisotropic wave propagation , 2003 .

[15]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[16]  S. Osher,et al.  A LEVEL SET METHOD FOR THE COMPUTATION OF MULTIVALUED SOLUTIONS TO QUASI-LINEAR HYPERBOLIC PDES AND HAMILTON-JACOBI EQUATIONS , 2003 .

[17]  G. Lambaré,et al.  Can we quantitatively image complex structures with rays , 2000 .

[18]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[19]  Håvar Gjøystdal,et al.  Traveltime and amplitude estimation using wavefront construction , 1993 .

[20]  Stanley Osher,et al.  Geometric Optics in a Phase-Space-Based Level Set and Eulerian Framework , 2002 .

[21]  P. Lions,et al.  Two approximations of solutions of Hamilton-Jacobi equations , 1984 .

[22]  William H. Press,et al.  Numerical recipes in C , 2002 .

[23]  Håvar Gjøystdal,et al.  Estimation of multivalued arrivals in 3D models using wavefront construction—Part I , 1996 .

[24]  George Em Karniadakis,et al.  A NEW TRIANGULAR AND TETRAHEDRAL BASIS FOR HIGH-ORDER (HP) FINITE ELEMENT METHODS , 1995 .

[25]  F. Reitich,et al.  State-of-the-Art, Trends, and Directions in Computational Electromagnetics , 2004 .

[26]  K. Bathe.  Finite Element Procedures , 1995 .

[27]  Jeff Hughes,et al.  Xpatch 4: the next generation in high frequency electromagnetic modeling and simulation software , 2000, Record of the IEEE 2000 International Radar Conference [Cat. No. 00CH37037].

[28]  S. Wandzurat,et al.  Symmetric quadrature rules on a triangle , 2003 .

[29]  Fernando Reitich,et al.  High-order RKDG Methods for Computational Electromagnetics , 2005, J. Sci. Comput..