Optimized High-Order Derivative and Dissipation Operators Satisfying Summation by Parts, and Applications in Three-dimensional Multi-block Evolutions

We construct optimized high-order finite differencing operators which satisfy summation by parts. Since these operators are not uniquely defined, we consider several optimization criteria: minimizing the bandwidth, the truncation error on the boundary points, the spectral radius, or a combination of these. We examine in detail a set of operators that are up to tenth order accurate in the interior, and we surprisingly find that a combination of these optimizations can improve the operators’ spectral radius and accuracy by orders of magnitude in certain cases. We also construct high-order dissipation operators that are compatible with these new finite difference operators and which are semi-definite with respect to the appropriate summation by parts scalar product. We test the stability and accuracy of these new difference and dissipation operators by evolving a three-dimensional scalar wave equation on a spherical domain consisting of seven blocks, each discretized with a structured grid, and connected through penalty boundary conditions. In particular, we find that the constructed dissipation operators are effective in suppressing instabilities that are sometimes otherwise present in the restricted full norm case.

[1]  Bertil Gustafsson,et al.  On the implementation of boundary conditions for the method of lines , 1998 .

[2]  Paolo Secchi,et al.  Well-posedness of characteristic symmetric hyperbolic systems , 1996 .

[3]  Paolo Secchi,et al.  The initial-boundary value problem for linear symmetric hyperbolic systems with characteristic boundary of constant multiplicity , 1996, Differential and Integral Equations.

[4]  Helmut Friedrich.  Conformal Einstein evolution , 2002 .

[5]  Magnus Svärd,et al.  On Coordinate Transformations for Summation-by-Parts Operators , 2004, J. Sci. Comput..

[6]  D. Gottlieb,et al.  Spectral methods for hyperbolic problems , 2001 .

[7]  Magnus Svärd,et al.  Steady-State Computations Using Summation-by-Parts Operators , 2005, J. Sci. Comput..

[8]  Pelle Olsson.  Supplement to Summation by Parts, Projections, and Stability. I , 1995 .

[9]  P. Olsson.  Summation by parts, projections, and stability. II , 1995 .

[10]  B. Schutz,et al.  Gravitational wave astronomy , 1999, gr-qc/9911034.

[11]  Magnus Svärd,et al.  Stable and Accurate Artificial Dissipation , 2004, J. Sci. Comput..

[12]  D. Gottlieb,et al.  A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy , 1999 .

[13]  Ken Mattsson,et al.  Boundary Procedures for Summation-by-Parts Operators , 2003, J. Sci. Comput..

[14]  John Shalf,et al.  The Cactus Framework and Toolkit: Design and Applications , 2002, VECPAR.

[15]  Manuel Tiglio,et al.  Cauchy-perturbative matching reexamined: Tests in spherical symmetry , 2005, gr-qc/0511163.

[16]  Anthony Skjellum,et al.  A High-Performance, Portable Implementation of the MPI Message Passing Interface Standard , 1996, Parallel Comput..

[17]  Greg Burns,et al.  LAM: An Open Cluster Environment for MPI , 2002 .

[18]  E. Tadmor.  SPECTRAL METHODS FOR HYPERBOLIC PROBLEMS , 1994 .

[19]  E. Lusk,et al.  Installation guide to mpich, a portable implementation of MPI , 1996 .

[20]  Oscar Reula,et al.  Multi-block simulations in general relativity: high-order discretizations, numerical stability and applications , 2005, Classical and Quantum Gravity.

[21]  William Gropp,et al.  User's Guide for mpich, a Portable Implementation of MPI Version 1.2.2 , 1996 .

[22]  E. Schnetter,et al.  A multi-block infrastructure for three-dimensional time-dependent numerical relativity , 2006, gr-qc/0602104.

[23]  H. Kreiss,et al.  Finite Element and Finite Difference Methods for Hyperbolic Partial Differential Equations , 1974 .

[24]  Jan Nordström,et al.  High-order finite difference methods, multidimensional linear problems, and curvilinear coordinates , 2001 .

[25]  Andrew Lumsdaine,et al.  A Component Architecture for LAM/MPI , 2003, PVM/MPI.

[26]  Barbara Zubik-Kowal,et al.  Spectral Versus Pseudospectral Solutions of the Wave Equation by Waveform Relaxation Methods , 2004, J. Sci. Comput..

[27]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[28]  Jeffrey Rauch,et al.  Symmetric positive systems with boundary characteristic of constant multiplicity , 1985 .

[29]  Jack Dongarra,et al.  Vector and Parallel Processing — VECPAR 2000 , 2001, Lecture Notes in Computer Science.

[30]  Scott H. Hawley,et al.  Evolutions in 3D numerical relativity using fixed mesh refinement , 2003, gr-qc/0310042.

[31]  William Gropp,et al.  Users guide for mpich, a portable implementation of MPI , 1996 .

[32]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .