On Finite Difference Schemes for the 3-d Wave Equation Using Non-cartesian Grids
نویسندگان
چکیده
In this paper, we investigate finite difference schemes for the 3-D wave equation using 27-point stencils on the cubic lattice, a 13-point stencil on the face-centered cubic (FCC) lattice, and a 9-point stencil on the body-centered cubic (BCC) lattice. The tiling of the wavenumber space for nonCartesian grids is considered in order to analyse numerical dispersion. Schemes are compared for computational efficiency in terms of minimising numerical wave speed error. It is shown that the 13-point scheme on the FCC lattice is more computationally efficient than 27-point schemes on the cubic lattice when less than 8% error in the wave speed is desired.
منابع مشابه
Supra-convergence of Linear Equations on Irregular Cartesian Grids
We prove the second order convergence of a class of finite difference schemes for linear heat equations, and wave equations on irregular grids. Numerical examples and convergence studies are provided to demonstrate our theoretical results.
متن کاملA supra-convergent finite difference scheme for the variable coefficient Poisson equation on non-graded grids
We introduce a method for solving the variable coefficient Poisson equation on non-graded Cartesian grids that yields second order accuracy for the solutions and their gradients. We employ quadtree (in 2D) and octree (in 3D) data structures as an efficient means to represent the Cartesian grid, allowing for constraint-free grid generation. The schemes take advantage of sampling the solution at ...
متن کاملNumerical studies of non-local hyperbolic partial differential equations using collocation methods
The non-local hyperbolic partial differential equations have many applications in sciences and engineering. A collocation finite element approach based on exponential cubic B-spline and quintic B-spline are presented for the numerical solution of the wave equation subject to nonlocal boundary condition. Von Neumann stability analysis is used to analyze the proposed methods. The efficiency, accu...
متن کاملHexagonal vs. Rectilinear Grids for Explicit Finite Difference Schemes for the Two-dimensional Wave Equation
Finite difference schemes for the 2-D wave equation operating on hexagonal grids and the accompanying numerical dispersion properties have received little attention in comparison to schemes operating on rectilinear grids. This paper considers the hexagonal tiling of the wavenumber plane in order to show that the hexagonal grid is a more natural choice to emulate the isotropy of the Laplacian op...
متن کاملThe Method of Difference Potentials for the Helmholtz Equation Using Compact High Order Schemes
The method of difference potentials was originally proposed by Ryaben’kii and can be interpreted as a generalized discrete version of the method of Calderon’s operators in the theory of partial differential equations. It has a number of important advantages; it easily handles curvilinear boundaries, variable coefficients, and non-standard boundary conditions while keeping the complexity at the ...
متن کامل