pith. sign in

arxiv: 2607.01207 · v1 · pith:VYVA6RX2new · submitted 2026-07-01 · ⚛️ physics.flu-dyn · astro-ph.CO· astro-ph.GA

No evidence of vorticity production from initially irrotational turbulent gravitational collapse

Axel Brandenburg , Evangelia Ntormousi , Jennifer Schober This is my paper

Pith reviewed 2026-07-02 05:04 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn astro-ph.COastro-ph.GA
keywords gravitational collapsevorticity productionturbulencedirect numerical simulationbarotropic equation of stateirrotational flowfluid dynamicsviscous effects
0
0 comments X

The pith

Gravitational collapse produces no vorticity from initially irrotational turbulence in direct simulations.

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

The paper tests whether gravitational collapse can generate vortical motions when the initial turbulent flow has no vorticity. Direct numerical simulations are performed with a barotropic equation of state and no magnetic fields, so that the only possible source of vorticity is viscosity. The simulations show that observed vorticity traces entirely to the initial conditions and does not increase due to the collapse flow itself. This matters for understanding whether collapse can seed the vortical turbulence required for small-scale dynamo action that would amplify magnetic fields during events such as star formation.

Core claim

In direct numerical simulations of gravitational collapse with an initially irrotational turbulent velocity field, a barotropic equation of state, and no magnetic fields, vorticity production occurs only through the initial turbulence and viscous effects, with no measurable contribution from the collapse flow within the parameter space accessible to the numerical resolution.

What carries the argument

Direct numerical simulations that enforce parallel pressure and density gradients via the barotropic equation of state, thereby restricting vorticity sources exclusively to viscosity and allowing isolation of any collapse-induced generation.

If this is right

  • Vorticity in collapsing regions originates from pre-collapse turbulence rather than the collapse dynamics.
  • Small-scale dynamo action during collapse requires initial vortical motions and is not seeded by the collapse itself.
  • Barotropic models of collapse do not artificially generate vorticity through the flow geometry.
  • Turbulence remains predominantly irrotational throughout the collapse when started from irrotational conditions.

Where Pith is reading between the lines

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

  • Realistic collapses may need non-barotropic thermodynamics or magnetic fields to produce additional vorticity sources not captured here.
  • Focus on setting realistic initial vorticity levels may be more important than modeling collapse-induced generation in astrophysical simulations.
  • The result applies only within current resolution limits; unresolved small-scale effects could behave differently at higher resolution.

Load-bearing premise

The numerical resolution is sufficient to capture any vorticity generation mechanism that would operate during real gravitational collapse.

What would settle it

A higher-resolution simulation under identical initial conditions and barotropic setup that shows a net increase in vorticity attributable to the collapse flow rather than initial conditions would falsify the central claim.

Figures

Figures reproduced from arXiv: 2607.01207 by Axel Brandenburg, Evangelia Ntormousi, Jennifer Schober.

Figure 2
Figure 2. Figure 2: Mesh Reynolds number for Runs A–G (shown in the same colors and line styles as in [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: (a) Time series of the Mach number, Ma = urms/cs0, for different values of ν, along with the nondimen￾sional kinetic energy dissipation, ǫK/c3 s0k0, as well as the time dependence of (b) h(∇ · u) 2 i/c2 s0k 2 0 and h(∇ × u) 2 i/c2 s0k 2 0, for Runs A (black), C (blue), D (green), F (orange), and G (red). The black and green dashed lines denote Runs B and E, which are lower resolution versions of Runs A and… view at source ↗
Figure 3
Figure 3. Figure 3: (b), where we have marked as a red line the average value during the collapse phase (1.3 ≤ t ≤ 2.7), as well as the early maximum at t ≈ 0.6. We refer to these two values as Rgen1 and Rgen2, respectively. In [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Scaling of the ratios Rgen1, Rgen2, Rdyn1, and Rdyn2, and the early peak values of the terms Tgen0, Tdis0, and Tdyn0 with ν. Runs B and E are lower resolution results of Runs A and D, respectively, and are shown as open symbols. The red line is proportional to ν −1 , which is the approximate scaling found for Rdyn. The approximate ν +0.4 scaling for Rdyn2 is more uncertain. Rdyn and Tdyn. Thus, poor resolu… view at source ↗
Figure 5
Figure 5. Figure 5: Compensated spectra of u, and uncompensated spectra of ln ρ and ω, for Run A with ν = 0.01 at t = 0, 10−4 , and 0.1 (all dotted lines), as well as 0.6 (solid black), 0.9 (blue), 1.3 (green), 1.6 (orange), and 2.5 (red). Note the k 4 subinertial range spectrum in ω, in analogy to a similar behavior of the magnetic field in MHD. The approximate k −11/3 spectrum is unrelated to a Kolmogorov spectrum and appea… view at source ↗
Figure 6
Figure 6. Figure 6: Similar to [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Scalings of (∇ · u)rms and ωrms with uini at t = 0.5 [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Dependence of ∆ ln B vs. PrMReω for Runs 15– 19 and 32–34 from Brandenburg & Ntormousi (2025), show￾ing the critical value of PrMReω being around 300. (Haugen et al. 2004; Elias-L´opez et al. 2023, 2024). In [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Similar to [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
read the original abstract

Gravitational collapse creates large amounts of kinetic energy that could potentially seed turbulence. If such turbulence were also suitable to initiate dynamo action, the resulting magnetic field would further modify the dynamics, especially on small length scales. However, a small-scale dynamo requires vortical turbulence, while the collapse produces mainly irrotational motions, which may not be efficient for dynamo action. Here, we study the efficiency of vorticity production during a turbulent collapse. We use a barotropic equation of state, where pressure and density gradients are parallel, and no magnetic field, so that vorticity can only be produced by viscosity. Using direct numerical simulations of gravitational collapse, we show that, for the parameter space accessible to our numerical resolution, this effect is related to the initial irrotational turbulence and is not a consequence of the collapse flow.

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

0 major / 2 minor

Summary. The manuscript reports direct numerical simulations of gravitational collapse starting from initially irrotational turbulent conditions. Using a barotropic equation of state and no magnetic fields (so that vorticity can arise only from viscous effects), the authors find no evidence that the collapse flow itself generates vorticity; any vorticity present is instead attributable to the initial turbulence, within the parameter space accessible to the employed numerical resolution.

Significance. If the result holds under the stated scoping, it establishes a controlled negative finding that purely gravitational collapse of irrotational flow does not efficiently source vortical motions. This baseline is useful for assessing the conditions required for small-scale dynamo action in astrophysical collapse problems such as star formation, and it highlights the dominant role of initial conditions over the collapse dynamics in this simplified setup.

minor comments (2)
  1. [Abstract and § Conclusions] The abstract and introduction clearly scope the negative result to the accessible numerical resolution and the barotropic, non-magnetized setup; this scoping should be repeated explicitly in the conclusions to avoid over-generalization.
  2. [Numerical methods] A brief resolution study or convergence test (e.g., comparing vorticity spectra or enstrophy evolution across at least two grid resolutions) would strengthen the claim that the absence of collapse-induced vorticity is not a numerical artifact.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary of our work and for recognizing the value of this controlled negative result. The recommendation for minor revision is noted. No specific major comments were provided in the report, so we have no points requiring response or revision at this time.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper reports results from direct numerical simulations of gravitational collapse under a barotropic EOS with no magnetic fields. The central finding—that vorticity production is attributable to initial irrotational turbulence rather than the collapse flow—is scoped explicitly to the accessible numerical resolution and is presented as an empirical observation from the runs, not as a derivation or prediction that reduces to fitted parameters or self-citations by construction. No load-bearing equations, uniqueness theorems, or ansatzes are invoked that would create the enumerated circularity patterns. The result is therefore self-contained against external benchmarks within its stated limits.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the assumptions that vorticity production is limited to viscous effects under the chosen equation of state and that the accessible numerical resolution captures all relevant physics.

free parameters (2)
  • numerical resolution
    Limits the parameter space for which the no-vorticity-production conclusion holds.
  • viscosity coefficient
    Determines the only allowed source of vorticity in the model.
axioms (2)
  • domain assumption Pressure is a function of density only (barotropic equation of state).
    Invoked to ensure pressure and density gradients remain parallel, eliminating baroclinic vorticity production.
  • domain assumption No magnetic field is present.
    Removes Lorentz-force contributions to vorticity.

pith-pipeline@v0.9.1-grok · 5675 in / 1338 out tokens · 18424 ms · 2026-07-02T05:04:00.711667+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

42 extracted references · 39 canonical work pages · 2 internal anchors

  1. [1]

    Bonnor , W. B. 1956, , 116, 351, 10.1093/mnras/116.3.351

    work page doi:10.1093/mnras/116.3.351 1956

  2. [2]

    2003, in Advances in Nonlinear Dynamics, ed

    Brandenburg , A. 2003, in Advances in Nonlinear Dynamics, ed. A. Ferriz-Mas & M. N \'u \ n ez , 269, 10.1201/9780203493137.ch9

    work page doi:10.1201/9780203493137.ch9 2003

  3. [3]

    2022, , 513, 2136, 10.1093/mnras/stac982

    Brandenburg , A., & Ntormousi , E. 2022, , 513, 2136, 10.1093/mnras/stac982

    work page doi:10.1093/mnras/stac982 2022

  4. [4]

    2025, , 990, 223, 10.3847/1538-4357/adf725

    ---. 2025, , 990, 223, 10.3847/1538-4357/adf725

    work page doi:10.3847/1538-4357/adf725 2025

  5. [5]

    2023, , 518, 6367, 10.1093/mnras/stac3555

    Brandenburg , A., Rogachevskii , I., & Schober , J. 2023, , 518, 6367, 10.1093/mnras/stac3555

    work page doi:10.1093/mnras/stac3555 2023

  6. [6]

    2025, , 983, 105, 10.3847/1538-4357/adbe38

    Brandenburg , A., & Scannapieco , E. 2025, , 983, 105, 10.3847/1538-4357/adbe38

    work page doi:10.3847/1538-4357/adbe38 2025

  7. [7]

    2001, , 379, 1153, 10.1051/0004-6361:20011400

    Brandenburg , A., & von Rekowski , B. 2001, , 379, 1153, 10.1051/0004-6361:20011400

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

  8. [8]

    2011, , 528, A145, 10.1051/0004-6361/201015661

    Del Sordo , F., & Brandenburg , A. 2011, , 528, A145, 10.1051/0004-6361/201015661

    work page doi:10.1051/0004-6361/201015661 2011

  9. [9]

    1955, , 37, 217

    Ebert , R. 1955, , 37, 217

    1955

  10. [10]

    2023, , 677, A46, 10.1051/0004-6361/202346696

    Elias-L \'o pez , A., Del Sordo , F., & Vigan \`o , D. 2023, , 677, A46, 10.1051/0004-6361/202346696

    work page doi:10.1051/0004-6361/202346696 2023

  11. [11]

    2024, , 690, A77, 10.1051/0004-6361/202450398

    ---. 2024, , 690, A77, 10.1051/0004-6361/202450398

    work page doi:10.1051/0004-6361/202450398 2024

  12. [12]

    2003, , 68, 016311, 10.1103/PhysRevE.68.016311

    Elperin , T., Kleeorin , N., & Rogachevskii , I. 2003, , 68, 016311, 10.1103/PhysRevE.68.016311

    work page doi:10.1103/physreve.68.016311 2003

  13. [13]

    1994, PhFl, 6, 1411, 10.1063/1.868255

    Falkovich , G. 1994, PhFl, 6, 1411, 10.1063/1.868255

    work page doi:10.1063/1.868255 1994

  14. [14]

    2011, , 107, 114504, 10.1103/PhysRevLett.107.114504

    Federrath , C., Chabrier , G., Schober , J., et al. 2011, , 107, 114504, 10.1103/PhysRevLett.107.114504

    work page doi:10.1103/physrevlett.107.114504 2011

  15. [15]

    2007 , month =

    Field , G. B., Blackman , E. G., & Keto , E. R. 2008, , 385, 181, 10.1111/j.1365-2966.2007.12609.x

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

  16. [16]

    S., & Sulem , P

    Frisch , U., She , Z. S., & Sulem , P. L. 1987, PhyD, 28, 382, 10.1016/0167-2789(87)90026-1

    work page doi:10.1016/0167-2789(87)90026-1 1987

  17. [17]

    Haugen , N. E. L., Brandenburg , A., & Mee , A. J. 2004, , 353, 947, 10.1111/j.1365-2966.2004.08127.x

    work page doi:10.1111/j.1365-2966.2004.08127.x 2004

  18. [18]

    2021, , 655, A3, 10.1051/0004-6361/202141650

    Hennebelle , P. 2021, , 655, A3, 10.1051/0004-6361/202141650

    work page doi:10.1051/0004-6361/202141650 2021

  19. [19]

    J., Mitra , D., & Brandenburg , A

    K \"a pyl \"a , P. J., Mitra , D., & Brandenburg , A. 2009, , 79, 016302, 10.1103/PhysRevE.79.016302

    work page doi:10.1103/physreve.79.016302 2009

  20. [20]

    P., Ruzmaikin , A

    Kazantsev , A. P., Ruzmaikin , A. A., & Sokolov , D. D. 1985, ZhETF, 61, 285

    1985

  21. [21]

    S., & Hennebelle , P

    Klessen , R. S., & Hennebelle , P. 2010, , 520, A17, 10.1051/0004-6361/200913780

    work page doi:10.1051/0004-6361/200913780 2010

  22. [22]

    1974, AN, 295, 93, 10.1002/asna.19742950205

    Krause , F., & R \"u diger , G. 1974, AN, 295, 93, 10.1002/asna.19742950205

    work page doi:10.1002/asna.19742950205 1974

  23. [23]

    Larson , R. B. 1969, , 145, 271, 10.1093/mnras/145.3.271

    work page doi:10.1093/mnras/145.3.271 1969

  24. [24]

    2003, Reports on Progress in Physics, 66, 1651, 10.1088/0034-4885/66/10/R03

    ---. 2003, Reports on Progress in Physics, 66, 1651, 10.1088/0034-4885/66/10/R03

    work page doi:10.1088/0034-4885/66/10/r03 2003

  25. [25]

    Levina , G. V. 2019, in JPhCS, Vol. 1336, JPhCS (IOP), 012007, 10.1088/1742-6596/1336/1/012007

    work page doi:10.1088/1742-6596/1336/1/012007 2019

  26. [26]

    2019, RSPSA, 475, 20180591, 10.1098/rspa.2018.0591

    Martins Afonso , M., Mitra , D., & Vincenzi , D. 2019, RSPSA, 475, 20180591, 10.1098/rspa.2018.0591

    work page doi:10.1098/rspa.2018.0591 2019

  27. [27]

    Malicious User Experience Design Research for Cybersecurity

    McKee , C. F., & Ostriker , E. C. 2007, , 45, 565, 10.1146/annurev.astro.45.051806.110602

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1146/annurev.astro.45.051806.110602 2007

  28. [28]

    2006 , month =

    Mee , A. J., & Brandenburg , A. 2006, , 370, 415, 10.1111/j.1365-2966.2006.10476.x

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

  29. [29]

    S., Sagdeev , R

    Moiseev , S. S., Sagdeev , R. Z., Tur , A. V., Khomenko , G. A., & Yanovskii , V. V. 1983, JETP, 58, 1149

    1983

  30. [30]

    2026, , 136, 091201, 10.1103/fp1v-xrr5

    Muhammed Irshad , P., Bhat , P., Subramanian , K., & Shukurov , A. 2026, , 136, 091201, 10.1103/fp1v-xrr5

    work page doi:10.1103/fp1v-xrr5 2026

  31. [31]

    1995, , 455, 536, 10.1086/176603

    Passot , T., Vazquez-Semadeni , E., & Pouquet , A. 1995, , 455, 536, 10.1086/176603

    work page doi:10.1086/176603 1995

  32. [32]

    2021, JOSS, 6, 2807, 10.21105/joss.02807

    Pencil Code Collaboration , Brandenburg , A., Johansen , A., et al. 2021, JOSS, 6, 2807, 10.21105/joss.02807

    work page doi:10.21105/joss.02807 2021

  33. [33]

    Penston , M. V. 1969, , 144, 425, 10.1093/mnras/144.4.425

    work page doi:10.1093/mnras/144.4.425 1969

  34. [34]

    H., Jones , T

    Porter , D. H., Jones , T. W., & Ryu , D. 2015, , 810, 93, 10.1088/0004-637X/810/2/93

    work page doi:10.1088/0004-637x/810/2/93 2015

  35. [35]

    Imprints of primordial magnetic fields in gravitational collapse during early structure formation

    Schober , J., Abramson , M., Mandal , S., Mtchedlidze , S., & Kahniashvili , T. 2026, arXiv e-prints, arXiv:2602.23263, 10.48550/arXiv.2602.23263

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2602.23263 2026

  36. [36]

    2023, , 2023, 042, 10.1088/1475-7516/2023/12/042

    Sharma , R., Dahl , J., Brandenburg , A., & Hindmarsh , M. 2023, , 2023, 042, 10.1088/1475-7516/2023/12/042

    work page doi:10.1088/1475-7516/2023/12/042 2023

  37. [37]

    Shu , F. H. 1977, , 214, 488, 10.1086/155274

    work page doi:10.1086/155274 1977

  38. [38]

    L., She , Z

    Sulem , P. L., She , Z. S., Scholl , H., & Frisch , U. 1989, JFM, 205, 341, 10.1017/S0022112089002065

    work page doi:10.1017/s0022112089002065 1989

  39. [39]

    Sur , S., Federrath , C., Schleicher , D. R. G., Banerjee , R., & Klessen , R. S. 2012, , 423, 3148, 10.1111/j.1365-2966.2012.21100.x

    work page doi:10.1111/j.1365-2966.2012.21100.x 2012

  40. [40]

    Sur , S., Schleicher , D. R. G., Banerjee , R., Federrath , C., & Klessen , R. S. 2010, , 721, L134, 10.1088/2041-8205/721/2/L134

    work page doi:10.1088/2041-8205/721/2/l134 2010

  41. [41]

    2020 a , , 890, 157, 10.3847/1538-4357/ab6e63

    Xu , S., & Lazarian , A. 2020 a , , 890, 157, 10.3847/1538-4357/ab6e63

    work page doi:10.3847/1538-4357/ab6e63 2020

  42. [42]

    2020 b , , 899, 115, 10.3847/1538-4357/aba7ba

    ---. 2020 b , , 899, 115, 10.3847/1538-4357/aba7ba

    work page doi:10.3847/1538-4357/aba7ba 2020