Solving wave propagation problems within heterogeneous media has been of great interest and has a wide range of applications in physics and engineering. The design of numerical methods for such general wave propagation problems is challenging because the energy conserving property has to be incorporated in the numerical algorithms in order to minimize the phase or shape errors after long time integration. In this paper, we focus on multi-dimensional wave problems and consider linear second-order wave equation in heterogeneous media. We develop and analyze an LDG method, in which numerical fluxes are carefully designed to maintain the energy preserving property and accuracy. Compatible high order energy conserving time integrators are also proposed. The optimal error estimates and the energy conserving property are proved for the semi-discrete methods. Our numerical experiments demonstrate optimal rates of convergence, and show that the errors of the numerical solutions do not grow significantly in time due to the energy conserving property.

Peskin's Immersed Boundary (IB) method is one of the most popular numerical methods for many years and has been applied to problems in mathematical biology, fluid mechanics, material sciences, and many other areas. Peskiness IB method is associated with discrete delta functions. It is believed that the IB method is first order accurate in the $L^{\infty}$ norm. But almost no rigorous proof could be found in the literature until recently [Mori, Comm. Pure. Appl. Math: 61:2008] in which the author showed that the velocity is indeed first order accurate for the Stokes equations with a periodic boundary condition. In this paper, we show first order convergence with a $\log h$ factor of the IB method for elliptic interface problems with Dirichlet boundary conditions.

In this paper, we analyze the Lax-Wendroff discontinuous Galerkin (LWDG) method for solving linear conservation laws. The method was originally proposed by Guo et al., where they applied local discontinuous Galerkin (LDG) techniques to approximate high order spatial derivatives in the Lax-Wendroff time discretization. We show that, under the standard CFL condition $\tau\leq \lambda h$ (where $\tau$ and $h$ are the time step and the maximum element length respectively and $\lambda>0$ is a constant) and uniform or non-increasing time steps, the second order schemes with piecewise linear elements and the third order schemes with arbitrary piecewise polynomial elements are stable in the $L^2$ norm. The specific type of stability may differ with different choices of numerical fluxes. Our stability analysis includes multidimensional problems with divergence-free coefficients. Besides solving the equation itself, the LWDG method also gives approximations to its time derivative simultaneously. We obtain optimal error estimates for both the solution $u$ and its first order time derivative $u_t$ in one dimension, and numerical examples are given to validate our analysis.

No abstract uploaded!

No abstract uploaded!