{"id":"1152c314-f2bb-40b8-9a04-69c41e54510a","arxiv_id":"2507.18210","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The paper identifies the zero-order consistency residue, amplified by background pressure, as the root cause of resolution-insensitive numerical damping in conservative SPH simulations of channel and free-surface flows.","lead":"A widely used mesh-free fluid simulation method quietly loses flow speed in long channels and suppresses wave motion in deep tanks. The cause is a residual force that grows when the flow sits under high background pressure, and the best existing correction only partially fixes it.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The free-surface leg of the root-cause claim is not isolated: the standing-wave energy decay is measured with Eq. (8)'s pressure limiter active, and the residue is only shown to be dissipative by assuming its direction; direct residue work is never computed.","rationale":"The paper identifies a real and practically important phenomenon, and the channel-flow evidence (periodic BC removes damping, channel-length and outlet-pressure trends, resolution insensitivity, partial RKGC correction) is coherent. The central claim, however, is that the same residue is the common root cause in free-surface flows, and that leg rests on the standing-wave energy decay. My concern is that this evidence is not causally decomposed: the limiter in Eq. (8) is active precisely in the inviscid standing-wave tests, and no diagnostic computes the work done by the residue. Also, the sign argument from Ref. [29] does not by itself guarantee net dissipation in an oscillating flow; one needs the actual force-velocity correlation. Neither issue is fatal: both are testable by post-processing the released SPHinXsys code, and the paper's dimensional argument and parametric trends remain plausible. The poor fit at H_w = 3.0 (R^2 = 0.2951 in Table 3) further weakens the quantitative water-depth claim. The appropriate disposition is therefore the same conditional acceptance the reader recommended, with the added requirement that the residue power be verified directly. I agree with the reader that this is an addressable gap rather than a fundamental flaw.","tokens_in":17641,"tokens_out":7502,"duration_ms":86187,"concrete_test":"Re-run the §4.2.1 standing-wave cases (A = 0.1 and A = 0.3, H_w = 1.0) with a diagnostic that evaluates, at every time step, the residue acceleration R_i = -(2 p_i/rho_i) sum_j grad W_ij V_j and the instantaneous residue power P_res(t) = sum_i m_i v_i dot R_i. Compare the time integral of P_res with the measured decrease of total mechanical energy (Fig. 4), and also correlate R_i with v_i over several wave periods. If integral P_res dt reproduces the reported energy decay, the residue is the active damping mechanism; if it does not, the standing-wave data are contaminated by the Eq. (8) pressure limiter or by other numerical dissipation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4's standing-wave experiments are the only quantitative evidence that the zero-order residue produces energy decay in free-surface flows, but they do not isolate the residue. The paper itself states in §2.1.1 that the dissipation limiter beta_ij = min(eta max(v_ij dot e_ij, 0), c0) in Eq. (8) 'is activated only in inviscid flow cases, such as the standing wave.' The measured dE/dt is therefore the sum of the residue contribution and this Riemann-solver dissipation, both growing with wave amplitude; no run with beta_ij = 0 or with a separate estimate of the limiter's power is reported. Separately, Eqs. (16)-(18) convert the residue into a dissipation rate by multiplying the residue force by v_i, but this is only a dissipation rate if sum_j grad W_ij V_j is systematically anti-aligned with v_i. That sign is imported from Ref. [29] and is not verified for the oscillating standing-wave flow, where v_i changes sign each half-period; net decay requires a sustained correlation between the residue force and velocity, not merely a fixed directional bias in the kernel sum. Until the residue power is computed directly from the simulation and shown to match the observed energy loss, the free-surface attribution remains underdetermined. The channel-flow evidence is less exposed because beta_ij is inactive there, but the paper's headline claim is that one residue explains both scenarios.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper identifies the \"zero-order consistency residue\" (Eq. 13: 2 p_i Σ_j ∇W_ij V_j) in the conservative SPH pressure-gradient discretization as the common root cause of non-physical numerical damping in two flow classes: gravity-driven free-surface flows (studied via inviscid standing waves) and pressure-driven internal flows (laminar and turbulent channels, plus the FDA nozzle). It argues that background pressure amplifies the residue, presents sensitivity studies of energy decay or velocity loss to wave amplitude, water depth, channel length, outlet pressure, and resolution, and evaluates the reverse kernel gradient correction (RKGC) as a partial remedy. The paper claims the damping is resolution-insensitive and only partially removable by RKGC, with residual velocity losses up to 17.47% in the long channel and 33.7% in the 3D FDA nozzle without correction.","tokens_in":17804,"tokens_out":6052,"duration_ms":63199,"significance":"If the root-cause attribution is correct, the paper offers a unified explanation for a long-observed excessive-dissipation problem in conservative SPH and gives practical guidance on background-pressure selection and correction schemes. The study is systematic, releases its code through the SPHinXsys repository, and its empirical trends (periodic BC removes damping, lower outlet pressure reduces loss, deeper water and larger amplitude increase decay, RKGC mitigates but does not eliminate) are internally consistent. However, the central causal claim hinges on isolating the residue contribution from other numerical dissipation, and that isolation is incomplete for the free-surface leg of the study.","major_comments":[{"comment":"The standing-wave energy-decay experiments in §4.2 are performed while the dissipation limiter β_ij = min(η max(v_ij · e_ij, 0), c_0) is active, because the paper states that this term is activated only in inviscid flow cases such as the standing wave. The measured dE*/dt is therefore the sum of the residue contribution and the Riemann-solver limiter dissipation, and no run with β_ij = 0 or separate estimate of the limiter's power is reported. Since both contributions grow with wave amplitude, the observed exponential decay cannot be uniquely attributed to the zero-order consistency residue. Please add a control simulation with the limiter disabled or directly compute the residue power and show that it accounts for the measured energy loss.","section":"§2.1.1, Eq. (8)"},{"comment":"The dissipative character of the residue is justified by the statement, imported from Ref. [29], that 'the direction of kernel gradient summation is generally opposite to the flow direction.' This directional bias is not verified for the oscillating standing-wave flow, where the velocity changes sign every half-period; a fixed bias cannot by itself produce net decay unless it is correlated with the velocity over the oscillation cycle. The paper never computes the work done by the residue term in the simulations, so the free-surface attribution remains underdetermined and the claim that a single residue explains both the free-surface and channel-flow observations is not fully supported.","section":"§3, after Eq. (13)"},{"comment":"The dimensional analysis concluding that the dimensionless energy dissipation rate is independent of water depth relies on the hydrostatic scaling p_i ∝ ρ g H_w. In the standing-wave test the pressure also contains a dynamic component proportional to ρ g A, so p_i/(ρ g H_w) = 1 + O(A/H_w), which is not exactly constant across the tested depths. The fitted slopes in Table 3 vary by roughly 10%, so the claim that the similarity is 'well explained' by Eq. (18) needs either a more careful pressure decomposition or an explicit demonstration that the dynamic-pressure correction is negligible in the fitted time interval.","section":"§4.2.2, Eqs. (16)-(18)"},{"comment":"The abstract's phrase 'resolution-insensitive' is not what Fig. 19 shows: the maximum centerline velocity loss increases from 14.67% at Nf = 20 to 17.47% at Nf = 80. The defensible statement is that the damping is not reduced by spatial refinement, but the paper should explain why the loss slightly grows with resolution, since an increasing error with Nf is not the behavior one would naively expect from a pure discretization error and is itself informative about the nature of the residue.","section":"§5.2.4, Fig. 19"}],"minor_comments":[{"comment":"In the opening sentence of Section 5, 'we firstly analysis this problem' should read 'we first analyze this problem.'","section":"§5"},{"comment":"The phrase 'we need to time the non-dimensionalized velocity' should be 'multiply by the non-dimensionalized velocity.'","section":"§4.2.2, after Eq. (17)"},{"comment":"In the Conclusion, 'The non-dimensionless energy dissipation rate' should be 'the dimensionless energy dissipation rate.'","section":"§7"},{"comment":"The label 'Inlet buffer' in Fig. 11 is not defined or discussed in the text; please clarify what region this refers to and why it is shown.","section":"Fig. 11"},{"comment":"For Hw = 3.0 the reported R² is 0.2951, which is too low to support the claim in §4.2.2 that the energy-decay curves are 'well consistent with the logarithmic law.' Please discuss this outlier explicitly rather than attributing it only to numerical oscillation.","section":"Table 3"},{"comment":"The paper says the RKGC technique is 'introduced'; since it is adopted from Ref. [34], 'applied' or 'evaluated' would be more accurate.","section":"Abstract and §1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript leans on several same-group references for load-bearing premises: the sign of the kernel-gradient summation (Ref. [29]), the boundedness of the consistency error under TVF (Ref. [31]), and the RKGC implementation (Ref. [34]). I would encourage the editor to ensure those premises are independently validated or, failing that, to require the authors to provide more direct in-paper evidence. The central gap is the missing isolation of the residue contribution in the free-surface tests; if the authors can supply the control experiment or a direct residue-power measurement, the root-cause claim would become much stronger. This is fixable within the scope of the paper, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nWhat you should know: the channel-flow half of this paper is strong and genuinely new. Long VIPO channels lose centerline velocity to a resolution-insensitive damping that periodic BCs eliminate, outlet pressure modulates, and RKGC only partly fixes. I know of no prior report of this channel-length-dependent loss, and the turbulence case shows the same mechanism eating TKE near the inlet. The FDA nozzle numbers (33.7% velocity loss without correction, 2.1% with) are a practical punchline. The dimensional argument for water-depth-independent dimensionless dissipation rates is neat and testable. Credit where due: the paper ships code, reports measured trends rather than fits to the model, and is honest about RKGC's limits.\n\nThe soft spot is the free-surface leg. The standing-wave energy decay is measured with the Eq. (8) pressure-limiter active, and the paper itself says that limiter operates only in inviscid flows. So the measured dE/dt mixes residue contribution with Riemann-solver dissipation. The residue's dissipative direction is imported from Colagrossi et al. and never verified in the oscillating flow, where velocity reverses each half-period; net decay needs a sustained correlation, not just a fixed kernel-sum bias. The paper does not compute residue power directly from the simulation. That is the load-bearing gap. Separately, the Hw = 3.0 fit has R2 = 0.30, yet the paper claims close agreement across depths; that claim needs qualification.\n\nThese are addressable, not fatal. The channel-flow evidence stands on its own: periodic BC removes the effect, outlet pressure trends match the hypothesis, resolution insensitivity shows it is not a discretization error, and RKGC partial recovery is consistent. The paper's central mechanism is plausible and the engineering relevance is clear.\n\nWho should read it: SPH method developers and CFD practitioners using conservative SPH in internal flows or wave tanks. General CFD readers can skip the method detail but should know the long-channel caveat.\n\nRecommendation: yes, send to peer review. The channel-flow result alone justifies referee time. Ask for direct residue-work measurement and a limiter-free run in the standing-wave tests, or at least a bound on the limiter's contribution. That is a revision request, not a rejection.","headline":"The channel-flow part is a solid, reproducible finding, but the free-surface attribution is underdetermined until the residue power is computed directly and the limiter dissipation is isolated.","tokens_in":18439,"tokens_out":1492,"would_cite":true,"duration_ms":19636,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"One residual term in conservative SPH pressure gradients is proposed as the common cause of spurious damping in waves and channels.","keywords":["smoothed particle hydrodynamics","zero-order consistency","kernel gradient correction","background pressure","numerical dissipation","free-surface flow","channel flow","nozzle benchmark"],"falsifier":"Run the inviscid standing wave case with a pressure-gradient discretization that is exactly zero-order consistent, for example a difference-based or fully corrected gradient, while keeping the same Riemann solver, particle distribution, and wave parameters; if the exponential energy decay remains at the same rate, the decay is not caused by the zero-order consistency residue.","tokens_in":17309,"feed_emoji":"🌊","tokens_out":11567,"duration_ms":105561,"temperature":0.7,"pith_summary":"The paper proposes that a single quantity, the zero-order consistency residue, is the common root cause of the non-physical numerical damping seen in conservative smoothed particle hydrodynamics (SPH) simulations of both gravity-driven free-surface flows and pressure-driven channel flows. The residue is the term $2p_i \\sum_j \\nabla W_{ij} V_j$ left over when the conservative pairwise-average pressure gradient fails to reproduce a constant pressure field, and it acts as a spurious dissipative force whose strength grows with local pressure. In numerical tests, the damping grows with water depth and initial wave amplitude in an inviscid standing wave, and with channel length and outlet pressure in long channels, while remaining insensitive to spatial resolution. The paper also shows that the reverse kernel gradient correction reduces but cannot fully remove the effect, leaving measurable velocity loss in the engineering nozzle benchmark. If correct, the result reframes a known accuracy complaint as a consistency defect with a specific mechanism, and it warns that grid refinement alone will not cure the problem.","feed_headline":"SPH's zero-order residue damps waves and channel flows","feed_subtitle":"A single kernel-gradient term grows with depth and channel length; spatial refinement cannot remove it.","key_machinery":"The load-bearing object is the zero-order consistency residue, the second term on the right of Eq. (13), $2p_i \\sum_j \\nabla W_{ij} V_j$. It arises because a pairwise-average pressure gradient cannot vanish on a constant field when the particle distribution is irregular; the kernel-gradient sum is the zero-order consistency error, and its multiplication by the local pressure $p_i$ makes deep water, long channels, and high outlet pressure the dangerous regimes. The paper uses this identity in two arguments: a dimensional analysis of the residue-induced acceleration that predicts the dimensionless energy-dissipation rate is independent of water depth, and a balance condition at the channel mid-section, $1/L - \\sum_j \\nabla W_{ij} V_j = 0$, which separates regimes where the residue dominates the driving pressure gradient. The reverse kernel gradient correction (RKGC), a momentum-conserving rebuild of kernel gradients, is the counter-mechanism used to test the identification.","core_discovery":"The central claim is that the zero-order gradient inconsistency of the strictly conservative SPH formulation is not just an accuracy defect but a source of physical-looking damping. When the conservative gradient is used to reproduce a constant pressure field, an irreducible residual term appears, $(\\nabla p)_i = \\sum_j W_{ij} V_j (\\nabla p)_i + 2p_i \\sum_j \\nabla W_{ij} V_j$, and the second term, named the zero-order consistency residue, dissipates the flow because the kernel-gradient summation runs opposite to the flow direction. In gravity-driven free surfaces, hydrostatic pressure makes the residue grow with water depth; in pressure-driven channels, the inlet pressure required by long domains does the same. The paper reports exponential energy decay in the standing wave, velocity loss that grows with channel length and outlet pressure but not with resolution, and a normalized energy-decay rate that is independent of water depth, as the depth scaling of pressure cancels in the dimensional analysis. The reverse kernel gradient correction (RKGC) mitigates the residue but leaves about 3.5% velocity loss in the longest tested laminar channel and 2.1% in the 3D nozzle benchmark, where the uncorrected simulation loses 33.7% of the axis velocity.","pith_inferences":["A direct testable extension would recompute the inviscid standing wave with a strictly zero-order-consistent pressure gradient and the same Riemann solver; if the exponential decay persists, the attribution to the residue would need revision.","The dimensional result implies a scaling prediction not tested in the paper: at fixed wave shape, the normalized decay rate should scale with the dynamic pressure and the smoothing-length-normalized consistency error, so varying the sound speed or smoothing length should move the decay rate in a predictable way.","Because the velocity loss is resolution-insensitive, standard convergence studies can look satisfactory while the inlet-region momentum balance is wrong; reporting the residue term as a routine diagnostic would expose the defect in future SPH work.","The same mechanism should appear in other conservative particle methods built on pairwise-average gradients and run with large background pressure, making cross-method comparison a way to see whether this is an SPH-specific issue or a general particle-approximation property."],"forward_implications":["In the inviscid standing wave, residue-induced energy decay follows a clear exponential law whose rate increases with initial wave amplitude and drops to nearly zero when the amplitude is zero.","In the normalized standing-wave variables, the energy-decay rate is essentially the same across water depths, which the dimensional analysis explains by the cancellation of depth in the hydrostatic and residue pressure scalings.","For velocity-inlet and pressure-outlet channels, maximum centerline velocity loss grows with channel length and outlet pressure, is absent under periodic boundary conditions, and gets slightly worse rather than better with refinement (14.67% to 17.47% as the number of particles across the channel rises from 20 to 80).","Kernel gradient correction reduces but does not eliminate the residue: about 3.5% velocity loss remains at length-to-height 60, turbulent kinetic energy under-prediction near the inlet is mitigated, and the 3D nozzle benchmark drops from a 33.7% uncorrected axis-velocity loss to 2.1% with correction.","In high-background-pressure engineering geometries, the residue cannot be neglected, so the paper advises caution with conservative SPH and suggests simplified physical models when the background pressure is unavoidable."],"supporting_citations":[{"why":"Supplies the rigorous derivation of the residue expression in Eq. (13) and the statement that the kernel-gradient sum opposes the flow.","marker":"[29]"},{"why":"First identifies the residual term that appears when the conservative gradient reproduces a constant pressure field.","marker":"[28]"},{"why":"Introduces the transport velocity formulation used to bound the zero-order consistency error without adding background pressure.","marker":"[10]"},{"why":"Provides the consistency-driven advection scheme the paper relies on to keep the consistency error bounded.","marker":"[31]"},{"why":"Defines the low-dissipation Riemann solver whose discretization is the conservative SPH scheme under study.","marker":"[23]"},{"why":"Supplies the reverse kernel gradient correction (RKGC) tested here as a partial remedy.","marker":"[34]"},{"why":"Provides the WCSPH-RANS formulation and wall function used for the turbulent long-channel simulation.","marker":"[21]"},{"why":"Earlier observation that gravity-wave dissipation increases with water depth, which the paper explains through the residue.","marker":"[15]"},{"why":"Supplies the reference simulation data for the nozzle benchmark used in the 3D comparison.","marker":"[35]"},{"why":"Supplies the experimental velocity data for the nozzle benchmark that the uncorrected and RKGC-corrected runs are measured against.","marker":"[36]"}],"fun_headline_variants":["SPH zero-order residue damps channels and waves","Zero-order residue in conservative SPH causes flow damping","SPH's residue damping is worse with background pressure","SPH's hidden residue damps flow, and fixes are limited","Ignoring SPH's zero-order residue drains channel flow"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The attribution of all the measured standing-wave energy decay to the zero-order consistency residue requires that the dissipation produced by the Riemann solver's pressure limiter, which is active in these inviscid flows, is negligible compared with the residue's effect.","fun_headline_variants_meta":{"raw":{"variants":["SPH zero-order residue damps channels and waves","Zero-order residue in conservative SPH causes flow damping","SPH's residue damping is worse with background pressure","SPH's hidden residue damps flow, and fixes are limited","Ignoring SPH's zero-order residue drains channel flow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001607,"raw_usage":{"total_tokens":6470,"prompt_tokens":1086,"completion_tokens":5384,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":702,"completion_tokens_details":{"reasoning_tokens":5304}},"tokens_in":702,"tokens_out":5384,"duration_ms":36816,"temperature":1.0,"reasoning_tokens":5304,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:37:57.045761+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the inviscid standing wave case with a pressure-gradient discretization that is exactly zero-order consistent, for example a difference-based or fully corrected gradient, while keeping the same Riemann solver, particle distribution, and wave parameters; if the exponential energy decay remains at the same rate, the decay is not caused by the zero-order consistency residue.","supporting_citations":[{"cited_title":"Colagrossi, B","cited_arxiv_id":null,"evidence_quote":"Supplies the rigorous derivation of the residue expression in Eq. (13) and the statement that the kernel-gradient sum opposes the flow."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"First identifies the residual term that appears when the conservative gradient reproduces a constant pressure field."},{"cited_title":"Adami, X","cited_arxiv_id":null,"evidence_quote":"Introduces the transport velocity formulation used to bound the zero-order consistency error without adding background pressure."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the consistency-driven advection scheme the paper relies on to keep the consistency error bounded."},{"cited_title":"Zhang, X","cited_arxiv_id":null,"evidence_quote":"Defines the low-dissipation Riemann solver whose discretization is the conservative SPH scheme under study."},{"cited_title":"Zhang, N","cited_arxiv_id":null,"evidence_quote":"Supplies the reverse kernel gradient correction (RKGC) tested here as a partial remedy."},{"cited_title":"A weakly compressible SPH method for RANS simulation of wall-bounded turbulent flows","cited_arxiv_id":"2501.18397","evidence_quote":"Provides the WCSPH-RANS formulation and wall function used for the turbulent long-channel simulation."},{"cited_title":"Colagrossi, A","cited_arxiv_id":null,"evidence_quote":"Earlier observation that gravity-wave dissipation increases with water depth, which the paper explains through the residue."},{"cited_title":"Huang, R","cited_arxiv_id":null,"evidence_quote":"Supplies the reference simulation data for the nozzle benchmark used in the 3D comparison."},{"cited_title":"Hariharan, M","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental velocity data for the nozzle benchmark that the uncorrected and RKGC-corrected runs are measured against."}],"review_version":1}