REVIEW 4 major objections 7 minor 1 cited by
The efficient implementation of transport velocity formulation
T0 review · 4 major / 7 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read SPH clumping fixed by smoothing-length transport correction
desk verdict Useful SPH limiter and broad benchmarks, but Eq. (14)'s coefficient is off by a factor of 2.19 and the central correction is effectively calibrated, not derived. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the limiter $\beta$ of Eq. (15), a switch built on the scalar disorder measure $h^2 |\sum_j \nabla W_{ij} V_j|^2$, where $h=h_{\min}$ in variable-resolution regions. When the measure is above the threshold $5\times10^{-4}$, $\beta = \min(C \cdot \text{measure}, 1)$ with $C=10^3$; otherwise $\beta=0$. This multiplier sits in front of the smoothing-length-scaled transport displacement $\Delta \mathbf{e}_r = -0.2\, \beta\, h_{\min}^2 \sum_j \nabla W_{ij} V_j$, so the same formula handles both uniform and disordered particle clouds without needing a background pressure.
What would settle it
Run a low-velocity variable-resolution simulation, such as a lid-driven cavity with a fine region embedded in a coarse region and maximum velocity below $10^{-3}$, and check whether the limiter stays off in the uniform coarse region while particles still develop clumps or an energy plateau appears. If overcorrection or clumping persists outside the tested parameter range, the threshold and gain are not general.
Extended reading notes
Core claim
On the paper's own terms, the transport-velocity correction need not depend on a global background pressure. The authors derive, from the background-pressure formulation of Adami et al. under dual-criteria time stepping, that the correction per advection step is a displacement proportional to $h^2 \sum_j \nabla W_{ij} V_j$ with coefficient $-0.2$. They then bound this correction with a limiter $\beta$ that is zero unless $h^2|\sum_j \nabla W_{ij} V_j|^2$ exceeds $5\times10^{-4}$, and saturates at 1 with gain $C=10^3$. The central discovery is that this disorder-triggered, smoothing-length-scaled limiter removes the overcorrection visible in low-velocity Taylor-Green flow while preserving the stabilising effect in disordered or variable-resolution regions.
Load-bearing premise
The limiter's fixed threshold and gain constants, tuned on the test cases, are assumed to separate reliably the particle distributions that need transport correction from those that do not, across all velocities and resolutions.
Editorial extensions
If this is right
- Low-velocity flows no longer show the energy decay plateau caused by the original background-pressure formulation.
- The correction scales with smoothing length, so variable-resolution SPH uses the same formula with $h=h_{\min}$ in refined regions.
- The limiter allows transport velocity to be used with dual-criteria time stepping without the sound-speed-dependent background pressure.
- Accuracy benchmarks for cavity flow, FSI, and the FDA nozzle match reference data, so the change does not sacrifice solution quality.
- Two- and three-dimensional tests show the same stabilising effect, so the method generalises beyond planar flows.
Reading between the lines
- The disorder measure could be made dimensionless and adaptive, replacing the fixed $5\times10^{-4}$ threshold with a value tied to kernel shape, which would remove the case-dependent tuning.
- The limiter is likely sensitive to the Wendland kernel and $h=1.3dp$ used here; kernels with different derivative behaviour may need recalibration of $C$ and the threshold.
- If the measure is local and cheap, it could also gate particle-shifting or artificial-viscosity corrections in other SPH variants, not just transport velocity.
- A direct comparison against the particle-shifting technique on the same variable-resolution FSI cases would clarify which strategy better preserves free-surface volume.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes an improved transport-velocity formulation for weakly-compressible SPH. In Section 3, the authors replace the background-pressure-dependent correction of Adami et al. with a smoothing-length-scaled correction, Eq. (16), and introduce a limiter, Eq. (15), that suppresses the correction for well-ordered particle distributions to avoid overcorrection in low-velocity flows. The formulation is claimed to extend the transport-velocity approach to variable-resolution simulations and to prevent particle clumping and void formation without the overcorrection of the original method. Section 4 validates the method on a Taylor-Green vortex, 2D and 3D lid-driven cavity flows, flow-induced vibration of an elastic beam, multi-resolution flow around a cylinder, and the 3D FDA nozzle, with comparisons to references and experiments.
Significance. The practical importance of a transport-velocity correction that scales with smoothing length rather than background pressure is clear: it could make the method applicable to multi-resolution SPH and address the known overcorrection at low velocities. The paper has several strengths: the benchmark suite is broad, including 2D and 3D cases, fluid-structure interaction, variable resolution, and an experimental device; the SPHinXsys code is publicly accessible, which supports reproducibility; and the proposed limiter targets a real limitation of the original formulation. However, the central coefficient in Eq. (16) is not derived from the stated constants in Eq. (14), and the limiter constants are calibrated from numerical cases without a sensitivity study. These issues place the core quantitative claim in need of revision before the result can be relied upon.
major comments (4)
- [Section 3, Eq. (14)] The coefficient ξ = -0.2 is not obtained from the stated formulas. Substituting Eq. (13) with α = 7.0 and Eq. (9) with CFL_ad = 0.25 into Eq. (12), and using m_j = ρ_j V_j with ρ_i ≈ ρ_0, gives Δe_r = -α CFL_ad^2 h^2 Σ_j ∇W_ij V_j = -7 × (0.25)^2 h^2 Σ_j ∇W_ij V_j = -0.4375 h^2 Σ_j ∇W_ij V_j. The printed value -0.2 is smaller by a factor of 2.1875. Since all validations in Section 4 use the coefficient -0.2 of Eq. (16), the derivation in Eq. (14) does not support the implemented correction; an independent re-implementation following the printed constants would apply a much larger correction and would not reproduce the reported results. The authors should either correct the derivation with consistent constants or explicitly state that ξ is a calibrated parameter. This issue is load-bearing because Eq. (16) is the central claim of the paper.
- [Section 3, Eq. (15)] The limiter constants C = 10^3 and the threshold 5 × 10^-4 are introduced 'according to the numerical cases'. Since the limiter is added specifically to remove the overcorrection observed with the original formulation, part of the reported improvement is built into these fitted parameters. The manuscript provides no sensitivity study with respect to C or the threshold, and no scale analysis showing that these values are robust across resolutions, Reynolds numbers, or variable-resolution ratios. Without such evidence, the general-accuracy claim for low-velocity and multi-resolution flows is not fully supported. At minimum, the authors should report results for nearby values of C and the threshold, or derive the constants from dimensional or error considerations.
- [Section 4.2] The statement that the lid-driven cavity results show 'approximate second-order convergence' is not backed by a convergence table or error norms. The figure shows velocity profiles at three resolutions, but no quantitative errors are reported. Please add a table of L1/L2 errors and observed orders for the velocity profiles, or soften the convergence claim.
- [Sections 3 and 4.4] The formulation uses h_min in Eqs. (14)-(16) for variable-resolution flows, but the manuscript does not justify why the minimum smoothing length is the correct length scale for the correction across a resolution interface, nor does it test sensitivity to the resolution ratio or to the location of the interface. Since the central novelty includes variable-resolution applicability, this assumption should be either derived or tested, for example by varying the refinement ratio and measuring differences in the solution.
minor comments (7)
- [Section 1] The word 'stablity' in the introduction should be 'stability'.
- [Section 2.3, Eq. (8)] The text refers to ρ* as the density prior to reinitialization, but Eq. (8) does not define ρ* or show how it enters the summation; please clarify the notation.
- [Section 2.2, Eqs. (4) and (5)] The notation p_L/p_R in Eq. (5) is inconsistent with p_l/p_r in Eq. (4); please unify.
- [Section 4.5] The phrase 'as depicted in the left panel of Figure 12 (left panel)' repeats 'left panel'; remove the duplicate.
- [Section 5] The word 'sceneries' in the conclusion should be 'scenarios'.
- [Table 1] The column header 'Amplitude in y direction/(D)' is unclear; please state explicitly that the amplitude is normalized by D.
- [Figure 13] The caption refers to panels (a) and (b), but the panels are not labeled in the figure; please add labels.
Circularity Check
The limiter constants and the -0.2 coefficient are fitted inputs presented as derived; the claimed derivation in Eq. (14) is arithmetically inconsistent with the printed alpha and CFL_ad values.
-
fitted input called prediction
[Section 3, Eq. (15) and Eq. (16)]
"beta = { min(C h^2 |Σ_j ∇W_ij V_j|^2, 1), h^2 |Σ_j ∇W_ij V_j|^2 > 5 × 10^-4; 0, otherwise } (15) with C = 10^3 according to the numerical cases."
The limiter's threshold (5×10^-4) and gain (C=10^3) are explicitly chosen 'according to the numerical cases', i.e. fitted to the low-velocity and variable-resolution behaviors that the paper then reports as validation. The claim that Eq. (16) avoids overcorrection in the Taylor-Green case is therefore a restatement of the tuning target rather than an independent prediction; the constants encode the desired outcome. Independent benchmark cases (cavity, FSI, FDA) keep the rest of the validation meaningful, but this specific improvement is built into the parameters.
-
other
[Section 3, Eq. (14); Eq. (16); constants in Sec. 2.4 and Eq. (13)]
"Upon inserting Eq. (13) into Eq. (12), the displacement correction ... is formulated as ∆e_r = 1/2 a_t(∆t_ad)^2 = ... = ξ h^2 Σ_j ∇W_ij V_j, (14) with the coefficient ξ = −0.2 accordingly."
The stated inputs are α=7.0 (Eq. 13) and CFL_ad=0.25 (Sec. 2.4). Direct substitution gives ξ = -α CFL_ad^2 = -7 × 0.25^2 = -0.4375, not -0.2. The printed 'ξ = -0.2 accordingly' does not follow from the displayed equations, so the central coefficient in Eq. (16) is an additional fitted value rather than a derived consequence. The central correction formula thus presupposes the number it claims to obtain; an outside implementation using the stated constants would apply a correction 2.19 times larger and could not reproduce the results.
full rationale
The paper is not globally circular: the lid-driven cavity, FSI, flow-around-cylinder, and FDA-nozzle results are checked against external references (Ghia et al., Turek-Hron, experiment/PIV), and the transport-velocity idea is taken from the independent Adami et al. formulation. However, two load-bearing pieces of the claimed improvement are inputs in disguise. First, the limiter constants in Eq. (15) are selected 'according to the numerical cases' and then the same low-velocity/variable-resolution behavior is presented as validation, which is fitted-input-called-prediction circularity for the overcorrection claim. Second, the derivation of Eq. (14) is arithmetically inconsistent with the printed alpha=7.0 and CFL_ad=0.25: the algebra yields xi=-0.4375, so the xi=-0.2 in Eq. (16) is not 'accordingly' derived but is a calibrated/ansatz constant. These issues concern the central correction formula, but because the method still passes external benchmarks and the code is open-source, the paper is only partially circular rather than wholly reducible to its inputs. Score 6.
Assumptions & free parameters
free parameters (4)
- Transport correction coefficient xi =
-0.2
- Limiter gain C =
10^3
- Limiter threshold =
5x10^-4
- Background-pressure coefficient alpha =
7.0
assumptions (4)
- domain assumption Weakly-compressible approximation with artificial EOS p = c^2(rho - rho_0) and c = 10Umax keeps density variation below 1%.
- domain assumption The Riemann-based SPH discretization with the dissipation limiter beta (eta=3) is a faithful discretization of the governing equations.
- domain assumption Density reinitialization via Eq. (8) removes density errors without biasing the transport-velocity result.
- ad hoc to paper Using the minimum smoothing length h_min in the transport correction is valid for variable-resolution flows.
Cite this review
Pith. "Pith review of The efficient implementation of transport velocity formulation." pith.science (2026). https://pith.science/paper/TXP4ZJKZ
@misc{pith2026241113992,
author = {Pith},
title = {Pith review of: The efficient implementation of transport velocity formulation},
year = {2026},
howpublished = {\url{https://pith.science/paper/TXP4ZJKZ}},
note = {Machine review of arXiv:2411.13992}
}
read the original abstract
The standard smoothed particle hydrodynamics (SPH) method suffers from tensile instability, resulting in particle clumping and void regions under negative pressure conditions. In this study, we extend the transport-velocity formulation of Adami et al. (2013) \cite{adami2013transport} in the weakly-compressible SPH (WCSPH) framework to address this long-standing issue. Rather than relying on background pressure, our modified and improved transport-velocity correction scales directly to the smoothing length, making it suitable for variable-resolution flows. Additionally, we introduce a limiter to the new formulation to prevent overcorrection, especially for flow with small velocities. These modifications enhance the general applicability of the transport velocity in fluid dynamics. Numerical tests involving low-velocity and variable-resolution cases demonstrate that the new formulation offers a general and accurate solution for multi-physics SPH simulations.
Figures
Figures from the paper (14 more)
Forward citations
Cited by 1 Pith paper
-
A weakly compressible SPH method for RANS simulation of wall-bounded turbulent flows
A weakly compressible SPH method with adaptive dissipation, limited particle transport, and a constant-y-plus wall treatment achieves convergent k-epsilon RANS simulations of wall-bounded turbulent channel flows.
Reference graph
Works this paper leans on
-
[1]
Adami, X
S. Adami, X. Hu, N. A. Adams, A transport-velocity formulation for smoothed particle hydrodynamics, Journal of Computational Physics 241 (2013) 292–307
2013
-
[2]
A numerical approach to the testing of the fission hypothesis, author=Lucy, Leon B, The astronomical journal 82 (1977) 1013–1024
work page 1977
-
[3]
R. A. Gingold, J. J. Monaghan, Smoothed particle hydrodynamics: theory and application to non-spherical stars, Monthly notices of the royal astronomical society 181 (3) (1977) 375–389
work page 1977
-
[4]
J. J. Monaghan, SPH compressible turbulence, Monthly Notices of the Royal Astronomical Society 335 (3) (2002) 843–852. 26
work page 2002
-
[5]
W. C. Welton, Two-dimensional PDF/SPH simulations of compressible turbu- lent flows, Journal of Computational Physics 139 (2) (1998) 410–443
work page 1998
-
[6]
J. J. Monaghan, Simulating free surface flows with SPH, Journal of computa- tional physics 110 (2) (1994) 399–406
work page 1994
-
[7]
L. D. Libersky, A. G. Petschek, T. C. Carney, J. R. Hipp, F. A. Allahdadi, High strain Lagrangian hydrodynamics: a three-dimensional SPH code for dynamic material response, Journal of computational physics 109 (1) (1993) 67–75
work page 1993
-
[8]
P. Randles, L. D. Libersky, Smoothed particle hydrodynamics: some recent improvements and applications, Computer methods in applied mechanics and engineering 139 (1-4) (1996) 375–408
work page 1996
Show all 50 references
-
[9]
S. M. Longshaw, B. D. Rogers, Automotive fuel cell sloshing under temporally and spatially varying high acceleration using GPU-based Smoothed Particle Hydrodynamics (SPH), Advances in Engineering Software 83 (2015) 31–44
2015
-
[10]
Z. Wang, C. Zhang, O. J. Haidn, X. Hu, An Eulerian SPH method with WENO reconstruction for compressible and incompressible flows, Journal of Hydrody- namics 35 (2) (2023) 210–221
2023
-
[11]
Khayyer, H
A. Khayyer, H. Gotoh, Y. Shimizu, On systematic development of FSI solvers in the context of particle methods, Journal of Hydrodynamics 34 (3) (2022) 395–407
2022
-
[12]
J. J. Monaghan, Smoothed particle hydrodynamics, Reports on progress in physics 68 (8) (2005) 1703. 27
2005
-
[13]
Schuessler, D
I. Schuessler, D. Schmitt, Comments on smoothed particle hydrodynamics, As- tronomy and Astrophysics, vol. 97, no. 2, Apr. 1981, p. 373-379. 97 (1981) 373–379
1981
-
[14]
G. R. Johnson, S. R. Beissel, Normalized smoothing functions for SPH im- pact computations, International Journal for Numerical Methods in Engineering 39 (16) (1996) 2725–2741
1996
-
[15]
Colagrossi, M
A. Colagrossi, M. Landrini, Numerical simulation of interfacial flows by smoothed particle hydrodynamics, Journal of computational physics 191 (2) (2003) 448–475
2003
-
[16]
Monaghan, On the problem of penetration in particle methods, Journal of Computational physics 82 (1) (1989) 1–15
J. Monaghan, On the problem of penetration in particle methods, Journal of Computational physics 82 (1) (1989) 1–15
1989
-
[17]
R. Xu, P. Stansby, D. Laurence, Accuracy and stability in incompressible SPH (ISPH) based on the projection method and a new approach, Journal of com- putational Physics 228 (18) (2009) 6703–6725
2009
-
[18]
X. Xu, P. Yu, A technique to remove the tensile instability in weakly compress- ible SPH, Computational Mechanics 62 (5) (2018) 963–990
2018
-
[19]
S. J. Lind, R. Xu, P. K. Stansby, B. D. Rogers, Incompressible smoothed particle hydrodynamics for free-surface flows: A generalised diffusion-based algorithm for stability and validations for impulsive flows and propagating waves, Journal of Computational Physics 231 (4) (201...
2012
-
[20]
Antuono, A
M. Antuono, A. Colagrossi, S. Marrone, D. Molteni, Free-surface flows solved by means of SPH schemes with numerical diffusive terms, Computer Physics Communications 181 (3) (2010) 532–549. 28
2010
-
[21]
P.-N. Sun, A. Colagrossi, S. Marrone, A.-M. Zhang, The δplus-SPH model: Simple procedures for a further improvement of the SPH scheme, Computer Methods in Applied Mechanics and Engineering 315 (2017) 25–49
2017
-
[22]
Vacondio, B
R. Vacondio, B. D. Rogers, P. K. Stansby, P. Mignosa, J. Feldman, Variable res- olution for SPH: a dynamic particle coalescing and splitting scheme, Computer Methods in Applied Mechanics and Engineering 256 (2013) 132–148
2013
-
[23]
Lyu, P.-N
H.-G. Lyu, P.-N. Sun, Further enhancement of the particle shifting technique: Towards better volume conservation and particle distribution in SPH simula- tions of violent free-surface flows, Applied Mathematical Modelling 101 (2022) 214–238
2022
-
[24]
Zhang, X
C. Zhang, X. Hu, N. A. Adams, A weakly compressible SPH method based on a low-dissipation Riemann solver, Journal of Computational Physics 335 (2017) 605–620
2017
-
[25]
Litvinov, X
S. Litvinov, X. Hu, N. A. Adams, Towards consistence and convergence of con- servative SPH approximations, Journal of Computational Physics 301 (2015) 394–401
2015
-
[26]
Zhang, X
C. Zhang, X. Y. Hu, N. A. Adams, A generalized transport-velocity formulation for smoothed particle hydrodynamics, Journal of Computational Physics 337 (2017) 216–232
2017
-
[27]
Zhang, M
C. Zhang, M. Rezavand, X. Hu, Dual-criteria time stepping for weakly com- pressible smoothed particle hydrodynamics, Journal of Computational Physics 404 (2020) 109135. 29
2020
-
[28]
Zhang, M
C. Zhang, M. Rezavand, Y. Zhu, Y. Yu, D. Wu, W. Zhang, J. Wang, X. Hu, SPHinXsys: An open-source multi-physics and multi-resolution library based on smoothed particle hydrodynamics, Computer Physics Communications 267 (2021) 108066
2021
-
[29]
E. F. Toro, Riemann solvers and numerical methods for fluid dynamics: a prac- tical introduction, Springer Science & Business Media, 2013
2013
-
[30]
X. Y. Hu, N. A. Adams, A multi-phase SPH method for macroscopic and meso- scopic flows, Journal of Computational Physics 213 (2) (2006) 844–861
2006
-
[31]
J. J. Monaghan, Smoothed particle hydrodynamics, Annual review of astronomy and astrophysics 30 (1) (1992) 543–574
1992
-
[32]
Zhang, W
S. Zhang, W. Zhang, C. Zhang, X. Hu, A Lagrangian free-stream boundary condition for weakly compressible smoothed particle hydrodynamics, Journal of Computational Physics 490 (2023) 112303
2023
-
[33]
Wendland, Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree, Advances in computational Mathematics 4 (1995) 389–396
H. Wendland, Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree, Advances in computational Mathematics 4 (1995) 389–396
1995
-
[34]
G. I. Taylor, A. E. Green, Mechanism of the production of small eddies from large ones, Proceedings of the Royal Society of London. Series A-Mathematical and Physical Sciences 158 (895) (1937) 499–521
1937
-
[35]
U. Ghia, K. N. Ghia, C. Shin, High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method, Journal of computational physics 48 (3) (1982) 387–411. 30
1982
-
[36]
Turek, J
S. Turek, J. Hron, Proposal for numerical benchmarking of fluid-structure in- teraction between an elastic object and laminar incompressible flow, Springer, 2006
2006
-
[37]
L. Han, X. Hu, SPH modeling of fluid-structure interaction, Journal of Hydro- dynamics 30 (2018) 62–69
2018
-
[38]
Zhang, M
C. Zhang, M. Rezavand, X. Hu, A multi-resolution SPH method for fluid- structure interactions, Journal of Computational Physics 429 (2021) 110028
2021
-
[39]
Bhardwaj, R
R. Bhardwaj, R. Mittal, Benchmarking a coupled immersed-boundary-finite- element solver for large-scale flow-induced deformation, AIAA journal 50 (7) (2012) 1638–1642
2012
-
[40]
F.-B. Tian, H. Dai, H. Luo, J. F. Doyle, B. Rousseau, Fluid–structure interac- tion involving large deformations: 3D simulations and applications to biological systems, Journal of computational physics 258 (2014) 451–469
2014
-
[41]
F. M. White, J. Majdalani, Viscous fluid flow, Vol. 3, McGraw-Hill New York, 2006
2006
-
[42]
Chiu, R.-K
P.-H. Chiu, R.-K. Lin, T. W. Sheu, A differentially interpolated direct forcing immersed boundary method for predicting incompressible Navier–Stokes equa- tions in time-varying complex geometries, Journal of Computational Physics 229 (12) (2010) 4476–4500
2010
-
[43]
D.-V. Le, B. C. Khoo, J. Peraire, An immersed interface method for viscous incompressible flows involving rigid and flexible boundaries, Journal of Compu- tational Physics 220 (1) (2006) 109–138. 31
2006
-
[44]
Brehm, C
C. Brehm, C. Hader, H. F. Fasel, A locally stabilized immersed boundary method for the compressible Navier–Stokes equations, Journal of Computational Physics 295 (2015) 475–504
2015
-
[45]
Russell, Z
D. Russell, Z. J. Wang, A Cartesian grid method for modeling multiple moving objects in 2D incompressible viscous flow, Journal of Computational Physics 191 (1) (2003) 177–205
2003
-
[46]
H. C. Ku, R. S. Hirsh, T. D. Taylor, A pseudospectral method for solution of the three-dimensional incompressible Navier-Stokes equations, Journal of Com- putational Physics 70 (2) (1987) 439–462
1987
-
[47]
Komaizi, A
D. Komaizi, A. Niknam, Hybrid finite-volume–finite-element scheme for 3D sim- ulation of thermal plasma arc configuration, AIP Advances 14 (1) (2024)
2024
-
[48]
S. F. Stewart, E. G. Paterson, G. W. Burgreen, P. Hariharan, M. Giarra, V. Reddy, S. W. Day, K. B. Manning, S. Deutsch, M. R. Berman, et al., Assess- ment of CFD performance in simulations of an idealized medical device: results of FDA’s first computational interlaboratory stu...
2012
-
[49]
Simulation of the FDA nozzle benchmark: A lattice Boltzmann study, au- thor=Huang, Feng and No¨ el, Romain and Berg, Philipp and Hosseini, Seyed Ali, Computer Methods and Programs in Biomedicine 221 (2022) 106863
2022
-
[50]
Hariharan, M
P. Hariharan, M. Giarra, V. Reddy, S. Day, K. Manning, S. Deutsch, S. Stew- art, M. Myers, M. Berman, G. Burgreen, et al., Experimental particle image velocimetry protocol and results database for validating computational fluid dy- 32 namic simulations of the fda benchmark noz...
2011
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.