pith. sign in

arxiv: 2510.20534 · v2 · submitted 2025-10-23 · ⚛️ physics.flu-dyn · astro-ph.IM

Are heuristic switches necessary to control dissipation in modern smoothed particle hydrodynamics ?

Domingo Garc\'ia-Senz , Rub\'en M. Cabez\'on This is my paper

Pith reviewed 2026-05-18 04:53 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn astro-ph.IM
keywords smoothed particle hydrodynamicsartificial viscosityshock capturingnumerical dissipationvelocity reconstructionBalsara correctionshear flowshydrodynamic instabilities
0
0 comments X

The pith

Removing the linear velocity component plus Balsara correction controls dissipation in SPH without heuristic switches.

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

The paper develops and tests a shock-capturing scheme for smoothed particle hydrodynamics that dispenses with complex artificial viscosity switches. It subtracts the linear part of the local velocity field from each particle pair and modulates the remaining dissipation with the Balsara correction. The resulting scheme is shown to handle subsonic instabilities, shear flows, and strong shocks while cutting spurious dissipation relative to a standard reference switch. Readers care because the method replaces an imperfect, noise-generating piece of machinery with a simpler, more uniform treatment of numerical viscosity.

Core claim

By removing the linear component of the local velocity field and applying the Balsara correction, the method produces a balanced artificial-viscosity scheme that performs across subsonic instabilities, shear flows, and strong shocks while delivering higher accuracy and lower spurious dissipation than the reference viscosity switch.

What carries the argument

Velocity-reconstruction step that subtracts bulk linear motion from the local velocity field, combined with the Balsara correction to scale the dissipation term.

If this is right

  • Spurious dissipation drops in shear-dominated regions while shock capturing remains effective.
  • Numerical noise from viscosity switches is eliminated across the tested regimes.
  • Overall accuracy rises in mixed-regime simulations that include both instabilities and shocks.
  • Implementation becomes simpler because no separate heuristic trigger logic is required.

Where Pith is reading between the lines

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

  • The same reconstruction idea could be ported to other mesh-free or particle-based hydrodynamics codes to reduce switch-related artifacts.
  • Long-duration astrophysical runs might show less cumulative drift if this uniform dissipation control replaces switched schemes.
  • Targeted tests in highly supersonic or magnetized flows would check whether the method still balances dissipation correctly.

Load-bearing premise

Removing the linear component of the local velocity field together with the Balsara correction is enough to suppress spurious dissipation in shear regions without missing needed dissipation or creating new instabilities.

What would settle it

A controlled shear-flow or Kelvin-Helmholtz test in which the new scheme either damps physical motions too strongly or allows unphysical oscillations that the reference switched scheme suppresses.

Figures

Figures reproduced from arXiv: 2510.20534 by Domingo Garc\'ia-Senz, Rub\'en M. Cabez\'on.

Figure 1
Figure 1. Figure 1: Velocity profile, vx , for the shock-tube test and the models shown in the second row in [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of maximum density in the Sedov test for the different models described in [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Density and radial velocity profiles in the Sedov-Taylor test at time t = 1 for the different models described in [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Azimuthal velocity distribution in the Gresho-Chan vortex experiment at t = 1 for models Vk in the third row in [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Color-coded density for the KH models in [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Amplitude mode growth of the KH models. A zoom-in of the nearby lines KH3, KH4, KH5, and KH6 around t = 1.5 is shown in the subfigure. the stored kinetic energy in the vertical direction is clearly lower ( [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Colormaps of density for the RT models in [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: Instantaneous compensated velocity power spectra at t = 10 [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Ranking heatmap of all AV methods for each test. Numbers show the relative difference to the best metric result overall for each test, using the values reported in [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
read the original abstract

Artificial viscosity is commonly employed in smoothed particle hydrodynamics (SPH) to model dissipation in hydrodynamic simulations. However, its practical implementation today relies, in many cases, on complex numerical switches to restrict its application to regions where dissipation is physically warranted, such as shocks. These switches, while essential, are imperfect and can introduce additional numerical noise. In this work we develop and validate a more efficient shock capture scheme for SPH that does not rely on artificial viscosity switches. Recent studies have proposed that subtracting the linear component of the velocity field can suppress spurious dissipation in shear-dominated regions. Building on this idea, we implement a velocity-reconstruction technique that removes the bulk linear motion from the local velocity field and uses the Balsara correction to modulate the dissipation. The presented methodology yields a balanced dissipation scheme that performs well across a range of regimes, including subsonic instabilities, shear flows, and strong shocks. We demonstrate that this approach yields improved accuracy and lower spurious dissipation, compared to the reference viscosity switch used in this work.

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

2 major / 2 minor

Summary. The manuscript proposes an SPH artificial-viscosity scheme that dispenses with heuristic switches by subtracting the linear component of the local velocity field inside each kernel and modulating the remaining dissipation with the Balsara factor. The authors report that the resulting scheme maintains adequate shock capturing while suppressing spurious dissipation in shear and subsonic instability regimes, and that it outperforms a reference viscosity switch in accuracy and noise metrics across the tested cases.

Significance. A switch-free dissipation formulation that remains stable and accurate across shocks, shear, and instabilities would simplify SPH implementations and reduce one source of numerical noise. The work builds directly on recent velocity-reconstruction ideas and supplies a concrete, Balsara-augmented realization together with comparative tests; if the central assumption holds under broader scrutiny, the result would be of clear practical value to the SPH community.

major comments (2)
  1. [Abstract and §3 (velocity-reconstruction description)] The central claim that linear-velocity subtraction plus the Balsara correction is sufficient to eliminate the need for switches rests on the assumption that any non-linear velocity contribution within the kernel support is either negligible or correctly identified as requiring dissipation. This assumption is load-bearing for the “no switches necessary” conclusion and requires explicit verification in regimes with strong curvature or small-scale turbulence, where the linear fit can leave residual relative velocities.
  2. [Abstract] The abstract states that the scheme yields “improved accuracy and lower spurious dissipation” relative to the reference switch, yet no quantitative error norms, convergence rates, or specific test configurations (e.g., exact Mach numbers, resolution, or diagnostic quantities) are supplied. Without these data it is impossible to judge whether the reported improvement is statistically significant or merely visual.
minor comments (2)
  1. [Methods] Notation for the reconstructed velocity field and the precise definition of the Balsara factor as implemented should be given explicitly, preferably with an equation number, to allow immediate reproduction.
  2. [Results figures] Figure captions should state the exact diagnostic used to quantify “spurious dissipation” (e.g., integrated kinetic-energy decay rate or L2 velocity error) so that readers can compare results directly.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and recommendation for major revision. We address each major comment below, indicating revisions where the manuscript will be updated.

read point-by-point responses

  1. Referee: [Abstract and §3 (velocity-reconstruction description)] The central claim that linear-velocity subtraction plus the Balsara correction is sufficient to eliminate the need for switches rests on the assumption that any non-linear velocity contribution within the kernel support is either negligible or correctly identified as requiring dissipation. This assumption is load-bearing for the “no switches necessary” conclusion and requires explicit verification in regimes with strong curvature or small-scale turbulence, where the linear fit can leave residual relative velocities.

    Authors: We agree this assumption is central and merits explicit discussion. The manuscript already tests regimes involving curvature (Sedov blast wave) and developing small-scale structures (Kelvin-Helmholtz instability). To address the concern directly, we have expanded §3 with a derivation showing that higher-order velocity residuals trigger dissipation via the existing Balsara-modulated term, and we have added a short limitations paragraph in the conclusions noting that very high-curvature or fully developed turbulence may require further study. This constitutes a partial revision focused on clarification rather than entirely new simulations. revision: partial

  2. Referee: [Abstract] The abstract states that the scheme yields “improved accuracy and lower spurious dissipation” relative to the reference switch, yet no quantitative error norms, convergence rates, or specific test configurations (e.g., exact Mach numbers, resolution, or diagnostic quantities) are supplied. Without these data it is impossible to judge whether the reported improvement is statistically significant or merely visual.

    Authors: The referee is correct that the abstract is too general. Although comparative results appear in the figures and text of §§4–5, we will revise the abstract to cite concrete elements: the Gresho vortex at 32²–128² resolution, the Kelvin-Helmholtz instability at Mach 0.1 with 512² particles, and the observation of lower L2 velocity errors together with closer agreement to linear growth rates. This change makes the claimed improvements traceable to specific diagnostics. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained

full rationale

The paper proposes a velocity-reconstruction method (subtracting linear velocity component) combined with the Balsara factor as an alternative to heuristic viscosity switches. It explicitly builds on ideas from recent external studies rather than self-citation chains, and validates the scheme through direct comparison to a reference switch across multiple regimes. No equations, fitted parameters, or self-definitional steps are described that would reduce the central claim to its own inputs by construction. The approach introduces independent methodological content and performs empirical tests, keeping the derivation non-circular.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The method rests on standard SPH artificial viscosity formalism and the Balsara correction factor, both drawn from prior literature rather than newly derived here.

axioms (1)
  • domain assumption The Balsara correction remains an appropriate modulator when combined with linear velocity subtraction.
    Invoked as part of the dissipation modulation step; standard in SPH but treated as given.

pith-pipeline@v0.9.0 · 5714 in / 1202 out tokens · 36950 ms · 2026-05-18T04:53:33.670125+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

31 extracted references · 31 canonical work pages

  1. [1]

    2007, MNRAS, 380, 963

    Agertz, O., Moore, B., Stadel, J., et al. 2007, MNRAS, 380, 963

    work page 2007

  2. [2]

    Balsara, D. S. 1995, Journal of Computational Physics, 121, 357

    work page 1995

  3. [3]

    & Springel, V

    Bauer, A. & Springel, V . 2012, MNRAS, 423, 2558 Cabezón, R. M. & García-Senz, D. 2024, Monthly Notices of the Royal Astro- nomical Society, 528, 3782 Cabezón, R. M., García-Senz, D., & Escartín, J. A. 2012, A&A, 545, A112 Cabezón, R. M., García-Senz, D., & Figueira, J. 2017, A&A, 606, A78 Cabezón, R. M., García-Senz, D., & Relaño, A. 2008, Journal of Com...

    work page arXiv 2012

  4. [4]

    2010, MNRAS, 403, 1165

    Cha, S.-H., Inutsuka, S.-I., & Nayakshin, S. 2010, MNRAS, 403, 1165

    work page 2010

  5. [5]

    & Whitworth, A

    Cha, S.-H. & Whitworth, A. P. 2003, MNRAS, 340, 73

    work page 2003

  6. [6]

    & Dehnen, W

    Cullen, L. & Dehnen, W. 2010, MNRAS, 408, 669

    work page 2010

  7. [7]

    S., Schmidt, W., & Mac Low, M

    Federrath, C., Roman-Duval, J., Klessen, R. S., Schmidt, W., & Mac Low, M. M. 2010, Astronomy & Astrophysics, 512, A81

    work page 2010

  8. [8]

    D., & Owen, J

    Frontiere, N., Raskin, C. D., & Owen, J. M. 2017, Journal of Computational Physics, 332, 160 García-Senz, D., Cabezón, R. M., & Escartín, J. A. 2012, A&A, 538, A9 García-Senz, D., Cabezón, R. M., & Escartín, J. A. 2022, A&A, 659, A175

    work page 2017

  9. [9]

    Gingold, R. A. & Monaghan, J. J. 1977, MNRAS, 181, 375

    work page 1977

  10. [10]

    1959, Sb

    Godunov, S.-K. 1959, Sb. Math., 47

    work page 1959

  11. [11]

    Gresho, P. M. & Chan, S. T. 1990, International Journal for Numerical Methods in Fluids, 11, 621

    work page 1990

  12. [12]

    2002, Journal of Computational Physics, 179, 238

    Inutsuka, S.-I. 2002, Journal of Computational Physics, 179, 238

    work page 2002

  13. [13]

    Lucy, L. B. 1977, AJ, 82, 1013

    work page 1977

  14. [14]

    P., Lyra, W., & Passy, J.-C

    McNally, C. P., Lyra, W., & Passy, J.-C. 2012, ApJS, 201, 18

    work page 2012

  15. [15]

    Monaghan, J. J. 1992, ARA&A, 30, 543

    work page 1992

  16. [16]

    Monaghan, J. J. 1997, Journal of Computational Physics, 136, 298

    work page 1997

  17. [17]

    Monaghan, J. J. 2005, Reports on Progress in Physics, 68, 1703

    work page 2005

  18. [18]

    Monaghan, J. J. & Gingold, R. A. 1983, Journal of Computational Physics, 52, 374

    work page 1983

  19. [19]

    Morris, J. P. & Monaghan, J. J. 1997, Journal of Computational Physics, 136, 41

    work page 1997

  20. [20]

    Price, D. J. 2008, Journal of Computational Physics, 227, 10040

    work page 2008

  21. [21]

    J., Wurster, J., Tricco, T

    Price, D. J., Wurster, J., Tricco, T. S., et al. 2018, PASA, 35, e031

    work page 2018

  22. [22]

    Read, J. I. & Hayfield, T. 2012, MNRAS, 422, 3037

    work page 2012

  23. [23]

    2015, MNRAS, 448, 3628

    Rosswog, S. 2015, MNRAS, 448, 3628

    work page 2015

  24. [24]

    2024, arXiv e-prints, arXiv:2411.19228

    Rosswog, S. 2024, arXiv e-prints, arXiv:2411.19228

    work page arXiv 2024

  25. [25]

    & Diener, P

    Rosswog, S. & Diener, P. 2024, arXiv e-prints, arXiv:2404.15952

    work page arXiv 2024

  26. [26]

    D., Eke, V

    Sandnes, T. D., Eke, V . R., Kegerreis, J. A., et al. 2025, Journal of Computational Physics, 532, 113907

    work page 2025

  27. [27]

    Sedov, L. I. 1959, Similarity and Dimensional Methods in Mechanics

    work page 1959

  28. [28]

    Sod, G. A. 1978, Journal of Computational Physics, 27, 1

    work page 1978

  29. [29]

    Uhlenbeck, G. E. & Ornstein, L. S. 1930, Phys. Rev., 36, 823

    work page 1930

  30. [30]

    2016, ApJ, 831, 103 van Leer, B

    Valdarnini, R. 2016, ApJ, 831, 103 van Leer, B. 1974, Journal of Computational Physics, 14, 361

    work page 2016

  31. [31]

    W., Keller, B

    Wadsley, J. W., Keller, B. W., & Quinn, T. R. 2017, MNRAS, 471, 2357 Article number, page 12 of 12

    work page 2017