Hybrid Simulations of Supersonic Shear Flows: II) Cosmic Ray Viscosity
Pith reviewed 2026-05-10 15:45 UTC · model grok-4.3
The pith
Cosmic ray particles introduce an effective viscosity in shear flows by carrying momentum across layers when their gyroradii are smaller than the shear scale.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Particles with large gyroradii act as long-range messengers that promote momentum exchange between layers, hence introducing a form of cosmic ray viscosity. Even when not energetically dominant, increasing the cosmic ray energy density generally enhances momentum transfer, provided that their gyroradii are smaller than the shear lengthscale.
What carries the argument
Cosmic ray viscosity, the process by which nonthermal particles with gyroradii comparable to or smaller than the shear scale transport momentum across layers in Kelvin-Helmholtz unstable flows.
If this is right
- Shear dissipation accelerates as cosmic ray energy density rises under the gyroradius condition, shortening the time for the flow to become turbulent.
- Energy from the initial shear is partitioned more toward particle heating and acceleration rather than remaining in bulk motion.
- Magnetic field amplification grows in tandem with the enhanced momentum mixing driven by the cosmic rays.
- The maximum energies reached by accelerated particles increase because the effective viscosity sustains stronger turbulence longer.
Where Pith is reading between the lines
- In astrophysical settings such as relativistic jets or galactic outflows, cosmic ray viscosity could regulate mixing rates even at modest cosmic ray pressures.
- Extending these results to three dimensions might reveal whether the viscosity alters the cascade of turbulent energy to smaller scales.
- Laboratory laser-plasma experiments with controlled shear and injected energetic particles could test the predicted dependence of dissipation on gyroradius.
Load-bearing premise
The hybrid kinetic-ion fluid-electron treatment of two-dimensional sinusoidal shear flows accurately represents the dominant momentum-exchange physics when gyroradii are smaller than the shear scale.
What would settle it
A direct measurement in a controlled plasma experiment or high-resolution observation showing that momentum transfer rate across the shear interface fails to increase with cosmic ray density when gyroradii satisfy the scale condition.
Figures
read the original abstract
In this second paper in a series dedicated to characterizing shear layers via 2D hybrid (kinetic ions -- fluid electrons) simulations, we study the dynamical role of nonthermal particles (cosmic rays, CRs), either spontaneously generated or pre-existing. We initialize Kolmogorov-type sinusoidal velocity shear flows unstable to the Kelvin--Helmholtz instability, which evolve nonlinearly into turbulence. Particles with large gyroradii act as long-range messengers that promote momentum exchange between layers, hence introducing a form of CR viscosity. Even when not energetically dominant, increasing the CR energy density generally enhances momentum transfer, provided that their gyroradii are smaller than the shear lengthscale. We consider flows ranging from subsonic to supersonic and assess the rate of shear dissipation, the partition of the initial kinetic energy among heating, thermal ion acceleration, CR reacceleration, and magnetic-field amplification, and the maximum energy attained by accelerated particles.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports 2D hybrid (kinetic ions, fluid electrons) simulations of Kelvin-Helmholtz-unstable Kolmogorov-type sinusoidal shear flows, with cosmic rays either pre-existing or spontaneously generated. It claims that particles with large gyroradii act as long-range messengers that enhance momentum exchange between layers (introducing a form of CR viscosity), that this enhancement occurs even when CRs are not energetically dominant provided gyroradii remain smaller than the shear lengthscale, and that the study quantifies shear dissipation rates, energy partitioning among heating, ion/CR acceleration and magnetic amplification, plus maximum particle energies, across subsonic-to-supersonic Mach numbers and varying CR energy densities.
Significance. If the quantitative results hold, the work supplies a first-principles demonstration of how nonthermal particles can mediate momentum transport in turbulent shear layers without being energetically dominant, which is relevant to astrophysical environments such as jets, supernova remnants and galactic winds. The parameter survey across Mach number and CR energy fraction, together with the hybrid kinetic treatment of ions, constitutes a concrete advance over purely fluid models.
major comments (2)
- The manuscript provides insufficient detail on numerical resolution, the precise method used to initialize or inject CR particles, and any convergence tests performed with respect to grid size or particle number. These omissions directly affect the reliability of the reported momentum-transfer enhancements and energy-partition fractions described in the abstract and results sections.
- All simulations are performed in 2D. Because 2D turbulence lacks vortex stretching, supports an inverse energy cascade, and yields different mixing and dissipation statistics than 3D turbulence, the manner in which CRs sample and transport momentum across shear layers may differ; the central claim that the gyroradius-to-shear-scale ratio controls CR viscosity therefore requires explicit discussion or supporting 3D tests to establish generality.
minor comments (1)
- The abstract refers to 'Kolmogorov-type sinusoidal velocity shear flows' without a concise definition or reference to the first paper in the series; a short clarifying phrase would improve accessibility.
Simulated Author's Rebuttal
We thank the referee for their constructive comments and positive evaluation of the work's significance. We address each major comment below, indicating where revisions will be made to improve clarity and completeness.
read point-by-point responses
-
Referee: The manuscript provides insufficient detail on numerical resolution, the precise method used to initialize or inject CR particles, and any convergence tests performed with respect to grid size or particle number. These omissions directly affect the reliability of the reported momentum-transfer enhancements and energy-partition fractions described in the abstract and results sections.
Authors: We agree that more explicit numerical details are required for reproducibility and to substantiate the robustness of the momentum-transfer and energy-partition results. In the revised manuscript we will expand the numerical methods section to specify the grid resolution (cells per shear scale length), the number of particles per cell for both ions and CRs, the precise initialization of the Kolmogorov-type sinusoidal shear profile and the CR distribution function (including how pre-existing CRs are sampled to achieve target energy densities and gyroradii), and the outcomes of convergence tests performed by varying grid size and particle number. These additions will demonstrate that the reported CR-viscosity enhancements remain stable within the parameter ranges explored. revision: yes
-
Referee: All simulations are performed in 2D. Because 2D turbulence lacks vortex stretching, supports an inverse energy cascade, and yields different mixing and dissipation statistics than 3D turbulence, the manner in which CRs sample and transport momentum across shear layers may differ; the central claim that the gyroradius-to-shear-scale ratio controls CR viscosity therefore requires explicit discussion or supporting 3D tests to establish generality.
Authors: We recognize the intrinsic limitations of two-dimensional turbulence, including the absence of vortex stretching and the inverse energy cascade, which can alter mixing and dissipation relative to three dimensions. Our study deliberately employs 2D hybrid simulations to isolate the long-range momentum-transport role of CRs in a controlled setting, consistent with the preceding paper in the series. In the revised manuscript we will insert a new subsection discussing the 2D approximation, arguing that the gyroradius-to-shear-scale ratio remains a controlling parameter because it governs the ability of particles to traverse shear layers irrespective of the details of the turbulent cascade; we will also qualify that quantitative transport rates may differ in 3D and outline the need for future three-dimensional extensions. Performing new 3D runs lies outside the scope of the present revision. revision: partial
Circularity Check
Minor self-citation to prior series paper; central results from direct hybrid simulations
full rationale
The paper presents results from 2D hybrid kinetic-ion/fluid-electron simulations of Kelvin-Helmholtz unstable shear flows. The central claim that large-gyroradius CRs act as long-range messengers enhancing momentum exchange (CR viscosity) when r_g < shear scale follows directly from the particle trajectories and momentum transport observed in the runs. No analytic derivation reduces to fitted parameters or self-defined quantities. The only self-reference is to the preceding paper in the series for setup details, which is not load-bearing for the viscosity conclusion. The work is therefore self-contained against external benchmarks and receives only the minimal score for a normal series citation.
Axiom & Free-Parameter Ledger
free parameters (2)
- CR energy density fraction
- Gyroradius to shear lengthscale ratio
axioms (2)
- domain assumption Hybrid approximation (kinetic ions, fluid electrons) captures the essential ion-scale dynamics
- domain assumption 2D geometry is sufficient to capture the Kelvin-Helmholtz evolution and CR transport
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Particles with large gyroradii act as long-range messengers that promote momentum exchange between layers, hence introducing a form of CR viscosity... provided that their gyroradii are smaller than the shear lengthscale.
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We use hybrid particle-in-cell (PIC) simulations with kinetic ions and fluid electrons...
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...
-
[2]
" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...
-
[3]
thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...
work page 2017
-
[4]
Armitage , P. J. 2011, title Dynamics of Protoplanetary Disks , , 49, 195, 10.1146/annurev-astro-081710-102521
-
[5]
Bai , X.-N., Caprioli , D., Sironi , L., & Spitkovsky , A. 2015, title Magnetohydrodynamic-PIC Method for Coupling Cosmic Rays with a Thermal Plasma: Application to Non-relativistic Shocks , , 809, 55, 10.1088/0004-637X/809/1/55
-
[6]
Balbus, S. A., & Hawley, J. F. 1998, title Instability, turbulence, and enhanced transport in accretion disks, Rev. Mod. Phys., 70, 1, 10.1103/RevModPhys.70.1
-
[7]
Bell , A. R. 1978, title The acceleration of cosmic rays in shock fronts. I , MNRAS, 182, 147. https://ui.adsabs.harvard.edu/abs/1978MNRAS.182..147B/abstract
work page 1978
-
[8]
Belyaev, M. A., & Rafikov, R. R. 2012, title SUPERSONIC SHEAR INSTABILITIES IN ASTROPHYSICAL BOUNDARY LAYERS, The Astrophysical Journal, 752, 115, 10.1088/0004-637x/752/2/115
-
[9]
Braginskii , S. I. 1965, title Transport Processes in a Plasma , Reviews of Plasma Physics, 1, 205
work page 1965
-
[10]
2014, ApJ, 783, 91, doi: 10.1088/0004-637X/783/2/91
Caprioli , D., & Spitkovsky , A. 2014, title Simulations of Ion Acceleration at Non-relativistic Shocks: I. Acceleration Efficiency , , 783, 91, 10.1088/0004-637X/783/2/91
-
[11]
Diffusive Shock Re-Acceleration
Caprioli, D., Zhang, H., & Spitkovsky, A. 2018, title Diffusive Shock Re-Acceleration, . 1801.01510
work page Pith review arXiv 2018
-
[12]
1961, Hydrodynamic and hydromagnetic stability
Chandrasekhar , S. 1961, Hydrodynamic and hydromagnetic stability
work page 1961
-
[13]
Chow, A., Rowan, M. E., Sironi, L., et al. 2023, title Linear analysis of the Kelvin-Helmholtz instability in relativistic magnetized symmetric flows, Monthly Notices of the Royal Astronomical Society, 524, 90, 10.1093/mnras/stad1833
-
[14]
S., Califano, F., & Pegoraro, F
De Camillis, S., Cerri, S. S., Califano, F., & Pegoraro, F. 2016, title Pressure anisotropy generation in a magnetized plasma configuration with a shear flow velocity, Plasma Physics and Controlled Fusion, 58, 045007, 10.1088/0741-3335/58/4/045007
-
[15]
Earl , J. A., Jokipii , J. R., & Morfill , G. 1988, title Cosmic-Ray Viscosity , , 331, L91, 10.1086/185242
-
[16]
Fermi , E. 1954, title Galactic Magnetic Fields and the Origin of Cosmic Radiation. , , 119, 1, 10.1086/145789
-
[17]
Frank , A., Jones , T. W., Ryu , D., & Gaalaas , J. B. 1996, title The Magnetohydrodynamic Kelvin-Helmholtz Instability: A Two-dimensional Numerical Study , , 460, 777, 10.1086/177009
-
[18]
Fraser, A. E., Kaminski, A. K., & Oishi, J. S. 2025, Nonmodal growth and optimal perturbations in magnetohydrodynamic shear flows, 10.48550/arXiv.2510.03141
-
[19]
Gonz \'a lez, A. G., & Gratton, J. 1994, title The Kelvin--Helmholtz instability in a compressible plasma: the role of the orientation of the magnetic field with respect to the flow, Journal of Plasma Physics, 51, 43, 10.1017/s0022377800017384
-
[20]
Haggerty , C. C., & Caprioli , D. 2019, title dHybridR: A Hybrid Particle-in-cell Code Including Relativistic Ion Dynamics , , 887, 165, 10.3847/1538-4357/ab58c8
-
[21]
Henri, P., Cerri, S. S., Califano, F., et al. 2013, title Nonlinear evolution of the magnetized Kelvin-Helmholtz instability: From fluid to kinetic modeling, Physics of Plasmas, 20, 10.1063/1.4826214
-
[22]
Hillas , A. M. 1984, title The Origin of Ultra-High-Energy Cosmic Rays , Ann. Rev. of A&A, 22, 425, 10.1146/annurev.aa.22.090184.002233
-
[23]
Hunt, J. C. R. 1966, title On the stability of parallel flows with parallel magnetic fields, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences, 293, 342, 10.1098/rspa.1966.0177
-
[24]
Kunz , M. W., Schekochihin , A. A., Cowley , S. C., et al. 2011, title A thermally stable heating mechanism for the intracluster medium: turbulence, magnetic fields and plasma instabilities , , 410, 2446, 10.1111/j.1365-2966.2010.17621.x
-
[25]
Kunz , M. W., Schekochihin , A. A., & Stone , J. M. 2014, title Firehose and Mirror Instabilities in a Collisionless Shearing Plasma , , 112, 205003, 10.1103/PhysRevLett.112.205003
-
[26]
Liang , N., & Caprioli , D. 2025, title Hybrid Simulations of Supersonic Shear Flows: I) Particle Acceleration , arXiv e-prints, arXiv:2511.10739, 10.48550/arXiv.2511.10739
-
[27]
Liu , M., Ruszkowski , M., Zweibel , E., et al. 2025, title Particle Acceleration in Magnetized Shear-Driven Turbulence , arXiv e-prints, arXiv:2512.12720, 10.48550/arXiv.2512.12720
-
[28]
Liu , W., Goodman , J., & Ji , H. 2006, title Simulations of Magnetorotational Instability in a Magnetized Couette Flow , , 643, 306, 10.1086/501495
-
[29]
Meshalkin, L., & Sinai, I. 1961, title Investigation of the Stability of a Stationary Solution of a System of Equations for the Plane Movement of an Incompressible Viscous Liquid, Journal of Applied Mathematics and Mechanics, 25, 1700, 10.1016/0021-8928(62)90149-1
-
[30]
Miura , A. 1997, title Compressible magnetohydrodynamic Kelvin-Helmholtz instability with vortex pairing in the two-dimensional transverse configuration , Physics of Plasmas, 4, 2871, 10.1063/1.872419
-
[31]
Oparina, E., & Troshkin, O. 2004, title Stability of Kolmogorov Flow in a Channel with Rigid Walls, Doklady Physics, 49, 583, 10.1134/1.1815419
-
[32]
Rieger, F. M. 2019, title An Introduction to Particle Acceleration in Shearing Flows, Galaxies, 7, 10.3390/galaxies7030078
-
[33]
Rieger , F. M., & Duffy , P. 2004, title Shear Acceleration in Relativistic Astrophysical Jets , , 617, 155, 10.1086/425167
-
[34]
Rieger , F. M., & Duffy , P. 2019, title Particle Acceleration in Shearing Flows: Efficiencies and Limits , , 886, L26, 10.3847/2041-8213/ab563f
-
[35]
Rosin, M. S., Schekochihin, A. A., Rincon, F., & Cowley, S. C. 2011, title A non-linear theory of the parallel firehose and gyrothermal instabilities in a weakly collisional plasma, Monthly Notices of the Royal Astronomical Society, 413, 7, 10.1111/j.1365-2966.2010.17931.x
-
[36]
Sandoval, A., Riquelme, M., Spitkovsky, A., & Bacchini, F. 2024, title Particle-in-cell simulations of the magnetorotational instability in stratified shearing boxes, Monthly Notices of the Royal Astronomical Society, 530, 1866, 10.1093/mnras/stae959
-
[37]
Schekochihin , A. A. 2022, title MHD turbulence: a biased review , Journal of Plasma Physics, 88, 155880501, 10.1017/S0022377822000721
-
[38]
Schekochihin, A. A., Cowley, S. C., Rincon, F., & Rosin, M. S. 2010, title Magnetofluid dynamics of magnetized cosmic plasma: firehose and gyrothermal instabilities, Monthly Notices of the Royal Astronomical Society, 405, 291, 10.1111/j.1365-2966.2010.16493.x
-
[39]
Simionescu, A., ZuHone, J., Zhuravleva, I., et al. 2019, title Constraining Gas Motions in the Intra-Cluster Medium, , 215, 24, 10.1007/s11214-019-0590-1
-
[40]
Sommeria, J., & Moreau, R. 1982, title Why, How, and When, MHD Turbulence Becomes Two-Dimensional, Journal of Fluid Mechanics, 118, 507, 10.1017/S0022112082001177
-
[41]
1992, title Instabilities in Two-Dimensional Spatially Periodic Flows
Thess, A. 1992, title Instabilities in Two-Dimensional Spatially Periodic Flows. Part I: Kolmogorov Flow, Physics of Fluids A: Fluid Dynamics (1989-1993), 4, 1385, 10.1063/1.858415
-
[42]
M., Al-Nussirat , S., Mostafavi , P., et al
Webb , G. M., Al-Nussirat , S., Mostafavi , P., et al. 2019, title Particle Acceleration by Cosmic Ray Viscosity in Radio-jet Shear Flows , , 881, 123, 10.3847/1538-4357/ab2fca
-
[43]
Webb , G. M., Barghouty , A. F., Hu , Q., & le Roux , J. A. 2018, title Particle Acceleration Due to Cosmic-ray Viscosity and Fluid Shear in Astrophysical Jets , , 855, 31, 10.3847/1538-4357/aaae6c
-
[44]
Zuhone , J. A., & Roediger , E. 2016, title Cold fronts: probes of plasma astrophysics in galaxy clusters , Journal of Plasma Physics, 82, 535820301, 10.1017/S0022377816000544
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.