{"id":"a33d960c-4507-4034-958f-4074ed98d6cb","arxiv_id":"1908.03873","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The paper presents new discontinuous Galerkin schemes for the Ostrovsky-Vakhnenko equation, with proofs of energy stability or Hamiltonian conservation, error estimates, and numerical tests for shock, peakon, cuspon, and loop solutions.","lead":"The authors construct and analyze discontinuous Galerkin numerical schemes for the Ostrovsky-Vakhnenko equation, a model for shallow-water and relaxing-medium waves. The schemes are designed to preserve or dissipate conserved wave quantities and to handle singular waveforms such as peakons, cuspons, and loop solitons.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Proposition 2.1's L2-stability proof contains an algebraic sign error in the entropy dissipation term (2.16), so the paper's foundational stability claim is not established as printed.","rationale":"The reader's conditional verdict remains appropriate. The sign error in (2.16) attacks the core stability proposition, but it is a localized algebraic mistake: replacing the incorrect second term with gamma/2*((v^+)^2 - (v^-)^2) repairs the displayed identity and restores the nonnegativity argument. Numerical experiments in Section 4 also support the stability claim. The reader's identified weakest assumption, the unclosed a priori bound (2.27), is a genuine gap in the error estimates, particularly for k=0 where the derived rate h^{1/2} does not improve on the assumed L2 bound. However, this is a standard bootstrap issue and affects only the convergence theorems, not the primary stability or conservation assertions. Both concerns justify conditional acceptance pending correction but do not warrant changing the reader's verdict.","tokens_in":18422,"tokens_out":24613,"duration_ms":257852,"concrete_test":"Evaluate the two sides of (2.16) with gamma=1, v_h^-=1, v_h^+=3, u_h^-=0, u_h^+=2, alpha=2, so that vhat=v^- and fhat=-1. The printed left-hand side equals -11/3, while the claimed right-hand side equals 16/3. If the corrected term -gamma*vhat*[[v]] + gamma/2*((v^+)^2 - (v^-)^2) - fhat*[[u]] + [[F(u)]] is used, both sides equal 16/3 and Theta is nonnegative. This single algebraic check settles whether the stability proof in the manuscript is valid as written.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central advertised property, L2 stability of the energy-stable DG scheme, is not proved as written. In Proposition 2.1, equation (2.16) defines the interface dissipation as Theta = -gamma*vhat*[[v]] - gamma*(1/2(v^-)^2 + 1/2(v^+)^2) - fhat*[[u]] + [[F(u)]], and then asserts that this equals gamma*(-vhat + {{v}})*[[v]] + (f(xi) - {{f}})*[[u]] + 1/2*alpha*[[u]]^2. These two expressions differ by a sign. Summing the cell identity (2.14) yields the correct interface contribution -gamma*vhat*[[v]] + gamma/2*((v^+)^2 - (v^-)^2) - fhat*[[u]] + [[F(u)]], not the printed second term. With the printed formula, Theta can be negative even in an admissible case: gamma=1, vhat=v^-=1, v^+=3, u^-=0, u^+=2, alpha=2 gives printed Theta = -11/3 while the asserted right-hand side is 16/3. Nonnegativity of the printed Theta therefore does not follow, so the proof of Proposition 2.1 and the foundation of the scheme's energy stability need correction. This is more load-bearing than the error-estimate bootstrap in Section 2.2.4, because it concerns the primary stability claim itself. The bootstrap gap is also real, especially for k=0, but it is secondary.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops discontinuous Galerkin discretizations for the Ostrovsky-Vakhnenko equation (1.1). It proposes two semi-discrete energy-stable DG schemes (Section 2.2), one Hamiltonian-conservative DG scheme (Section 2.3), and proves, respectively, L2 stability, Hamiltonian conservation, and suboptimal h^{k+1/2} error estimates for the energy-stable schemes. For singular solutions, the equation is mapped through a hodograph transformation to a coupled dispersionless system, for which two additional DG schemes are constructed. Numerical experiments cover smooth, shock, peakon, cuspon, and loop-soliton solutions, with convergence tables. The advertised theoretical results are plausible and the numerics are informative, but several load-bearing algebraic steps in the proofs are incorrect as printed.","tokens_in":18781,"tokens_out":17625,"duration_ms":174899,"significance":"If the proofs are repaired, the paper would provide a useful family of structure-preserving DG methods for an equation with peakon, cuspon, and loop solutions. The hodograph-based indirect approach is interesting and the numerical experiments (Tables 4.1-4.5) show clean, credible convergence behavior. The paper also gives implementation details for the schemes, which is valuable for reproducibility. However, the central stability and conservation proofs contain algebraic errors as printed, so the theoretical contribution currently requires substantial revision.","major_comments":[{"comment":"The displayed equality in (2.16) is algebraically false. The correct interface contribution obtained by summing (2.14) over adjacent cells is -gamma vhat[[v_h]] + gamma/2((v_h^+)^2 - (v_h^-)^2) - fhat[[u_h]] + [[F(u_h)]], which is the second expression in (2.16), not the first. With the first expression as printed, Theta need not be nonnegative; for gamma=1, vhat=v^-=1, v^+=3, u^-=0, u^+=2, alpha=2 one gets Theta=-11/3 while the claimed right-hand side is 16/3. Separately, the assertion that (f(xi)-{{f}})[[u]] >= 0 follows from monotonicity is false (for f(s)=s^2/2, u^-=-1, u^+=1 the product is -2/3); the alpha/2[[u]]^2 term is essential for dissipativity. Proposition 2.1 needs a corrected cell identity and a correct proof of nonnegativity of the full dissipation term.","section":"Section 2.2.4, assumption (2.27)"},{"comment":"Theorems 2.3 and 2.4 are stated as unconditional convergence results for k>=0, but both proofs rely on the a priori bound ||u-u_h||_{L2(I)} <= h stated in (2.27), which is never derived; no bootstrap argument is supplied. As printed, the error estimates are conditional on this proximity assumption and do not follow from the scheme alone. The case k=0 also needs attention because the special projections P^+- in (2.23)-(2.24) are not defined for piecewise constants. Please either prove the bound, complete the bootstrap induction, or state the theorems with the assumption and justify its admissibility.","section":"Section 2.3, Proposition 2.5, Eq. (2.59)"},{"comment":"The displayed expression for Theta in (2.59) is asserted to vanish, but this is not true in general. The first four terms cancel by the product-jump identity [[ab]]={{a}}[[b]]+{{b}}[[a]], leaving {{v_h}}[[v_h]] - 1/2(v_h^-)^2 + 1/2(v_h^+)^2 = (v_h^+)^2 - (v_h^-)^2, which is nonzero when v_h is discontinuous (as it generally is in this DG scheme). Thus the algebra in the proof of Hamiltonian conservation, as printed, does not establish (2.60). The derivation in (2.54)-(2.59) should be redone carefully, and any missing boundary cancellations must be identified.","section":"Section 2.3, Proposition 2.5, Eq. (2.59)"}],"minor_comments":[{"comment":"The sentence 'The choice of vhat (2.6)' should refer to Eq. (2.5), since the Lax-Friedrichs flux is defined in (2.6).","section":"Section 2.2.3, after Eq. (2.16)"},{"comment":"The source term is written as 'f = cos 2(x+t)'; this should be printed as cos(2(x+t)) to avoid confusion with the square of the cosine.","section":"Example 4.1"},{"comment":"The identity (2.37) for the bilinear form is asserted without derivation; the projection identities and the treatment of the cross terms (xi_u, xi_v) should be displayed, since this is a central step in the error estimate.","section":"Equation (2.37)"},{"comment":"The introduction promises that the paper is concluded in Section 5, but the concluding remarks appear as Section 4.1 before the references; the numbering should be corrected.","section":"Organization"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the first DG treatment of the Ostrovsky-Vakhnenko equation, and the main numerical claims hold up. The paper builds energy-stable and Hamiltonian-conservative DG schemes in the standard LDG mold, extends the hodograph-transform trick from the authors' short-pulse work to handle cuspons and loops, and backs it with clean convergence tables. I would send it to a referee.\n\nThe catch is in the central stability proof. In Proposition 2.1, equation (2.16) asserts an identity for Theta that is wrong as printed: the first expression differs from the second by a term of size gamma (v^+)^2. Plug in gamma=1, v^-=1, v^+=3, u^-=0, u^+=2, alpha=2 and you get -11/3 versus 16/3. The second expression is the correct interface contribution--derive it from summing (2.12)-(2.13) and you get exactly gamma(-vhat+{{v}})[[v]] + ... --so the stability argument itself is recoverable with a one-line sign correction. But as published, the proof of the paper's advertised property is not valid. That is a mandatory fix, not a stylistic quibble.\n\nSecond soft spot: the error estimates in Theorems 2.3-2.4 use the standard a priori assumption ||u-u_h||<=h to control nonlinear terms, but the stated h^{k+1/2} bound only closes that assumption for k>=1. For k=0 the bootstrap fails, so the theorem as stated for all k>=0 is not actually proved. This is a real gap, though it mirrors a fair amount of the DG literature and the numerics in Section 4 are for P1 and P2.\n\nWhat the paper does well: the numerical experiments are honest and cover smooth, shock, peakon, cuspon, and loop solutions; the convergence rates match the claims; the flux choices are standard; the self-citations are to the authors' own related short-pulse paper, which is relevant, not padding. No code or data is provided, and there is no comparison with the existing finite-difference schemes of Coclite-Ridder-Risebro, which would have strengthened the story.\n\nBottom line: the paper is a useful, competent extension with a fixable flaw in the central proof, and one genuine k=0 gap in the error analysis. It deserves serious peer review. I'd accept it as a referee candidate and require those two corrections before publication.","headline":"The paper is the first DG treatment of the Ostrovsky-Vakhnenko equation with clean numerics for peakons, cuspons, loops and shocks, but the central L2-stability proof contains a misprinted identity that needs a sign correction before the claim is trustworthy.","tokens_in":19260,"tokens_out":10098,"would_cite":true,"duration_ms":90175,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M60","65M12","35Q53","35Q51"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that DG discretizations of the Ostrovsky-Vakhnenko equation can be made $L^2$ energy stable or exactly Hamiltonian-conserving, with convergence rates $O(h^{k+1/2})$ for smooth solutions and a hodograph-transformation…","keywords":["Ostrovsky-Vakhnenko equation","discontinuous Galerkin method","Lax-Friedrichs flux","energy stability","Hamiltonian conservation","hodograph transformation","coupled dispersionless system","soliton solutions"],"falsifier":"Run Scheme 1 on the smooth periodic sine test of Section 4.1 on a sequence of meshes and monitor the semi-discrete energy derivative: any positive value of $\\frac{d}{dt}\\|u_h\\|_{L^2}^2$, or an $L^2$ convergence order consistently below $k+\\frac12$, would contradict Proposition 2.1 or Theorem 2.3. Alternatively, evaluate the cell entropy-production term $\\Theta$ in (2.16) on random admissible DG states; a single negative value refutes the nonnegativity claim that the stability proof relies on.","tokens_in":18228,"feed_emoji":"🌊","tokens_out":11715,"duration_ms":105460,"temperature":0.7,"pith_summary":"The paper builds discontinuous Galerkin (DG) spatial discretizations for the Ostrovsky-Vakhnenko equation, a nonlinear wave equation used to model high-frequency waves in relaxing media and a limit of rotating shallow-water models. It claims two families of schemes: energy-stable schemes that drive the $L^2$ norm of the numerical solution downward, and a Hamiltonian-conservative scheme that preserves the invariant $H=\\int_I(-\\frac16u^3+\\frac12v^2)\\,dx$ in the semi-discrete setting. For smooth solutions, both energy-stable schemes are proved to converge at order $O(h^{k+1/2})$ in $L^2$ for degree-$k$ polynomials. To handle singular solutions such as cuspon and loop solitons, the paper uses a hodograph transformation to the coupled dispersionless system and constructs two further DG schemes there. If these claims hold, the methods offer structure-preserving long-time integration and a practical computational route to singular soliton profiles.","feed_headline":"DG schemes keep Ostrovsky-Vakhnenko waves stable and energy-exact","feed_subtitle":"Three DG schemes: L2-stable, Hamiltonian-conserving, and a coordinate-change route to cuspon and loop waves.","key_machinery":"The load-bearing mechanism is the first-order reformulation of the equation together with the numerical-flux design at cell interfaces. The Lax-Friedrichs flux $\\hat f(u)=\\frac12(f(u^+)+f(u^-)-\\alpha(u^+-u^-))$ introduces controlled dissipation; combined with the sign-dependent upwind flux $\\hat v_h=v_h^-$ for $\\gamma>0$ or $v_h^+$ for $\\gamma<0$, it makes the cell entropy-production terms nonnegative, yielding $L^2$ stability. The Hamiltonian scheme instead uses central fluxes $\\hat w_h=\\{\\{w_h\\}\\}$ and $\\hat v_h=\\{\\{v_h\\}\\}$, which make the discrete Hamiltonian telescope exactly under summation. The error analysis follows the projection-based technique with the $L^2$ projection and the special projections $P^{\\pm}$, using the a priori closeness assumption (2.27) to control nonlinear flux terms. For singular solutions, the hodograph transformation recasts the OV equation as the coupled dispersionless system $q_s=u_y$, $u_{ys}+\\gamma q u+c(1-q)=0$, and DG discretizations of this system provide the cuspon and loop profiles.","core_discovery":"The central discovery is that the Ostrovsky-Vakhnenko equation $(u_t+uu_x)_x+\\gamma u=0$, rewritten as the first-order system $u_t+(\\frac12u^2)_x+\\gamma v=0$, $v_x=u$, admits DG discretizations with deliberately chosen numerical fluxes that inherit its two invariants. Scheme 1 and Scheme 2 use the Lax-Friedrichs flux for $f(u)=\\frac12u^2$ and an upwind flux for $v$ chosen by the sign of $\\gamma$; both satisfy $\\frac{d}{dt}\\|u_h\\|_{L^2}^2\\le 0$. Scheme 3 is a Hamiltonian-conservative variant with central fluxes $\\hat w_h=\\{\\{w_h\\}\\}$ and $\\hat v_h=\\{\\{v_h\\}\\}$ that makes $\\frac{d}{dt}\\int_I(-\\frac16u_h^3+v_h^2)\\,dx=0$. For smooth periodic solutions, Theorems 2.3 and 2.4 give the suboptimal error bound $\\|u-u_h\\|_{L^2}\\le Ch^{k+1/2}$ for both energy-stable schemes. Finally, the hodograph transformation $dx=\\rho^{-1}dy+u\\,ds$, $dt=ds$ maps singular solutions to the coupled dispersionless system, and DG schemes for that system reproduce cuspon and loop soliton solutions.","pith_inferences":["Inference: the convergence proof rests on the unproved bootstrap assumption (2.27); completing it would require showing the discrete solution remains within $O(h)$ of the exact solution, for instance by induction on time using the stability estimate.","Inference: because the Hamiltonian scheme uses non-dissipative central fluxes, it will likely not control shock formation, so the energy-stable schemes remain the appropriate choice for shock computations.","Inference: the hodograph route should extend to other short-wave and peakon equations with bi-Hamiltonian structure, provided the transformation stays nonsingular; wave-breaking events where $\\rho$ vanishes would be the natural limitation.","Inference: a direct numerical monitor of the quantity $\\|u-u_h\\|_{L^2}/h$ on the smooth test case would quickly reveal whether the a priori assumption is actually satisfied and whether the stated rates are sharp."],"forward_implications":["Smooth solutions of the OV equation can be advanced by degree-$k$ DG schemes with a guaranteed $L^2$ convergence rate of $O(h^{k+1/2})$, and numerical tests indicate the full $O(h^{k+1})$ order in many cases.","The energy-stable schemes keep the discrete $L^2$ energy nonincreasing in the semi-discrete limit, so they are reliable candidates for computations involving shocks and peakons where oscillation control matters.","The Hamiltonian-conservative scheme preserves the discrete Hamiltonian exactly in space, which should suppress long-time phase drift for smooth and peakon solutions.","Cuspon and loop soliton solutions, which are hard for direct methods because the solution graph folds, become computable through the hodograph-transformed coupled dispersionless system, with the integration variant giving one extra order of accuracy for $u$.","The requirement that $v$ satisfy a boundary condition or zero-mean constraint is built into the schemes, so the discrete methods are consistent with the two standard well-posedness settings for the OV equation."],"supporting_citations":[{"why":"supplies the bi-Hamiltonian structure that motivates the conserved energy and Hamiltonian used in the scheme design","marker":"[2]"},{"why":"justifies the boundary and zero-mean constraints on v and provides the entropy-solution context","marker":"[7]"},{"why":"companion finite-difference convergence analysis for the Ostrovsky-Hunter equation that the constraints build on","marker":"[22]"},{"why":"gives the cuspon and loop soliton solutions via hodograph transformations","marker":"[8]"},{"why":"provides the two-component generalization and integrable semi-discrete analogue used for the transformed system","marker":"[9]"},{"why":"supplies the template for DG methods via hodograph transformations for short-pulse-type equations","marker":"[41]"},{"why":"provides the nonlinear-flux error-estimate technique used in Theorems 2.3 and 2.4","marker":"[38]"},{"why":"extends the error-estimate technique to semi-discrete LDG methods for nonlinear convection-diffusion and KdV equations","marker":"[32]"},{"why":"supplies the projection and interpolation estimates used throughout the error analysis","marker":"[4]"}],"fun_headline_variants":["DG method tames Ostrovsky-Vakhnenko shocks and solitons","Energy-stable DG schemes for Ostrovsky-Vakhnenko equation","Three DG schemes for OV waves: stable, conservative, singular","Hodograph route to DG schemes for cuspon and loop solitons","DG discretizations preserve Ostrovsky-Vakhnenko invariants"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The convergence-rate theorems assume that the discrete solution stays within distance $h$ of the exact solution in $L^2$ for small $h$, and the paper does not prove this bootstrap condition, so the rates are conditional on a closeness property that is not derived from the scheme.","fun_headline_variants_meta":{"raw":{"variants":["DG method tames Ostrovsky-Vakhnenko shocks and solitons","Energy-stable DG schemes for Ostrovsky-Vakhnenko equation","Three DG schemes for OV waves: stable, conservative, singular","Hodograph route to DG schemes for cuspon and loop solitons","DG discretizations preserve Ostrovsky-Vakhnenko invariants"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000666,"raw_usage":{"total_tokens":3075,"prompt_tokens":1020,"completion_tokens":2055,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":636,"completion_tokens_details":{"reasoning_tokens":1957}},"tokens_in":636,"tokens_out":2055,"duration_ms":13450,"temperature":1.0,"reasoning_tokens":1957,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:01:11.061309+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run Scheme 1 on the smooth periodic sine test of Section 4.1 on a sequence of meshes and monitor the semi-discrete energy derivative: any positive value of $\\frac{d}{dt}\\|u_h\\|_{L^2}^2$, or an $L^2$ convergence order consistently below $k+\\frac12$, would contradict Proposition 2.1 or Theorem 2.3. Alternatively, evaluate the cell entropy-production term $\\Theta$ in (2.16) on random admissible DG states; a single negative value refutes the nonnegativity claim that the stability proof relies on.","supporting_citations":[{"cited_title":"Brunelli and S","cited_arxiv_id":null,"evidence_quote":"supplies the bi-Hamiltonian structure that motivates the conserved energy and Hamiltonian used in the scheme design"},{"cited_title":"Coclite, J","cited_arxiv_id":null,"evidence_quote":"justifies the boundary and zero-mean constraints on v and provides the entropy-solution context"},{"cited_title":"Ridder and A.M","cited_arxiv_id":null,"evidence_quote":"companion finite-difference convergence analysis for the Ostrovsky-Hunter equation that the constraints build on"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the cuspon and loop soliton solutions via hodograph transformations"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the two-component generalization and integrable semi-discrete analogue used for the transformed system"},{"cited_title":"Discontinuous Galerkin methods for short pulse type equations via hodograph transformations","cited_arxiv_id":"1907.07842","evidence_quote":"supplies the template for DG methods via hodograph transformations for short-pulse-type equations"},{"cited_title":"Zhang and C.-W","cited_arxiv_id":null,"evidence_quote":"provides the nonlinear-flux error-estimate technique used in Theorems 2.3 and 2.4"},{"cited_title":"Xu and C.-W","cited_arxiv_id":null,"evidence_quote":"extends the error-estimate technique to semi-discrete LDG methods for nonlinear convection-diffusion and KdV equations"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the projection and interpolation estimates used throughout the error analysis"}],"review_version":1}