Accurate simulation of seismic waves is of critical importance in a variety of geophysical applications. Based on recent works on staggered discontinuous Galerkin methods, we have developed a new method for the simulations of seismic waves, which has energy conservation and extremely low grid dispersion, so that it naturally provided accurate numerical simulations of wave propagation useful for geophysical applications and was a generalization of classical staggered-grid finite-difference methods. Moreover, it could handle with ease irregular surface topography and discontinuities in the subsurface models. Our new method discretized the velocity and the stress tensor on this staggered grid, with continuity imposed on different parts of the mesh. The symmetry of the stress tensor was enforced by the Lagrange multiplier technique. The resulting method was an explicit scheme, requiring the solutions of a

No abstract uploaded!

In this paper, we generalize the point integral method to solve the general elliptic PDEs with variable coefficients and corresponding eigenvalue problems with Neumann, Robin and Dirichlet boundary conditions on point cloud. The main idea is using integral equations to approximate the original PDEs. The integral equations are easy to discretize on the point cloud. The truncation error of the integral approximation is analyzed. Numerical examples are presented to demonstrate that PIM is an effective method to solve the elliptic PDEs with smooth coefficients on point cloud.

We develop a method searching for quasi-conformal maps from point-cloud surfaces to Euclidean plane. The maps can be got by directly solving Beltrami equation with proper boundary condition or solving optimization problem of the variational formu- lation. The point integral methd (PIM) is used for discretization on point clouds. Numerical experiments suggest that the method converges as the density of points in- creases.

The goal of this paper is to prove the well-posedness for the governing equations which are used for cylindrical cloaking simulation. A new time-domain finite element scheme is developed to solve the governing equations. Numerical results demonstrating the cloaking phenomenon with the cylindrical cloak are presented. We finally extend the analysis and simulation to an elliptical cloak model.