{"id":"d2a75e72-2e53-4875-a2f1-8956ddb4b658","arxiv_id":"2510.17500","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A Roe-type finite volume scheme for a two-dimensional system of nonlocal conservation laws with obstacles converges to the unique entropy solution.","lead":"This paper proves that a numerical scheme for simulating several classes of particles on a conveyor belt with walls is convergent. It extends earlier scalar results to multi-class systems and matches microscopic particle simulations qualitatively.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 2.2 proves well-posedness only for a smoothed Heaviside; no limit as ||H'||∞→∞ is established, so the original discontinuous model (1.2) is not covered.","rationale":"I read the paper as an extension of [14] to systems with obstacles. The proof structure is standard and the estimates in Appendix A are consistent aside from minor typos (e.g., signs in A.4/A.5 and the circular K^c formula in Prop. 3.9). The most important limitation is not an internal error but the mismatch between the model (1.2) with Heaviside and the theorem's smooth-H assumption. All compactness estimates depend on L_H; hence the theorem is only for fixed smooth H. The reader's conditional verdict is appropriate: the paper is a solid analysis of the regularized system, but the abstract's claim about treating the discontinuity is overstated. Other flagged issues (omitted proof of Lemma 3.8, reliance on [9] to appear) are secondary: Lemma 3.8 is plausibly analogous to the cited scalar results, and [9] is a reference gap, not a demonstrated flaw. I therefore agree with the reader's weakest assumption and leave the verdict unchanged.","tokens_in":29871,"tokens_out":14333,"duration_ms":114589,"concrete_test":"Re-run the two-class test of Section 4 with the arctan steepness parameter in (4.1) increased from 50 to 100, 200, and 500 (or with H_δ(u)=1/2+(1/π)arctan(u/δ), δ→0), keeping the grid fixed and respecting (3.13). Compute L1 differences between successive density fields and the outflow curves (4.2). If the sequence does not converge as the smoothing steepens, or if max density grows without bound, the smoothed solutions do not have a discontinuous-H limit and Theorem 2.2 cannot be invoked for the original model. An analytic alternative: in a 1D analogue, attempt to prove a BV bound uniform in L_H; a counterexample with BV growing as L_H→∞ would settle the gap.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 2.2 is proved under Assumption (H), which replaces the Heaviside function in (1.2) by a smooth H with L_H = ||H'||∞. The original problem stated in (1.1)-(1.2) uses the discontinuous Heaviside. Every quantitative estimate depends on L_H: the CFL condition (3.13), the L∞ growth constant C_c^∞ in (3.12), and the BV constants K_1^c, K_2^c in (3.16)-(3.17). As the approximation approaches the true Heaviside, L_H → ∞, so Δt → 0 and the a-priori bounds blow up; no compactness argument passes the regularized solutions to a solution of the discontinuous model. Section 1 explicitly admits that the smoothing is needed for stability and that the maximal-density constraint is lost. Thus the central claim does not cover the model 'as originally formulated': the paper proves well-posedness and convergence for a family of regularized models, but provides no limit as the regularization vanishes. The numerical tests (4.1) also use a fixed arctan regularization rather than a sequence converging to H. This is a scope gap rather than an internal inconsistency of Theorem 2.2, but it is load-bearing for the abstract/conclusion claim that the discontinuity in the flux has been 'treated'.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a two-dimensional system of nonlocal conservation laws modeling heterogeneous material flow on conveyor belts, with N particle classes, obstacles encoded via an augmented density, and a dynamic velocity activated by a Heaviside function H. The main result, Theorem 2.2, asserts existence, uniqueness, and BV/L1 estimates for entropy weak solutions obtained as limits of approximate solutions generated by a Roe finite-volume scheme with dimensional splitting (Algorithm 3.1). The proofs rely on positivity, L1, L∞, BV-in-space-time, and discrete entropy estimates, plus a Lipschitz continuous dependence result imported from a companion paper [9]. The numerical section compares two configurations (small/large particle ordering) with microscopic simulations. The central claim is proven only under Assumption (H), which replaces the discontinuous Heaviside by a smooth approximation with bounded derivative L_H, and no limit is established as the smoothing vanishes; the paper explicitly acknowledges that the maximal density constraint is lost.","tokens_in":30229,"tokens_out":2586,"duration_ms":23757,"significance":"If the results hold as stated, the paper is a useful extension of the scalar nonlocal model in [14] to a system of multiple particle classes with obstacles, giving explicit a priori bounds and a finite-volume convergence framework. The technical estimates in Appendix A are elaborate and appear carefully derived. The numerical experiments provide a qualitative validation of the model's behavior. However, the scope of the mathematical theorem is narrower than the abstract and conclusion suggest: well-posedness is proved for a family of smoothed Heaviside models, not for the original discontinuous model (1.1)-(1.2). Two load-bearing ingredients are also deferred: the discrete entropy inequality (Lemma 3.8) is stated without proof, and uniqueness/stability (Proposition 3.9) is imported from an unpublished 'to appear' source. These gaps make the central claim conditional rather than fully established in the submitted manuscript.","major_comments":[{"comment":"The original model (1.1)-(1.2) uses the discontinuous Heaviside H, but Assumption (H) replaces it by a smooth approximation with derivative bounded by L_H. Theorem 2.2 is therefore proved only for the smoothed model. Every quantitative estimate depends on L_H: the CFL condition (3.13), the L∞ growth constant C_c^∞ in (3.12), the BV constants K_1^c, K_2^c in (3.16)-(3.17), and the Lipschitz constant in Proposition 3.9. As the smoothing approaches the true Heaviside, L_H → ∞, so Δt → 0 and the a priori bounds blow up; no compactness argument is given to pass to a solution of the discontinuous model. Section 1 explicitly states that the smoothing is needed for stability and that the maximal-density constraint is lost. Thus the central claim does not cover the model as originally formulated. This is a scope gap rather than an internal inconsistency, but it is load-bearing for the abstract's","section":"§2, Assumption (H); Theorem 2.2; §1, last paragraph"},{"comment":"The discrete entropy inequality is a load-bearing component: it is used to pass to the limit and identify the limit as an entropy solution. However, the proof is omitted with the comment that it is 'entirely analogous' to [2, Proposition 2.8] and [1, Lemma 2.8]. The setting here is not identical: the scheme uses dimensional splitting, a system of N equations, and a nonlocal flux with an additional convolution through H. A rigorous derivation, or at least a precise statement that the cited result applies verbatim to this setting (including the treatment of the extra sgn terms involving v_stat and J), should be provided. Without it, Theorem 2.2's existence claim is incomplete.","section":"§3.5, Lemma 3.8"},{"comment":"Uniqueness in Theorem 2.2 rests entirely on Proposition 3.9, whose proof invokes 'Theorem 2.4 of [9]', a self-authored paper listed as 'to appear' and hence not independently verifiable by the reader. The manuscript does not state the assumptions (ν) and (J) of that theorem nor reproduce its statement, so the reader cannot check that the present setting satisfies them. Since uniqueness is part of the central well-posedness claim, this is a substantive gap. The authors should either include a self-contained proof of the Lipschitz dependence (or a full statement of the external result with all hypotheses) or make the dependence on the unpublished work explicit in the abstract/introduction.","section":"§3.6, Proposition 3.9"},{"comment":"The numerical tests use a fixed arctan approximation of the Heaviside with parameter 50. This is a single regularized model, not a sequence of approximations converging to the discontinuous H. The numerical section therefore cannot serve as evidence for the behavior of the original model, nor does it test the range of validity of the estimates as L_H grows. If the scope of the paper is the regularized model, this is acceptable, but it should be stated clearly; if the intent is to approximate the original model, convergence in the regularization parameter should be examined numerically.","section":"§4.2, equation (4.1)"}],"minor_comments":[{"comment":"The statement says 'for all c∈{0,...,N}' but the proof fixes c∈{1,...,N}; the index set should be {1,...,N}.","section":"§3.3, Lemma 3.4"},{"comment":"The proof defines K_1^c and K_2^c in (3.31)-(3.32) and (3.33)-(3.34), while the theorem statement refers to (3.16)-(3.17). The constants are consistent only after the inequalities 'K_1^c, K_3^c ≤ K_1^c; K_2^c, K_4^c ≤ K_2^c' are used, but the notation is confusing. Please align the definitions or explain the relationship.","section":"§3.4, Proposition 3.5"},{"comment":"In the last inequality of (3.28), the term |ρ^c_{i,j}+ρ^c_{i-1,j}| appears; based on the preceding lines it should be |ρ^c_{i,j}-ρ^c_{i-1,j}|. Please correct the typo.","section":"§3.4, equation (3.28)"},{"comment":"Figure 2 shows time values t=1,2,3 while Figure 3 shows t=100,200,300. The reason for this difference is not explained; if the two cases use different time scales or parameter regimes, this should be stated in the text.","section":"§4, Figures 2 and 3"},{"comment":"The notation 'k' is used in the convolution sums (e.g., ∂_1η_c(x_{i+1/2}-k, y_j-ℓ)) where the summation index h was introduced; this is a minor typo that should be fixed.","section":"Appendix A, (A.7)-(A.11)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically substantial but the advertised result is narrower than claimed. The reliance on an unpublished companion theorem for uniqueness and the omitted proof of the discrete entropy inequality are typical refereeing concerns that can be addressed by including the missing details. The smoothing gap is a scope issue: the authors should either extend the analysis to the discontinuous Heaviside or rewrite the abstract and conclusion to state precisely that the theorem applies to a regularized family. With those changes, the paper would be a solid contribution to the numerical analysis of nonlocal systems."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real contribution here is the extension of the scalar Roe-scheme analysis from [14] to an N-class system in 2D, with obstacles handled through an augmented density. The coupled BV/L∞ estimates and the technical Appendix A are the substance, and they look largely correct. The Roe scheme is genuinely less diffusive than the Lax-Friedrichs approach in [1,9], so this has practical value for conveyor belt modeling. I also give credit for being explicit that the smoothing of H is needed and that the maximal density constraint is lost.\n\nThe main soft spot is the one the stress-test flagged: Theorem 2.2 is proved under Assumption (H), where the discontinuous Heaviside of the original model (1.2) is replaced by a smooth approximation. Every constant—CFL, L∞ growth, BV bounds—blows up as L_H → ∞, and there is no limit argument back to the original problem. The abstract's claim about treating \"the discontinuity in the flux function\" overstates what is actually proved. That is a scope gap, not an internal inconsistency, but it is load-bearing for the paper's central claim.\n\nThe second issue is that Lemma 3.8, the discrete entropy inequality that underpins the convergence proof, is stated with \"the proof is omitted.\" For a convergence paper this is the largest hole. It may well be analogous to [1,2], but a referee needs to see the argument, especially because the Roe flux with dimensional splitting and the boundary terms are exactly where the analogy can break.\n\nThird, Proposition 3.9 (Lipschitz dependence and uniqueness) leans on Theorem 2.4 of [9], a self-authored paper listed as \"to appear.\" That is not independently checkable. The proof here is only a verification of assumptions, so if [9] changes, this paper changes. Also, the constant K^c(t) is left in the abstract form e^{∫K^c(s)ds} without explicit dependence on the norms, which weakens the claimed quantitative estimate.\n\nMinor: numerical tests are qualitative, no code or data are shipped, and the regularization in the tests is a fixed arctan rather than a sequence converging to H. That limits what the numerics say about the original model.\n\nThe paper is worth refereeing. The estimates in Appendix A are elaborate and likely correct, and the system extension is a natural next step. But before acceptance I would want Lemma 3.8 proved or precisely referenced, the scope of the smoothing stated honestly in the abstract and conclusion, and the dependence on [9] resolved. I would not cite it in my own work until those pieces are in place, but I would bring it to a reading group: it is a good example of how far the standard compactness framework pushes, and where the limits are.","headline":"Careful extension of the scalar Roe scheme to nonlocal systems with obstacles, but the headline well-posedness claim covers only a smoothed Heaviside, and a key entropy estimate is asserted without proof.","tokens_in":30700,"tokens_out":1941,"would_cite":false,"duration_ms":19362,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35L65","65M12"],"pacs":[],"model":"deepseek-v4-flash","headline":"For a two-dimensional nonlocal system describing conveyor-belt flow of mixed particle classes, this paper proves existence, uniqueness, and Lipschitz stability of entropy solutions via a convergent Roe finite-volume scheme.","keywords":["nonlocal conservation laws","heterogeneous material flow","conveyor belts","Roe scheme","dimensional splitting","entropy solution","bounded variation","boundary conditions"],"falsifier":"Simulate a two-class conveyor-belt configuration with initial density locally exceeding r_max and take a sequence of Heaviside approximations with increasing steepness; if the approximate solutions fail to converge, or violate the predicted BV bounds, or the time step must shrink to zero, the central well-posedness result for the smoothed family would be contradicted.","tokens_in":29749,"feed_emoji":"📦","tokens_out":4159,"duration_ms":35249,"temperature":0.7,"pith_summary":"This paper extends a scalar nonlocal model of conveyor-belt particle transport to a system of N interacting classes in two space dimensions, with walls and obstacles folded into the nonlocal convolution. The central result is that, for bounded-variation initial data and a smooth approximation of the Heaviside activation in the velocity, the Cauchy problem has a unique entropy solution for every finite time horizon, obtained as the limit of approximate solutions built by a Roe scheme with dimensional splitting. The proof supplies explicit L1, L∞, BV and time-compactness bounds, and a Lipschitz continuous dependence on initial data. A sympathetic reader would care because this is the first well-posedness and convergence proof for the multi-class, obstacle-aware version of the model, and the numerical tests reproduce qualitative mixing effects seen in microscopic particle simulations.","feed_headline":"Multi-class conveyor-belt flow proven to have unique solutions","feed_subtitle":"A convergent Roe scheme with obstacles matches microscopic particle simulations.","key_machinery":"The core machinery is the modified Roe finite-volume scheme with dimensional splitting (Algorithm 3.1), combined with the representation of obstacles as an additional fictitious density class that enters the nonlocal convolution. The dynamic velocity uses a smooth approximation of the Heaviside function with Lipschitz constant L_H; this constant appears in the CFL condition, the L∞ and BV bounds, and the Lipschitz stability estimate, and all estimates degrade as L_H grows. The discrete entropy inequality of Lemma 3.8 upgrades consistency to convergence toward the unique entropy solution.","core_discovery":"The discovery is that the discontinuous Heaviside activation in the dynamic velocity can be replaced by a smooth approximation, and that the resulting regularized system admits a unique entropy weak solution whose discrete approximations are controlled by estimates that are explicit in the smoothing parameter. The convergence argument is carried by a Roe-type scheme with dimensional splitting: the nonlocal velocity is computed once per time step, the flux is split into static and dynamic parts, and the BV and entropy estimates are proved directly on the discrete level. The paper also proves that the solution map is Lipschitz in L1 with respect to initial data, and that the numerical scheme i","pith_inferences":["Because the well-posedness proof requires a smooth Heaviside approximation and the maximal-density constraint is lost, the theorem does not cover the original discontinuous model; a limit passage as L_H→∞ would be needed to close that gap.","The explicit dependence of the CFL and bounds on L_H suggests an adaptive time-stepping strategy: as the smoothing sharpens, the step size must shrink, which could make the scheme expensive for near-discontinuous activations.","The obstacle-as-fictitious-density construction is a reusable device: any region that should be impermeable can be encoded as a high-density class, potentially extending to multi-domain or moving-obstacle settings.","The same Roe-based framework might be transferred to other nonlocal multi-class models in two dimensions, such as pedestrian or traffic flows, where the flux discontinuity arises from a threshold activation."],"forward_implications":["Initial data in (L∞∩BV) produce a unique entropy solution for every T>0, with explicit L1 conservation, L∞ growth, and BV growth bounds.","The Roe scheme with dimensional splitting converges to the entropy solution, giving a practical numerical method for multi-class belt flow with obstacles.","The solution depends Lipschitz continuously on the initial data, so small measurement errors in the initial density lead to controlled errors at later times.","The model, when smoothed, reproduces class-dependent mixing and overtaking effects that match microscopic simulations, suggesting it captures the essential physics of size-heterogeneous cargo."],"fun_headline_variants":["Convergent Roe scheme for heterogeneous conveyor flow","Unique solutions proven for nonlocal conveyor belts","Smoothing regularization ensures unique conveyor flow","Nonlocal model matches particles via convergent Roe scheme","Convergent scheme for nonlocal conveyor material flow"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing assumption is that the Heaviside activation is smoothed into a Lipschitz function; every estimate in the paper scales with the size of its derivative, so if the original discontinuous Heaviside is used, the argument—and the theorem—no longer applies.","fun_headline_variants_meta":{"raw":{"variants":["Convergent Roe scheme for heterogeneous conveyor flow","Unique solutions proven for nonlocal conveyor belts","Smoothing regularization ensures unique conveyor flow","Nonlocal model matches particles via convergent Roe scheme","Convergent scheme for nonlocal conveyor material flow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000989,"raw_usage":{"total_tokens":3945,"prompt_tokens":576,"completion_tokens":3369,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":320,"completion_tokens_details":{"reasoning_tokens":3299}},"tokens_in":320,"tokens_out":3369,"duration_ms":20438,"temperature":1.0,"reasoning_tokens":3299,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T09:01:21.296074+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate a two-class conveyor-belt configuration with initial density locally exceeding r_max and take a sequence of Heaviside approximations with increasing steepness; if the approximate solutions fail to converge, or violate the predicted BV bounds, or the time step must shrink to zero, the central well-posedness result for the smoothed family would be contradicted.","supporting_citations":[],"review_version":1}