pith. sign in

arxiv: 1906.08660 · v1 · pith:KJAVNBUDnew · submitted 2019-06-20 · 🌌 astro-ph.SR · astro-ph.EP

Focusing of nonlinear eccentric waves in astrophysical discs

Elliot M. Lynch , Gordon I. Ogilvie This is my paper

Pith reviewed 2026-05-25 19:13 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.EP
keywords eccentric wavesastrophysical discsnonlinear wavesshort-wavelength limitaveraged LagrangianWKB approximationwave steepeningradial structure
0
0 comments X

The pith

Nonlinear eccentric waves in astrophysical discs steepen in nonlinearity and eccentricity but remain bounded by linear WKB solutions.

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

The paper develops a fully nonlinear approximation for the short-wavelength limit of eccentric waves in discs, using the averaged Lagrangian method to handle cases where small eccentricities still produce strong pressure gradients. It derives explicit conditions under which nonlinearity and eccentricity increase as the waves move through a radially structured background. The central result is that the nonlinear solution stays bounded by the corresponding linear WKB solution. This approach matters because it gives a tractable way to track highly nonlinear wave evolution without solving the complete hydrodynamic equations at every scale.

Core claim

We develop a fully nonlinear approximation to the short-wavelength limit of eccentric waves in astrophysical discs, based on the averaged Lagrangian method of Whitham (1965). In this limit there is a separation of scales between the rapidly varying eccentric wave and the background disc. Despite having small eccentricities, such rapidly varying waves can be highly nonlinear, potentially approaching orbital intersection, and this can result in strong pressure gradients in the disc. We derive conditions for the steepening of nonlinearity and eccentricity as the waves propagate in a radially structured disc in this short-wavelength limit and show that the behaviour of the solution can be bounde

What carries the argument

Averaged Lagrangian method, which averages over the rapid oscillations of the eccentric wave while treating the disc background as slowly varying to permit a fully nonlinear treatment.

If this is right

  • Waves can develop strong pressure gradients as they approach orbital intersection under the derived steepening conditions.
  • The nonlinear solution is always bounded above by the linear WKB solution in both nonlinearity and eccentricity.
  • Wave propagation and focusing in radially structured discs can be tracked using this scale-separated approximation.
  • The method applies directly when the background disc varies slowly compared with the wave.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The bounding result suggests linear WKB calculations can serve as a conservative upper limit for estimating nonlinear effects in disc models.
  • Similar averaged-Lagrangian treatments might extend to other short-wavelength waves in discs if scale separation is present.
  • The approach could help predict where pressure gradients become dynamically important without full nonlinear simulations.

Load-bearing premise

The short-wavelength limit permits a clean separation of scales between the rapidly varying eccentric wave and the slowly varying background disc.

What would settle it

A numerical hydrodynamical simulation of an eccentric wave propagating in a disc with known radial structure in which the nonlinear eccentricity exceeds the amplitude predicted by the linear WKB solution at the same location.

Figures

Figures reproduced from arXiv: 1906.08660 by Elliot M. Lynch, Gordon I. Ogilvie.

Figure 1
Figure 1. Figure 1: Geometric part of the Hamiltonian in the short￾wavelength limit. In each case the unimportant constant term has been subtracted off. Most forms of F of interest are brack￾eted by the behaviour of these functions. In particular we expect all should have the property of being monotonically increasing even functions which diverge at q = 1. The dashed lines are series expansion in q including terms up to O(q 2… view at source ↗
Figure 2
Figure 2. Figure 2: Phase space of the small scale nonlinear oscillator for F = F (iso) showing the Hamiltonian in terms of q and πq/k. and 2D isothermal discs, as the breathing mode is indepen￾dent of q owing to the vertical specific enthalpy gradient be￾ing independent of horizontal compression. The 3D isother￾mal disc will have an O(1) contribution to N 2 . It is useful to reformulate the Lagrangian viewing q as the coordi… view at source ↗
Figure 3
Figure 3. Figure 3: Variation of the eccentric mode on the short length￾scale (ϕ) for different maximum nonlinearities (q+) in an isother￾mal disc. Top (rescaled) eccentricity, bottom nonlinearity q. The eccentricity tends to a triangle wave and the nonlinearity to a square wave in the nonlinear limit q+ → 1. The averaged Lagrangian is given by hLi := 1 2π Z 2π 0 L dϕ = Z H ◦ ahLiˆ da , (37) where here the angle brackets deno… view at source ↗
Figure 4
Figure 4. Figure 4: shows how q+ varies on the disc length scale dependent on λ. As expected the solution is a power law at low q+ where linear theory is valid. As the nonlinearity in￾creases the solution turns over due to the increasingly strong pressure forces which are able to more effectively balance the precessional forces in the nonlinear regime so that q+ need not be as high as would be predicted in linear theory. The … view at source ↗
Figure 5
Figure 5. Figure 5: Eccentric amplitude on the disc scale for λ = −1, as the theory is linear in eccentricity, the absolute value of this is unimportant. Dashed lines correspond to power law solutions from linear theory. and µ. Unlike q the equations are linear in e and as such we are free to rescale our solutions by a constant value. As with q+ the solution is a power law in the linear regime and turns over at increasing non… view at source ↗
Figure 6
Figure 6. Figure 6: Same as [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: Same as [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
read the original abstract

We develop a fully nonlinear approximation to the short-wavelength limit of eccentric waves in astrophysical discs, based on the averaged Lagrangian method of Whitham (1965). In this limit there is a separation of scales between the rapidly varying eccentric wave and the background disc. Despite having small eccentricities, such rapidly varying waves can be highly nonlinear, potentially approaching orbital intersection, and this can result in strong pressure gradients in the disc. We derive conditions for the steepening of nonlinearity and eccentricity as the waves propagate in a radially structured disc in this short-wavelength limit and show that the behaviour of the solution can be bounded by the behaviour of the WKB solution to the linearised equations.

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 / 1 minor

Summary. The paper develops a fully nonlinear short-wavelength approximation for eccentric waves in astrophysical discs using Whitham's averaged Lagrangian method. In this limit, it derives conditions for the steepening of nonlinearity and eccentricity as waves propagate in a radially structured disc and shows that the nonlinear solution behavior can be bounded by the WKB solution to the linearised equations, despite small eccentricities allowing highly nonlinear regimes approaching orbital intersection.

Significance. If the central bounding result holds, this provides a valuable analytic framework for nonlinear wave dynamics in discs that extends beyond linear theory while retaining the short-wavelength separation. The application of the established Whitham averaged Lagrangian to this regime, yielding falsifiable steepening conditions without free parameters, is a strength that could inform models of disc evolution and wave focusing.

major comments (1)
  1. [Abstract] Abstract: the claim that 'the behaviour of the solution can be bounded by the behaviour of the WKB solution to the linearised equations' is load-bearing for the paper's main result, yet the scale-separation assumption underlying the averaged Lagrangian is not explicitly re-validated in the regime where eccentricity approaches orbital intersection and 'strong pressure gradients' develop (as noted in the abstract itself). If this assumption fails, both the steepening conditions and the bounding statement lose justification.
minor comments (1)
  1. [Abstract] The abstract mentions 'conditions for the steepening' but does not preview their explicit form or dependence on disc structure; adding a brief statement would improve clarity for readers.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their detailed review and for highlighting this key point about the foundational assumptions of our analysis. We address the major comment below and agree that explicit clarification is warranted.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that 'the behaviour of the solution can be bounded by the behaviour of the WKB solution to the linearised equations' is load-bearing for the paper's main result, yet the scale-separation assumption underlying the averaged Lagrangian is not explicitly re-validated in the regime where eccentricity approaches orbital intersection and 'strong pressure gradients' develop (as noted in the abstract itself). If this assumption fails, both the steepening conditions and the bounding statement lose justification.

    Authors: The short-wavelength assumption in the Whitham averaged Lagrangian requires only that the eccentric wave varies rapidly relative to the background disc structure (i.e., wavelength much smaller than radial scale height of the disc), independent of the local amplitude. Strong pressure gradients develop on the wave scale itself and are already incorporated into the nonlinear Lagrangian; they do not introduce new rapid variations on the background scale. The bounding by linear WKB solutions is derived within this fixed scale-separation framework and holds formally even as local eccentricity approaches orbital intersection. Nevertheless, we agree the manuscript would benefit from an explicit statement re-confirming the assumption's validity in the high-nonlinearity limit. We will add a short paragraph (likely in Section 2) deriving that the averaging procedure remains justified provided the wavelength condition is satisfied, with a brief estimate showing pressure-gradient length scales remain tied to the wave rather than the disc. revision: yes

Circularity Check

0 steps flagged

No circularity; applies external Whitham framework to new regime

full rationale

The derivation rests on the established averaged Lagrangian method of Whitham (1965), an external reference cited in the abstract. The paper derives steepening conditions and a bounding statement relative to linear WKB solutions within the short-wavelength scale-separation assumption. No quoted equations reduce a claimed prediction or result to a fitted parameter, self-definition, or self-citation chain by construction. The central claims remain independent of the inputs and are not forced by renaming or ansatz smuggling. This is the normal case of an honest non-finding.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption of scale separation in the short-wavelength limit that enables the averaged Lagrangian treatment; no free parameters or invented entities are mentioned.

axioms (1)
  • domain assumption Short-wavelength limit permits separation of scales between rapidly varying eccentric wave and background disc
    Stated explicitly in the abstract as the foundation for the nonlinear approximation.

pith-pipeline@v0.9.0 · 5638 in / 1094 out tokens · 28033 ms · 2026-05-25T19:13:25.392696+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

43 extracted references · 43 canonical work pages · 2 internal anchors

  1. [1]

    J., Ogilvie G

    Barker A. J., Ogilvie G. I., 2016, @doi [ ] 10.1093/mnras/stw580 , http://adsabs.harvard.edu/abs/2016MNRAS.458.3739B 458, 3739

    work page doi:10.1093/mnras/stw580 2016

  2. [2]

    W., Farihi J., Redfield S., Bachman S., Parsons S

    Cauley P. W., Farihi J., Redfield S., Bachman S., Parsons S. G., G \"a nsicke B. T., 2018, @doi [ ] 10.3847/2041-8213/aaa3d9 , http://adsabs.harvard.edu/abs/2018ApJ...852L..22C 852, L22

    work page doi:10.3847/2041-8213/aaa3d9 2018

  3. [3]

    H., Bate M

    D'Angelo G., Lubow S. H., Bate M. R., 2006, @doi [ ] 10.1086/508451 , http://adsabs.harvard.edu/abs/2006ApJ...652.1698D 652, 1698

    work page doi:10.1086/508451 2006

  4. [4]

    W., Latter H

    Dewberry J. W., Latter H. N., Ogilvie G. I., 2018, @doi [ ] 10.1093/mnras/sty385 , http://adsabs.harvard.edu/abs/2018MNRAS.476.4085D 476, 4085

    work page doi:10.1093/mnras/sty385 2018

  5. [5]

    W., Latter H

    Dewberry J. W., Latter H. N., Ogilvie G. I., 2019, @doi [ ] 10.1093/mnras/sty3061 , http://adsabs.harvard.edu/abs/2019MNRAS.483.1609D 483, 1609

    work page doi:10.1093/mnras/sty3061 2019

  6. [6]

    C., Chiang E., 2015, @doi [ ] 10.1088/0004-637X/812/2/94 , http://adsabs.harvard.edu/abs/2015ApJ...812...94D 812, 94

    Duffell P. C., Chiang E., 2015, @doi [ ] 10.1088/0004-637X/812/2/94 , http://adsabs.harvard.edu/abs/2015ApJ...812...94D 812, 94

    work page doi:10.1088/0004-637x/812/2/94 2015

  7. [7]

    Excitation of trapped oscillations in discs around black holes

    Ferreira B. T., Ogilvie G. I., 2008, in Axelsson M., ed., American Institute of Physics Conference Series Vol. 1054, American Institute of Physics Conference Series. pp 111--116 ( @eprint arXiv 0810.0116 ), @doi 10.1063/1.3002492

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1063/1.3002492 2008

  8. [8]

    T., Ogilvie G

    Ferreira B. T., Ogilvie G. I., 2009, @doi [ ] 10.1111/j.1365-2966.2008.14070.x , http://adsabs.harvard.edu/abs/2009MNRAS.392..428F 392, 428

    work page doi:10.1111/j.1365-2966.2008.14070.x 2009

  9. [9]

    Goldreich P., Sari R., 2003, @doi [ ] 10.1086/346202 , http://adsabs.harvard.edu/abs/2003ApJ...585.1024G 585, 1024

    work page doi:10.1086/346202 2003

  10. [10]

    Goodchild S., Ogilvie G., 2006, @doi [ ] 10.1111/j.1365-2966.2006.10197.x , http://adsabs.harvard.edu/abs/2006MNRAS.368.1123G 368, 1123

    work page doi:10.1111/j.1365-2966.2006.10197.x 2006

  11. [11]

    Kato S., 1990, , http://adsabs.harvard.edu/abs/1990PASJ...42...99K 42, 99

    work page 1990

  12. [12]

    Kato S., 2008, @doi [ ] 10.1093/pasj/60.1.111 , http://adsabs.harvard.edu/abs/2008PASJ...60..111K 60, 111

    work page doi:10.1093/pasj/60.1.111 2008

  13. [13]

    Kley W., Dirksen G., 2006, @doi [ ] 10.1051/0004-6361:20053914 , http://adsabs.harvard.edu/abs/2006A

    work page doi:10.1051/0004-6361:20053914 2006

  14. [14]

    Lee E., Goodman J., 1999, @doi [ ] 10.1046/j.1365-8711.1999.02764.x , http://adsabs.harvard.edu/abs/1999MNRAS.308..984L 308, 984

    work page doi:10.1046/j.1365-8711.1999.02764.x 1999

  15. [15]

    M., Lithwick Y., 2018, arXiv e-prints, http://adsabs.harvard.edu/abs/2018arXiv181111758L

    Lee W.-K., Dempsey A. M., Lithwick Y., 2018, arXiv e-prints, http://adsabs.harvard.edu/abs/2018arXiv181111758L

    work page 2018

  16. [16]

    M., Lithwick Y., 2019, arXiv e-prints, https://ui.adsabs.harvard.edu/abs/2019arXiv190605290L p

    Lee W.-K., Dempsey A. M., Lithwick Y., 2019, arXiv e-prints, https://ui.adsabs.harvard.edu/abs/2019arXiv190605290L p. arXiv:1906.05290

    work page arXiv 2019

  17. [17]

    H., 1991a, @doi [ ] 10.1086/170647 , http://adsabs.harvard.edu/abs/1991ApJ...381..259L 381, 259

    Lubow S. H., 1991a, @doi [ ] 10.1086/170647 , http://adsabs.harvard.edu/abs/1991ApJ...381..259L 381, 259

    work page doi:10.1086/170647

  18. [18]

    H., 1991b, @doi [ ] 10.1086/170648 , http://adsabs.harvard.edu/abs/1991ApJ...381..268L 381, 268

    Lubow S. H., 1991b, @doi [ ] 10.1086/170648 , http://adsabs.harvard.edu/abs/1991ApJ...381..268L 381, 268

    work page doi:10.1086/170648

  19. [19]

    I., Milosavljevi \'c M., 2008, @doi [ ] 10.1086/523869 , http://adsabs.harvard.edu/abs/2008ApJ...672...83M 672, 83

    MacFadyen A. I., Milosavljevi \'c M., 2008, @doi [ ] 10.1086/523869 , http://adsabs.harvard.edu/abs/2008ApJ...672...83M 672, 83

    work page doi:10.1086/523869 2008

  20. [20]

    J., et al., 2016, @doi [ ] 10.1093/mnras/stv2603 , http://adsabs.harvard.edu/abs/2016MNRAS.455.4467M 455, 4467

    Manser C. J., et al., 2016, @doi [ ] 10.1093/mnras/stv2603 , http://adsabs.harvard.edu/abs/2016MNRAS.455.4467M 455, 4467

    work page doi:10.1093/mnras/stv2603 2016

  21. [21]

    Fast and Slow Precession of Gaseous Debris Disks Around Planet-Accreting White Dwarfs

    Miranda R., Rafikov R. R., 2018, preprint, http://adsabs.harvard.edu/abs/2018arXiv180200464M ( @eprint arXiv 1802.00464 )

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  22. [22]

    J., Lai D., 2017, @doi [ ] 10.1093/mnras/stw3189 , http://adsabs.harvard.edu/abs/2017MNRAS.466.1170M 466, 1170

    Miranda R., Mu \ n oz D. J., Lai D., 2017, @doi [ ] 10.1093/mnras/stw3189 , http://adsabs.harvard.edu/abs/2017MNRAS.466.1170M 466, 1170

    work page doi:10.1093/mnras/stw3189 2017

  23. [23]

    I., 2001, @doi [ ] 10.1046/j.1365-8711.2001.04416.x , http://adsabs.harvard.edu/abs/2001MNRAS.325..231O 325, 231

    Ogilvie G. I., 2001, @doi [ ] 10.1046/j.1365-8711.2001.04416.x , http://adsabs.harvard.edu/abs/2001MNRAS.325..231O 325, 231

    work page doi:10.1046/j.1365-8711.2001.04416.x 2001

  24. [24]

    T., Ogilvie G

    Ogilvie G. I., 2008, @doi [ ] 10.1111/j.1365-2966.2008.13484.x , http://adsabs.harvard.edu/abs/2008MNRAS.388.1372O 388, 1372

    work page doi:10.1111/j.1365-2966.2008.13484.x 2008

  25. [25]

    I., Barker A

    Ogilvie G. I., Barker A. J., 2014, @doi [ ] 10.1093/mnras/stu1795 , http://adsabs.harvard.edu/abs/2014MNRAS.445.2621O 445, 2621

    work page doi:10.1093/mnras/stu1795 2014

  26. [26]

    I., Lynch E

    Ogilvie G. I., Lynch E. M., 2019, @doi [ ] 10.1093/mnras/sty3436 , https://ui.adsabs.harvard.edu/abs/2019MNRAS.483.4453O 483, 4453

    work page doi:10.1093/mnras/sty3436 2019

  27. [27]

    T., 1991, , http://adsabs.harvard.edu/abs/1991PASJ...43...75O 43, 75

    Okazaki A. T., 1991, , http://adsabs.harvard.edu/abs/1991PASJ...43...75O 43, 75

    work page 1991

  28. [28]

    T., 1997, , http://adsabs.harvard.edu/abs/1997A

    Okazaki A. T., 1997, , http://adsabs.harvard.edu/abs/1997A

    work page 1997

  29. [29]

    Papaloizou J. C. B., 2002, @doi [ ] 10.1051/0004-6361:20020490 , http://adsabs.harvard.edu/abs/2002A

    work page doi:10.1051/0004-6361:20020490 2002

  30. [30]

    Papaloizou J. C. B., Savonije G. J., 2006, @doi [ ] 10.1051/0004-6361:20065407 , http://adsabs.harvard.edu/abs/2006A

    work page doi:10.1051/0004-6361:20065407 2006

  31. [31]

    Papaloizou J. C. B., Nelson R. P., Masset F., 2001, @doi [ ] 10.1051/0004-6361:20000011 , http://adsabs.harvard.edu/abs/2001A

    work page doi:10.1051/0004-6361:20000011 2001

  32. [32]

    J., Lodato G., 2018, @doi [ ] 10.1093/mnras/stx3094 , http://adsabs.harvard.edu/abs/2018MNRAS.474.4460R 474, 4460

    Ragusa E., Rosotti G., Teyssandier J., Booth R., Clarke C. J., Lodato G., 2018, @doi [ ] 10.1093/mnras/stx3094 , http://adsabs.harvard.edu/abs/2018MNRAS.474.4460R 474, 4460

    work page doi:10.1093/mnras/stx3094 2018

  33. [33]

    Rice W. K. M., Armitage P. J., Hogg D. F., 2008, @doi [ ] 10.1111/j.1365-2966.2007.12817.x , http://adsabs.harvard.edu/abs/2008MNRAS.384.1242R 384, 1242

    work page doi:10.1111/j.1365-2966.2007.12817.x 2008

  34. [34]

    P., Booth R

    Rosotti G. P., Booth R. A., Clarke C. J., Teyssandier J., Facchini S., Mustill A. J., 2017, @doi [ ] 10.1093/mnrasl/slw184 , http://adsabs.harvard.edu/abs/2017MNRAS.464L.114R 464, L114

    work page doi:10.1093/mnrasl/slw184 2017

  35. [35]

    J., Heemskerk M

    Savonije G. J., Heemskerk M. H. M., 1993, , http://adsabs.harvard.edu/abs/1993A

    work page 1993

  36. [36]

    I., Sunyaev R

    Shakura N. I., Sunyaev R. A., 1973, , http://adsabs.harvard.edu/abs/1973A

    work page 1973

  37. [37]

    I., 2016, @doi [ ] 10.1093/mnras/stw521 , http://adsabs.harvard.edu/abs/2016MNRAS.458.3221T 458, 3221

    Teyssandier J., Ogilvie G. I., 2016, @doi [ ] 10.1093/mnras/stw521 , http://adsabs.harvard.edu/abs/2016MNRAS.458.3221T 458, 3221

    work page doi:10.1093/mnras/stw521 2016

  38. [38]

    I., 2017, @doi [ ] 10.1093/mnras/stx426 , http://adsabs.harvard.edu/abs/2017MNRAS.467.4577T 467, 4577

    Teyssandier J., Ogilvie G. I., 2017, @doi [ ] 10.1093/mnras/stx426 , http://adsabs.harvard.edu/abs/2017MNRAS.467.4577T 467, 4577

    work page doi:10.1093/mnras/stx426 2017

  39. [39]

    Thun D., Kley W., 2018, @doi [ ] 10.1051/0004-6361/201832804 , http://adsabs.harvard.edu/abs/2018A

    work page doi:10.1051/0004-6361/201832804 2018

  40. [40]

    Thun D., Kley W., Picogna G., 2017, @doi [ ] 10.1051/0004-6361/201730666 , http://adsabs.harvard.edu/abs/2017A

    work page doi:10.1051/0004-6361/201730666 2017

  41. [41]

    Tremaine S., 2001, @doi [ ] 10.1086/319398 , http://adsabs.harvard.edu/abs/2001AJ....121.1776T 121, 1776

    work page doi:10.1086/319398 2001

  42. [42]

    B., 1965, @doi [Journal of Fluid Mechanics] 10.1017/S0022112065000745 , http://adsabs.harvard.edu/abs/1965JFM....22..273W 22, 273

    Whitham G. B., 1965, @doi [Journal of Fluid Mechanics] 10.1017/S0022112065000745 , http://adsabs.harvard.edu/abs/1965JFM....22..273W 22, 273

    work page doi:10.1017/s0022112065000745 1965

  43. [43]

    write newline

    " write newline "" before.all 'output.state := FUNCTION fin.entry write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION new.sentence output.state after.block = 'skip output.state before.all = 'skip after.sentence 'output.state := if if FUNCTION not #0 #1 if FUNCTION and 'skip pop #0 if FUNCTION or pop #1...

    work page