{"id":"271ab6b9-14f5-4b98-ad84-daaaa209d291","arxiv_id":"2505.17026","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Numerical simulations of a coupled Burgers-KdV system show that solitons meeting a Burgers bore can be refracted, reflected, or fused, and the system locally returns to the integrable Gardner equation.","lead":"This paper studies a coupled Burgers-KdV system that models how a current interacts with waves, and uses simulations to show that waves meeting a current can be refracted, reflected, or fused. It matters because soliton fusion and reflection are rare behaviors that could be relevant for real-world wave-current interactions such as tidal bores.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Section 3's stated simulation equation omits the factor 3 in the Burgers nonlinearity, so the reported refraction/reflection threshold may be for a different PDE than system (1).","rationale":"The analytical parts of the paper—Proposition 2.1, the m=0 reduction to Gardner, and the exact compound soliton solutions—are internally coherent and do not depend on the disputed numerics. The load-bearing problem is with the evidential basis of the strongest numerical claim. The reader's weakest assumption correctly flags viscosity, resolution, and missing numerical details, but my concern is more specific and is visible directly in the manuscript text: the equation actually simulated is printed with u u_x rather than 3u u_x, and the viscosity coefficient is given as both 0.02 and 0.01. Because the refraction/reflection threshold and the classification outcomes are read off Figure 7, an equation mismatch of this kind would invalidate the central numerical conclusions. Without code or complete numerical specifications, the discrepancy cannot be resolved by inspection; a targeted rerun with both equation forms and both viscosity values is the decisive check. Since the analytical contributions remain sound and the issue is potentially fixable by corrected simulations, the conditional verdict should stand unchanged.","tokens_in":11874,"tokens_out":5503,"duration_ms":60789,"concrete_test":"Reproduce the Figure 7 parameter sweep using the viscous form of the actual system (1), u_t + 3u u_x = ηu_xx - v∂x(3v^2 + γv_xx), at η=0.01 and η=0.02, and compare with the printed simulation equation u_t + u u_x = ηu_xx - v∂x(3v^2 + γv_xx) at the same η values. Classify outcomes by x*(t)=argmax_x v and by the number of outgoing solitons. If the refraction/reflection boundary or the c*≈u0/5 slope changes by more than a few percent between the two equation forms, or between η=0.01 and η=0.02, the threshold claim is not established for system (1).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claims—refraction/reflection classification and the threshold c*≈u0/5—rest entirely on the Section 3 simulations, but the equation stated for those simulations does not match the system under study. System (1) and its derivation in Eq. (30) give u_t + 3u u_x = -v∂x(3v^2 + γv_xx). Section 3 instead adds viscosity 'as u_t + u u_x = ηu_xx - v∂x(3v^2 + γv_xx)', dropping the factor 3. This changes the Burgers shock structure and the bore-front speed, and since the threshold is defined in the bore's moving frame, a missing factor 3 can shift the classification boundary. The same section also sets η=0.02 while the Figure 7 summary caption says η=0.01, and no grid, domain, time-step, or convergence details are reported. If the code followed the printed u_t + u u_x equation, the numerical results are not evidence about system (1); if the code followed system (1), the manuscript misstates the simulated equation. Either way, the numerical support for the strongest claim is not reproducible as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Burgers–swept KdV system (1), a coupled model of a Burgers-type current and a KdV-type wave field. It gives two variational derivations, shows that the momentum variable m = u - v^2/2 is Lie-transported, and observes that when m = 0 the system reduces to the integrable Gardner equation. Exact compound soliton solutions, solitary-wave relations, and periodic solutions are presented. The main numerical claims are that the system spontaneously creates regions with m → 0, that a Burgers bore interacting with compound solitons leads to refraction, reflection, or fusion depending on parameters, and that a threshold c* ≈ u0/5 separates refraction from reflection.","tokens_in":12126,"tokens_out":3616,"duration_ms":38106,"significance":"The variational derivation and the travelling-wave analysis are coherent and provide a useful geometric perspective on a wave–current interaction system. The explicit reduction to the Gardner equation for m = 0 is a clean structural observation, and the proposed classification of bore–soliton interactions is physically interesting if the numerical results are reliable. However, the numerical section as written is not reproducible: the simulated equation differs from system (1) by a factor of 3 in the Burgers nonlinearity, the viscosity parameter is stated inconsistently, and no resolution or convergence information is given. The exact soliton solution is also presented as a special m = 0 reduction rather than as a solution of the fully coupled system, which should be stated more carefully. With the numerical evidence repaired, the paper could make a solid contribution; in its current form the central numerical claims are not established.","major_comments":[{"comment":"The simulation equation is stated as u_t + u u_x = η u_xx − v ∂x(3v^2 + γ v_xx), but the system under study in (1) and in Eq. (30) is u_t + 3u u_x = −v ∂x(3v^2 + γ v_xx). The missing factor 3 changes the Burgers shock structure and the bore-front speed. Since the refraction/reflection threshold is defined in the bore's moving frame, the reported classification may be for a PDE different from system (1). The authors must state exactly which equation was integrated and provide the corresponding code or a detailed numerical specification.","section":"Section 3, first paragraph"},{"comment":"The viscosity coefficient is set to η = 0.02 in the text of Section 3, while the caption of Figure 7 reports η = 0.01. In addition, no grid resolution, domain size, time step, or convergence tests are reported. The threshold c* ≈ u0/5 is read directly from these simulations, and without convergence and error quantification the numerical classification of refraction versus reflection is not supported as written.","section":"Section 3 and Figure 7"},{"comment":"The 'exact compound soliton' of Theorem 2.1 is obtained by invoking Proposition 2.1, i.e., by assuming m = 0 initially, which reduces the system to the Gardner equation. The solution is therefore a Gardner soliton expressed in the (u,v) variables, not a solution of the fully coupled system with general m. This should be stated explicitly as a special reduced solution. In addition, Lemma 2.1 with u0 = 0 gives the condition (u−c)^2(u − v^2/2) = 0; the text selects u = v^2/2 without discussing the u = c branch, which is excluded only by the boundary condition u → 0 at infinity. This exclusion should be justified.","section":"Theorem 2.1 and Lemma 2.1"}],"minor_comments":[{"comment":"Equation (15) appears to have a typo: integrating the travelling-wave equations gives (u(v) + 3v − c)v + γ v_xx = const, but the displayed equation omits the coefficient γ in front of v_xx.","section":"Remark 2.5 / Eq. (15)"},{"comment":"The word 'auxialary' in the first sentence of Appendix A.1 is a typo for 'auxiliary'.","section":"Appendix A.1"},{"comment":"The caption defines x*(t) = argmax_x v(t), but the text does not specify how the refraction/reflection outcome is classified from the time series, nor how the threshold c* is extracted from the dashed lines. A precise classification criterion would help reproducibility.","section":"Figure 7"}],"recommendation":"major_revision","confidential_remarks":"The numerical section is the weakest part of the manuscript: the stated simulated equation does not match system (1), the viscosity parameter is inconsistent between text and figure caption, and no numerical details are provided. These issues are fixable, but they are load-bearing for the paper's central claims, so a major revision is appropriate before the paper can be considered for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The main thing you should know: the analytical work is solid but the numerical centerpiece is currently not reproducible, because the equation printed for the simulations doesn't match the system studied. That needs to be fixed before anyone trusts the refraction/reflection threshold.\n\nWhat's genuinely new is the numerical phenomenology. The variational derivation of the Burgers–swept KdV system is clean, and the Hamiltonian/Lie-transport perspective is a nice way to see how m = u - v^2/2 is advected. The exact compound soliton and the periodic solutions are correct, though they are basically Gardner solutions in disguise once m=0, so the analytical novelty is modest. The new content is the observation that bores can refract, reflect, and fuse compound solitons, with an apparent threshold c* ≈ u0/5. If that holds, it's an interesting wave-current interaction result.\n\nThe soft spot is the numerics. Section 3 states the simulated equation as u_t + u u_x = η u_xx - v∂x(3v^2 + γ v_xx), but system (1) has 3u u_x. You can't just drop the 3 when regularizing a shock. If the code used the printed equation, the results are for a different PDE and the threshold doesn't apply to (1). If the code used (1), then the manuscript misstates the simulated equation. Either way, the numerical claim is not reproducible as written. On top of that, η is given as 0.02 in the text and 0.01 in the Figure 7 caption, and there are no grid, timestep, or convergence checks. So the central classification rests on runs we can't validate.\n\nThe reader's concern about m=0 is fair: the exact solution is a Gardner soliton rearranged, so the 'exact' part is less novel than it looks. But the paper's real claim is about the m≠0 dynamics, and that is independently observed.\n\nWho is this for? Someone working on KdV-type wave-current models, or on near-integrable systems where solitons re-emerge after interactions. They'll get a good variational derivation and an intriguing but as-yet-unverified numerical story.\n\nRecommendation: I'd send it to peer review, but with a request that the authors reconcile the simulated equation with system (1), fix the viscosity inconsistency, and provide enough numerical detail (resolution, convergence) to make the threshold believable. That's a substantial but fixable revision.","headline":"A useful variational derivation and exact solutions, but the simulation section's equation mismatch and lack of numerical details make the central threshold claim unsupported as written.","tokens_in":12648,"tokens_out":2943,"would_cite":false,"duration_ms":28289,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q53","35Q51","37K10","76B15"],"pacs":[],"model":"deepseek-v4-flash","headline":"A coupled Burgers–KdV system drives bore–soliton collisions back toward the integrable Gardner equation.","keywords":["Burgers-swept KdV system","Gardner equation","compound soliton","wave-current interaction","tidal bore","soliton refraction","soliton reflection","soliton fusion"],"falsifier":"A convergence study of the same overtaking collisions at successively smaller viscosity ($\\eta=0.02$, $0.01$, $0.005$) and higher resolution would settle the claim: if the refraction–reflection boundary shifts by more than the numerical error, or if the $m\\to0$ balance breaks down at longer times, the central claim fails. One could also check whether the speed and mass of the escaping fused soliton depend measurably on $\\eta$.","tokens_in":11692,"feed_emoji":"🌊","tokens_out":9951,"duration_ms":87858,"temperature":0.7,"pith_summary":"This paper claims that a shallow-water model coupling the Burgers and KdV equations, obtained by Doppler-shifting KdV into the frame of the Burgers flow, is governed by a mismatch momentum $m=u-\\tfrac12 v^2$. When $m=0$ everywhere, the system is exactly the integrable Gardner equation, and exact compound soliton solutions follow; when $m\\neq0$, numerical simulations show the mismatch relaxing toward zero in localized structures, so Gardner-type solitons emerge spontaneously. Simulations of a Burgers bore overtaking compound solitons identify an approximate threshold $c_*\\approx u_0/5$: slower solitons are refracted into the bore, faster ones are reflected as outgoing compound solitons that carry part of the bore away, and near the threshold two individually refracted solitons can fuse into one faster escaping soliton. The paper matters because it gives a concrete route from a non-integrable wave–current interaction back to an integrable equation, with collision outcomes organized by a single speed parameter.","feed_headline":"A speed threshold decides whether solitons pass a bore or bounce off","feed_subtitle":"Coupled Burgers–KdV simulations show solitons refracting into bores, reflecting off them, or fusing near a threshold.","key_machinery":"The load-bearing object is the momentum mismatch $m=u-\\tfrac12 v^2$, a 1-form density transported by the Burgers velocity $u$ through the Lie-transport equation $m_t+(m\\partial_x+\\partial_x m)u=0$. When $m=0$, the second equation becomes the integrable Gardner equation, so $m=0$ is an invariant integrable submanifold; the numerical claim is that localized structures drive $m$ back toward zero. The derivation uses a Doppler-shifted variational principle that couples the Burgers and KdV Lagrangians, yielding the Hamiltonian form and the coadjoint transport structure. The collision classification is carried by the approximate linear threshold $u_0\\approx 5c$, inferred from pseudo-spectral simulations of the system with viscosity $\\eta u_{xx}$ added to form a smooth bore front.","core_discovery":"The central discovery is that the Burgers–swept KdV system,\n$$u_t+3uu_x=-v\\partial_x($3v^{2}$+\\gamma v_{xx}),\\qquad v_t+6vv_x+\\gamma v_{xxx}=-\\partial_x(uv),$$\nis best read through the Lie-transported momentum $m=u-\\tfrac12 v^2$. If $m=0$ initially, the system reduces exactly to Gardner's equation for $v$, and the paper proves this reduction while giving closed-form compound soliton solutions and periodic travelling waves. Numerically, even with $m\\neq0$ initial data, the system self-organizes into localized regions where $m\\to0$, restoring Gardner dynamics and producing trains of compound solitons. In overtaking collisions, the simulations classify outcomes as refraction (the soliton slows, enters the bore, and leaves a wave train), reflection (the soliton exits the far side as a faster compound soliton carrying bore mass away), or fusion (two solitons merge near the bore front into one escaping soliton). The boundary between refraction and reflection is approximately $u_0\\approx 5c$, i.e. $c_*\\approx u_0/5$, in the small-amplitude regime.","pith_inferences":["Threshold derivation: I would expect the $c_*\\approx u_0/5$ line to follow from a momentum-flux balance at the bore front, so a matched-asymptotics or modulation-theory derivation should reproduce it independently of the numerical viscosity; if it cannot, the threshold is likely a numerical artifact.","Generality to other wave equations: the same Doppler-shifted coupling applied to other integrable equations (modified KdV, nonlinear Schrödinger) should produce analogous refraction–reflection–fusion classifications, providing a sharp test of whether the threshold behaviour is universal.","Integrability recurrence: the local return of $m$ to zero is reminiscent of recurrence of integrable behaviour in an infinite-dimensional system, making this model a concrete laboratory for the infinite-dimensional analogue of Fermi–Pasta–Ulam recurrence.","Experimental signature: the theory predicts that in a controlled flume with a current of height $u_0$, an incident soliton with speed below about $0.2u_0$ (in the bore frame) is absorbed, while one above it escapes the bore as a faster soliton."],"forward_implications":["Localized emergence: if the $m\\to0$ relaxation is generic, the asymptotic outcome of bore–soliton collisions is a train of compound Gardner solitons, opening the way to inverse-scattering methods for the final state.","Predictive threshold: the reported relation $c_*\\approx u_0/5$ means bore height and incident soliton speed alone decide, to leading order, whether an overtaking collision ends in absorption into the bore or escape ahead of it.","Fusion changes soliton number: the bore can merge two solitons that would each be refracted into one reflected soliton, so the outgoing count is not fixed by the incoming count; identical, evenly spaced inputs can yield different numbers of escaping solitons as bore height changes.","Closed-form benchmarks: the exact compound soliton and Jacobi-elliptic periodic solutions provide concrete initial data for testing simulations and reduced models of the interaction.","Exact integrable subcase: whenever $m=0$ initially, Gardner-equation soliton dynamics persist exactly inside the coupled system; the paper's compound solitons are Gardner solitons modulated by an interaction factor $w(\\xi)$."],"supporting_citations":[{"why":"Introduces the Burgers–swept KdV system whose collision dynamics this paper classifies.","marker":"Dombret et al. [2025]"},{"why":"Supplies the Lagrangian-reduction framework used to derive the coupled equations.","marker":"Holm et al. [2023]"},{"why":"Established the integrable Gardner equation, which is the $m=0$ limit of the system.","marker":"[Gardner et al., 1967]"},{"why":"Describes compound Gardner solitons in variable media, the asymptotic wave-train behaviour seen in the simulations.","marker":"[Gorshkov et al., 2012]"},{"why":"Provides the pseudo-spectral numerical method used for all simulations.","marker":"[Burns et al., 2020]"},{"why":"Underlies the Dirac–Frenkel variational coupling of the Burgers and KdV Lagrangians.","marker":"[Frenkel, 1934]"},{"why":"Motivates the Doppler-shifted dispersion relation used in the variational derivation.","marker":"Pizzo and Salmon [2021]"}],"fun_headline_variants":["Solitons near a bore: refract, reflect, or fuse at a speed threshold","Speed threshold decides soliton fate: refraction, reflection, or fusion","Compound solitons: refraction, reflection, and fusion in Burgers–KdV","At u0≈5c, solitons refract, reflect, or fuse — no middle ground","Soliton-bore interaction: pass, bounce, or merge, chosen by speed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The collision classification and threshold $c_*\\approx u_0/5$ are read from finite-time simulations that add viscosity $\\eta=0.02$ to the Burgers equation to smooth the bore front, so the load-bearing premise is that this viscosity and the chosen resolution faithfully represent the dynamics of the inviscid system.","fun_headline_variants_meta":{"raw":{"variants":["Solitons near a bore: refract, reflect, or fuse at a speed threshold","Speed threshold decides soliton fate: refraction, reflection, or fusion","Compound solitons: refraction, reflection, and fusion in Burgers–KdV","At u0≈5c, solitons refract, reflect, or fuse — no middle ground","Soliton-bore interaction: pass, bounce, or merge, chosen by speed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000865,"raw_usage":{"total_tokens":3735,"prompt_tokens":912,"completion_tokens":2823,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":2713}},"tokens_in":528,"tokens_out":2823,"duration_ms":22568,"temperature":1.0,"reasoning_tokens":2713,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:00:38.070056+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A convergence study of the same overtaking collisions at successively smaller viscosity ($\\eta=0.02$, $0.01$, $0.005$) and higher resolution would settle the claim: if the refraction–reflection boundary shifts by more than the numerical error, or if the $m\\to0$ balance breaks down at longer times, the central claim fails. One could also check whether the speed and mass of the escaping fused soliton depend measurably on $\\eta$.","supporting_citations":[{"cited_title":"Holm, Ruiao Hu, Oliver D","cited_arxiv_id":null,"evidence_quote":"Introduces the Burgers–swept KdV system whose collision dynamics this paper classifies."},{"cited_title":"Particle description of the interaction between wave packets and point vortices","cited_arxiv_id":null,"evidence_quote":"Motivates the Doppler-shifted dispersion relation used in the variational derivation."}],"review_version":1}