TY - JOUR

T1 - The method of difference potentials for the helmholtz equation using compact high order schemes

AU - Medvinsky, M.

AU - Tsynkov, S.

AU - Turkel, E.

N1 - Funding Information:
Acknowledgements Work supported by the US NSF under grant # DMS-0810963, US–Israel Binational Science Foundation (BSF) under grant # 2008094, US AFOSR under grant # FA9550-10-1-0092, and US ARO under grant # W911NF-11-1-0384.

PY - 2012/10

Y1 - 2012/10

N2 - 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 level of a finite-difference scheme on a regular structured grid. The method of difference potentials assembles the overall solution of the original boundary value problem by repeatedly solving an auxiliary problem. This auxiliary problem allows a considerable degree of flexibility in its formulation and can be chosen so that it is very efficient to solve. Compact finite difference schemes enable high order accuracy on small stencils at virtually no extra cost. The scheme attains consistency only on the solutions of the differential equation rather than on a wider class of sufficiently smooth functions. Unlike standard high order schemes, compact approximations require no additional boundary conditions beyond those needed for the differential equation itself. However, they exploit two stencilsone applies to the left-hand side of the equation and the other applies to the right-hand side of the equation. We shall show how to properly define and compute the difference potentials and boundary projections for compact schemes. The combination of the method of difference potentials and compact schemes yields an inexpensive numerical procedure that offers high order accuracy for non-conforming smooth curvilinear boundaries on regular grids. We demonstrate the capabilities of the resulting method by solving the inhomogeneous Helmholtz equation with a variable wavenumber with high order (4 and 6) accuracy on Cartesian grids for non-conforming boundaries such as circles and ellipses.

AB - 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 level of a finite-difference scheme on a regular structured grid. The method of difference potentials assembles the overall solution of the original boundary value problem by repeatedly solving an auxiliary problem. This auxiliary problem allows a considerable degree of flexibility in its formulation and can be chosen so that it is very efficient to solve. Compact finite difference schemes enable high order accuracy on small stencils at virtually no extra cost. The scheme attains consistency only on the solutions of the differential equation rather than on a wider class of sufficiently smooth functions. Unlike standard high order schemes, compact approximations require no additional boundary conditions beyond those needed for the differential equation itself. However, they exploit two stencilsone applies to the left-hand side of the equation and the other applies to the right-hand side of the equation. We shall show how to properly define and compute the difference potentials and boundary projections for compact schemes. The combination of the method of difference potentials and compact schemes yields an inexpensive numerical procedure that offers high order accuracy for non-conforming smooth curvilinear boundaries on regular grids. We demonstrate the capabilities of the resulting method by solving the inhomogeneous Helmholtz equation with a variable wavenumber with high order (4 and 6) accuracy on Cartesian grids for non-conforming boundaries such as circles and ellipses.

KW - Boundary projections

KW - Calderon's operators

KW - Compact differencing

KW - Curvilinear boundaries

KW - Difference potentials

KW - High order accuracy

KW - Regular grids

KW - Variable coefficients

UR - http://www.scopus.com/inward/record.url?scp=84865686540&partnerID=8YFLogxK

U2 - 10.1007/s10915-012-9602-y

DO - 10.1007/s10915-012-9602-y

M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???

AN - SCOPUS:84865686540

VL - 53

SP - 150

EP - 193

JO - Journal of Scientific Computing

JF - Journal of Scientific Computing

SN - 0885-7474

IS - 1

ER -