pith. sign in

arxiv: 2503.12805 · v3 · submitted 2025-03-17 · 🧮 math.NA · cs.NA

A fast Fourier spectral method for wave kinetic equation

Pith reviewed 2026-05-23 00:19 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords wave kinetic equationfast Fourier transformspectral methodwave turbulencespherical integraldouble convolutionBoltzmann operator
0
0 comments X

The pith

The wave kinetic equation can be solved by a fast Fourier spectral method that turns its nonlinear term into an FFT-evaluable double convolution.

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

The paper develops a fast Fourier spectral method for the wave kinetic equation. It reformulates the nonlinear operator as a spherical integral analogous to the Boltzmann operator. This allows conservation of mass and momentum to yield a double convolution in Fourier space that FFT can handle. The computational cost drops from O(N^{3d}) to O(M N^d log N). This matters because it enables numerical studies of wave turbulence that were previously impractical due to high dimensionality.

Core claim

The central claim is that the high-dimensional nonlinear wave kinetic operator can be reformulated as a spherical integral, and that the conservation of mass and momentum then produces a double convolution structure in Fourier space which the fast Fourier transform evaluates at much lower cost while preserving the original properties.

What carries the argument

The double convolution structure in Fourier space that arises from the spherical-integral reformulation of the wave kinetic operator and the conservation laws.

Load-bearing premise

The reformulation of the wave kinetic operator as a spherical integral is mathematically valid and preserves conservation properties exactly.

What would settle it

A side-by-side computation of the nonlinear term on a small discrete frequency set using both the original definition and the double-convolution FFT method, checking for agreement to high precision.

Figures

Figures reproduced from arXiv: 2503.12805 by Kunlun Qi, Lian Shen, Li Wang.

Figure 1
Figure 1. Figure 1: Stationary solution in 2D with S = 5, Nr = N = 128, Ns = Nsig = 12. 4.1.2. Time evolution test in 2D. Now, we perform numerical tests for the 2D equation (2.1). The time discretization is performed by the fourth-order Runge￾Kutta (RK4) methods. Example 2 (Isotropic case). We consider the following isotropic initial condition: (4.1) f 0 (k) = 1 3 (δu(|k|) + δu(|k| − 0.2)), where δu(k) is an approximated del… view at source ↗
Figure 2
Figure 2. Figure 2: Top-to-bottom view of the time evolution of the solution profile for the isotropic [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Discontinous initial condition f 0 dis in 2D for N = Nr = 64, Ns = Nsig = 12, S = 3. Example 4 (Non-isotropic case). In this example, we test a conjecture regarding the non-isotropic case. As discussed in the theoretical work [15, pp. 5-7] and demonstrated by the numerical experiment in Example 3, an isotropic initial condition can lead to a blow-up of the solution at the origin. It is conjectured that a n… view at source ↗
Figure 4
Figure 4. Figure 4: Evolution of mass and energy in 2D for N = Nr, Ns = Nsig = 12 and ∆t = 0.1. K1 K2 K3 K4 K N = 16 0.27s 0.17s 0.20s 0.18s 0.82s N = 32 2.11s 1.45s 1.60s 1.39s 6.55s N = 64 22.38s 14.67s 17.73s 16.21s 70.99s [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Time evolution of the solution profile for the non-isotropic initial condition [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Evolution of mass and energy in 3D for N = Nr, Ns = Nsig, S = 0.33 and ∆t = 0.1. cubic Schr¨odinger equation. The core idea of our approach is to reformulate the high-dimensional nonlinear wave kinetic operator as a spherical integral, drawing an analogy to the classical Boltzmann collision operator. Additionally, we extract a double-convolution structure from the numerical system by applying Fourier spect… view at source ↗
read the original abstract

The central object in wave turbulence theory is the wave kinetic equation (WKE), which is an evolution equation for wave action density and acts as the wave analog of the Boltzmann kinetic equations for particle interactions. Despite recent exciting progress in the theoretical aspects of the WKE, numerical developments have lagged behind. In this paper, we introduce a fast Fourier spectral method for solving the WKE. The key idea lies in reformulating the high-dimensional nonlinear wave kinetic operator as a spherical integral, analogous to the classical Boltzmann collision operator. The conservation of mass and momentum leads to a double convolution structure in Fourier space, which can be efficiently handled by the fast Fourier transform (FFT), reducing the computational cost from $O(N^{3d})$ to $O(M N^d \log N)$ with $N$-frequency nodes and $M \ll N^{2d-1}$ in $d$ dimensions. We demonstrate the accuracy and efficiency of the proposed method through several numerical tests in both 2D and 3D, revealing and conjecturing some interesting and unique features of this equation.

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

1 major / 2 minor

Summary. The paper claims to introduce a fast Fourier spectral method for the wave kinetic equation (WKE) by reformulating the high-dimensional nonlinear operator as a spherical integral analogous to the Boltzmann collision operator. Conservation of mass and momentum produces a double convolution structure in Fourier space that is evaluated via FFT, reducing cost from O(N^{3d}) to O(M N^d log N) with M ≪ N^{2d-1}. Accuracy and efficiency are illustrated by numerical tests in 2D and 3D, from which some features of the WKE are conjectured.

Significance. If the spherical-integral reformulation is exact, preserves the required conservation laws without hidden approximation, and the discretization with M points retains the convolution property, the method would enable feasible high-dimensional WKE simulations that are currently prohibitive. This would be a notable contribution to numerical wave turbulence, provided the complexity reduction is realized in practice.

major comments (1)
  1. [Abstract] Abstract (and the central derivation of the method): the efficiency claim O(M N^d log N) is load-bearing on the assertion that the spherical-integral reformulation yields an exact double convolution amenable to direct FFT evaluation. No error analysis, convergence rates, or explicit verification that the M-point spherical discretization preserves the convolution structure and conservation identities without residual error is supplied; numerical tests alone do not establish this.
minor comments (2)
  1. The numerical experiments should report quantitative error measures (e.g., relative L^2 errors versus reference solutions) and demonstrate observed convergence rates with respect to both N and M.
  2. Clarify the precise dispersion relation employed and how the resonance manifold is discretized on the sphere so that readers can assess whether the Boltzmann analogy holds exactly for the WKE.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback. We respond to the major comment below.

read point-by-point responses
  1. Referee: [Abstract] Abstract (and the central derivation of the method): the efficiency claim O(M N^d log N) is load-bearing on the assertion that the spherical-integral reformulation yields an exact double convolution amenable to direct FFT evaluation. No error analysis, convergence rates, or explicit verification that the M-point spherical discretization preserves the convolution structure and conservation identities without residual error is supplied; numerical tests alone do not establish this.

    Authors: The reformulation of the wave kinetic operator as a spherical integral follows exactly from the resonance condition and the conservation of mass and momentum; this yields an exact double-convolution structure in Fourier space that is evaluated by FFT. The M-point quadrature is a standard discretization of the sphere, and the numerical experiments demonstrate that the discrete operator preserves the conservation identities to machine precision. We nevertheless agree that an explicit error analysis of the quadrature and its effect on the convolution property would strengthen the presentation. We will add a dedicated subsection on the quadrature error and its impact on conservation in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No circularity: algorithmic reformulation stands on explicit conservation identities

full rationale

The paper's central contribution is a numerical algorithm that recasts the WKE collision operator as a spherical integral (via mass/momentum deltas) and then exploits the resulting double-convolution structure for FFT evaluation. This is a direct, constructive mapping whose validity is asserted from the conservation laws themselves; the complexity reduction O(N^{3d}) to O(M N^d log N) follows immediately from the algebraic form once the spherical measure is introduced. No parameter is fitted to data and then relabeled a prediction, no uniqueness theorem is imported from the authors' prior work, and no ansatz is smuggled via self-citation. The abstract and method description present the spherical-integral step as an exact equivalence that preserves the required convolution property without approximation, making the efficiency claim a straightforward consequence of the reformulation rather than a self-referential loop. Because the derivation chain is self-contained against the stated conservation identities and does not reduce any claimed result to its own inputs by construction, the circularity score is 0.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the spherical-integral reformulation and the exact emergence of the double-convolution structure from conservation laws; these are domain assumptions rather than new axioms or fitted parameters.

axioms (1)
  • domain assumption The nonlinear wave kinetic operator can be exactly rewritten as a spherical integral analogous to the classical Boltzmann collision operator.
    This equivalence is the foundational step that enables the subsequent Fourier-space convolution structure.

pith-pipeline@v0.9.0 · 5714 in / 1314 out tokens · 90842 ms · 2026-05-23T00:19:23.893109+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

45 extracted references · 45 canonical work pages

  1. [1]

    A. M. Balk and V. E. Zakharov , Stability of weak-turbulence Kolmogorov spectra , in Nonlinear waves and weak turbulence, vol. 182 of Amer. Math. Soc. Transl. Ser. 2, Amer. Math. Soc., Providence, RI, 1998, pp. 31–81. This manuscript is for review purposes only. 20 K. QI, L. SHEN AND L. WANG

  2. [2]

    J. W. Banks, T. Buckmaster, A. O. Korotkevich, G. Kovacic, and J. Shatah , Direct verification of the kinetic description of wave turbulence for finite-size systems dominated by interactions among groups of six waves , Phys. Rev. Lett., 129 (2022), pp. Paper No. 034101, 6

  3. [3]

    Buckmaster, P

    T. Buckmaster, P. Germain, Z. Hani, and J. Shatah , On the kinetic wave turbulence description for NLS , Quart. Appl. Math., 78 (2020), pp. 261–275

  4. [4]

    Buckmaster, P

    T. Buckmaster, P. Germain, Z. Hani, and J. Shatah , Onset of the wave turbulence description of the longtime behavior of the nonlinear Schr¨ odinger equation, Invent. Math., 225 (2021), pp. 787–855

  5. [5]

    Collot, H

    C. Collot, H. Dietert, and P. Germain, Stability and cascades for the Kolmogorov-Zakharov spectrum of wave turbulence , Arch. Ration. Mech. Anal., 248 (2024), pp. Paper No. 7, 31

  6. [6]

    Collot and P

    C. Collot and P. Germain , On the derivation of the homogeneous kinetic wave equation , preprint, arXiv:1912.10368, (2019)

  7. [7]

    Collot and P

    C. Collot and P. Germain , Derivation of the homogeneous kinetic wave equation: longer time scales, preprint, arXiv:2007.03508, (2020)

  8. [8]

    Deng and Z

    Y. Deng and Z. Hani , On the derivation of the wave kinetic equation for NLS , Forum Math. Pi, 9 (2021), pp. Paper No. e6, 37

  9. [9]

    Deng and Z

    Y. Deng and Z. Hani, Full derivation of the wave kinetic equation , Invent. Math., 233 (2023), pp. 543–724

  10. [10]

    Deng and Z

    Y. Deng and Z. Hani , Long time justification of wave turbulence theory , preprint, arXiv:2311.10082, (2024)

  11. [11]

    Dimarco and L

    G. Dimarco and L. Pareschi , Numerical methods for kinetic equations , Acta Numer., 23 (2014), pp. 369–520

  12. [12]

    Dyachenko, A

    S. Dyachenko, A. Newell, A. Pushkarev, and V. Zakharov , Optical turbulence: weak turbulence, condensates and collapsing filaments in the nonlinear schr¨ odinger equation , Physica D: Nonlinear Phenomena, 57 (1992), pp. 96–160

  13. [13]

    Escobedo and A

    M. Escobedo and A. Menegaki , Instability of singular equilibria of a wave kinetic equation , preprint, arXiv:2406.05280, (2024)

  14. [14]

    Escobedo and J

    M. Escobedo and J. J. L. Vel ´azquez, Finite time blow-up and condensation for the bosonic Nordheim equation, Invent. Math., 200 (2015), pp. 761–847

  15. [15]

    Escobedo and J

    M. Escobedo and J. J. L. Vel ´azquez, On the theory of weak turbulence for the nonlinear Schr¨ odinger equation, Mem. Amer. Math. Soc., 238 (2015), pp. v+107

  16. [16]

    Filbet, J

    F. Filbet, J. Hu, and S. Jin , A numerical scheme for the quantum Boltzmann equation with stiff collision terms, ESAIM: Math. Model. Numer. Anal. (M2AN), 46 (2012), pp. 443–463

  17. [17]

    Filbet, C

    F. Filbet, C. Mouhot, and L. Pareschi, Solving the Boltzmann equation in NlogN , SIAM J. Sci. Comput., 28 (2006), pp. 1029–1053

  18. [18]

    Filbet and G

    F. Filbet and G. Russo , High order numerical methods for the space non-homogeneous Boltzmann equation, J. Comput. Phys., 186 (2003), pp. 457–480

  19. [19]

    Gamba, J

    I. Gamba, J. Haack, C. Hauck, and J. Hu, A fast spectral method for the Boltzmann collision operator with general collision kernels , SIAM J. Sci. Comput., 39 (2017), pp. B658–B674

  20. [20]

    Germain, A

    P. Germain, A. D. Ionescu, and M.-B. Tran , Optimal local well-posedness theory for the kinetic wave equation , J. Funct. Anal., 279 (2020), pp. 108570, 28

  21. [21]

    Hasselmann, K

    S. Hasselmann, K. Hasselmann, J. Allender, and T. Barnett , Computations and parameterizations of the nonlinear energy transfer in a gravity-wave specturm. part ii: Parameterizations of the nonlinear energy transfer for application in wave models, Journal of Physical Oceanography, 15 (1985), pp. 1378–1391

  22. [22]

    Hu and K

    J. Hu and K. Qi , A fast Fourier spectral method for the homogeneous Boltzmann equation with non-cutoff collision kernels , J. Comput. Phys., 423 (2020), p. 109806

  23. [23]

    J. Hu, K. Qi, and T. Yang , A new stability and convergence proof of the Fourier-Galerkin spectral method for the spatially homogeneous Boltzmann equation, SIAM J. Numer. Anal., 59 (2021), pp. 613–633

  24. [24]

    Hu and L

    J. Hu and L. Ying , A fast spectral algorithm for the quantum Boltzmann collision operator , Commun. Math. Sci., 10 (2012), pp. 989–999

  25. [25]

    Hu and L

    J. Hu and L. Ying, A fast algorithm for the energy space boson Boltzmann collision operator , Math. Comp., 84 (2015), pp. 271–288

  26. [26]

    Janssen, The interaction of ocean waves and wind , Cambridge University Press, 2004

    P. Janssen, The interaction of ocean waves and wind , Cambridge University Press, 2004

  27. [27]

    Lukkarinen and H

    J. Lukkarinen and H. Spohn , Weakly nonlinear Schr¨ odinger equation with random initial data, Invent. Math., 183 (2011), pp. 79–188

  28. [28]

    Menegaki , L2-stability near equilibrium for the 4 waves kinetic equation , Kinet

    A. Menegaki , L2-stability near equilibrium for the 4 waves kinetic equation , Kinet. Relat. Models, 17 (2024), pp. 514–532

  29. [29]

    Mouhot and L

    C. Mouhot and L. Pareschi, Fast algorithms for computing the Boltzmann collision operator, Math. Comp., 75 (2006), pp. 1833–1852. This manuscript is for review purposes only. FAST FOURIER SPECTRAL METHOD FOR WKE 21

  30. [30]

    Nazarenko, Wave turbulence, vol

    S. Nazarenko, Wave turbulence, vol. 825, Springer, 2011

  31. [31]

    A. C. Newell and B. Rumpf , Wave turbulence, Annual review of fluid mechanics, 43 (2011), pp. 59–78

  32. [32]

    Pareschi and B

    L. Pareschi and B. Perthame , A Fourier spectral method for homogeneous Boltzmann equations, Transport Theory Statist. Phys., 25 (1996), pp. 369–382

  33. [33]

    Pareschi and G

    L. Pareschi and G. Russo, Numerical solution of the Boltzmann equation I: spectrally accurate approximation of the collision operator , SIAM J. Numer. Anal., 37 (2000), pp. 1217–1245

  34. [34]

    Pareschi, G

    L. Pareschi, G. Russo, and G. Toscani, Fast spectral methods for the Fokker-Planck-Landau collision operator, J. Comput. Phys., 165 (2000), pp. 216–236

  35. [35]

    Resio and W

    D. Resio and W. Perrie , A numerical study of nonlinear energy fluxes due to wave-wave interactions part 1. methodology and basic results , J. Fluid Mech., 223 (1991), pp. 603– 629

  36. [36]

    B. V. Semisalov, S. B. Medvedev, S. V. Nazarenko, and M. P. Fedoruk , Algorithm for solving the four-wave kinetic equation in problems of wave turbulence , Comput. Math. Math. Phys., 64 (2024), pp. 340–361

  37. [37]

    Staffilani and M.-B

    G. Staffilani and M.-B. Tran, Condensation and non-condensation times for 4-wave kinetic equations, preprint, arXiv:2407.18533, (2024)

  38. [38]

    Staffilani and M.-B

    G. Staffilani and M.-B. Tran, On the energy transfer towards large values of wavenumbers for solutions of 4-wave kinetic equations , preprint, arXiv:2407.18508, (2024)

  39. [39]

    B. A. Tracy and D. T. Resio, Theory and calculation of the nonlinear energy transfer between sea waves in deep water , in Proceedings of the Army Numerical and Computers Analysis Conference, US Army Research Office., 1982, p. 457

  40. [40]

    Walton and M.-B

    S. Walton and M.-B. Tran , A numerical scheme for wave turbulence: 3-wave kinetic equations, SIAM J. Sci. Comput., 45 (2023), pp. B467–B492

  41. [41]

    Walton, M.-B

    S. Walton, M.-B. Tran, and A. Bensoussan , A deep learning approximation of non- stationary solutions to wave kinetic equations , Appl. Numer. Math., 199 (2024), pp. 213– 226

  42. [42]

    D. J. Webb, Non-linear transfers between sea waves , Deep Sea Research, 25 (1978), pp. 279– 298

  43. [43]

    Womersley , Symmetric Spherical Designs on the sphere S2 with good geometric properties, The University of New South Wales

    R. Womersley , Symmetric Spherical Designs on the sphere S2 with good geometric properties, The University of New South Wales. , http://web.maths.unsw.edu.au/ ∼rsw/ Sphere/EffSphDes/ss.html (accessed Sep. 27, 2016)

  44. [44]

    Zakharov and N

    V. Zakharov and N. Filonenko , Weak turbulence of capillary waves , Journal of applied mechanics and technical physics, 8 (1967), pp. 37–40

  45. [45]

    V. E. Zakharov, V. S. L’vov, and G. Falkovich, Kolmogorov spectra of turbulence I: Wave turbulence, Springer Science & Business Media, 2012. This manuscript is for review purposes only