Pith. sign in

REVIEW 3 major objections 5 minor 32 references

On Two Conservative HDG Schemes for Nonlinear Klein-Gordon Equation

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read The paper proves optimal O(h^{k+1}) HDG convergence for the nonlinear Klein-Gordon equation, O(h^{k+2}) after post-processing, and a fully discrete scheme that conserves energy while keeping error O(h^{k+1} + (Δt)^2).

desk verdict The semidiscrete HDG analysis and energy conservation are solid and new, but the fully discrete error claims rest on a false lemma and an unstated fourth-order time regularity assumption. read the letter →

arxiv 2411.15572 v1 pith:XY7OBPUP submitted 2024-11-23 math.NA cs.NA

classification math.NAcs.NA MSC 65N1265N22
keywords HDGhybridizablediscontinuousGalerkinKlein-Gordonequationenergyconservationapriorierrorestimatessuperconvergencepost-processingnonlinearwave
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper analyzes two hybridizable discontinuous Galerkin (HDG) schemes for the nonlinear Klein-Gordon equation u_tt − Δu + f(u) = 0 with f(u) = u³ − u, a locally Lipschitz nonlinearity. It proves that the semidiscrete approximation converges at the optimal spatial rate O($h^{{k+1}}$) for both the displacement and the flux when piecewise polynomials of degree k are used, and that a cheap element-by-element post-processing of the displacement improves this to O($h^{{k+2}}$) for k ≥ 1. For the time discretization, a second-order conservative finite-difference scheme is shown to conserve a discrete energy exactly while keeping the overall error at O($h^{{k+1}}$ + (Δt)²); a non-conservative variant that is linear at each time step achieves the same rates. A second HDG variant, with displacement polynomials of degree k+1, gives O($h^{{k+2}}$) for the displacement. Numerical experiments with polynomial degrees 1 to 3 confirm the predicted convergence orders and show energy errors at machine precision.

What carries the argument

The load-bearing machinery is the HDG spatial discretization together with its projection-based error analysis. The method seeks a triple (u_h, q_h, û_h) of piecewise polynomials of degree k, where q_h approximates the flux ∇u, u_h the displacement, and û_h a single-valued trace on the mesh faces; the numerical flux is defined by q̂_h·n = q_h·n − τ(u_h − û_h). The HDG projection (Π_W, Π_V) is defined locally by L² orthogonality conditions and a trace condition on element faces, and it transfers the regularity of the exact solution into sharp O($h^{{k+1}}$) projection errors. The discrete energy identity supplies the boundedness that tames the locally Lipschitz nonlinear term, and for the fully discrete scheme the symmetric quotient F($U^{{n+1}}$,$U^{{n−1}}$) = (F($U^{{n+1}}$) − F($U^{{n−1}}$))/($U^{{n+1}}$ − $U^{{n−1}}$) with F(s) = (1−s²)²/4 makes the discrete energy exactly conserved.

What would settle it

Take the fully discrete conservative scheme on the unit square with the manufactured solution u(x,y,t) = $t^{{3.2}}$ sin(πx) sin(πy), whose fourth time derivative is not square-integrable near t = 0; if the observed temporal order of convergence is still 2 in Δt, the u_tttt hypothesis is not necessary, whereas a drop below 2 would show that the missing hypothesis is real. A complementary test is to solve the same problem on a domain with a re-entrant corner, where the spatial solution is only $H^{{k+1}}$; the postprocessed O($h^{{k+2}}$) rate should degrade if it depends on elliptic regularity.

Watch

Extended reading notes

Core claim

The paper establishes that the HDG framework, which rewrites the second-order wave equation as a first-order system in a displacement variable u_h and a flux variable q_h together with a single-valued trace û_h on element faces, can handle the locally Lipschitz nonlinearity f(u) = u³ − u without any global Lipschitz condition. The key step is the discrete energy identity E(t) = 1/2(‖u_ht‖² + ‖q_h‖² + ‖u_h − û_h‖²₊ + 2(F(u_h),1)), with F(s) = (1−s²)²/4, which is conserved exactly and yields the L^p bounds needed to control the nonlinear term. On the strength of this identity and the HDG projection, the paper proves optimal O($h^{{k+1}}$) error estimates for displacement and flux, superconvergence O($h^{{k+2}}$) after local post-processing for k ≥ 1, and a fully discrete error bound O($h^{{k+1}}$ + (Δt)²) for a conservative time-stepping scheme that preserves the discrete energy exactly. The analysis also covers a non-conservative linear-per-step scheme and a mixed-degree variant, and the claims are backed by numerical experiments.

Load-bearing premise

For the fully discrete error bound to hold, the exact solution must have a square-integrable fourth time derivative, a regularity condition that is used in the Taylor remainder estimates but not stated in the theorem; if the solution is less smooth in time, the (Δt)² part of the claimed convergence is not proven.

Editorial extensions

If this is right

  • The semidiscrete HDG approximation converges at rate O(h^{k+1}) in L² for both displacement and flux, and the locally post-processed displacement converges at rate O(h^{k+2}) for polynomial degree k ≥ 1.
  • The fully discrete conservative scheme conserves the discrete energy exactly, E^{n+1/2} = E^{3/2}, and has error O(h^{k+1} + (Δt)²).
  • The non-conservative variant, which solves only linear systems at each time step, attains the same optimal error rates.
  • The mixed-degree HDG variant, with displacement polynomials of degree k+1 and flux polynomials of degree k, yields O(h^{k+2}) for the displacement and O(h^{k+1}) for the flux.
  • Numerical experiments with polynomial degrees 1 through 3 match the predicted convergence orders and give energy errors at machine precision.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural extension is to apply the same symmetric-quotient energy-conserving time discretization to other Hamiltonian wave equations with a polynomial potential, such as the sine-Gordon equation; the proof structure would likely carry over as long as the potential gives a priori L^p bounds.
  • The gap between the stated regularity in the fully discrete theorem and the regularity actually used in the Taylor remainders suggests that the (Δt)² estimate could potentially be proved under weaker time regularity by adapting the paper's own nonstandard energy technique, and a numerical test with rough-in-time data would reveal whether the missing hypothesis is essential.
  • Because the discrete energy is an exact invariant, the conservative scheme is a candidate for long-time integration of the Klein-Gordon equation; one could test this by measuring phase and amplitude errors over very long intervals.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper develops and analyzes two hybridizable discontinuous Galerkin (HDG) schemes for the nonlinear Klein-Gordon equation u_tt - Δu + f(u) = 0 with f(u) = u^3 - u. The main claims are: optimal O(h^{k+1}) semidiscrete error for both displacement and flux, O(h^{k+2}) superconvergence for a local postprocessed displacement when k ≥ 1, discrete energy conservation for a fully discrete conservative scheme, fully discrete errors of order O(h^{k+1} + (Δt)^2) for both conservative and non-conservative time discretizations, and a variant HDG method with O(h^{k+2}) displacement error. The analysis uses HDG projection estimates, L^p bounds for the discrete solution to handle the locally Lipschitz nonlinearity, a Baker-type argument for the displacement, and a Brouwer fixed-point argument for discrete well-posedness. Numerical experiments for k = 1, 2, 3 and for several exact solutions are reported and confirm the predicted orders and energy conservation.

Significance. If the results are correct, the paper would be a useful contribution to HDG methods for nonlinear hyperbolic problems: it addresses a locally Lipschitz polynomial nonlinearity without imposing a global Lipschitz condition, proves discrete energy conservation, and provides postprocessing superconvergence together with fully discrete second-order time error estimates. The semidiscrete part is structurally substantial, and the L^p technique used to control the nonlinearity is a genuine strength. However, the fully discrete claims currently rest on an unstated fourth-order time regularity assumption and on a false algebraic inequality in Lemma 5.1, while the superconvergence proof for the variant in Section 6 is only sketched. No fitting of parameters or circularity was detected; the concerns are about technical correctness of key proof steps.

major comments (3)
  1. [Section 5.1, Lemma 5.1] The inequality stated in Lemma 5.1 and displayed in Eq. (5.8) is false as written. Taking Φ1 = Ψ1 = M, Φ2 = Ψ2 = 0 gives a left-hand side behaving like ‖M^3 - M‖^2, which is of order M^6 for large M, whereas the right-hand side is of order M^4 since each difference factor contributes M^2 and the bracketed factor contributes M^2. This lemma is used in the uniqueness proof at Eq. (5.16) and in the error estimate at Eq. (5.30) to reduce the nonlinear term to a difference squared in L^6 times a bounded factor. The lemma is repairable by putting a quartic power in the bracketed factor and then using the uniform L^6 bounds from Theorem 5.1, but as stated the estimate is not valid and the subsequent proofs relying on it lack support.
  2. [Theorem 5.3, Eq. (5.28)] The fully discrete error analysis in Theorem 5.3 needs an unstated regularity assumption on the exact solution. In Eq. (5.28), the Taylor remainder terms R*_1 and R*_2 involve u_tttt, and the displayed bound contains ∫_{t_{n-1}}^{t_{n+1}} ‖u_tttt(s)‖^2 ds. The theorem statement, however, imposes no regularity hypotheses on u beyond the semidiscrete context. The same expansion is used in Theorem 5.4. Without an explicit assumption such as u_tttt ∈ L^2(0,T;L^2(Ω)), the right-hand side of (5.28) is not controlled, so the claimed O((Δt)^2) rates in (5.25), (5.26), (5.39), and (5.40) are not proved. This is not a cosmetic omission: second-order accuracy for this implicit midpoint-type scheme genuinely depends on this fourth-order time smoothness, and the theorem must state it.
  3. [Section 6.2, proof of Eq. (6.5)] The proof of the O(h^{k+2}) superconvergence result (6.5) for the variant HDG method is only sketched. After the decomposition into Q1,...,Q8, several terms are bounded, but the text then says that the remaining terms are handled 'systematically' or 'similarly as done in subsection 4.2' without giving the actual estimates. In particular, the nonlinear terms Q4 and Q8 and the time-integrated φ terms require the local Lipschitz structure of f and are not straightforward copies of the earlier argument because the projection and error equations for the variant differ. Since the abstract explicitly advertises the error estimates for this variant, the proof needs to be completed rather than summarized.
minor comments (5)
  1. [Eq. (5.37c)] In the non-conservative scheme, the flux equation (5.37c) writes the test function as v instead of v_h; this should be corrected.
  2. [Proof of Theorem 5.1] In the proof of Theorem 5.1, the sentence 'By adding ... on both sides of (5.37c)' appears to refer to Eq. (5.5), not to (5.37c).
  3. [Lemma 3.1] In the proof of Lemma 3.1, the sentence 'an application of Theorem 2.1 implies' is misleading; the argument used is a Cauchy-Schwarz/Hölder bound, not the energy conservation theorem.
  4. [Theorem 6.1] The statement of Theorem 6.1 says 'u0, u1 ∈ L∞(0,T;H^{k+2})', but u0 and u1 are initial data, not time-dependent functions; the statement should read u0, u1 ∈ H^{k+2}.
  5. [Notation, Section 5] The time-difference notation is overloaded: both ∂_t and overbar-∂_t are introduced in Section 5.1, and the distinction between ∂_tϕ^n and ∂_tϕ^n is easy to miss; a summary table or consistent use of different symbols would improve readability.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular derivation: the error and energy results are proved from external HDG projection theory and standard energy arguments; self-citations are minor and not load-bearing.

full rationale

The paper's central claims—semidiscrete O(h^{k+1}) convergence, postprocessing O(h^{k+2}), discrete energy conservation, and fully discrete O(h^{k+1}+(Δt)^2) error estimates—are not obtained by assuming the quantities they purport to prove. The error analysis uses the standard decomposition e_u=θ_u−ϑ_u, with projection errors ϑ controlled by the external HDG projection results of [11] (Lemma 2.1) and with θ estimated by energy arguments (Lemmas 3.3–3.5, Theorems 3.1–3.2). The fully discrete analysis of Theorem 5.3 proceeds from the same θ-error equations and Taylor expansions; the cited works involving the authors, [24,29,32], are used for standard technique (Baker-type estimates and Taylor remainder bounds) and are not themselves the target results. No parameter is fitted to data, and no quantity is defined in terms of the quantity being predicted. The discrete energy E(t) in Theorem 2.1 is defined from the scheme variables and shown to be conserved by direct differentiation of the HDG equations, not by invoking conservation. The concerns raised by the skeptic are correctness issues rather than circularity: Eq. (5.28) requires u_tttt∈L²(0,T;L²) although Theorem 5.3 states no explicit regularity hypotheses on u, and Lemma 5.1's claimed inequality appears scale-inconsistent and false as written (e.g., Φ1=Ψ1=M, Φ2=Ψ2=0). These may require repair or additional hypotheses, but they do not make the derivation circular. The low nonzero score reflects only the presence of minor self-citations that are not load-bearing.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The central estimates rest on standard finite element machinery (HDG projection, elliptic regularity, duality) and on an unstated fourth-order time regularity assumption; no free parameters are fitted.

assumptions (5)
  • standard math The HDG projection (Pi_V, Pi_W) from [11] is well-defined and has the approximation properties in Lemma 2.1.
    Invoked throughout Section 3 without proof; the paper cites Appendix A of [11].
  • domain assumption On convex polygonal domains, the Poisson problem with L^q data has W^{2,q} regularity.
    Used in Lemma 3.1 to derive L^p bounds for the discrete solution; this requires Omega convex and is standard but unproved.
  • domain assumption The exact solution u is smooth enough to lie in H^{k+2} in space and W^{3,1} or W^{2,2} in time, and f(u) = u^3 - u is locally Lipschitz with u bounded in L^infty.
    Stated in the semidiscrete theorems and used for the mean-value estimates and the dual problem; the bounds fail if u is only H^{k+1}.
  • ad hoc to paper For the fully discrete scheme, u_tttt lies in L^2(0,T;L^2(Omega)) although Theorem 5.3 does not state this.
    The Taylor remainders R*_1 and R*_2 in Eq (5.28) are controlled with fourth time derivatives; without this hypothesis the (Delta t)^2 error bound is not proved.
  • domain assumption Proposition 4.1's H^2 bound for the time-integrated dual solution holds with constant independent of h.
    Assumed in the postprocessing superconvergence analysis; only a sketch of the energy argument is given.

how reviews work

0 comments
Cite this review

Pith. "Pith review of On Two Conservative HDG Schemes for Nonlinear Klein-Gordon Equation." pith.science (2026). https://pith.science/paper/XY7OBPUP

@misc{pith2026241115572,
  author       = {Pith},
  title        = {Pith review of: On Two Conservative HDG Schemes for Nonlinear Klein-Gordon Equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XY7OBPUP}},
  note         = {Machine review of arXiv:2411.15572}
}
abstract

In this article, a hybridizable discontinuous Galerkin (HDG) method is proposed and analyzed for the Klein-Gordon equation with local Lipschitz-type non-linearity. {\it A priori} error estimates are derived, and it is proved that approximations of the flux and the displacement converge with order $O(h^{k+1}),$ where $h$ is the discretizing parameter and $k$ is the degree of the piecewise polynomials to approximate both flux and displacement variables. After post-processing of the semi-discrete solution, it is shown that the post-processed solution converges with order $O(h^{k+2})$ for $k \geq 1.$ Moreover, a second-order conservative finite difference scheme is applied to discretize in time %second-order convergence in time. and it is proved that the discrete energy is conserved with optimal error estimates for the completely discrete method. %Since at each time step, one has to solve a nonlinear system of algebraic equations, To avoid solving a nonlinear system of algebraic equations at each time step, a non-conservative scheme is proposed, and its error analysis is also briefly established. Moreover, another variant of the HDG scheme is analyzed, and error estimates are established. Finally, some numerical experiments are conducted to confirm our theoretical findings.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

  1. [1]

    T. Achouri. Finite difference schemes for the two-dimensional se milinear wave equation. Numerical Methods for Partial Differential Equations , 35(1):200–221, 2019. 29

  2. [2]

    Baccouch

    M. Baccouch. Optimal error estimates of the local discontinuou s Galerkin method for the two- dimensional sine-Gordon equation on cartesian grids. International Journal of Numerical Analysis & Modeling, 16(3):436–462, 2019

  3. [3]

    G. A. Baker. Error estimates for finite element methods for sec ond order hyperbolic equations. SIAM journal on numerical analysis , 13(4):564–576, 1976

  4. [4]

    Brezzi, J

    F. Brezzi, J. Douglas, and L. D. Marini. Two families of mixed finite ele ments for second order elliptic problems. Numerische Mathematik , 47:217–235, 1985

  5. [5]

    G. Chen, B. Cockburn, J. Singler, and Y. Zhang. Superconverg ent interpolatory HDG methods for reaction diffusion equations. Journal of Scientific Computing , 81(3):2188–2212, 2019

  6. [6]

    Cockburn

    B. Cockburn. The hybridizable discontinuous Galerkin methods. I n Proceedings of the International Congress of Mathematicians 2010, (volume 4) Vol. I: Plenary Lectures and Ceremonies Vols. II–IV: Invited Lectures, pages 2749–2775. World Scientific, 2010

  7. [7]

    Cockburn

    B. Cockburn. Hybridizable discontinuous Galerkin methods for se cond-order elliptic problems: overview, a new result and open problems. Japan Journal of Industrial and Applied Mathematics , 40(3):1637–1676, 2023

  8. [8]

    Cockburn, Z

    B. Cockburn, Z. Fu, A. Hungria, L. Ji, M. A. Sanchez, and F.J. Sa yas. Stormer-Numerov HDG methods for acoustic waves. Journal of Scientific Computing , 75(2):597–624, 2018

Show all 32 references
  1. [9]

    Cockburn and J

    B. Cockburn and J. Gopalakrishnan. The derivation of hybridizab le discontinuous Galerkin methods for stokes flow. SIAM Journal on Numerical Analysis , 47(2):1092–1125, 2009

  2. [10]

    Cockburn, J

    B. Cockburn, J. Gopalakrishnan, and R. Lazarov. Unified hybr idization of discontinuous Galerkin, mixed, and continuous Galerkin methods for second order elliptic pro blems. SIAM Journal on Numerical Analysis, 47(2):1319–1365, 2009

  3. [11]

    Cockburn, J

    B. Cockburn, J. Gopalakrishnan, and F.-J. Sayas. A projectio n-based error analysis of HDG meth- ods. Mathematics of Computation , 79(271):1351–1367, 2010

  4. [12]

    Cockburn, W

    B. Cockburn, W. Qiu, and K. Shi. Conditions for superconverge nce of HDG methods for second- order elliptic problems. Mathematics of Computation , 81(279):1327–1353, 2012

  5. [13]

    Cockburn and V

    B. Cockburn and V. Quenneville-B´ elair. Uniform-in-time superc onvergence of the HDG methods for the acoustic wave equation. Mathematics of Computation , 83(285):65–85, 2014

  6. [14]

    Cockburn and K

    B. Cockburn and K. Shi. Conditions for superconvergence of H DG methods for Stokes flow. Math- ematics of Computation , 82(282):651–671, 2013

  7. [15]

    Cockburn, J

    B. Cockburn, J. R. Singler, and Y. Zhang. Interpolatory HDG m ethod for parabolic semilinear PDEs. Journal of Scientific Computing , 79(3):1777–1800, 2019

  8. [16]

    J. F. Currie, J. A. Krumhansl, A. R. Bishop, and S. E. Trullinger. Statistical mechanics of one- dimensional solitary-wave-bearing scalar fields: Exact results and ideal-gas phenomenology. Physical Review B , 22(2):477, 1980

  9. [17]

    D. B. Duncan. Sympletic finite difference approximations of the n onlinear Klein-Gordon equation. SIAM Journal on numerical analysis , 34(5):1742–1760, 1997

  10. [18]

    Jim´ enez and L

    S. Jim´ enez and L. V´ azquez. Analysis of four numerical schem es for a nonlinear Klein-Gordon equation. Applied Mathematics and Computation , 35(1):61–94, 1990

  11. [19]

    Johnson and V

    C. Johnson and V. Thom´ ee. Error estimates for some mixed fin ite element methods for parabolic type problems. RAIRO. Analyse num´ erique, 15(1):41–78, 1981

  12. [20]

    Kapitanski

    L. Kapitanski. Global and unique weak solutions of nonlinear wave equations. Mathematical Research Letters, 1(2):211–223, 1994

  13. [21]

    S. Kesavan. Topics in Functional Analysis and Applications . New Age International Pvt Limited Publishers, 2015. 30

  14. [22]

    R. M. Kirby, S. J. Sherwin, and B. Cockburn. To HDG or to CG: A c omparative study. Journal of Scientific Computing , 51(1):183–212, 2012

  15. [23]

    H. Kragh. Equation with the many fathers. The Klein–Gordon eq uation in 1926. American Journal of Physics , 52(11):1024–1033, 1984

  16. [24]

    A. K. Pani and J. Y. Yuan. Mixed finite element method for a stro ngly damped wave equation. Numerical Methods for Partial Differential Equations: An In ternational Journal , 17(2):105–119, 2001

  17. [25]

    W. Qiu, J. Shen, and K. Shi. An HDG method for linear elasticity with strong symmetric stresses. Mathematics of Computation , 87(309):69–93, 2018

  18. [26]

    Qiu and K

    W. Qiu and K. Shi. An HDG method for convection diffusion equation . Journal of Scientific Computing, 66(1):346–357, 2016

  19. [27]

    S´ anchez, T

    N. S´ anchez, T. S´ anchez-Vizuet, and M. Solano. Error analy sis of an unfitted HDG method for a class of non-linear elliptic problems. Journal of Scientific Computing , 90(3):1–28, 2022

  20. [28]

    S´ anchez, T

    N. S´ anchez, T. S´ anchez-Vizuet, and M. E. Solano. A priori a nd a posteriori error analysis of an unfitted HDG method for semi-linear elliptic problems. Numerische Mathematik , 148(4):919–958, 2021

  21. [29]

    R. K. Sinha and A. K. Pani. The effect of spatial quadrature on fi nite element galerkin approxima- tions to hyperbolic integro-differential equations. Numerical functional analysis and optimization , 19(9-10):1129–1153, 1998

  22. [30]

    Tourigny

    Y. Tourigny. Product approximation for nonlinear Klein-Gordon equations. IMA journal of Nu- merical Analysis, 10(3):449–462, 1990

  23. [31]

    Y. Wong, Q. Chang, and L. Gong. An initial-boundary value proble m of a nonlinear Klein-Gordon equation. Applied Mathematics and Computation , 84(1):77–93, 1997

  24. [32]

    Yadav and A

    S. Yadav and A. K. Pani. Superconvergent discontinuous Galer kin methods for nonlinear parabolic initial and boundary value problems. Journal of Numerical Mathematics , 27(3):183–202, 2019. 31

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.