Optimal grids for anisotropic problems

Spectral convergence of optimal grids for anisotropic problems is both numerically observed and explained. For elliptic problems, the gridding algorithm is reduced to a Stieltjes rational approximation on an interval of a line in the complex plane instead of the real axis as in the isotropic case. We show rigorously why this occurs for a semi-infinite and bounded interval. We then extend the gridding algorithm to hyperbolic problems on bounded domains. For the propagative modes, the problem is reduced to a rational approximation on an interval of the negative real semiaxis, similarly to in the isotropic case. For the wave problem we present numerical examples in 2-D anisotropic media.

[1]  W. Mccrea,et al.  Boundary Conditions for the Wave Equation , 1933 .

[2]  A. Majda,et al.  Radiation boundary conditions for acoustic and elastic wave calculations , 1979 .

[3]  John Sylvester,et al.  The Dirichlet to Neumann map and applications , 1989 .

[4]  Alfio Quarteroni,et al.  Theory and Application of Steklov-Poincaré Operators for Boundary-Value Problems , 1991 .

[5]  Sofia Davydycheva,et al.  A Finite Difference Scheme for Elliptic Equations with Rough Coefficients Using a Cartesian Grid Nonconforming to Interfaces , 1999 .

[6]  Vladimir Druskin,et al.  Gaussian Spectral Rules for the Three-Point Second Differences: I. A Two-Point Positive Definite Problem in a Semi-Infinite Domain , 1999, SIAM J. Numer. Anal..

[7]  John L. Tassoulas,et al.  CONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION , 2000 .

[8]  Vladimir Druskin,et al.  Application of the Difference Gaussian Rules to Solution of Hyperbolic Problems , 2000 .

[9]  Vladimir Druskin,et al.  Optimal finite difference grids and rational approximations of the square root I. Elliptic problems , 2000 .

[10]  Liliana Borcea,et al.  INVERSE PROBLEMS PII: S0266-5611(02)33630-X Optimal finite difference grids for direct and inverse Sturm–Liouville problems , 2002 .

[11]  Martin J. Gander,et al.  Optimized Schwarz Methods without Overlap for the Helmholtz Equation , 2002, SIAM J. Sci. Comput..

[12]  Vladimir Druskin,et al.  Application of the difference Gaussian rules to solution of hyperbolic problems: global expansion , 2002 .

[13]  Vladimir Druskin,et al.  Three-point finite-difference schemes, Padé and the spectral Galerkin method. I. One-sided impedance approximation , 2001, Math. Comput..

[14]  Murthy N. Guddati,et al.  On Optimal Finite-Difference Approximation of PML , 2003, SIAM J. Numer. Anal..

[15]  Sofia Davydycheva,et al.  An efficient finite‐difference scheme for electromagnetic logging in 3D anisotropic inhomogeneous media , 2003 .

[16]  Vladimir Druskin,et al.  Gaussian spectral rules for second order finite-difference schemes , 2000, Numerical Algorithms.