An embedded boundary method for the wave equation with discontinuous coefficients

Journal Article · Mon Sep 26 04:00:00 EDT 2005 · SIAM Journal of Scientific Computing, vol. 28, n/a, December 5, 2006, pp. 2054-2074

OSTI ID:899432

Kreiss, H O; Petersson, N A

A second order accurate embedded boundary method for the two-dimensional wave equation with discontinuous wave propagation speed is described. The wave equation is discretized on a Cartesian grid with constant grid size and the interface (across which the wave speed is discontinuous) is allowed to intersect the mesh in an arbitrary fashion. By using ghost points on either side of the interface, previous embedded boundary techniques for the Neumann and Dirichlet problems are generalized to satisfy the jump conditions across the interface to second order accuracy. The resulting discretization of the jump conditions has the desirable property that each ghost point can be updated independently of all other ghost points, resulting in a fully explicit time-integration method. Numerical examples are given where the method is used to study electro-magnetic scattering of a plane wave by a dielectric cylinder. The numerical solutions are evaluated against the analytical solution due to Mie, and point-wise second order accuracy is confirmed.

• Research Organization: Lawrence Livermore National Laboratory (LLNL), Livermore, CA
• Sponsoring Organization: USDOE
• DOE Contract Number: W-7405-ENG-48
• OSTI ID: 899432
• Report Number(s): UCRL-JRNL-215702
• Journal Information: SIAM Journal of Scientific Computing, vol. 28, n/a, December 5, 2006, pp. 2054-2074, Journal Name: SIAM Journal of Scientific Computing, vol. 28, n/a, December 5, 2006, pp. 2054-2074
• Country of Publication: United States
• Language: English

Related Subjects

36 MATERIALS SCIENCE
42 ENGINEERING
58 GEOSCIENCES
99 GENERAL AND MISCELLANEOUS
ACCURACY
ANALYTICAL SOLUTION
DIELECTRIC MATERIALS
DIRICHLET PROBLEM
NUMERICAL SOLUTION
SCATTERING
VELOCITY
WAVE EQUATIONS
WAVE PROPAGATION