Two-Dimensional Fourier Continuation and Applications
From MaRDI portal
Publication:5075686
DOI10.1137/20M1373189zbMATH Open1490.35109arXiv2010.03901OpenAlexW3092632390MaRDI QIDQ5075686FDOQ5075686
Publication date: 11 May 2022
Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)
Abstract: 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.
Full work available at URL: https://arxiv.org/abs/2010.03901
Trigonometric approximation (42A10) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Wave equation (35L05) Trigonometric interpolation (42A15)
Cites Work
- Title not available (Why is that?)
- Title not available (Why is that?)
- Title not available (Why is that?)
- Title not available (Why is that?)
- Linear integral equations
- Partial differential equations. 4th ed
- On the Fourier Extension of Nonperiodic Functions
- A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas
- A comparison of numerical algorithms for Fourier extension of the first, second, and third kinds
- High-order unconditionally stable FC-AD solvers for general smooth domains. I: Basic elements
- Accurate, high-order representation of complex three-dimensional surfaces via Fourier continuation analysis
- Rapidly convergent two-dimensional quasi-periodic Green function throughout the spectrum-including wood anomalies
- A spectral FC solver for the compressible Navier-Stokes equations in general domains. I: Explicit time-stepping
- High-order unconditionally stable FC-AD solvers for general smooth domains. II: Elliptic, parabolic and hyperbolic PDEs; theoretical considerations
- On the numerical stability of Fourier extensions
- A Fourth-Order Accurate Embedded Boundary Method for the Wave Equation
- Immersed boundary smooth extension: a high-order method for solving PDE on arbitrary smooth domains using Fourier spectral methods
- An FC-based spectral solver for elastodynamic problems in general three-dimensional domains
- An adaptive fast multipole accelerated Poisson solver for complex geometries
- Partition of unity extension of functions on complex domains
- Function Approximation on Arbitrary Domains Using Fourier Extension Frames
- FFT-based high order central difference schemes for three-dimensional Poisson's equation with various types of boundary conditions
- A boundary integral algorithm for the Laplace Dirichlet-Neumann mixed eigenvalue problem
- 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
Cited In (8)
- A fast, high-order scheme for evaluating volume potentials on complex 2D geometries via area-to-line integral conversion and domain mappings
- A fictitious domain spectral method for solving the Helmholtz equation in exterior domains
- A double Fourier sphere method for \(d\)-dimensional manifolds
- Rapid Evaluation of Newtonian Potentials on Planar Domains
- Numerical analytic continuation
- An FMM Accelerated Poisson Solver for Complicated Geometries in the Plane Using Function Extension
- Spectrally accurate solutions to inhomogeneous elliptic PDE in smooth geometries using function intension
- Title not available (Why is that?)
This page was built for publication: Two-Dimensional Fourier Continuation and Applications
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5075686)
