Two-dimensional Fourier Continuation and applications

This paper presents a "two-dimensional Fourier Continuation" method (2D-FC) for construction of bi-periodic extensions of smooth non-periodic functions defined over general two-dimensional smooth domains. The approach can be directly generalized to domains of any given dimensionality, and even to non-smooth domains, but such generalizations are not considered here. The 2D-FC extensions are produced in a two-step procedure. In the first step the one-dimensional Fourier Continuation method is applied along a discrete set of outward boundary-normal directions to produce, along such directions, continuations that vanish outside a narrow interval beyond the boundary. Thus, the first step of the algorithm produces "blending-to-zero along normals" for the given function values. In the second step, the extended function values are evaluated on an underlying Cartesian grid by means of an efficient, high-order boundary-normal interpolation scheme. A Fourier Continuation expansion of the given function can then be obtained by a direct application of the two-dimensional FFT algorithm. Algorithms of arbitrarily high order of accuracy can be obtained by this method. The usefulness and performance of the proposed two-dimensional Fourier Continuation method are illustrated with applications to the Poisson equation and the time-domain wave equation within a bounded domain. As part of these examples the novel "Fourier Forwarding" solver is introduced which, propagating plane waves as they would in free space and relying on certain boundary corrections, can solve the time-domain wave equation and other hyperbolic partial differential equations within general domains at computing costs that grow sublinearly with the size of the spatial discretization.

[1]  Erik Lehto,et al.  Partition of unity extension of functions on complex domains , 2017, J. Comput. Phys..

[2]  Robert D. Guy,et al.  Immersed boundary smooth extension: A high-order method for solving PDE on arbitrary smooth domains using Fourier spectral methods , 2015, J. Comput. Phys..

[3]  Antoine J. Cerfon,et al.  An adaptive fast multipole accelerated Poisson solver for complex geometries , 2016, J. Comput. Phys..

[4]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[5]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[6]  Daan Huybrechs,et al.  On the Fourier Extension of Nonperiodic Functions , 2010, SIAM J. Numer. Anal..

[7]  Mark Lyon,et al.  High-order unconditionally stable FC-AD solvers for general smooth domains II. Elliptic, parabolic and hyperbolic PDEs; theoretical considerations , 2010, J. Comput. Phys..

[8]  Mark Lyon,et al.  High-order unconditionally stable FC-AD solvers for general smooth domains I. Basic elements , 2010, J. Comput. Phys..

[9]  Oscar P. Bruno,et al.  An FC-based spectral solver for elastodynamic problems in general three-dimensional domains , 2016, J. Comput. Phys..

[10]  William H. Press,et al.  Numerical Recipes 3rd Edition: The Art of Scientific Computing , 2007 .

[11]  Nilima Nigam,et al.  A boundary integral algorithm for the Laplace Dirichlet-Neumann mixed eigenvalue problem , 2014, J. Comput. Phys..

[12]  Jean-Luc Vay,et al.  PPPS-2013: Topic 1.2: A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas , 2013, 2013 Abstracts IEEE International Conference on Plasma Science (ICOPS).

[13]  Kareem T. Elgindy.  A high-order embedded domain method combining a Predictor-Corrector-Fourier-Continuation-Gram method with an integral Fourier pseudospectral collocation method for solving linear partial differential equations in complex domains , 2019, J. Comput. Appl. Math..

[14]  Nathan Albin,et al.  A spectral FC solver for the compressible Navier-Stokes equations in general domains I: Explicit time-stepping , 2011, J. Comput. Phys..

[15]  Shan Zhao,et al.  FFT-based high order central difference schemes for three-dimensional Poisson's equation with various types of boundary conditions , 2020, J. Comput. Phys..

[16]  Robert H. Halstead,et al.  Matrix Computations , 2011, Encyclopedia of Parallel Computing.

[17]  Oscar P. Bruno,et al.  Rapidly convergent two-dimensional quasi-periodic Green function throughout the spectrum - including Wood anomalies , 2014, J. Comput. Phys..

[18]  Daan Huybrechs,et al.  On the Numerical Stability of Fourier Extensions , 2012, Foundations of Computational Mathematics.

[19]  R. Kanwal.  Linear Integral Equations , 1925, Nature.

[20]  N. Anders Petersson,et al.  A Fourth-Order Accurate Embedded Boundary Method for the Wave Equation , 2012, SIAM J. Sci. Comput..