Skip to content
DraconiticResearch
Search Sources API
201,683 studies31 sourcesUpdated 38 min ago
PreprintCC BY 4.0arXiv

Lattice Green's functions for Poisson problems with two unbounded and one periodic direction

  • Gilles Poncelet
  • Philippe Chatelain
Open at sourcePDFDOICite

Abstract

Fourier-based Poisson solvers treat periodic directions spectrally and unbounded ones by the Hockney--Eastwood convolution with a free-space kernel. With two unbounded and one periodic direction, the setting of vortex tubes, wakes and jets, the kernel used so far by these solvers combines the two-dimensional kernel on the zero periodic mode with the symbol of the operator sampled on the doubled grid on every other mode. We show that the second ingredient is not exact: the error is the field of the periodic images of the source and independent of the mesh size. While the period Lz is not much longer than the unbounded extent Lx, it is of order exp(−2πd/Lz), with d the distance to the nearest image; beyond, it grows like (Lz/Lx)2. The error is thus small on cubic domains, but larger than the solution itself on elongated ones. We remove the approximation by extending the integral representation of the exact second-order kernel of Hou and Colonius: the screened kernel of each mode is evaluated in closed form along one unbounded direction and by Gauss--Legendre quadrature along the other, for the high-order dimension-split and Mehrstellen stencils, at a one-off cost of about one solve. The kernels, implemented in the open-source library \flups{}, are verified by comparison with independent quadratures, by a padding test sensitive to the boundary condition enforced, and by convergence at the order of each stencil on an axially modulated tube source, on domains up to $64$ times longer than wide and up to 2.7×1011 unknowns.

Related

  1. Formalizing Fefferman's Official Periodic Navier-Stokes Problem (B) in Lean 4: Global Linear Strain Exclusion on T3, Biot-Savart Direction Depletion Identities, and Conditional Reductions via Beale-Kato-Majda
    Benjamin Frohman
    Zenodo (CERN European Organization for Nuclear Research)PreprintLicence not stated5 citedGeneral mathematicsOpenAlex
  2. An Upper Bound for the Mean Speed of Transition Fronts and Unbounded Front Widths for Fisher KPP Equations in Almost Periodic Media
    Xing Liang, Linfeng Xu, Qi Zhou +1
    arXivPreprintCC BY 4.0Full textAnalysis & PDEsarXiv
  3. Last-Iterate Convergence of Stochastic Gradient Methods for Unbounded Problems
    Kevin Kurian Thomas Vaidyan, Ahmet Alacaoglu, Michael P. Friedlander
    arXivPreprintarXiv licenceFull textOptimisation & controlStatisticsarXiv
  4. Numerical Construction of Quasi-Periodic Solutions for Nonlinear PDEs: II. Unbounded Perturbation
    Mingwei Fu, Bin Shi
    arXivPreprintCC BY-NC-SA 4.0Full textNumerical analysis & scientific computingarXiv
  5. A Localized Core-Tail Fourier-Laguerre Frame Method with Adaptive Frequency Detection for Unbounded-Domain Problems
    Chenyang Wang, Zhenyu Zhao, Tinggang Zhao
    arXivPreprintarXiv licenceFull textNumerical analysis & scientific computingarXiv

Cite

BibTeXRIS
BibTeX
@misc{poncelet2026lattice,
  title = {{Lattice Green's functions for Poisson problems with two unbounded and one periodic direction}},
  author = {Gilles Poncelet and Philippe Chatelain},
  howpublished = {arXiv},
  year = {2026},
  month = nov,
  doi = {10.48550/arxiv.2611.03458},
  eprint = {2611.03458},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2611.03458},
  note = {CC BY 4.0},
}
RIS
TY  - UNPB
TI  - Lattice Green's functions for Poisson problems with two unbounded and one periodic direction
AU  - Gilles Poncelet
AU  - Philippe Chatelain
PY  - 2026
DA  - 2026/11/02
PB  - arXiv
DO  - 10.48550/arxiv.2611.03458
UR  - https://arxiv.org/abs/2611.03458
AB  - Fourier-based Poisson solvers treat periodic directions spectrally and unbounded ones by the Hockney--Eastwood convolution with a free-space kernel. With two unbounded and one periodic direction, the setting of vortex tubes, wakes and jets, the kernel used so far by these solvers combines the two-dimensional kernel on the zero periodic mode with the symbol of the operator sampled on the doubled grid on every other mode. We show that the second ingredient is not exact: the error is the field of the periodic images of the source and independent of the mesh size. While the period L_z is not much longer than the unbounded extent L_x, it is of order exp(−2πd/L_z), with d the distance to the nearest image; beyond, it grows like (L_z/L_x)^2. The error is thus small on cubic domains, but larger than the solution itself on elongated ones. We remove the approximation by extending the integral representation of the exact second-order kernel of Hou and Colonius: the screened kernel of each mode is evaluated in closed form along one unbounded direction and by Gauss--Legendre quadrature along the other, for the high-order dimension-split and Mehrstellen stencils, at a one-off cost of about one solve. The kernels, implemented in the open-source library \flups{}, are verified by comparison with independent quadratures, by a padding test sensitive to the boundary condition enforced, and by convergence at the order of each stencil on an axially modulated tube source, on domains up to $64$ times longer than wide and up to 2.7×10^11 unknowns.
ER  -