pith. sign in

arxiv: 2606.25642 · v1 · pith:PCAQH6ROnew · submitted 2026-06-24 · 🧮 math.NA · cs.NA

A hybrid C⁰-interior penalty method for the nematic Helmholtz--Korteweg equation

Tim van Beeck This is my paper

Pith reviewed 2026-06-25 20:35 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords nematic Helmholtz-Korteweg equationhybrid interior penalty methodC0 finite elementsstability analysisconvergence ratesanisotropic PDEacoustic wave propagationliquid crystals
0
0 comments X

The pith

The C0-hybrid interior penalty method for the nematic Helmholtz-Korteweg equation is stable for any polynomial degree at least two when anisotropy is sufficiently small.

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

This paper develops a C0-hybrid interior penalty finite element method to solve the nematic Helmholtz-Korteweg equation, which augments the classical Helmholtz equation with two anisotropic fourth-order terms to model time-harmonic acoustic waves in nematic Korteweg fluids. The discretization uses hybrid spaces and interior penalty terms to handle the higher-order derivatives without requiring global C1 continuity. The authors prove stability for polynomial degrees greater than or equal to two, independent of spatial dimension, by refining earlier analysis with the Cordes condition, and they establish convergence of discrete solutions to the continuous solution under minimal regularity. They also derive convergence rates when the solution has extra smoothness. This flexibility makes the method practical for three-dimensional problems and domains with curved boundaries.

Core claim

The authors introduce a C0-hybrid interior penalty discretization of the nematic Helmholtz--Korteweg equation and prove that it is stable for polynomial degrees greater than or equal to two, independent of the spatial dimension, when the anisotropy is sufficiently small. They further show convergence to the continuous solution under minimal regularity and derive rates for smoother solutions. The discretization provides greater flexibility than C1-conforming methods for applications in three dimensions and on curved domains.

What carries the argument

The C0-hybrid interior penalty discretization, which combines hybrid finite element spaces with interior penalty terms to weakly enforce continuity across element interfaces for the fourth-order anisotropic terms.

If this is right

  • Stability holds for every polynomial degree at least two, independent of spatial dimension, when anisotropy stays below the threshold.
  • Discrete solutions converge to the continuous solution under only minimal regularity assumptions.
  • Convergence rates are obtained whenever the continuous solution possesses additional regularity.
  • The method applies directly to three-dimensional problems and curved domains without needing globally C1 elements.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach could extend to other fourth-order anisotropic equations in materials modeling where C1 elements are impractical to construct.
  • It may support efficient adaptive refinement on complex geometries that arise in liquid-crystal device simulations.
  • The Cordes-based analysis might generalize to time-dependent or nonlinear variants of the same equation family.

Load-bearing premise

The anisotropy parameter must remain below the threshold needed for the Cordes condition to produce the stability estimate.

What would settle it

Numerical evidence of instability or failure of convergence for polynomial degree two in three dimensions at a small positive anisotropy value would contradict the claimed stability.

read the original abstract

The nematic Helmholtz--Korteweg equation models the propagation of time-harmonic acoustic waves in nematic Korteweg fluids, such as nematic liquid crystals. The PDE augments the classical Helmholtz equation with two additional fourth-order terms, one of which is anisotropic in the direction of the nematic field. We refine the previous continuous analysis of Farrell et al. (2025) by using the Cordes condition and present a $C^0$-hybrid interior penalty discretization. The proposed discretization offers greater flexibility than $C^1$-conforming methods and is well-suited for applications in three dimensions and on curved domains. We prove stability of the method for any polynomial degree greater than or equal to two, independent of the spatial dimension, provided that the anisotropy is sufficiently small. Further, we show that the sequence of discrete solutions converges to the continuous solution under minimal regularity assumptions and derive convergence rates if the continuous solution has additional regularity. Finally, we illustrate the capabilities of the method through numerical examples.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript introduces a C⁰-hybrid interior penalty finite element discretization for the nematic Helmholtz–Korteweg equation, which augments the Helmholtz equation with anisotropic fourth-order terms. Building on prior continuous analysis, it invokes the Cordes condition to establish stability of the method for any polynomial degree p ≥ 2, independent of spatial dimension, when the anisotropy parameter is sufficiently small. Discrete solutions are shown to converge to the continuous solution under minimal regularity assumptions, with convergence rates derived under additional regularity; the claims are illustrated by numerical examples.

Significance. If the stability and convergence results hold with the claimed uniformity, the work supplies a flexible, implementable discretization well-suited to three-dimensional problems and curved domains, offering a practical advance over C¹-conforming methods for modeling time-harmonic waves in nematic fluids. The Cordes-based analysis and the dimension-/degree-independent stability (when anisotropy is controlled) would constitute a substantive contribution to the numerical analysis of fourth-order anisotropic PDEs.

major comments (2)
  1. [Abstract / stability analysis] Abstract and stability theorem (presumably §4 or §5): the central claim that stability holds for any p ≥ 2 and any spatial dimension rests on the Cordes condition producing coercivity when anisotropy is below an unspecified threshold. No explicit form of this threshold is supplied, nor is it shown that the resulting coercivity constants remain independent of dimension and degree; without these details the uniformity assertion cannot be verified and is load-bearing for the main result.
  2. [Convergence section] Convergence analysis (presumably §5): the passage from stability to convergence under minimal regularity is asserted, but the argument must explicitly confirm that the hybrid penalty terms and the anisotropic operator together preserve the Cordes-derived estimates without introducing dimension- or degree-dependent factors; the current abstract-level statement leaves this step unexamined.
minor comments (2)
  1. [Introduction / preliminaries] Notation for the anisotropy parameter and the precise statement of the Cordes condition should be introduced early and used consistently throughout the analysis sections.
  2. [Numerical results] Numerical examples would benefit from a table or figure explicitly reporting the observed convergence rates alongside the theoretical predictions for different polynomial degrees and anisotropy values.

Simulated Author's Rebuttal

2 responses · 0 unresolved

Thank you for the referee's thorough review and constructive feedback on our manuscript. We are pleased that the referee recognizes the potential significance of our C0-hybrid interior penalty method for the nematic Helmholtz-Korteweg equation. Below, we provide point-by-point responses to the major comments. We will incorporate revisions to address the points raised regarding explicit thresholds and detailed convergence arguments.

read point-by-point responses

  1. Referee: [Abstract / stability analysis] Abstract and stability theorem (presumably §4 or §5): the central claim that stability holds for any p ≥ 2 and any spatial dimension rests on the Cordes condition producing coercivity when anisotropy is below an unspecified threshold. No explicit form of this threshold is supplied, nor is it shown that the resulting coercivity constants remain independent of dimension and degree; without these details the uniformity assertion cannot be verified and is load-bearing for the main result.

    Authors: We thank the referee for highlighting this important point. Upon review, we agree that an explicit form of the anisotropy threshold would strengthen the presentation and allow verification of the uniformity. In the revised manuscript, we will derive and state an explicit bound on the anisotropy parameter in terms of the Cordes constant, the penalty parameters, and other relevant quantities. Additionally, we will include a remark or appendix tracing the dependence of the coercivity constant through the proof to confirm its independence from the polynomial degree p ≥ 2 and the spatial dimension. This addresses the load-bearing nature of the uniformity claim. revision: yes

  2. Referee: [Convergence section] Convergence analysis (presumably §5): the passage from stability to convergence under minimal regularity is asserted, but the argument must explicitly confirm that the hybrid penalty terms and the anisotropic operator together preserve the Cordes-derived estimates without introducing dimension- or degree-dependent factors; the current abstract-level statement leaves this step unexamined.

    Authors: We acknowledge that the convergence section would benefit from a more explicit verification that the hybrid penalty terms and the anisotropic terms preserve the key estimates from the stability analysis. In the revision, we will expand the relevant section to include a detailed outline of how the consistency and stability estimates are combined to obtain convergence under minimal regularity assumptions. We will explicitly note that no dimension- or degree-dependent factors are introduced beyond those already controlled in the stability result, thereby confirming the passage from stability to convergence. revision: yes

Circularity Check

0 steps flagged

No circularity: stability and convergence rest on external Cordes condition and standard analysis

full rationale

The paper's central claims are a stability proof for the C0-hybrid IP discretization (any p≥2, any dimension, anisotropy below a threshold) and convergence under minimal regularity. These are derived using the Cordes condition on the anisotropic fourth-order operator, which is an external, standard tool from PDE theory (not defined or fitted within the paper). The abstract explicitly refines prior continuous analysis by Farrell et al. (2025) via this condition, with no self-citations, no fitted parameters renamed as predictions, and no equations that reduce to inputs by construction. The anisotropy threshold is a hypothesis for the coercivity estimate to hold, not a self-definitional loop. The derivation chain is therefore self-contained against external mathematical benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The stability proof depends on the Cordes condition applying to the anisotropic fourth-order term; no free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption The Cordes condition holds for the coefficients of the anisotropic term
    Invoked to obtain the stability estimate independent of dimension.

pith-pipeline@v0.9.1-grok · 5706 in / 1183 out tokens · 23334 ms · 2026-06-25T20:35:43.501820+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

65 extracted references · 15 canonical work pages

  1. [1]

    Selinger, J.V.: Introduction to the Theory of Soft Matter: from Ideal Gases to Liquid Crystals vol. 408. Springer, (2016)

    2016

  2. [2]

    Molecular Crystals and Liquid Crystals647(1), 1–27 (2017)

    Ball, J.M.: Mathematics and liquid crystals. Molecular Crystals and Liquid Crystals647(1), 1–27 (2017)

    2017

  3. [3]

    Chapman and Hall/CRC, New York (2018)

    Virga, E.G.: Variational Theories for Liquid Crystals. Chapman and Hall/CRC, New York (2018)

    2018

  4. [4]

    Acta Numer

    Wang, W., Zhang, L., Zhang, P.: Modelling and computation of liquid crystals. Acta Numerica 30, 765–851 (2021) https://doi.org/10.1017/S0962492921000088

    work page doi:10.1017/s0962492921000088 2021

  5. [5]

    Lord, A.E., Labes, M.M.: Anisotropic Ultrasonic Properties of a Nematic Liquid Crystal. Phys. Rev. Lett.25, 570–572 (1970) https://doi.org/10.1103/PhysRevLett.25.570

    work page doi:10.1103/physrevlett.25.570 1970

  6. [6]

    Physical Review E66(5 Pt 1), 051708–051708 (2002)

    Selinger, J.V., Spector, M.S., Greanya, V.A., Weslowski, B.T., Shenoy, D.K., Shashidhar, R.: Acoustic realignment of nematic liquid crystals. Physical Review E66(5 Pt 1), 051708–051708 (2002)

    2002

  7. [7]

    Springer, (2012)

    Sonnet, A.M., Virga, E.G.: Dissipative Ordered Fluids: Theories for Liquid Crystals. Springer, (2012)

    2012

  8. [8]

    Physical Review E111(3), 035413 (2025)

    Farrell, P.E., Zerbinati, U.: Time-harmonic waves in korteweg and nematic-korteweg fluids. Physical Review E111(3), 035413 (2025)

    2025

  9. [9]

    https://arxiv.org/abs/2503.10771

    Farrell, P.E., van Beeck, T., Zerbinati, U.: Analysis and numerical analysis of the Helmholtz- Korteweg equation (2025). https://arxiv.org/abs/2503.10771

    Pith/arXiv arXiv 2025

  10. [10]

    Physical Review E80(3 Pt 1), 031705– 031705 (2009)

    Virga, E.G.: Variational theory for nematoacoustics. Physical Review E80(3 Pt 1), 031705– 031705 (2009)

    2009

  11. [11]

    Mathematische Annalen131(3), 278–312 (1956)

    Cordes, H.O.: ¨Uber die erste randwertaufgabe bei quasilinearen differentialgleichungen zweiter ordnung in mehr als zwei variablen. Mathematische Annalen131(3), 278–312 (1956)

    1956

  12. [12]

    SIAM Journal on Numerical Analysis51(4), 2088–2106 (2013)

    Smears, I., S¨ uli, E.: Discontinuous galerkin finite element approximation of nondivergence form elliptic equations with cordes coefficients. SIAM Journal on Numerical Analysis51(4), 2088–2106 (2013)

    2088

  13. [13]

    Acta Numerica26, 137–303 (2017)

    Neilan, M., Salgado, A.J., Zhang, W.: Numerical analysis of strongly nonlinear pdes. Acta Numerica26, 137–303 (2017)

    2017

  14. [14]

    The Aeronautical Journal72(692), 701–709 (1968)

    Argyris, J.H., Fried, I., Scharpf, D.W.: The TUBA family of plate elements for the matrix displacement method. The Aeronautical Journal72(692), 701–709 (1968)

    1968

  15. [15]

    Clough, R.W., Tocher, J.L.: Finite element stiffness matrices for analysis of plate bending, 515– 546 (1965) 27

    1965

  16. [16]

    Journal of Approximation Theory7(4), 334–351 (1973)

    ˇZen´ ıˇ sek, A.: Polynomial approximation on tetrahedrons in the finite element method. Journal of Approximation Theory7(4), 334–351 (1973)

    1973

  17. [17]

    Computer methods in applied mechanics and engineering196(13-16), 1851–1863 (2007)

    S¨ uli, E., Mozolevski, I.: hp-version interior penalty DGFEMs for the biharmonic equation. Computer methods in applied mechanics and engineering196(13-16), 1851–1863 (2007)

    2007

  18. [18]

    Journal of Scientific Computing30(3), 465–491 (2007)

    Mozolevski, I., S¨ uli, E., B¨ osing, P.R.: hp-version a priori error analysis of interior penalty discon- tinuous Galerkin finite element approximations to the biharmonic equation. Journal of Scientific Computing30(3), 465–491 (2007)

    2007

  19. [19]

    Computational Methods in Applied Mathematics3(4), 596–607 (2003)

    Mozolevski, I., S¨ uli, E.: A priori error analysis for the hp-version of the discontinuous Galerkin finite element method for the biharmonic equation. Computational Methods in Applied Mathematics3(4), 596–607 (2003)

    2003

  20. [20]

    International Journal of Numerical Analysis and Modeling16(5), 825–846 (2019)

    Dong, Z.: Discontinuous Galerkin methods for the biharmonic problem on polygonal and poly- hedral meshes. International Journal of Numerical Analysis and Modeling16(5), 825–846 (2019)

    2019

  21. [21]

    SIAM Journal on Numerical Analysis60(5), 2626–2656 (2022)

    Dong, Z., Ern, A.: Hybrid high-order and weak Galerkin methods for the biharmonic problem. SIAM Journal on Numerical Analysis60(5), 2626–2656 (2022)

    2022

  22. [22]

    https: //arxiv.org/abs/2504.16608

    Liang, Y., Tran, N.T.: A hybrid high-order method for the biharmonic problem (2025). https: //arxiv.org/abs/2504.16608

    Pith/arXiv arXiv 2025

  23. [23]

    Computer Methods in Applied Mechanics and Engineering191(34), 3669–3750 (2002) https://doi.org/10.1016/S0045-7825(02)00286-4

    Engel, G., Garikipati, K., Hughes, T.J.R., Larson, M.G., Mazzei, L., Taylor, R.L.: Continu- ous/discontinuous finite element approximations of fourth-order elliptic problems in structural and continuum mechanics with applications to thin beams and plates, and strain gradient elas- ticity. Computer Methods in Applied Mechanics and Engineering191(34), 3669–...

    work page doi:10.1016/s0045-7825(02)00286-4 2002

  24. [24]

    Journal of Scientific Computing22(1), 83–118 (2005)

    Brenner, S.C., Sung, L.-Y.:C 0 interior penalty methods for fourth order elliptic boundary value problems on polygonal domains. Journal of Scientific Computing22(1), 83–118 (2005)

    2005

  25. [25]

    Brenner, S.C.:C 0 Interior Penalty Methods, pp. 79–147. Springer, Berlin, Heidelberg (2012). https://doi.org/10.1007/978-3-642-23914-4 2

    work page doi:10.1007/978-3-642-23914-4 2012

  26. [26]

    SIAM journal on numerical analysis49(2), 869–892 (2011)

    Brenner, S.C., Neilan, M.: AC 0 interior penalty method for a fourth order elliptic singular perturbation problem. SIAM journal on numerical analysis49(2), 869–892 (2011)

    2011

  27. [27]

    SIAM Journal on Numerical Analysis47(2), 1319–1365 (2009)

    Cockburn, B., Gopalakrishnan, J., Lazarov, R.: Unified hybridization of discontinuous Galerkin, mixed, and continuous Galerkin methods for second order elliptic problems. SIAM Journal on Numerical Analysis47(2), 1319–1365 (2009)

    2009

  28. [28]

    Master’s thesis, RWTH Aachen (May 2010)

    Lehrenfeld, C.: Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems. Master’s thesis, RWTH Aachen (May 2010). https://doi.org/10.25625/O4VBYH

    work page doi:10.25625/o4vbyh 2010

  29. [29]

    Computer Methods in Applied Mechanics and Engineering283, 1–21 (2015) 28

    Di Pietro, D.A., Ern, A.: A hybrid high-order locking-free method for linear elasticity on general meshes. Computer Methods in Applied Mechanics and Engineering283, 1–21 (2015) 28

    2015

  30. [30]

    Computer Methods in Applied Mechanics and Engineering307, 339–361 (2016)

    Lehrenfeld, C., Sch¨ oberl, J.: High order exactly divergence-free hybrid discontinuous Galerkin methods for unsteady incompressible flows. Computer Methods in Applied Mechanics and Engineering307, 339–361 (2016)

    2016

  31. [31]

    IMA Journal of Numerical Analysis44(1), 24–57 (2024)

    Dong, Z., Ern, A.:C 0-hybrid high-order methods for biharmonic problems. IMA Journal of Numerical Analysis44(1), 24–57 (2024)

    2024

  32. [32]

    SIAM Journal on Scientific Computing43(6), 3912–3937 (2021)

    Fu, G.: Uniform Auxiliary Space Preconditioning for HDG Methods for Elliptic Operators with a Parameter Dependent Low Order Term. SIAM Journal on Scientific Computing43(6), 3912–3937 (2021)

    2021

  33. [33]

    Halla, M., Lehrenfeld, C., Stocker, P.: A new T-compatibility condition and its application to the discretization of the damped time-harmonic Galbrun’s equation (2022)

    2022

  34. [34]

    https: //doi.org/10.25625/KGHQWV

    van Beeck, T.: On Stable Discontinuous Galerkin Discretizations for Galbrun’s Equation. https: //doi.org/10.25625/KGHQWV

    work page doi:10.25625/kghqwv

  35. [35]

    ESAIM: Mathematical Modelling and Numerical Analysis60(3), 1135–1175 (2026)

    Halla, M., Lehrenfeld, C., van Beeck, T.: Hybrid discontinuous galerkin discretizations for the damped time-harmonic galbrun’s equation. ESAIM: Mathematical Modelling and Numerical Analysis60(3), 1135–1175 (2026)

    2026

  36. [36]

    CRC Press, (2019)

    Sayas, F.J., Brown, T.S., Hassell, M.E.: Variational Techniques for Elliptic Partial Differential Equations: Theoretical Tools and Advanced Applications. CRC Press, (2019)

    2019

  37. [37]

    Unified Transform for Boundary Value Problems: Applications and Advances, 93–159 (2014)

    Spence, E.A.: Chapter 6: Overview of variational formulations for linear elliptic pdes. Unified Transform for Boundary Value Problems: Applications and Advances, 93–159 (2014)

    2014

  38. [38]

    Computers & Mathematics with Applications27(12), 69–84 (1994)

    Demkowicz, L.: Asymptotic convergence in finite and boundary element methods: Part 1: Theoretical results. Computers & Mathematics with Applications27(12), 69–84 (1994)

    1994

  39. [39]

    Springer series in Computational Mathematics, vol

    Sauter, S.A., Schwab, C.: Boundary Element Methods. Springer series in Computational Mathematics, vol. 39. Springer, Berlin (2011). https://doi.org/10.1007/978-3-540-68093-2

    work page doi:10.1007/978-3-540-68093-2 2011

  40. [40]

    Computer Methods in Applied Mechanics and Engineering317, 318–340 (2017) https://doi.org/10.1016/j.cma.2016

    Bespalov, A., Haberl, A., Praetorius, D.: Adaptive FEM with coarse initial mesh guarantees optimal convergence rates for compactly perturbed elliptic problems. Computer Methods in Applied Mechanics and Engineering317, 318–340 (2017) https://doi.org/10.1016/j.cma.2016. 12.014

    work page doi:10.1016/j.cma.2016 2017

  41. [41]

    Halla, M.: Convergence analysis of nonconforming h (div)-finite elements for the damped time- harmonic galbrun’s equation: M. halla. Numerische Mathematik, 1–31 (2026)

    2026

  42. [42]

    Journal of Soviet Mathematics15(6), 675–705 (1981)

    Vainikko, G.: Regular convergence of operators and approximate solution of equations. Journal of Soviet Mathematics15(6), 675–705 (1981)

    1981

  43. [43]

    Di Pietro, D., Ern, A.: Mathematical Aspects of Discontinuous Galerkin Methods vol. 69. Heidelberg (2012). https://doi.org/10.1007/978-3-642-22980-0

    work page doi:10.1007/978-3-642-22980-0 2012

  44. [44]

    29 Mathematics of computation28(128), 959–962 (1974)

    Schatz, A.H.: An observation concerning ritz-galerkin methods with indefinite bilinear forms. 29 Mathematics of computation28(128), 959–962 (1974)

    1974

  45. [45]

    eprint:https: / / epubs

    Grisvard, P.: Elliptic Problems in Nonsmooth Domains. Society for Industrial and Applied Mathematics, (2011). https://doi.org/10.1137/1.9781611972030

    work page doi:10.1137/1.9781611972030 2011

  46. [46]

    Springer, New York (2004)

    Renardy, M., Rogers, R.C.: An Introduction to Partial Differential Equations. Springer, New York (2004)

    2004

  47. [47]

    Maugeri, A., Palagachev, D.K., Softova, L.G.: Elliptic and Parabolic Equations with Discontin- uous Coefficients vol. 109. Wiley Verlag GmbH & Co., (2000)

    2000

  48. [48]

    arXiv preprint arXiv:2601.14548 (2026)

    Sprekeler, T.: A cordes framework for stationary fokker–planck–kolmogorov equations. arXiv preprint arXiv:2601.14548 (2026)

    arXiv 2026

  49. [49]

    Kato, T.: Perturbation Theory for Linear Operators vol. 132. Springer, Berlin (1966)

    1966

  50. [50]

    Journal of mathematical analysis and applications441(1), 235–258 (2016)

    Cuenin, J.-C., Tretter, C.: Non-symmetric perturbations of self-adjoint operators. Journal of mathematical analysis and applications441(1), 235–258 (2016)

    2016

  51. [51]

    Springer, Heidelberg (2021)

    Ern, A., Guermond, J.-L.: Finite Elements II: Galerkin Approximation, Elliptic and Mixed PDEs. Springer, Heidelberg (2021)

    2021

  52. [52]

    Journal of Differential Equations266(6), 2869– 2923 (2019) https://doi.org/10.1016/j.jde.2018.08.048

    Graham, I.G., Pembery, O.R., Spence, E.A.: The helmholtz equation in heterogeneous media: A priori bounds, well-posedness, and resonances. Journal of Differential Equations266(6), 2869– 2923 (2019) https://doi.org/10.1016/j.jde.2018.08.048

    work page doi:10.1016/j.jde.2018.08.048 2019

  53. [53]

    Computer Methods in Applied Mechanics and Engineering193(23-26), 2565–2580 (2004)

    Dawson, C., Sun, S., Wheeler, M.F.: Compatible algorithms for coupled flow and transport. Computer Methods in Applied Mechanics and Engineering193(23-26), 2565–2580 (2004)

    2004

  54. [54]

    and Guermond, J.-L

    Ern, A., Guermond, J.-L.: Finite Elements I: Approximation and Interpolation. Springer, Heidelberg (2021). https://doi.org/10.1007/978-3-030-56341-7

    work page doi:10.1007/978-3-030-56341-7 2021

  55. [55]

    IMA Journal of Numerical Analysis30(3), 777–798 (2010)

    Brenner, S.C., Gudi, T., Sung, L.-y.: An a posteriori error estimator for a quadratic C 0-interior penalty method for the biharmonic problem. IMA Journal of Numerical Analysis30(3), 777–798 (2010)

    2010

  56. [56]

    IMA journal of numerical analysis 31(1), 281–298 (2011)

    Georgoulis, E.H., Houston, P., Virtanen, J.: An a posteriori error indicator for discontinuous galerkin approximations of fourth-order elliptic problems. IMA journal of numerical analysis 31(1), 281–298 (2011)

    2011

  57. [57]

    Calcolo56(4), 44 (2019)

    Brenner, S.C., Sung, L.-Y.: Virtual enriching operators. Calcolo56(4), 44 (2019)

    2019

  58. [58]

    Numerische Mathematik102(2), 231–255 (2005)

    Brenner, S.C., Wang, K.: Two-level additive Schwarz preconditioners for C 0 interior penalty methods. Numerische Mathematik102(2), 231–255 (2005)

    2005

  59. [59]

    Comptes Rendus

    Burman, E., Ern, A.: Discontinuous Galerkin approximation with discrete variational principle for the nonlinear Laplacian. Comptes Rendus. Math´ ematique346(17-18), 1013–1016 (2008) 30

    2008

  60. [60]

    IMA Journal of Numerical Analysis29(4), 827–855 (2009) https://doi.org/10.1093/imanum/drn038

    Buffa, A., Ortner, C.: Compact embeddings of broken Sobolev spaces and applications. IMA Journal of Numerical Analysis29(4), 827–855 (2009) https://doi.org/10.1093/imanum/drn038

    work page doi:10.1093/imanum/drn038 2009

  61. [61]

    Mathematics of Computation 79(271), 1303–1330 (2010)

    Di Pietro, D.A., Ern, A.: Discrete functional analysis tools for discontinuous Galerkin methods with application to the incompressible Navier–Stokes equations. Mathematics of Computation 79(271), 1303–1330 (2010)

    2010

  62. [62]

    Comput- ing and Visualization in Science1(1), 41–52 (1997) https://doi.org/10.1007/s007910050004

    Sch¨ oberl, J.: Netgen an advancing front 2d/3d-mesh generator based on abstract rules. Comput- ing and Visualization in Science1(1), 41–52 (1997) https://doi.org/10.1007/s007910050004

    work page doi:10.1007/s007910050004 1997

  63. [63]

    Sch¨ oberl, J.: C++ 11 implementation of finite elements in NGSolve (2014)

    2014

  64. [64]

    Zenodo (2026)

    van Beeck, T.: Replication data for: A hybrid C0-interior penalty method for the nematic Helmholtz–Korteweg equation. Zenodo (2026). https://doi.org/10.5281/zenodo.20761015

    work page doi:10.5281/zenodo.20761015 2026

  65. [65]

    Ciarlet, P.:T-coercivity: a practical tool for the study of variational formulations in Hilbert spaces (2025) 31

    2025