Superposition of Multi-Valued Solutions in High Frequency Wave Dynamics

The weakly coupled WKB system captures high frequency wave dynamics in many applications. For such a system a level set method framework has been recently developed to compute multi-valued solutions to the Hamilton-Jacobi equation and evaluate position density accordingly. In this paper we propose two approaches for computing multi-valued quantities related to density, momentum as well as energy. Within this level set framework we show that physical observables evaluated in Jin et al. (J. Comput. Phys. 210(2):497–518, [2005]; J. Comput. Phys. 205(1):222–241, [2005]) are simply the superposition of their multi-valued correspondents. A series of numerical tests is performed to compute multi-valued quantities and validate the established superposition properties.

[1]  John D. Towers.  Two methods for discretizing a delta function supported on a level set , 2007, J. Comput. Phys..

[2]  S. Osher,et al.  Regular Article: A PDE-Based Fast Local Level Set Method , 1999 .

[3]  S. Osher,et al.  Multi-Valued Solution and Level Set Methods in Computational High Frequency Wave Propagation , 2006 .

[4]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[5]  Jacques-Louis Lions,et al.  Nonlinear partial differential equations and their applications , 1998 .

[6]  Li-Tien Cheng.  An Efficient Level Set Method for Constructing Wavefronts in Three Space Dimensions , 2006 .

[7]  S. Osher,et al.  A PDE-Based Fast Local Level Set Method 1 , 1998 .

[8]  Huijiang Zhao.  FIRST ORDER QUASILINEAR EQUATIONS IN SEVERAL INDEPENDENT VARIABLES WITH SINGULAR INITIAL DATA Lp(p , 1996 .

[9]  Hailiang Liu,et al.  Computing multi-valued velocity and electric fields for 1D Euler--Poisson equations , 2007 .

[10]  Wang Hai-bing,et al.  High-order essentially non-oscillatory schemes for Hamilton-Jacobi equations , 2006 .

[11]  Stanley Osher,et al.  Computing multi-valued physical observables for the high frequency limit of symmetric hyperbolic systems , 2005 .

[12]  Y. Brenier,et al.  A kinetic formulation for multi-branch entropy solutions of scalar conservation laws , 1998 .

[13]  Ami Harten,et al.  Preliminary results on the extension of eno schemes to two-dimensional problems , 1987 .

[14]  A. Harten.  ENO schemes with subcell resolution , 1989 .

[15]  R. Newcomb.  VISCOSITY SOLUTIONS OF HAMILTON-JACOBI EQUATIONS , 2010 .

[16]  L. Gosse.  Using K-Branch Entropy Solutions for Multivalued Geometric Optics Computations , 2002 .

[17]  B. Engquist,et al.  Multi-phase computations in geometrical optics , 1996 .

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

[19]  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 .

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

[21]  Pierre Charrier,et al.  Nonlinear Hyperbolic Problems , 1987 .

[22]  Hailiang Liu,et al.  A field-space-based level set method for computing multi-valued solutions to 1D Euler-Poisson equations , 2007, J. Comput. Phys..

[23]  B. Engquist,et al.  Discretization of Dirac delta functions in level set methods , 2005 .

[24]  R. LeVeque,et al.  Analysis of a one-dimensional model for the immersed boundary method , 1992 .

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

[26]  L. Chambers.  Linear and Nonlinear Waves , 2000, The Mathematical Gazette.

[27]  Chohong Min.  Local level set method in high dimension and codimension , 2004 .

[28]  Li-Tien Cheng,et al.  Efficient level set methods for constructing wavefronts in three spatial dimensions , 2007, J. Comput. Phys..

[29]  赵会江.  FIRST ORDER QUASILINEAR EQUATIONS IN SEVERAL INDEPENDENT VARIABLES WITH SINGULAR INITIAL DATA L~p（P＜∞） , 1996 .

[30]  Timothy J. Barth,et al.  High-order methods for computational physics , 1999 .

[31]  Laurent Gosse,et al.  TWO MOMENT SYSTEMS FOR COMPUTING MULTIPHASE SEMICLASSICAL LIMITS OF THE SCHRÖDINGER EQUATION , 2003 .

[32]  Shi Jin,et al.  Multi-phase computations of the semiclassical limit of the Schrödinger equation and related problems: Whitham vs Wigner , 2003 .

[33]  Chi-Wang Shu,et al.  High Order ENO and WENO Schemes for Computational Fluid Dynamics , 1999 .

[34]  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 .

[35]  Stanley Osher,et al.  Computing multivalued physical observables for the semiclassical limit of the Schrödinger equation , 2005 .

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

[37]  S. Osher,et al.  High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations , 1990 .

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