Evidence on the incompatibility of smoothed particle hydrodynamics and eddy viscosity models for large eddy simulations
Pith reviewed 2026-05-19 10:25 UTC · model grok-4.3
The pith
SPH methods function as Lagrangian large eddy simulations whose overlapping particles generate implicit subfilter stresses that clash with added eddy viscosity models.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
SPH methods operate intrinsically as Lagrangian Large Eddy Simulations for turbulent flows with strongly overlapping discretization elements. These overlapping elements in combination with numerical errors cause a significant amount of implicit subfilter stresses. In the Taylor-Green flow at Re=10^4 these stresses are relevant where turbulent fluctuations are created. Adding eddy viscosity models degrades the turbulent transition process due to the non-locality of these methods.
What carries the argument
Non-locality of SPH particle interactions together with implicit subfilter stresses arising from overlapping discretization elements and numerical errors.
If this is right
- Turbulent transition in SPH remains difficult because implicit subfilter stresses already act where fluctuations originate.
- Direct addition of eddy viscosity models to SPH is expected to interfere with rather than supplement the built-in stresses.
- New closure strategies for SPH must account for the non-local spreading of modeled stresses across overlapping particles.
- The challenge of capturing turbulence in current SPH codes is explained by the intrinsic LES-like character of the method itself.
Where Pith is reading between the lines
- The same non-locality argument may apply to other meshless or particle methods that use overlapping supports for derivative approximation.
- Turbulence models based on integral or non-local operators could be more compatible with SPH than classical local eddy-viscosity closures.
- Redesigning the base SPH discretization to reduce particle overlap might lessen the implicit stresses and allow standard models to work better.
- High-Reynolds-number applications of SPH may require entirely different modeling approaches that treat the method as an already-filtered Lagrangian system.
Load-bearing premise
The observed degradation of turbulent transition when eddy viscosity is added stems specifically from the non-locality of SPH rather than from other numerical artifacts, discretization choices, or test-case specifics.
What would settle it
Repeating the Taylor-Green vortex simulation at Re=10^4 with a strictly local discretization method or with SPH kernel support reduced enough to eliminate significant overlap, then checking whether adding the same eddy viscosity model still degrades the transition.
read the original abstract
In this work, we will present evidence for the incompatibility of Smoothed Particle Hydrodynamics (SPH) methods and eddy viscosity models. Taking a coarse-graining perspective, we physically argue that SPH methods operate intrinsically as Lagrangian Large Eddy Simulations (LES) for turbulent flows with strongly overlapping discretization elements. However, these overlapping elements in combination with numerical errors cause a significant amount of implicit subfilter stresses (SFS). Considering a Taylor-Green flow at $Re=10^4$, the SFS will be shown to be relevant where turbulent fluctuations are created, explaining why turbulent flows are challenging even for current SPH methods. Although one might hope to mitigate the implicit SFS using eddy viscosity models, we show a degradation of the turbulent transition process, which is rooted in the non-locality of these methods.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript argues that SPH methods intrinsically function as Lagrangian LES for turbulent flows because of strongly overlapping discretization elements. These overlaps, combined with numerical errors, generate significant implicit subfilter stresses (SFS). The authors support this with a coarse-graining physical argument and then demonstrate, for a Taylor-Green vortex at Re=10^4, that adding conventional eddy-viscosity models degrades the turbulent transition process; they attribute the degradation to the non-local character of SPH.
Significance. If the central claim is confirmed, the result would be relevant for the development of turbulence closures in Lagrangian particle methods. It would indicate that standard eddy-viscosity models are not readily portable to SPH and would motivate the search for closures that respect the non-local, implicit-filtering nature of the discretization. The work also supplies a concrete numerical example (Taylor-Green at Re=10^4) that can serve as a benchmark for future SPH turbulence studies.
major comments (2)
- [Taylor-Green flow at Re=10^4] Taylor-Green flow section (Re=10^4 results): the degradation of transition when an eddy-viscosity term is added is presented as evidence of incompatibility rooted in non-locality. However, the manuscript does not report auxiliary simulations that vary the kernel support radius, the projection of the viscous term onto particles, or the filter consistency; without such controls it is difficult to exclude the possibility that the observed effect arises from a specific discretization choice rather than from intrinsic non-locality of SPH.
- [Introduction / coarse-graining perspective] Physical argument for implicit SFS: the claim that overlapping elements plus numerical errors produce a 'significant amount' of implicit subfilter stress is central to the incompatibility thesis, yet the manuscript provides only a high-level coarse-graining sketch. A quantitative estimate (e.g., an order-of-magnitude calculation of the implicit stress tensor or a comparison against an explicit filter) would make the argument load-bearing rather than suggestive.
minor comments (2)
- [Abstract] The abstract states that 'turbulent flows are challenging even for current SPH methods' but does not cite prior quantitative studies of SPH transition or decay rates; adding one or two key references would place the new results in context.
- [Physical argument] Notation for the implicit SFS tensor is introduced without an explicit definition or relation to the standard filtered Navier-Stokes equations; a short equation block would improve clarity.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments on our manuscript. We address each major comment point by point below and indicate where revisions will be made to strengthen the presentation.
read point-by-point responses
-
Referee: [Taylor-Green flow at Re=10^4] Taylor-Green flow section (Re=10^4 results): the degradation of transition when an eddy-viscosity term is added is presented as evidence of incompatibility rooted in non-locality. However, the manuscript does not report auxiliary simulations that vary the kernel support radius, the projection of the viscous term onto particles, or the filter consistency; without such controls it is difficult to exclude the possibility that the observed effect arises from a specific discretization choice rather than from intrinsic non-locality of SPH.
Authors: We agree that additional controls would help isolate the role of non-locality. The non-locality is intrinsic to SPH because it stems directly from the overlapping kernel supports that define the discretization, rather than from any single implementation detail. To address the concern, we will add auxiliary simulations in the revised manuscript that vary the kernel support radius while keeping other parameters fixed, demonstrating that the degradation of transition persists. We will also clarify that the viscous term projection and filter consistency follow standard SPH practice for the weakly compressible formulation used here; a short discussion will be added to explain why these choices do not alter the fundamental incompatibility arising from the Lagrangian particle representation. revision: yes
-
Referee: [Introduction / coarse-graining perspective] Physical argument for implicit SFS: the claim that overlapping elements plus numerical errors produce a 'significant amount' of implicit subfilter stress is central to the incompatibility thesis, yet the manuscript provides only a high-level coarse-graining sketch. A quantitative estimate (e.g., an order-of-magnitude calculation of the implicit stress tensor or a comparison against an explicit filter) would make the argument load-bearing rather than suggestive.
Authors: We acknowledge that the coarse-graining argument is presented at a conceptual level in the current manuscript. To make the claim more quantitative, we will revise the introduction to include an order-of-magnitude estimate of the implicit subfilter stress. This estimate will be based on the typical overlap volume of neighboring kernels and the magnitude of numerical truncation errors at the resolutions employed, with a direct comparison to the subfilter stresses that would arise from an explicit filter of comparable width. The addition will be supported by references to existing analyses of SPH as an implicit filter. revision: yes
Circularity Check
No significant circularity; physical argument and simulation evidence are independent of inputs
full rationale
The paper advances a physical argument that overlapping SPH discretization elements plus numerical errors produce implicit subfilter stresses, then reports Taylor-Green Re=10^4 results showing degraded transition upon addition of eddy viscosity. This chain does not reduce any prediction to a fitted parameter by construction, nor does it rely on self-citation load-bearing or imported uniqueness theorems. The attribution of degradation to non-locality is an interpretive claim supported by the observed simulation outcomes rather than a tautological redefinition of the inputs. The derivation remains self-contained against external benchmarks and does not exhibit any of the enumerated circularity patterns.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption SPH discretization elements overlap strongly and numerical errors produce implicit subfilter stresses that act like unresolved turbulence scales.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.