Pith. sign in

REVIEW 1 major objections 3 minor 3 cited by

A New Approach to Direct Discretization of Wave Kinetic Equations with Application to a Nonlinear Schrodinger System in 2D

T0 review · 1 major / 3 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A new direct discretization of wave kinetic equations approximates resonant manifolds piecewise linearly and evaluates collision integrals by midpoint quadrature, achieving second-order accuracy.

desk verdict A promising general scheme for 3-wave kinetic equations, but the printed normal in Eq. (16b) is geometrically wrong and must be corrected before the method is reproducible. read the letter →

arxiv 2509.03432 v1 pith:O7N4JNTD submitted 2025-09-03 math.NA cs.NA

J. W. Banks , J. Shatah This is my paper · ORCID
classification math.NAcs.NA MSC 35Q5565M0665R20
keywords wavekineticequationresonantmanifoldcollisionintegralmidpointquadraturenonlinearSchrödingermethodoflinesself-convergenceturbulence
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper introduces a general numerical method for wave kinetic equations (WKEs), the statistical equations that describe how averaged wave spectra evolve through resonant wave interactions. The method approximates the resonant manifold piecewise linearly and then evaluates the collision integral by midpoint quadrature, avoiding analytic parametrization. For model collision integrals in 1D and 2D it demonstrates second-order accuracy; for the full 2D nonlinear Schrödinger WKE, self-convergence studies show near-second-order rates. It then compares WKE solutions with 180-member ensemble averages of the underlying NLS equations and finds agreement for both isotropic and anisotropic initial data.

What carries the argument

The resonant manifold carried by the collision integral's Dirac delta is the central object. The scheme replaces it with the zero level set of a cell-wise bilinear interpolant Lg, treats the manifold piecewise as line segments, and evaluates the integral as Lf(sample) times segment length divided by the interpolated normal derivative of g. This converts a singular delta integral into a sum of ordinary quadrature weights, which is what allows the WKE right-hand side to be evaluated without analytic parametrization.

What would settle it

Run the 2D model integrals (Tests I–V) on grids where resonant contours cross cells at 45 degrees, using exactly the printed Eq. (16b) normal: if max-norm errors stall or diverge instead of decreasing at second order, the printed scheme is not the one delivering the reported convergence. Running with a corrected normal should reproduce the reported rates.

Watch

Extended reading notes

Core claim

The central claim is that the right-hand side of a wave kinetic equation—a singular integral over a resonant manifold forced by Dirac deltas—can be discretized directly and generically by (i) bilinearly interpolating the resonance function g on each grid cell, (ii) tracing the zero-level segment of that interpolant as a piecewise-linear approximate manifold, and (iii) applying midpoint quadrature with a bilinearly interpolated numerator and gradient. In the 2D setting this yields Eq. (17), the midpoint approximation that defines the scheme. The paper demonstrates the scheme on six model collision integrals, obtains exact second-order convergence, and then embeds it in a method-of-lines time

Load-bearing premise

The reported second-order accuracy rests on the denominator in the quadrature weight being the true normal derivative of the resonance function at the manifold segment; as printed, the normal in Eq. (16b) is actually a tangent for generic segments, so unless the implementation uses a corrected normal the accuracy claim does not follow from the equations in the paper.

Editorial extensions

If this is right

  • Collision integrals can be evaluated at second order on uniform Cartesian grids without analytic parametrization of resonance conditions.
  • Because the method only needs grid evaluations of the resonance function, its gradient, and the integrand, it extends naturally to other WKEs with quadratic nonlinearity.
  • The WKE and ensemble-averaged NLS solutions agree at kinetic time τ=1 for both isotropic and anisotropic top-hat initial data, supporting the kinetic description for this model.
  • Smooth initial data yield near-second-order self-convergence rates; non-smooth top-hat data remain close to second order.
  • Higher-order quadrature, implicit time stepping, and application to capillary or gravity wave models are natural extensions within the same framework.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same piecewise-linear manifold reconstruction could be combined with adaptive mesh refinement near resonance curves, an option the paper does not explore.
  • The printed normal formula in Eq. (16b) is not a normal for segments at 45 degrees to the grid axes, so reproducing the reported second-order rates requires an unstated correction to that formula in the implementation.
  • The method's grid-based generality suggests it could be applied to WKEs whose resonant manifolds are not circular and for which analytic parametrizations are harder to construct.
  • The agreement seen with 180 ensemble members suggests that finite-size and discretization effects have fallen below the kinetic description's error at L=64 with 1001×1001 spatial points; larger ensembles would sharpen this conclusion.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

1 major / 3 minor

Summary. The manuscript proposes a general numerical method for direct discretization of the right-hand side of wave kinetic equations (WKEs). The collision integral is evaluated by approximating the resonant manifold as a piecewise linear curve: the resonance function g is interpolated bilinearly on each cell, the zero-level segments are located, and the integral over each segment is approximated by midpoint quadrature using the segment length and the projected gradient of g. The method is developed for a 2D nonlinear Schrödinger system and its WKE (3). In 1D and 2D model problems with exact solutions, the scheme is reported to converge at second order (Figures 2, 5, 7). For the full WKE, Richardson self-convergence studies give rates near 2 for smooth Gaussian data and roughly 1.6 for top-hat data (Table 2). Finally, WKE solutions are compared with 180-member ensemble averages of the underlying NLS dynamics at kinetic time τ=1, with good agreement for sufficiently large domains and spatial resolutions.

Significance. If the method is correct, it is a useful contribution: it avoids analytic parametrization of the resonant manifold and special polar grids, and it directly discretizes the collision integral in a way that is applicable to general 3-wave WKEs. The paper supports the claim with exact-solution model problems, self-convergence studies, and an independent comparison against ensemble NLS simulations, which is a substantial validation effort. The principal reservation is that the printed definition of the normal vector in Eq. (16b) is not a normal, so the scheme actually implemented and validated cannot be the scheme described by Eqs. (16)-(17). Since Eq. (17) is explicitly identified as the defining second-order midpoint approximation, this is a load-bearing error that must be fixed before the manuscript is acceptable.

major comments (1)
  1. [Eq. (16b), Eq. (17)] The vector n defined in Eq. (16b) is not a unit normal to the segment from (ξ_a,η_a) to (ξ_b,η_b). A true normal is, e.g., (-(η_b-η_a), ξ_b-ξ_a)/ΔK. With the printed definition, n·v = 2(ξ_b-ξ_a)(η_b-η_a)/ΔK for the segment vector v = (ξ_b-ξ_a, η_b-η_a). This dot product is nonzero generically and is largest for segments at 45° to the axes, where the printed n is parallel to the segment. For a chord of a circular resonant manifold at 45°, n·∇g ≈ 0, so the denominator |L∇_n g| in Eq. (17) vanishes or is spuriously small. The reported second-order convergence in Tests I–V and in the full WKE runs therefore cannot have been produced by Eq. (16b) as printed: either the implementation used a true normal and the paper contains an unstated correction to its defining formula, or the printed scheme would have failed at generic 45° chords. This is not cosmetic, because Eq. (17) is called 'the secon
minor comments (3)
  1. [Eq. (5), Eq. (7b)] The interpolation denominators are written as Δk, but the interpolation is in the variable ξ over [ξ_J, ξ_{J+1}]. Unless the ξ-grid spacing is intentionally identical to the k-grid spacing, the denominator should be Δξ = ξ_{J+1}-ξ_J. Please clarify the notation.
  2. [Figure 5 caption/legend] The legend appears to list 'Test II' twice; presumably the second occurrence should be 'Test III'. This makes the convergence plot harder to read.
  3. [Section 4.4] The NLS ensemble comparisons are presented and discussed only through visual inspection of slices and 2D colormaps. Since the claim is 'reasonable agreement', a quantitative measure (e.g., relative error in a chosen norm as a function of L and M) would strengthen the conclusion, though this is not required for the main numerical-method claim.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the accuracy claims are validated against exact analytic solutions and direct NLS ensemble averages, not against the scheme's own inputs.

full rationale

The derivation chain is self-contained and externally anchored. The WKE (3) is derived from the NLS model (2) in Appendix A using random-phase assumptions; this derivation is independent of the numerical scheme. The numerical method in Section 3 is defined by piecewise-linear/bilinear manifold approximation and midpoint quadrature (Eqs. 5-17), and its claimed second-order accuracy is tested against exact analytic solutions: Eq. (10e) for the 1D models, Table 1 plus exact results for the 2D tests, and Eq. (19f) for the general 2D RHS model. No parameter is fitted to any of these exact answers: the scheme's quadrature coefficients are fixed, and the only tuned constants (Lambda_S=2.25 in Eq. 24a and P=3/4 in Eq. 27a) are explicit experimental choices that do not encode the target outputs. The full WKE evolution is assessed by Richardson self-convergence (Section 4.3) and, crucially, by comparison with 180-member ensemble averages of direct NLS simulations (Section 4.4), an independent benchmark outside the fitted or derived quantities. Self-citations appear only as methodological support: [16] for the kinetic-limit ensemble procedure, and [27,28] for the Richardson extrapolation technique; neither is load-bearing for the mathematical correctness of the new discretization. The possible inconsistency in Eq. (16b) regarding the normal vector is a numerical correctness/implementation concern, not a circularity of the kind where an output is equivalent to an input by construction. Thus there is no exhibitable circular step, and the appropriate score is 0.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new physical entities and fits no parameters to its validation targets. The scheme's free parameters are experimental settings (safety factor, initial-data shapes, limit exponent). The load-bearing domain assumptions are the random-phase closure and the kinetic limit, both standard in wave turbulence, plus smoothness and truncation assumptions on the numerical side.

free parameters (4)
  • Time-step safety factor Lambda_S = 2.25
    Hand-chosen stability margin in Delta_tau = Lambda_S h / lambda_h (Eq. 24a); the paper states the precise value has little effect.
  • Kinetic-limit exponent P = 3/4
    Sets epsilon = L^(-P) and final time tF = L^(2P) for the NLS ensembles (Eq. 27); follows the setup of reference [16].
  • Top-hat radius r0 = 2*pi
    Initial condition scale for WKE and NLS runs (Eq. 25b); the text says it is chosen without real intent.
  • Gaussian amplitude A = 5
    Initial-condition amplitude for the Gaussian WKE runs (Eq. 25a); the text says it is chosen without real intent.
assumptions (5)
  • domain assumption Random phase approximation: <a_k>=0, <a_k a-bar_k>=delta_{k,k'} n_k, all odd and mixed cubic moments vanish (Eq. A.14)
    Standard wave-turbulence closure; enters the WKE derivation in Appendix A and is the statistical assumption behind the NLS-to-WKE comparison.
  • domain assumption Kinetic limit: L to infinity with epsilon to 0 and tau = t epsilon^2 = O(1) (Eqs. A.32, A.33)
    The WKE is valid in this distinguished limit; used in Section 4.4 to justify the finite-box NLS ensemble protocol.
  • standard math Dispersion positivity omega_k = 1 + ||k||^2 > 0 forces the 'funny' matching terms delta(k) delta(Omega) to vanish (Eq. A.32)
    Terms involving omega_0 = 1 or repeated frequencies cannot satisfy the resonance deltas, so the extra terms drop out.
  • domain assumption The spectrum n_k is smooth enough that bilinear interpolation of g, f, and the normal derivative is O(h^2) accurate at cell interiors (Section 3.2)
    Underlies the claimed second-order convergence; no error analysis is given for degenerate cells or manifolds tangent to the grid.
  • domain assumption Finite-domain truncation k in [-20,20]^2 with negligible boundary effects (Section 4.1)
    The WKE and its collision integral are restricted to a finite box; the paper checks boundary behavior only on a 2D model problem (Figure 6).

how reviews work

0 comments
Cite this review

Pith. "Pith review of A New Approach to Direct Discretization of Wave Kinetic Equations with Application to a Nonlinear Schrodinger System in 2D." pith.science (2026). https://pith.science/paper/O7N4JNTD

@misc{pith2026250903432,
  author       = {Pith},
  title        = {Pith review of: A New Approach to Direct Discretization of Wave Kinetic Equations with Application to a Nonlinear Schrodinger System in 2D},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/O7N4JNTD}},
  note         = {Machine review of arXiv:2509.03432}
}
read the original abstract

Wave Kinetic Equations (WKEs) are often used to describe the evolution of ensemble averaged wave amplitudes for nonlinear wave systems. In the present manuscript we describe a new approach to direct numerical simulation of solutions to WKEs. This new method relies on a piecewise polynomial approximation of the resonant manifold, followed by numerical quadrature of the collision integral. The approach is general in nature, and is discussed in detail here for a particular nonlinear Schrodinger model in 2 spatial dimensions. Detailed convergence studies demonstrate 2nd-order accuracy for model collision integrals, and self-convergence studies for the WKE show near 2nd-order rates. Furthermore, comparison of the WKE approximation to ensemble averages of the NLS illustrate the efficacy of the method and the validity of the WKE, for both isotropic and an-isotropic solutions.

Figures

Figures reproduced from arXiv: 2509.03432 by the authors.

Figure 1
Figure 1. Plot of exact soln and approximation (left), and the error in the approximation (right) for Test i. [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Max-norm convergence results showing error vs. grid size for Tests i and ii. Clear 2nd order convergence [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Plot of exact soln and approximation (left) and the error in the approximation (right) for Test ii. The [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (12 more)
Figure 4
Figure 4. Figure 4: Plot of approximate resonant manifolds for 2D Test I, V, and VI. Note that for case VI, the resonant [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: Convergence vs. grid size for 2D tests independent of ( [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]
Figure 6
Figure 6. Figure 6: Plot of computed solution (left), and error (right) for 2D model of Section 3.2.2 with ( [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]
Figure 7
Figure 7. Figure 7: Convergence vs. grid size for 2D model of Section 3.2.2 with ( [PITH_FULL_IMAGE:figures/full_fig_p012_7.png]
Figure 8
Figure 8. Figure 8: Plot of wk = log10(nk) at τ = 1 for initial Gaussian initial condition with A = 5 (left), and top-hat initial condition with r0 = 2π (right). Both cases are anisotropic and use k (x) 0 = 5 and k (y) 0 = 3. 4.3. Self-Convergence Studies Given initial conditions of the f…
Figure 9
Figure 9. Figure 9: Plot of estimated max-norm error in wk at τ = 1 for the anisotropic Gaussian (left), and anisotropic top-hat (right). In both cases, the maximum error occurs in the domain interior, and the error appears to be well-behaved throughout the domain. For the top-hat, one ca…
Figure 10
Figure 10. Figure 10: At left is the solution of the WKE at τ = 1. At center is a 180 member ensemble for the NLS with L = 32 and Mx = My = 301. At right is a 180 member ensemble for the NLS with L = 64 and Mx = My = 1001. All results are plotted using log10, which is also indicated by the…
Figure 11
Figure 11. Figure 11: Line slices along k (y) = 0 comparing the WKE solution to NLS ensemble averages. At left the NLS runs use L = 32 and Mx = My = 301. At center the NLS runs use L = 64 and Mx = My = 601. At right the NLS runs use L = 64 and Mx = My = 1001. the second, the length of the …
Figure 12
Figure 12. Figure 12: Comparison of the WKE solution for k (y) = 0 to the NLS solutions plotted vs. radius k for every k-space point. At left the NLS runs use L = 32 and Mx = My = 301. At center the NLS runs use L = 64 and Mx = My = 601. At right the NLS runs use L = 64 and Mx = My = 1001.…
Figure 13
Figure 13. Figure 13: At left is the solution of the WKE at τ = 1 for the anisotropic top-hat. At center is a 180 member ensemble for the NLS with L = 32 and Mx = My = 301. At right is a 180 member ensemble for the NLS with L = 64 and Mx = My = 1001. All results are plotted using log10, wh…
Figure 14
Figure 14. Figure 14: Line slices along k (y) = 0 comparing the WKE solution to NLS ensemble averages. The general agreement is quite good for wave numbers in roughly k (x) ∈ [−8, 12], with improvement for both increasing L and x-grid resolution. -20 -10 0 10 20 k 10-15 10-10 10-5 100 |n| …
Figure 15
Figure 15. Figure 15: Line slices along k (y) = 2 comparing the WKE solution to NLS ensemble averages. The general agreement remains quite good for wave numbers in roughly k (x) ∈ [−8, 12], with improvement for both increasing L and x-grid resolution. 5. Conclusions In the present manuscri…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Entropy Structures and Long-Time Relaxation for 3-Wave Kinetic Equations

    math.AP 2026-05 unverdicted novelty 8.0 of 10

    New entropy structures from one-sided balance conditions on interaction weights yield global weak L1_loc solutions to 3-wave kinetic equations and prove their local relaxation to zero equilibrium.

  2. Spectral Algorithms for 3-Wave Kinetic and $C_{12}$ Quantum Boltzmann Equations with General Resonance Manifolds in $\mathbb{R}^d$

    math.NA 2026-08 conditional novelty 6.0 of 10

    A fast FFT-based spectral method computes 3-wave kinetic and C12 quantum Boltzmann collision operators in O(L(2N)^{2d} log(2N)) operations and shows apparent energy cascades in non-radial 3D simulations.

  3. Global time-analytic strong solutions for a class of 3-wave kinetic equations

    math-ph 2026-07 accept novelty 6.0 of 10

    For radial nonnegative data, a regularized 3-wave kinetic equation has a unique global strong solution that stays nonnegative and is time-analytic.

Reference graph

Works this paper leans on

36 extracted references · 35 canonical work pages · cited by 3 Pith papers

  1. [1]

    M. J. Ablowitz, Nonlinear Dispersive Waves: Asymptotic Analysis and Solitons, Cambridge Texts in Applied Mathematics, Cambridge University Press, 2011

    work page 2011

  2. [2]

    G. B. Whitham, Linear and Nonlinear Waves, Wiley-Interscience, New York, 1974

    work page 1974

  3. [3]

    Lannes, The Water Waves Problem: Mathematical Analysis and Asymptotics, American Mathematical Society, Philadelphia, 2013

    D. Lannes, The Water Waves Problem: Mathematical Analysis and Asymptotics, American Mathematical Society, Philadelphia, 2013

    work page 2013

  4. [4]

    Galtier, Physics of Wave Turbulence, Cambridge University Press, Cambridge, 2022

    S. Galtier, Physics of Wave Turbulence, Cambridge University Press, Cambridge, 2022

    work page 2022

  5. [5]

    Agrawal, Fiber-Optic Communication Systems, Wiley, Hoboken, 2019

    G. Agrawal, Fiber-Optic Communication Systems, Wiley, Hoboken, 2019

    work page 2019

  6. [6]

    Drazin, R

    P. Drazin, R. Johnson, Solitons: An Introduction, Cambridge University Press, Cambridge, 1989

    work page 1989

  7. [7]

    G. J. Komen, L. Cavaleri, M. Donelan, K. Hasselmann, S. Hasselmann, P. A. E. M. Janssen, Dynamics and modelling of ocean waves, Vol. 532, Cambridge University Press, Cambridge, 1994

    work page 1994

  8. [8]

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

    work page 2004

Show all 36 references
  1. [9]

    P. A. E. M. Janssen, M. Onorato, The intermediate water depth limit of the zakharov equation and consequences for wave prediction, Journal of Physical Oceanography 37 (10) (2007) 2389– 2400

    2007

  2. [10]

    Bodineau, I

    T. Bodineau, I. Gallagher, L. Saint-Raymond, S. Simonella, Long-time derivation at equilib- rium of the fluctuating boltzmann equation, Ann. Probab. 52 (1) (2024) 217–295

    2024

  3. [11]

    L. Boltzmann, Weitere studien ¨ uber das warme gleichgenicht unter gasmolak¨ ulen (further studies on the thermal equilibrium of gas molecules), Sitzungsberichte Akademie der Wis- senschaften 66 (1872) 275–370

  4. [12]

    Y. Deng, Z. Hani, X. Ma, Long time derivation of the boltzmann equation from hard sphere dynamics (2025)

    2025

  5. [13]

    Nazarenko, Wave Turbulence (Lecture Notes in Physics, Vol

    S. Nazarenko, Wave Turbulence (Lecture Notes in Physics, Vol. 825), Springer, Heidelberg, 2011

    2011

  6. [14]

    V. E. Zakharov, V. S. Lvov, G. Falkovich, Kolmogorov Spectra of Turbulence I: Wave Turbu- lence, Springer, Heidelberg, 1992

    1992

  7. [15]

    Buckmaster, P

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

    2021

  8. [16]

    J. W. Banks, T. Buckmaster, A. Korotkevich, G. Kovacic, 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) 034101

    2022

  9. [17]

    Y. Deng, Z. Hani, Full derivation of the wave kinetic equation, Inventiones mathematicae 233 (2) (2023) 543–724. 20

    2023

  10. [18]

    Pan, Understanding of weak turbulence of capillary waves, Ph.D

    Y. Pan, Understanding of weak turbulence of capillary waves, Ph.D. thesis, Massachusetts Institute of Technology (2017)

    2017

  11. [19]

    Webb, Non-linear transfers between sea waves, Deep Sea Research 25 (3) (1978) 279–298

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

    1978

  12. [20]

    Resio, W

    D. Resio, 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) 603–629

    1991

  13. [21]

    Krstulovic, V

    G. Krstulovic, V. Labarre, Wavkins.jl : an efficient and modular julia software for solving wave kinetic equations (2025). URL https://arxiv.org/abs/2504.00252

    2025 arXiv

  14. [22]

    Walton, M.-B

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

    2023

  15. [23]

    Walton, M.-B

    S. Walton, M.-B. Tran, Numerical schemes for 3-wave kinetic equations: A complete treatment of the collision operator, J. Comput. Phys. 538 (2025) 114147

    2025

  16. [24]

    Engquist, A.-K

    B. Engquist, A.-K. Tornberg, R. Tsai, Discretization of dirac delta functions in level set meth- ods, J. Comput. Phys. 207 (1) (2005) 28–51

    2005

  17. [25]

    U. M. Ascher, L. R. Petzold, Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations, SIAM, Philadelphia, 1998

    1998

  18. [26]

    W. D. Henshaw, D. W. Schwendeman, Parallel computation of three-dimensional flows using overlapping grids with adaptive mesh refinement, J. Comput. Phys. 227 (16) (2008) 7469–7502

    2008

  19. [27]

    J. W. Banks, W. D. Henshaw, J. N. Shadid, An evaluation of the FCT method for high-speed flows on structured overlapping grids, J. Comput. Phys. 228 (15) (2009) 5349–5369

    2009

  20. [28]

    J. W. Banks, T. D. Aslam, Richardson extrapolation for linearly degenerate discontinuities, J. Sci. Comput. 57 (2013) 1–18

    2013

  21. [29]

    J. W. Banks, W. D. Henshaw, On the coefficients in finite difference series expansions of derivatives, SIAM J. Numer. Anal

  22. [30]

    Fornberg, A Practical Guide to Pseudospectral Methods, Cambridge University Press, Cambridge, 1996

    B. Fornberg, A Practical Guide to Pseudospectral Methods, Cambridge University Press, Cambridge, 1996

    1996

  23. [31]

    Fornberg, Classroom note: Calculation of weights in finite difference formulas, SIAM Rev

    B. Fornberg, Classroom note: Calculation of weights in finite difference formulas, SIAM Rev. 40 (3) (1998) 685–691. Appendix A. Derivation of WKE for NLS model Here we present a derivation of the WKE (3) from the governing equations of motion (2), which are repeated here as ι∂...

    1998

  24. [33]

    + 2α0 1 2 α4 2 3 (s + ιΩ0 12)(s + ιΩ0 12 + Ω42 3 ) − 2(ιµ)2 X 1,2,3,4 α2 0 1 α34 0 (s + ιΩ20 1 )(s + ιΩ20 1 + ιΩ34 0 ) + 2α2 0 1 α3 4 0 (s + ιΩ20 1 )(s + ιΩ20 1 + Ω3 40) + 2(ιµ)2 X 1,2,3,4 α20 1 α 1 34 (s + ιΩ20 1 )(s + ιΩ20 1 + ιΩ1

  25. [34]

    + 2α0 1 2α4 1 3 (s + ιΩ0 12)(s + ιΩ0 12 + ιΩ41 3 ) + (ιµ)2 X 1,2,3,4 α0 1 2 α 2 34 (s + ιΩ0 12)(s + ιΩ0 12 + Ω2

  26. [35]

    dt+ c.c

    + 2α20 1 α4 1 3 (s + ιΩ20 1 )(s + ιΩ20 1 + ιΩ41 3 ) − 2(ιµ)2 X 1,2,3,4 α 2 0 1 α34 2 (s + ιΩ20 1 )(s + ιΩ20 1 + Ω34 2 ) + 2α 2 0 1 α3 4 2 (s + ιΩ20 1 )(s + ιΩ20 1 + Ω3 42) + Z ∞ 0 . . . dt+ c.c.. (A.12) By symmetry the 2nd and 3rd terms can be combined to yield sL[|ak|2](s) − ...

  27. [36]

    + 2α0 1 2α4 1 3 (s + ιΩ0 12)(s + ιΩ0 12 + ιΩ41 3 ) + 2µ2 X 1,2,3,4 α2 0 1 α34 0 (s + ιΩ20 1 )(s + ιΩ20 1 + ιΩ34 0 ) + 2α2 0 1 α3 4 0 (s + ιΩ20 1 )(s + ιΩ20 1 + Ω3 40) − 2µ2 X 1,2,3,4 α20 1 α 1 34 (s + ιΩ20 1 )(s + ιΩ20 1 + ιΩ1

  28. [37]

    dt+ c.c

    + 2α20 1 α4 1 3 (s + ιΩ20 1 )(s + ιΩ20 1 + ιΩ41 3 ) + 2µ2 X 1,2,3,4 α 2 0 1 α34 2 (s + ιΩ20 1 )(s + ιΩ20 1 + Ω34 2 ) + 2α 2 0 1 α3 4 2 (s + ιΩ20 1 )(s + ιΩ20 1 + Ω3 42) + Z ∞ 0 . . . dt+ c.c.. (A.13) As indicated, the idea will now be to take the expectation over assumed rando...

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.