{"id":"167b286f-1ddb-4c31-9707-67211664a136","arxiv_id":"2507.01736","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The authors introduce OF-EDG, a discontinuous Galerkin scheme with solution-dependent damping and a jump penalty, and prove stability plus a priori error estimates for it.","lead":"This paper combines two established numerical methods to build a high-accuracy scheme for wave equations whose solutions develop sharp jumps. It adds automatic damping and a jump penalty to an energy-based discontinuous Galerkin method, keeping solutions clean at discontinuities without sacrificing smooth accuracy.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Nonlinear nonsmooth tests (Examples 4.4, 4.5, 4.7, 4.8) violate the G>0 assumption behind Remark 4.1's energy estimate; the advertised semilinear oscillation-free claim currently rests on unreported numerics rather than on the paper's theory, and a positive-G analogue would settle it.","rationale":"The reader's conditional verdict is appropriate. I checked the linear stability proof: the flux, penalty, and damping contributions telescope to nonpositive jump squares, so Theorem 2.1 is internally consistent. The a priori error estimate is plausible for smooth solutions and its suboptimality statements match the reported rates. The single load-bearing defect is the semilinear claim. Remark 4.1 sets a sufficient condition G>0 for the inherited energy estimate; all nonlinear discontinuous tests instead use G≤0, and the 2D tests additionally use χ=0 with no analysis. Since the abstract's selling point includes robustness for nonsmooth solutions with nonlinear sources, this is not a peripheral detail: the oscillation-free plots for Examples 4.4-4.5 and 4.7-4.8 are the only evidence in that regime, and they are not covered by the paper's theorem. The proposed positive-G rerun is a minimal, decisive check of whether the behavior is reproducible under the stated hypothesis. If it is, the paper needs only an honest limitation statement; if it is not, the semilinear robustness claim should be withdrawn or re-proved. Either way the current conditional verdict stands.","tokens_in":22441,"tokens_out":10058,"duration_ms":124826,"concrete_test":"Run the discontinuous initial-data experiment of Example 4.5 with the same p=2, q=1 mesh, time step, and OF-EDG parameters, but replace the source by g(u)=-4u^3 so that G(u)=u^4>0 and the test lies inside the hypothesis of Remark 4.1. Compare the computed profile at t=0.25 with the current Figure 4.7 result: if oscillations appear, the semilinear nonsmooth claim fails even under the paper's own assumptions; if the profile remains oscillation-free, the advertised behavior is reproducible inside the proven regime and the negative-G examples should be explicitly reclassified as outside the theory.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing gap is the mismatch between the semilinear stability hypothesis and the discontinuous nonlinear experiments. Remark 4.1 states that the χ=1 extension inherits an energy estimate from [3] only when G(u):=-∫_0^u g(z)dz>0, and the discrete energy (4.56) contains G(u_h). Every nonsmooth nonlinear test violates this: Example 4.4 with g(u)=160 sin u gives G=160(cos u-1)≤0; Example 4.5 with g(u)=4u^3 gives G=-u^4≤0; Examples 4.7-4.8 use the same forms in 2D and, moreover, scheme (4.55) is run with χ=0, for which no stability statement is given. The paper never flags these violations. Hence the central claim that OF-EDG is robust for nonsmooth solutions of wave equations 'with nonlinear source terms' has no theoretical support in exactly the regime advertised; it rests on visual agreement with a CTCS reference, with no code, no S-flux parameter, and no parameter-tuning record. The linear stability theorem (2.1) and the smooth nonlinear convergence test (Example 4.2, where G=1-cos u≥0) are not affected.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes an oscillation-free energy-based discontinuous Galerkin (OF-EDG) method for the second-order wave equation. The spatial discretization augments the EDG scheme of Appelö and Hagstrom with an interior-penalty jump term and with projection-based damping terms whose coefficients depend on jumps of the numerical solution. The authors prove a semi-discrete energy-stability estimate for general fluxes in one and two dimensions, derive a priori error estimates under H^{p+1}/H^{q+1} regularity, and present numerical experiments for smooth and discontinuous solutions, including semilinear sine-Gordon and Klein-Gordon source terms. The paper claims optimal convergence for smooth solutions and oscillation-free behavior for nonsmooth solutions, with no fitted parameters used to produce the reported rates.","tokens_in":22755,"tokens_out":8179,"duration_ms":90544,"significance":"If the claims hold, the method is a practically attractive high-order DG scheme for wave propagation with discontinuous data: the stability proof is clean and self-contained, the damping and penalty design is clearly motivated, and the smooth convergence tests against exact solutions support the analysis. The numerical evidence for oscillation suppression is visually convincing. However, the paper's central robustness claim for nonlinear source terms currently rests on examples that violate the stated semilinear stability hypothesis, and the error-analysis proof contains a regularity gap; these issues make the overall significance conditional until repaired.","major_comments":[{"comment":"The semilinear stability statement in Remark 4.1 is invoked for numerical tests whose source terms do not satisfy its hypothesis G(u)>0. For g(u)=160 sin(u) one has G(u)=160(cos(u)-1)<=0, and for g(u)=4u^3 one has G(u)=-u^4<=0; thus the energy estimate (4.56) is not available for Examples 4.4, 4.5, 4.7, or 4.8. Moreover, the two-dimensional versions are run with chi=0, for which no stability statement is given anywhere in the paper. Since the abstract advertises robustness for nonsmooth solutions with nonlinear source terms, the paper should either add discontinuous semilinear tests satisfying G>0, extend the analysis to chi=0, or explicitly state that the nonlinear nonsmooth results are empirical only.","section":"Remark 4.1, Eq. (4.56), Examples 4.4-4.5 and 4.7-4.8"},{"comment":"The proof of Theorem 2.2 contains a step that is not justified under the stated H^{p+1} regularity. In (2.27) the authors bound ||(u_x-P_{l-1}u_x)||_{L2(I_j)} first by h^{max(1,l)}|u|_{H^{max(1,l)+1}(I_j)}, then by a quantity involving |∂^{max(1,l)+1}u|∞, and finally drop that L∞ seminorm, writing the result as O(h^{max(1,l)+1/2}) with no factor containing u. Without an additional W^{p+1,∞} assumption or a careful elementwise Sobolev-embedding argument with constants under control, the chain does not follow from H^{p+1} regularity. The same issue appears in the two-dimensional estimate (3.51). Because this chain feeds directly into (2.31) and hence into the Gronwall argument, the a priori error bound as stated is not fully proven.","section":"Section 2.4, Eqs. (2.27)-(2.31); Section 3.3, Eq. (3.51)"}],"minor_comments":[{"comment":"The abstract contains the typo 'apriori error estimates'; it should be 'a priori error estimates'. Also, in the Introduction, 'This rest of the paper' should be 'The rest of the paper'.","section":"Abstract and Section 1"},{"comment":"In the third displayed equation of (4.55), the term 'qP' should be a summation symbol or explicitly written as a sum over l=0,...,q; as printed it is not readable.","section":"Remark 4.1, Eq. (4.55)"},{"comment":"The initial data are written as 'u(t,0)=...', but they are initial conditions in space and should be 'u(x,0)=...'.","section":"Examples 4.4 and 4.5"},{"comment":"The comparison with the CTCS reference is only visual; no quantitative error or mesh-convergence data is reported for the nonsmooth nonlinear tests, and the choice of damping/penalty parameters for those runs is not described. Reporting such details would make the robustness claims easier to reproduce and assess.","section":"Examples 4.6-4.8 and Figure 4.6"},{"comment":"Theorem 2.2 is stated in the energy norm, while the convergence tests report the L2 error of u_h; the indirect relationship between the two is only mentioned in Remark 2.2 and is not used to compare theory with the observed rates.","section":"Section 2.4 and Section 4.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for a numerical analysis journal and the linear stability/error-analysis core is worth publishing after revision. The most serious issue is the mismatch between the semilinear stability hypothesis G>0 and the discontinuous nonlinear experiments; this should be fixed before acceptance, either by adding compliant examples or by providing analysis for chi=0. The regularity gap in the error proof is also fixable within the manuscript's scope. I do not see grounds for rejection, provided the authors address these points."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The combination is new and largely works, but the semilinear claims for nonsmooth data are running ahead of the theory. In short: OF-EDG is a sensible merge of EDG with OF damping plus a jump-based penalty, and the penalty is not decoration—it fixes a genuine degeneracy, since plain EDG cannot even evolve piecewise-constant states whose discontinuities sit on cell interfaces. The stability theorem (2.1) is clean, covers all flux choices and all nonnegative parameters, and the proof is transparent. The error analysis is rougher: the chain (2.27)-(2.31) hides pointwise seminorm constants and effectively assumes more than stated H^{p+1} regularity, but the resulting rates match the numerical results, including on nonuniform meshes. The smooth nonlinear test (Sine-Gordon breather) is a nice touch.\n\nThe soft spot is the nonlinear discontinuous experiments. Remark 4.1 inherits an energy estimate from [3] only when G(u)>0. Example 4.4 uses g(u)=160 sin u, giving G=160(cos u−1)≤0; Example 4.5 uses g(u)=4u^3, giving G=−u^4≤0; the 2D examples use the same forms and, worse, run with χ=0, for which no stability statement is given. The paper never flags this mismatch. The smooth convergence test is unaffected (G=1−cos u≥0), but the advertised oscillation-free behavior for semilinear problems with discontinuities is currently supported only by the plots against a CTCS reference. That is not fatal—the numerics may well be right—but the claim needs to be either backed by a positive-G analogue or explicitly downgraded to empirical.\n\nReproducibility is also thinner than it should be: no code, no reported S-flux parameter, and Remark 4.2 admits the χ=1 scheme needs 'careful adjustment' of damping and penalty parameters for discontinuous cases. That makes the nonsmooth results harder to trust than the smooth ones.\n\nWho is this for? Anyone working on DG for wave propagation, especially with shocks or discontinuous data in acoustics or seismics. It is an incremental but solid contribution with a theorem you can verify in an afternoon. It deserves peer review, not desk rejection, but the referees should push for a resolution of the G>0 gap and for code/data.","headline":"Useful new OF-EDG scheme with a clean linear stability proof, but the nonlinear nonsmooth claims run ahead of the theory.","tokens_in":23276,"tokens_out":2844,"would_cite":true,"duration_ms":28875,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M60","65M12","65M15","35L05"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper develops an oscillation-free energy-based discontinuous Galerkin method for the second-order wave equation, proves energy stability and a priori error estimates for smooth solutions, and gives numerical evidence that the method…","keywords":["discontinuous Galerkin method","wave equation","nonsmooth solution","oscillation-free","energy stability","a priori error estimate","semilinear wave equation","SSP-RK3"],"falsifier":"Compute the discrete energy $E_h(t)$ from (4.56) while running Example 4.5 with $g(u)=4u^3$ and piecewise constant initial data; because $G(u)=-u^4<0$ the hypothesis of the energy estimate fails, so any measured growth of $E_h$ would show the nonlinear examples are outside the theorem's coverage, and a comparison spike in the OF-EDG profile against a fine CTCS reference would refute the oscillation-free claim for that test.","tokens_in":22224,"feed_emoji":"🌊","tokens_out":7907,"duration_ms":85806,"temperature":0.7,"pith_summary":"The paper aims to give the energy-based discontinuous Galerkin method for second-order wave equations the ability to handle discontinuous solutions. It adds two mechanisms to the standard EDG scheme: cell-wise damping terms, scaled by jumps of derivatives, that activate near discontinuities, and a jump penalty on the solution itself that prevents piecewise-constant initial data from stalling. On a periodic interval and on Cartesian meshes in two dimensions the resulting OF-EDG scheme is shown to dissipate the discrete energy $\\int ((\\partial_x u_h)^2 + v_h^2)\\,dx$ for every choice of the flux parameters with nonnegative coefficients, and to satisfy an a priori error estimate in the energy norm. Numerical tests exhibit optimal convergence rates for smooth solutions and oscillation-free profiles for discontinuous linear and semilinear waves. If correct, the method is a high-order DG option for nonsmooth wave propagation that needs no limiter and keeps its stability proof.","feed_headline":"Oscillation-free DG method gains stability proof for nonsmooth waves","feed_subtitle":"Adaptive damping plus a jump penalty tames wave discontinuities while preserving high-order accuracy.","key_machinery":"The workhorse is the semi-discrete OF-EDG scheme (2.11), written in first-order form with $v_h$ approximating $u_t$, using the general EDG flux family that contains the central, alternating, and Sommerfeld fluxes as special choices. Two nonlinear dissipation layers are added per cell: oscillation-free damping terms with coefficients $\\sigma_j^l$ and $\\tilde\\sigma_j^l$ proportional to jumps of the $l$-th derivatives of $u_h$ and $v_h$ at the two cell interfaces, projected through the local $L^2$ projection $P_{l-1}$, and an interior-penalty term with strength $c/h^2$ acting on the jump of $u_h$ itself. The damping terms vanish for smooth cells and switch on near discontinuities, while the penalty term provides the force that piecewise-constant solutions require. The energy identity obtained from testing with $(u_h,v_h)$ reduces all interface and damping contributions to nonpositive squares, which is exactly what produces the stability inequality and the error estimate.","core_discovery":"The central claim is that combining oscillation-free damping with an interior-penalty term inside the EDG framework yields a scheme that is simultaneously energy-stable, high-order accurate, and effectively non-oscillatory for nonsmooth wave solutions. Stability holds semidiscretely for any nonnegative penalty parameter $c$ and any nonnegative flux coefficients $\\tau,\\beta$, because every added term contributes a nonpositive jump or projection residual to the rate of change of $E_h$. The error analysis gives the energy error bound of order $h^{2\\gamma}$ with exponent $\\gamma=\\min(p',q')$, where $p'$ is $p-1$ or $p-1/2$ depending on whether $\\tau=0$ and $q'$ is $q$ or $q+1/2$ depending on whether $\\beta=0$, under the polynomial-degree restriction $p-2\\le q\\le p$. The authors further report that in practice the alternating and Sommerfeld fluxes reach the optimal rate $p+1$ for displacement, while the central flux reaches it for odd polynomial degrees only. For discontinuous data, the penalty term moves piecewise-constant values that pure damping would leave untouched, and the damping term removes the spurious oscillations the penalty alone would create.","pith_inferences":["The paper leaves the discrete-in-time analysis untouched; the reported SSP-RK3 step sizes are empirical matchings of the spatial accuracy, so a fully discrete stability proof and CFL condition are a natural next step rather than a proven property.","Because the damping coefficients are built from jumps of derivatives, they double as a discontinuity indicator; one could plausibly use their magnitude to drive h-refinement or p-adaptivity, an application the authors do not mention.","For nonlinear sources with $G(u)\\le0$, such as the two discontinuous tests in Examples 4.4 and 4.5, the stability theorem does not apply; a plausible working hypothesis is that dissipation from the OF terms, not the semilinear energy, is what keeps those runs bounded, which could be tested by switching the damping and penalty terms off."],"forward_implications":["Stability of the semidiscrete scheme holds for every member of the flux family with $c,\\tau,\\beta\\ge0$, so users can choose an energy-conserving flux such as the alternating flux or a dissipative Sommerfeld flux without losing the energy bound.","For smooth solutions the energy error is $O(h^{2\\gamma})$ with $\\gamma=\\min(p',q')$, and the reported experiments show optimal $p+1$ convergence for the alternating and Sommerfeld fluxes and for the central flux with odd $p$.","For piecewise-constant initial data whose jumps lie on cell interfaces, the damping terms alone do nothing and the penalty term is what drives the scheme to the correct solution; the full OF-EDG combination then suppresses the oscillations near the jumps.","The same stability and error analysis extends to two-dimensional problems on Cartesian meshes, and the numerical examples show the oscillation-free behavior persists there.","For semilinear sources the scheme with $\\chi=1$ inherits an energy-stable treatment under $G(u)>0$, and the two-dimensional runs use $\\chi=0$ with empirically stable results."],"supporting_citations":[{"why":"Supplies the EDG formulation in first-order form with the general flux family that this paper modifies.","marker":"[1]"},{"why":"Provides the semilinear EDG energy estimate and the $G(u)>0$ condition invoked in Remark 4.1 for the nonlinear extension.","marker":"[3]"},{"why":"Contributes the oscillation-free damping mechanism with projection-based coefficients that the scheme adapts to the wave equation.","marker":"[19]"},{"why":"Supplies the interior-penalty idea and its $h$-scaling, used here to move piecewise-constant solutions.","marker":"[14]"},{"why":"Gives the discontinuous semilinear wave test problems that Examples 4.4 and 4.5 are inspired by and compare against.","marker":"[6]"},{"why":"Supports the oscillation-free approach for hyperbolic systems that the present method extends to second-order wave equations.","marker":"[16]"},{"why":"Supplies the interpolation and inverse inequalities that convert projection errors into the stated powers of $h$ in the error estimate.","marker":"[8]"}],"fun_headline_variants":["Damping plus penalty makes EDG stable for nonsmooth waves","Stable EDG for nonsmooth waves: damping plus penalty","EDG with interior penalty and damping yields non-oscillatory waves","High-order EDG stays oscillation-free for discontinuous waves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"For the nonlinear wave tests the stability statement is inherited from a semilinear energy argument that assumes $G(u)=-\\int_0^u g(z)\\,dz>0$, but the test sources $g(u)=160\\sin u$ and $g(u)=4u^3$ violate that condition, and the paper does not acknowledge the mismatch, leaving those oscillation-free results resting on numerical observation rather than on the proved estimate.","fun_headline_variants_meta":{"raw":{"variants":["Damping plus penalty makes EDG stable for nonsmooth waves","Stable EDG for nonsmooth waves: damping plus penalty","EDG with interior penalty and damping yields non-oscillatory waves","High-order EDG stays oscillation-free for discontinuous waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001009,"raw_usage":{"total_tokens":4236,"prompt_tokens":885,"completion_tokens":3351,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":501,"completion_tokens_details":{"reasoning_tokens":3280}},"tokens_in":501,"tokens_out":3351,"duration_ms":26015,"temperature":1.0,"reasoning_tokens":3280,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:45:11.525182+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the discrete energy $E_h(t)$ from (4.56) while running Example 4.5 with $g(u)=4u^3$ and piecewise constant initial data; because $G(u)=-u^4<0$ the hypothesis of the energy estimate fails, so any measured growth of $E_h$ would show the nonlinear examples are outside the theorem's coverage, and a comparison spike in the OF-EDG profile against a fine CTCS reference would refute the oscillation-free claim for that test.","supporting_citations":[{"cited_title":"Appel¨ o and T","cited_arxiv_id":null,"evidence_quote":"Supplies the EDG formulation in first-order form with the general flux family that this paper modifies."},{"cited_title":"Appel¨ o, T","cited_arxiv_id":null,"evidence_quote":"Provides the semilinear EDG energy estimate and the $G(u)>0$ condition invoked in Remark 4.1 for the nonlinear extension."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contributes the oscillation-free damping mechanism with projection-based coefficients that the scheme adapts to the wave equation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the interior-penalty idea and its $h$-scaling, used here to move piecewise-constant solutions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the discontinuous semilinear wave test problems that Examples 4.4 and 4.5 are inspired by and compare against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supports the oscillation-free approach for hyperbolic systems that the present method extends to second-order wave equations."}],"review_version":1}