pith. sign in

arxiv: 2606.06247 · v1 · pith:3G6NPE5Mnew · submitted 2026-06-04 · ⚛️ physics.flu-dyn

A high-order Fourier Continuation (FC)-based spectral incompressible Smoothed Particle Hydrodynamics (ISPH) scheme for general boundary conditions in wall-bounded domains

Meixuan Lin , Georgios Fourtakas , Benedict D.Rogers This is my paper

Pith reviewed 2026-06-27 23:32 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords Fourier Continuationspectral ISPHwall-bounded flowsincompressible SPHhigh-order convergencevortex dynamicsprojection methodFFT discretization
0
0 comments X

The pith

Incorporating Fourier Continuation into spectral ISPH enables high-order simulation of wall-bounded incompressible flows with general boundary conditions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops a Fourier Continuation method to extend a spectral smoothed particle hydrodynamics scheme from periodic to wall-bounded domains. Velocity and pressure are extended using polynomial-based continuations so that the domain becomes both periodic and sufficiently smooth. The SPH discretization then proceeds in frequency space via the convolution theorem and fast Fourier transforms. Combined with a projection time integrator and spectral pressure Poisson solver, the resulting scheme is tested on classical benchmarks. The central result is that the approach delivers high-order convergence while capturing vortex structures in flows subject to non-periodic wall boundaries.

Core claim

A polynomial-based Fourier Continuation algorithm is applied to the velocity and pressure fields to render the computational domain periodic and C^p smooth; the spectral ISPH discretization is then performed in frequency space by FFT, and the combination with projection-based time integration and a spectral PPE solver produces high-order convergence together with accurate capture of complex vortex dynamics for wall-bounded incompressible flows under general boundary conditions.

What carries the argument

Polynomial-based Fourier Continuation (FC) extension of velocity and pressure, which enforces periodicity and C^p smoothness across the domain for subsequent FFT-based frequency-space SPH discretization regardless of boundary condition type.

Load-bearing premise

The polynomial-based Fourier continuation extension of velocity and pressure preserves the spectral accuracy and stability of the subsequent SPH discretization in frequency space without introducing new truncation or aliasing errors at the artificial periodic interfaces.

What would settle it

A sequence of lid-driven cavity simulations at successively doubled particle resolutions in which the measured L2 error norms fail to decrease at the claimed high-order rate after the FC extension is applied.

Figures

Figures reproduced from arXiv: 2606.06247 by Benedict D.Rogers, Georgios Fourtakas, Meixuan Lin.

Figure 1
Figure 1. Figure 1: Zero padding of the SPH kernel gradient ∂W ∂x , shown as the centre-row cross￾section: (a) original kernel gradient; (b) circularly shifted, kernel wrapped to the corners for circular convolution via FFT. WG8(q) = αG8  4 − 6q 2 + 2q 4 − q 6 6  e −q 2 (13) WG10(q) = αG10  5 − 10q 2 + 4q 4 − 2q 6 6 + q 8 24 e −q 2 (14) where q = r/h is the normalised distance, h is the smoothing length and is set to 8dx … view at source ↗
Figure 2
Figure 2. Figure 2: Real part of the Gaussian kernel function Re [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Resolving power of different orders of the Gaussian kernel functions for the [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Polynomial fitting and extrapolation of the Fourier continuation on the example [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Blending function σ(t) and FC-extended area of the example function. 4.4. Incorporation of the Neumann Boundary Condition For the wall-bounded flows simulated in this work, pressure satisfies Neu￾mann boundary conditions, therefore modifications are needed for the afore￾mentioned FC algorithm to account for the pressure gradient at the bound￾aries ∂P ∂n = 0. The derivative of the fitted polynomial p(ξ) is … view at source ↗
Figure 6
Figure 6. Figure 6: Fourier continuation of a non-periodic function [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The computational procedure of one time step by the Fourier continuation-based [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Convergence rate of gradient and Laplacian operators by the FC-based spectral [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Convergence rate of gradient operator by the FC-based spectral SPH scheme [PITH_FULL_IMAGE:figures/full_fig_p018_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Velocity profile of the Gaussian vortex leaving the domain at different time. [PITH_FULL_IMAGE:figures/full_fig_p019_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Convergence rate of u in the Gaussian vortex convection case. Under these conditions, the steady-state velocity profile is governed by the analytical solution shown in Equation (43): u(y) = Fx 2ν y(1.0 − y). (43) All simulations are initialised with the analytical velocity profile at a Reynolds number of Re=100 and all the cases below are run for t =1 s at which point the L2 norm of the velocities field h… view at source ↗
Figure 12
Figure 12. Figure 12: The L2 evolution of the velocities in the Poiseuille flow case. 6.4. Plane Couette flow The previous test case validates the capability of the scheme to simulate the homogenous no-slip boundary condition. In order to demonstrate that the FC scheme can be applied generally, the plane Couette flow is used as a test case, where the boundary conditions are u = 0 at y = 0 (stationary wall) and u = 1 m/s at y =… view at source ↗
Figure 13
Figure 13. Figure 13: Pressure profile at t = 1.0s of the Poiseuille flow case. 0.005 0.01 0.015 0.02 dx 10-10 10-9 10-8 10-7 10-6 10-5 10-4 L 2 error of u FC-spectral SPH 4th order reference 5th order reference [PITH_FULL_IMAGE:figures/full_fig_p021_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Convergence rate at t = 1 s of the Poiseuille flow case. and is advanced in time until the velocity field u converges to the analytical steady-state solution of plane Couette flow u(y) = y. Before reaching steady state, velocity profile u is compared to the transient analytical solution at 21 [PITH_FULL_IMAGE:figures/full_fig_p021_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Boundary condition of the plane Couette flow. [PITH_FULL_IMAGE:figures/full_fig_p022_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Comparison of the analytical transient velocity solution and the velocity ob [PITH_FULL_IMAGE:figures/full_fig_p023_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Vorticity contours of the vortex dipole at [PITH_FULL_IMAGE:figures/full_fig_p025_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Comparison of the TKE and the enstrophy for the FC-based spectral ISPH [PITH_FULL_IMAGE:figures/full_fig_p026_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: Vorticity contours of the vortex dipole at [PITH_FULL_IMAGE:figures/full_fig_p028_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: Vorticity at y = −1 at different t from the FC-spectral ISPH scheme and the Chebyshev pseudospectral solver in [30]. 29 [PITH_FULL_IMAGE:figures/full_fig_p029_20.png] view at source ↗
read the original abstract

In this paper, a high-order Fourier Continuation (FC) algorithm is introduced into the spectral smoothed particle hydrodynamics (SPH) scheme to simulate the wall-bounded incompressible flows. This work aims to extend the spectral ISPH scheme towards the high-order simulation of flows with non-periodic wall boundary conditions. Herein, a polynomial-based Fourier continuation technique is applied to the velocity and pressure to make the domain both periodic and Cp smooth. The spatial SPH discretisation is performed subsequently in the frequency space on the FC-extended domain by building upon the convolution theorem using fast Fourier transform (FFT). The incorporation of Neumann boundary conditions is straightforward, and more generally, the FC method enforces periodicity across the domain regardless of the boundary condition type. The convergence order, additional computational cost, and implementation technique of the FC method are also discussed. Combined with a projection-based time integration scheme and a spectral PPE solver, the FC-based spectral ISPH framework is validated against several classical CFD benchmarks. The principal finding of this work is that the incorporation of FC techniques enables the spectral ISPH scheme to simulate wall-bounded flows with high-order convergence, and accurately capturing complex vortex dynamics. This work therefore represents a step towards a fully high-order spectral Lagrangian SPH solver with complex geometries

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript introduces a high-order Fourier Continuation (FC)-based spectral incompressible Smoothed Particle Hydrodynamics (ISPH) scheme for wall-bounded incompressible flows with general boundary conditions. A polynomial-based FC technique is applied to velocity and pressure fields to enforce periodicity and C^p smoothness on the domain; spatial discretization then proceeds in frequency space via the convolution theorem and FFT. The scheme is combined with a projection-based time integrator and a spectral pressure Poisson equation (PPE) solver. Validation is performed on classical CFD benchmarks, with the central claim that FC incorporation enables high-order convergence while accurately capturing complex vortex dynamics in non-periodic domains.

Significance. If the quantitative high-order convergence claim holds after verification of the FC extension step, the work would constitute a meaningful advance toward fully spectral Lagrangian solvers for wall-bounded flows, extending the applicability of FFT-based SPH operators beyond periodic domains without loss of accuracy order. The approach of using FC to enable periodicity regardless of boundary-condition type is a clear technical strength if interface artifacts are shown to be controlled.

major comments (2)
  1. [Abstract] Abstract: the assertion of 'high-order convergence' and 'validation on classical benchmarks' is unsupported by any quantitative convergence rates, L2 or L-infinity error tables, or explicit verification that post-FC SPH operators retain spectral accuracy; without these data the central claim that polynomial FC preserves the order cannot be assessed.
  2. [FC method and discretization] FC extension and frequency-space discretization sections: the load-bearing assumption that the polynomial-based FC step (applied prior to FFT convolution) introduces neither truncation nor aliasing errors at the artificial periodic interfaces is stated but not accompanied by an error analysis or numerical test isolating the extension operator; any such artifact would directly degrade the global convergence order for wall-bounded cases independent of the projection or PPE components.
minor comments (2)
  1. [Implementation technique discussion] Clarify the precise polynomial degree and continuity order C^p chosen for the FC extension and state whether these parameters are held fixed across all reported benchmarks.
  2. [Convergence order and cost discussion] The additional computational cost of the FC step is mentioned but not quantified relative to the baseline spectral ISPH; a table or scaling plot would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough and constructive review. The two major comments identify areas where additional quantitative support and isolated testing would strengthen the manuscript. We address each point below and have revised the manuscript to incorporate the requested data and analysis.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the assertion of 'high-order convergence' and 'validation on classical benchmarks' is unsupported by any quantitative convergence rates, L2 or L-infinity error tables, or explicit verification that post-FC SPH operators retain spectral accuracy; without these data the central claim that polynomial FC preserves the order cannot be assessed.

    Authors: We agree that the abstract claims require explicit quantitative backing. While the original manuscript discusses the convergence order of the FC method, it does not present tabulated L2/L∞ errors or rates for the post-FC operators. In the revised version we have added a dedicated results subsection containing convergence tables (L2 and L∞ norms versus particle spacing) for velocity and pressure on the benchmark problems, together with a direct comparison of spectral accuracy before and after the FC extension step. revision: yes

  2. Referee: [FC method and discretization] FC extension and frequency-space discretization sections: the load-bearing assumption that the polynomial-based FC step (applied prior to FFT convolution) introduces neither truncation nor aliasing errors at the artificial periodic interfaces is stated but not accompanied by an error analysis or numerical test isolating the extension operator; any such artifact would directly degrade the global convergence order for wall-bounded cases independent of the projection or PPE components.

    Authors: This observation is correct. The original text asserts the smoothness and periodicity properties of the polynomial FC but does not isolate the extension operator with a dedicated error study. We have added a new subsection that applies the FC procedure to known analytic functions with non-periodic boundary data, computes the pointwise and L2 errors in the extended periodic domain, and verifies that interface truncation/aliasing remain below the level of the subsequent spectral truncation, thereby confirming that the global convergence order is preserved. revision: yes

Circularity Check

0 steps flagged

No significant circularity; method introduction and benchmark validation are independent of inputs

full rationale

The paper presents a new numerical technique combining Fourier Continuation with spectral ISPH for wall-bounded flows, followed by validation on classical CFD benchmarks. No equations, parameters, or quantities are fitted to a data subset and then relabeled as predictions of closely related quantities. No self-citations are invoked as load-bearing uniqueness theorems or to smuggle in ansatzes. The convergence claims rest on numerical experiments rather than reducing by construction to the method's own definitions or prior self-referential results. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that the FC extension can be treated as an exact periodicization step whose only effect is to enable FFT operations; no free parameters, new physical entities, or ad-hoc axioms are named in the abstract.

axioms (1)
  • standard math Convolution theorem holds for the SPH kernel after the FC-extended fields are made periodic and C^p smooth.
    Invoked when spatial discretization is performed in frequency space via FFT.

pith-pipeline@v0.9.1-grok · 5778 in / 1290 out tokens · 20223 ms · 2026-06-27T23:32:47.278020+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 5 canonical work pages

  1. [1]

    D. J. Price, Smoothed particle hydrodynamics and magnetohydrody- namics, Journal of Computational Physics 231 (3) (2012) 759–794

    2012

  2. [2]

    D. J. Price, J. Wurster, T. S. Tricco, C. Nixon, S. Toupin, A. Pet- titt, C. Chan, D. Mentiplay, G. Laibe, S. Glover, et al., Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics, Publications of the Astronomical Society of Australia 35 (2018). doi:10.1017/pasa.2018.25

    work page doi:10.1017/pasa.2018.25 2018

  3. [3]

    Monaghan, A

    J. Monaghan, A. Kocharyan, Sph simulation of multi-phase flow, Com- puter Physics Communications 87 (1) (1995) 225–235

    1995

  4. [4]

    X. Hu, N. Adams, A multi-phase sph method for macroscopic and meso- scopic flows, Journal of Computational Physics 213 (2) (2006) 844–861

    2006

  5. [5]

    Z. Chen, Z. Zong, M. Liu, L. Zou, H. Li, C. Shu, An sph model for multiphase flows with complex interfaces and large density differences, Journal of Computational Physics 283 (2015) 169–188

    2015

  6. [6]

    Vacondio, C

    R. Vacondio, C. Altomare, M. D. Leffe, X. Hu, D. L. Touzé, S. Lind, J.-C.Marongiu, S.Marrone, B.D.Rogers, A.Souto-Iglesias, Grandchal- lenges for smoothed particle hydrodynamics numerical schemes, Com- putational Particle Mechanics 8 (2021) 575–588

    2021

  7. [7]

    Moxey, C

    D. Moxey, C. D. Cantwell, Y. Bao, A. Cassinelli, G. Castiglioni, S. Chun, E. Juda, E. Kazemi, K. Lackhove, J. Marcon, G. Mengaldo, D. Serson, M. Turner, H. Xu, J. Peiró, R. M. Kirby, S. J. Sherwin, Nektar++: Enhancing the capability and application of high-fidelity spectral/hp el- ement methods, Computer Physics Communications 249 (2020) 107110

    2020

  8. [8]

    S. J. Lind, B. D. Rogers, P. K. Stansby, Review of smoothed particle hydrodynamics: towards converged lagrangian flow modelling, Proceed- ings of the Royal Society A: Mathematical, Physical and Engineering Sciences 476 (2241) (2020) 20190801

    2020

  9. [9]

    Meng, P.-N

    Z.-F. Meng, P.-N. Sun, P.-P. Wang, B. C. Khoo, A.-M. Zhang, High-order SPH: A review of the method and applications, Archives of Computational Methods in Engineering 33 (2025) 1409–1444. doi:10.1007/s11831-025-10346-0. 32

    work page doi:10.1007/s11831-025-10346-0 2025

  10. [10]

    Sibilla, An algorithm to improve consistency in smoothed particle hydrodynamics, Computers & Fluids 118 (2015) 148–158

    S. Sibilla, An algorithm to improve consistency in smoothed particle hydrodynamics, Computers & Fluids 118 (2015) 148–158

    2015

  11. [11]

    Nasar, G

    A. Nasar, G. Fourtakas, S. Lind, J. King, B. Rogers, P. Stansby, High- order consistent sph with the pressure projection method in 2-d and 3-d, Journal of Computational Physics 444 (2021) 110563

    2021

  12. [12]

    S. Lind, P. Stansby, High-order eulerian incompressible smoothed par- ticle hydrodynamics with transition to lagrangian free-surface motion, Journal of Computational Physics 326 (2016) 290–311

    2016

  13. [13]

    Nasar, G

    A. Nasar, G. Fourtakas, S. Lind, B. Rogers, P. Stansby, J. King, High- order velocity and pressure wall boundary conditions in eulerian incom- pressible sph, Journal of Computational Physics 434 (2021) 109793

    2021

  14. [14]

    Z. Wang, B. Zhang, O. J. Haidn, X. Hu, A fourth-order kernel for im- proving numerical accuracy and stability in eulerian sph for fluids and total lagrangian sph for solids, Journal of Computational Physics 519 (2024) 113385

    2024

  15. [15]

    T. Gao, T. Liang, L. Fu, A new smoothed particle hydrodynamics method based on high-order moving-least-square targeted essentially non-oscillatory scheme for compressible flows, Journal of Computational Physics 489 (2023) 112270

    2023

  16. [16]

    Meng, P.-N

    Z.-F. Meng, P.-N. Sun, Y. Xu, P.-P. Wang, A.-M. Zhang, High-order eulerian sph scheme through w/teno reconstruction based on primitive variables for simulating incompressible flows, Computer Methods in Ap- plied Mechanics and Engineering 427 (2024) 117065

    2024

  17. [17]

    Vergnaud, G

    A. Vergnaud, G. Oger, D. Le Touzé, Investigations on a high order sph scheme using weno reconstruction, Journal of Computational Physics 477 (2023) 111889

    2023

  18. [18]

    M. Lin, G. Fourtakas, B. D. Rogers, A novel high-order spectral in- compressible smoothed particle hydrodynamics (isph) scheme with an immersed boundary method (ibm), Journal of Computational Physics 540 (2025) 114264

    2025

  19. [19]

    L. N. Trefethen, Spectral Methods in MATLAB, Society for Industrial and Applied Mathematics, 2000. 33

    2000

  20. [20]

    O. P. Bruno, M. Lyon, High-order unconditionally stable fc-ad solvers for general smooth domains i. basic elements, Journal of Computational Physics 229 (6) (2010) 2009–2033

    2010

  21. [21]

    Albin, O

    N. Albin, O. P. Bruno, A spectral FC solver for the compress- ible Navier–Stokes equations in general domains I: Explicit time- stepping, Journal of Computational Physics 230 (16) (2011) 6248–6270. doi:10.1016/j.jcp.2011.04.023

    work page doi:10.1016/j.jcp.2011.04.023 2011

  22. [22]

    M. Lyon, O. P. Bruno, High-order unconditionally stable fc-ad solvers for general smooth domains ii. elliptic, parabolic and hyperbolic pdes; theoretical considerations, Journal of Computational Physics 229 (9) (2010) 3358–3381

    2010

  23. [23]

    S. J. Cummins, M. Rudman, An SPH Projection Method, Journal of Computational Physics 152 (2) (1999) 584–607

    1999

  24. [24]

    A. J. Chorin, Numerical solution of the navier-stokes equations, Math- ematics of Computation 22 (104) (1968) 745–762

    1968

  25. [25]

    Fontana, O

    M. Fontana, O. P. Bruno, P. D. Mininni, P. Dmitruk, Fourier con- tinuation method for incompressible fluids with boundaries, Computer Physics Communications 256 (2020) 107482

    2020

  26. [26]

    Lyon, A fast algorithm for Fourier continuation, SIAM Journal on Scientific Computing 33 (6) (2011) 3241–3260

    M. Lyon, A fast algorithm for Fourier continuation, SIAM Journal on Scientific Computing 33 (6) (2011) 3241–3260. doi:10.1137/11082436X

    work page doi:10.1137/11082436x 2011

  27. [27]

    V.Fuka, Poisfft–afreeparallelfastpoissonsolver, AppliedMathematics and Computation 267 (2015) 356–364

    2015

  28. [28]

    J. P. Boyd, Chebyshev and Fourier Spectral Methods, 2nd Edition, Dover Publications, Mineola, New York, 2001

    2001

  29. [29]

    X. Wang, J. Wang, X. Wang, C. Yu, A pseudo-spectral fourier collocation method for inhomogeneous elliptical inclusions with partial differential equations, Mathematics 10 (3) (2022) 296. doi:10.3390/math10030296

    work page doi:10.3390/math10030296 2022

  30. [30]

    Clercx, C.-H

    H. Clercx, C.-H. Bruneau, The normal and oblique collision of a dipole with a no-slip boundary, Computers & Fluids 35 (3) (2006) 245–279. 34

    2006

  31. [31]

    Keetels, U

    G. Keetels, U. D’Ortona, W. Kramer, H. Clercx, K. Schneider, G. van Heijst, Fourier spectral and wavelet solvers for the incompressible navier–stokes equations with volume-penalization: Convergence of a dipole–wall collision, Journal of Computational Physics 227 (2) (2007) 919–945

    2007

  32. [32]

    Laizet, E

    S. Laizet, E. Lamballais, High-order compact schemes for incompress- ible flows: A simple and efficient method with quasi-spectral accuracy, Journal of Computational Physics 228 (16) (2009) 5989–6015

    2009

  33. [33]

    O. P. Bruno, J. Paul, Two-dimensional fourier continuation and applica- tions, SIAM Journal on Scientific Computing 44 (2) (2022) A964–A992. 35

    2022