pith. sign in

arxiv: 2604.10600 · v1 · submitted 2026-04-12 · 🧮 math.NA · cs.NA

Nonconforming hp-FE/BE coupling on unstructured meshes based on Nitsche's method

Alexey Chernov , Peter Hansbo , Erik Marc Schetzke This is my paper

Pith reviewed 2026-05-10 16:13 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Nitsche's methodhp finite elementsboundary elementsnon-matching meshesnonconforming couplinga priori error estimatesstabilization parameterdomain decomposition
0
0 comments X

The pith

Nitsche's method yields a stable hp finite-element and boundary-element coupling on non-matching unstructured meshes.

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

The paper constructs a nonconforming coupling between hp finite-element and boundary-element approximations that permits completely independent meshes and polynomial degrees in each subdomain. Nitsche's method is used to enforce transmission conditions weakly at the interface, producing a symmetric positive-definite discrete problem that avoids any inf-sup condition. An explicit positive lower bound for the stabilization parameter is derived to guarantee coercivity, after which standard a priori error estimates are proved in the energy norm and in L2 for both quasi-uniform and geometrically refined hp meshes. The analysis is carried out in detail for a bounded domain and extends with only minor changes to the case of an unbounded boundary-element subdomain or to decompositions with more than two subregions.

Core claim

A positive-definite hp-FE/BE formulation on non-conforming meshes is obtained by adding consistency and penalty terms from Nitsche's method to the interface; the resulting bilinear form is coercive once the stabilization function exceeds an explicitly computable threshold that depends on the local mesh sizes and polynomial degrees, and optimal a priori error bounds then follow for both analytic and singular solutions under quasi-uniform or geometrically graded hp refinements.

What carries the argument

Nitsche interface terms that weakly enforce continuity and flux balance across non-matching FE/BE boundaries while supplying an explicit stabilization threshold.

If this is right

  • The same Nitsche formulation remains stable and convergent when applied to pure finite-element or pure boundary-element problems.
  • The analysis extends directly to decompositions involving more than two subdomains.
  • A priori error estimates hold for both quasi-uniform meshes and geometrically refined hp meshes near singularities.
  • The method applies with similar arguments to exterior problems that contain an unbounded boundary-element subdomain.

Where Pith is reading between the lines

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

  • The explicit threshold supplies a practical, computable rule for choosing the stabilization parameter from local mesh and degree data without trial-and-error tuning.
  • Independent hp-refinement across subdomains becomes feasible without the need to enforce mesh conformity at artificial interfaces.
  • Positive-definiteness of the resulting system may reduce the cost of iterative solvers compared with saddle-point formulations that arise from mortar-type couplings.

Load-bearing premise

The stabilization parameter must exceed an explicitly derived positive threshold that depends on the local mesh sizes and polynomial degrees at the interface.

What would settle it

Numerical computation of the smallest eigenvalue of the discrete system on a fixed mesh and fixed polynomial degrees when the stabilization parameter is deliberately set below the paper's explicit threshold.

read the original abstract

We construct and analyse a $hp$-FE/BE coupling on non-matching meshes, based on Nitsche's method. Both the mesh size and the polynomial degree are changed to improve accuracy. Nitsche's method leads to a positive definite formulation, thus, unlike the mortar method, it does not require the Babu\v{s}ka-Brezzi condition for stability. The method is stable provided the stabilization function is larger than a certain threshold. We obtain an explicit bound for the threshold and derive a priori error estimates. Our analysis can be easily extended to the pure FE or the pure BE decomposition as well as to the case of more than two subdomains. The problem in a bounded domain is considered in detail, but the case of an unbounded BE subdomain and a bounded FE subdomain follows with similar arguments. We develop convergence analysis and provide numerical examples for quasi-uniform as well as geometrically refined $hp$ discretisations in both subdomains with analytic and singular solutions.

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

1 major / 2 minor

Summary. The paper constructs and analyzes a nonconforming hp-FE/BE coupling method on non-matching unstructured meshes using Nitsche's method. It derives an explicit threshold for the stabilization parameter to ensure stability (positive-definiteness without Babuška-Brezzi), obtains a priori error estimates in suitable norms, and extends the analysis to pure FE/BE cases, multiple subdomains, and unbounded domains. Numerical examples are provided for quasi-uniform and geometrically refined hp-discretizations with analytic and singular solutions.

Significance. If the explicit stability threshold holds with the claimed p-dependence, the method supplies a practical, positive-definite alternative to mortar couplings for hp-FE/BE problems on non-matching meshes, with a priori estimates that support both quasi-uniform and geometric refinement. The explicit bound and extension to unbounded domains are notable strengths for implementation.

major comments (1)
  1. [§4] §4 (stability analysis), around the derivation of the threshold for the Nitsche stabilization function: the explicit bound relies on hp-trace inequalities applied to non-matching interface terms after integration by parts. The constants in these inequalities grow at least linearly with p (and potentially worse on unstructured non-matching meshes); if the paper's threshold does not absorb the precise p-scaling of every such constant, coercivity can fail for large p even when the stated threshold is respected. This is load-bearing for the positive-definiteness claim and all subsequent a priori estimates.
minor comments (2)
  1. [Introduction] The abstract and introduction claim the analysis 'can be easily extended' to pure FE/BE and multiple subdomains; a brief outline of the necessary modifications would strengthen the presentation.
  2. [Numerical examples] Numerical examples section: clarify the precise definition of the stabilization function used in the experiments (e.g., whether it is taken exactly at the derived threshold or a multiple) to allow direct comparison with the theory.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comment on the stability analysis. The positive assessment of the method's practical advantages for hp-FE/BE coupling on non-matching meshes is appreciated. We address the single major comment point-by-point below, with no standing objections.

read point-by-point responses

  1. Referee: [§4] §4 (stability analysis), around the derivation of the threshold for the Nitsche stabilization function: the explicit bound relies on hp-trace inequalities applied to non-matching interface terms after integration by parts. The constants in these inequalities grow at least linearly with p (and potentially worse on unstructured non-matching meshes); if the paper's threshold does not absorb the precise p-scaling of every such constant, coercivity can fail for large p even when the stated threshold is respected. This is load-bearing for the positive-definiteness claim and all subsequent a priori estimates.

    Authors: We thank the referee for this observation on the p-dependence. In Section 4, the threshold for the stabilization function is derived in the proof of the coercivity result (Theorem 4.3) by applying the hp-trace inequalities of Lemma 4.2 to the interface terms obtained after integration by parts. These inequalities are the standard ones for hp-FEM on shape-regular meshes (with constants scaling as O(p^2) for the relevant L^2 trace terms in 2D), and the non-matching character of the interface is accounted for by the Nitsche terms without introducing additional p-growth. The explicit threshold (displayed after Equation (4.12)) is chosen strictly larger than the sum of all such p-dependent contributions from both subdomains. Under the shape-regularity assumptions stated in Section 2, the constants remain uniform even for unstructured meshes. Consequently, the positive-definiteness holds for any polynomial degree p whenever the threshold is respected, and the subsequent a priori estimates follow directly. No revision to the analysis is required, though we are happy to add a clarifying remark on the absorption of the p-scaling if the editor deems it useful. revision: no

Circularity Check

0 steps flagged

No circularity: stability threshold and error estimates derived from standard variational arguments and trace inequalities

full rationale

The derivation establishes coercivity of the Nitsche bilinear form by bounding interface terms via trace inequalities and choosing the stabilization parameter above an explicit threshold derived from those bounds. A priori estimates then follow from the resulting coercivity, consistency, and approximation properties. No step reduces a claimed result to a fitted input, self-citation chain, or definitional renaming; the analysis is self-contained against external Sobolev-space benchmarks and does not invoke prior author-specific uniqueness theorems as load-bearing premises.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard mathematical assumptions from finite element theory rather than new postulates.

axioms (2)
  • standard math Standard trace and inverse inequalities hold on the interface for the chosen polynomial degrees
    Invoked to derive the explicit stability threshold
  • domain assumption The solution possesses sufficient Sobolev regularity for the a priori estimates
    Required for both analytic and singular solution cases

pith-pipeline@v0.9.0 · 5477 in / 1279 out tokens · 47910 ms · 2026-05-10T16:13:33.157734+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

37 extracted references · 37 canonical work pages

  1. [1]

    In: Asymptotic and Numerical Methods for Partial Differential Equations with Critical Parameters, pp

    Bernardi, C., Maday, Y., Patera, A.T.: Domain decomposition by the mortar element method. In: Asymptotic and Numerical Methods for Partial Differential Equations with Critical Parameters, pp. 269–286. Springer, ??? (1993)

    work page 1993

  2. [2]

    Abhandlungen aus dem Mathematischen Seminar der Universit¨ at Hamburg36, 9–15 (1971)

    Nitsche, J.A.: ¨Uber ein Variationsprinzip zur L¨ osung von Dirichlet-Problemen bei Verwendung von Teilr¨ aumen, die keinen Randbedingungen unterworfen sind. Abhandlungen aus dem Mathematischen Seminar der Universit¨ at Hamburg36, 9–15 (1971)

    work page 1971

  3. [3]

    ESAIM: Mathematical Modelling and Numerical Analysis37(2), 209–225 (2003)

    Becker, R., Hansbo, P., Stenberg, R.: A finite element method for domain decomposition with non-matching grids. ESAIM: Mathematical Modelling and Numerical Analysis37(2), 209–225 (2003)

    work page 2003

  4. [4]

    IMA Journal of Numerical Analysis23(2), 331–358 (2003) https://doi.org/10.1093/imanum/23.2.331

    Heinrich, B., Nicaise, S.: The Nitsche mortar finite element method for trans- mission problems with singularities. IMA Journal of Numerical Analysis23(2), 331–358 (2003) https://doi.org/10.1093/imanum/23.2.331

    work page doi:10.1093/imanum/23.2.331 2003

  5. [5]

    IMA Journal of Numerical Analysis26(3), 541–562 (2006) https://doi.org/ 10.1093/imanum/dri046 33

    Heinrich, B., Jung, B.: The Fourier-finite-element method with Nitsche mortar- ing. IMA Journal of Numerical Analysis26(3), 541–562 (2006) https://doi.org/ 10.1093/imanum/dri046 33

    work page doi:10.1093/imanum/dri046 2006

  6. [6]

    Numerical Algorithms91, 997–1019 (2022)

    Betcke, T., Bosy, M., Burman, E.: Hybrid coupling of finite element and boundary element methods using Nitsche’s method and the Calderon projection. Numerical Algorithms91, 997–1019 (2022)

    work page 2022

  7. [7]

    SIAM Journal on Scientific Computing41(3), 1357–1384 (2019)

    Betcke, T., Burman, E., Scroggs, M.W.: Boundary element methods with weakly imposed boundary conditions. SIAM Journal on Scientific Computing41(3), 1357–1384 (2019)

    work page 2019

  8. [8]

    Computing (5), 207–213 (1970)

    Babuˇ ska, I.: The finite element method for elliptic equations with discontinuous coefficients. Computing (5), 207–213 (1970)

    work page 1970

  9. [9]

    Springer, ??? (2008)

    Steinbach, O.: Numerical Approximation Methods for Elliptic Boundary Value Problems: Finite and Boundary Elements. Springer, ??? (2008)

    work page 2008

  10. [10]

    ESAIM: Mathematical Modelling and Numerical Analysis29(7), 779– 817 (1995)

    Carstensen, C., Stephan, E.P.: Adaptive coupling of boundary elements and finite elements. ESAIM: Mathematical Modelling and Numerical Analysis29(7), 779– 817 (1995)

    work page 1995

  11. [11]

    SIAM Journal on Mathematical Analysis19(3), 613–626 (1988) https: //doi.org/10.1137/0519043

    Costabel, M.: Boundary integral operators on Lipschitz domains: Elementary results. SIAM Journal on Mathematical Analysis19(3), 613–626 (1988) https: //doi.org/10.1137/0519043

    work page doi:10.1137/0519043 1988

  12. [12]

    Springer, ??? (2004)

    Sauter, S., Schwab, C.: Randelementmethoden: Analyse, Numerik und Implemen- tierung Schneller Algorithmen. Springer, ??? (2004)

    work page 2004

  13. [13]

    Stephan, E.P.: Coupling of Boundary Element Methods and Finite Element Meth- ods, pp. 1–40. John Wiley & Sons, Ltd, ??? (2017). https://doi.org/10.1002/ 9781119176817.ecm2008

    work page 2017

  14. [14]

    Oxford University Press, ??? (1998)

    Schwab, C.:p- andhp-finite Element Methods Theory and Applications in Solid and Fluid Mechanics. Oxford University Press, ??? (1998)

    work page 1998

  15. [15]

    Mathematical Methods in the Applied Sciences12(5), 413–427 (1990) https://doi.org/10.1002/ mma.1670120506

    Babuˇ ska, I., Guo, B.Q., Stephan, E.P.: On the exponential convergence of the h–pversion for boundary element Galerkin methods on polygons. Mathematical Methods in the Applied Sciences12(5), 413–427 (1990) https://doi.org/10.1002/ mma.1670120506

    work page 1990

  16. [16]

    ESAIM: Mathematical Modelling and Numerical Analysis37(3), 495–514 (2003) https://doi.org/10.1051/m2an:2003039

    Hansbo, A., Hansbo, P., Larson, M.G.: A finite element method on composite grids based on Nitsche’s method. ESAIM: Mathematical Modelling and Numerical Analysis37(3), 495–514 (2003) https://doi.org/10.1051/m2an:2003039

    work page doi:10.1051/m2an:2003039 2003

  17. [17]

    Computational Methods in Applied Mathe- matics1(4), 367–382 (2001) https://doi.org/10.2478/cmam-2001-0024

    Lazarov, R.D., Tomov, S.Z., Vassilevski, P.S.: Interior penalty discontinuous approximations of elliptic problems. Computational Methods in Applied Mathe- matics1(4), 367–382 (2001) https://doi.org/10.2478/cmam-2001-0024

    work page doi:10.2478/cmam-2001-0024 2001

  18. [18]

    Sch¨ otzau, D., Wihler, T.: Exponential convergence of mixedhp-DGFEM for Stokes flow in polygons. Numer. Math. (96), 339–361 (2003) https://doi.org/10. 34 1007/s00211-003-0478-5

    work page 2003

  19. [19]

    SIAM Journal on Numerical Analysis 39(6), 2133–2163 (2002) https://doi.org/10.1137/S0036142900374111

    Houston, P., Schwab, C., S¨ uli, E.: Discontinuoushp-finite element methods for advection-diffusion-reaction problems. SIAM Journal on Numerical Analysis 39(6), 2133–2163 (2002) https://doi.org/10.1137/S0036142900374111

    work page doi:10.1137/s0036142900374111 2002

  20. [20]

    Mathematics of Computation69(230), 521–546 (2000)

    Seshaiyer, P., Suri, M.: Uniformhpconvergence results for the mortar finite element method. Mathematics of Computation69(230), 521–546 (2000)

    work page 2000

  21. [21]

    Numerische Mathematik (84), 173–197 (1999)

    Belgacem, F.: The mortar finite element method with Lagrange multipliers. Numerische Mathematik (84), 173–197 (1999)

    work page 1999

  22. [22]

    Springer, (2001)

    Wohlmuth, B.I.: Discretization Methods and Iterative Solvers Based on Domain Decomposition. Springer, (2001)

    work page 2001

  23. [23]

    Calcolo (41), 115–137 (2004) https://doi.org/10.1007/ s10092-004-0087-4

    Fritz, A., H¨ ueber, S., Wohlmuth, B.: A comparison of mortar and Nitsche tech- niques for linear elasticity. Calcolo (41), 115–137 (2004) https://doi.org/10.1007/ s10092-004-0087-4

    work page 2004

  24. [24]

    Mathe- matical Models and Methods in Applied Sciences15(04), 575–622 (2005) https: //doi.org/10.1142/S0218202505000480

    Costabel, M., Dauge, M., Schwab, C.: Exponential convergence ofhp-fem for Maxwell equations with weighted regularization in polygonal domains. Mathe- matical Models and Methods in Applied Sciences15(04), 575–622 (2005) https: //doi.org/10.1142/S0218202505000480

    work page doi:10.1142/s0218202505000480 2005

  25. [25]

    Banach center publications15, 175–251 (1985) https://doi.org/10.4064/-15-1-175-251

    Costabel, M., Stephan, E.: Boundary integral equations for mixed boundary value problems in polygonal domains and galerkin approximation. Banach center publications15, 175–251 (1985) https://doi.org/10.4064/-15-1-175-251

    work page doi:10.4064/-15-1-175-251 1985

  26. [26]

    Computer Methods in Applied Mechanics and Engineering 192(25), 2765–2773 (2003) https://doi.org/10.1016/S0045-7825(03)00294-9

    Warburton, T., Hesthaven, J.S.: On the constants inhp-finite element trace inverse inequalities. Computer Methods in Applied Mechanics and Engineering 192(25), 2765–2773 (2003) https://doi.org/10.1016/S0045-7825(03)00294-9

    work page doi:10.1016/s0045-7825(03)00294-9 2003

  27. [27]

    ESAIM: Mathematical Modelling and Numerical Analysis21(2), 199–238 (1987) https://doi.org/10.1051/m2an/1987210201991

    Babuˇ ska, I., Suri, M.: Theh-pversion of the finite element method with quasiu- niform meshes. ESAIM: Mathematical Modelling and Numerical Analysis21(2), 199–238 (1987) https://doi.org/10.1051/m2an/1987210201991

    work page doi:10.1051/m2an/1987210201991 1987

  28. [28]

    Journal of Scientific Computing (42) (2009) https://doi.org/10.1007/s10915-009-9315-z

    Georgoulis, E.H., Hall, E., Melenk, J.M.: On the suboptimality of thep-version interior penalty discontinuous Galerkin method. Journal of Scientific Computing (42) (2009) https://doi.org/10.1007/s10915-009-9315-z

    work page doi:10.1007/s10915-009-9315-z 2009

  29. [29]

    Computers and Mathematics with Applications67(4) (2014) https://doi.org/10.1016/j.camwa.2013.12.016

    Melenk, J.M., Wurzer, T.: On the stability of the boundary trace of the poly- nomialL 2-projection on triangles and tetrahedra. Computers and Mathematics with Applications67(4) (2014) https://doi.org/10.1016/j.camwa.2013.12.016

    work page doi:10.1016/j.camwa.2013.12.016 2014

  30. [30]

    PhD thesis, TU Darmstadt (1989) 35

    Petersdorff, T.: Randwertprobleme der elastizit¨ atstheorie f¨ ur polyeder - sin- gularit¨ aten und approximation mit randelementmethoden. PhD thesis, TU Darmstadt (1989) 35

    work page 1989

  31. [31]

    ESAIM: Mathematical Modelling and Numerical Analysis25(6), 783–807 (1991)

    Stephan, E.P., Suri, M.: Theh-pversion of the boundary element method on polygonal domains with quasiuniform meshes. ESAIM: Mathematical Modelling and Numerical Analysis25(6), 783–807 (1991)

    work page 1991

  32. [32]

    Part II: The trace spaces and application to the boundary value problems with nonhomogeneous boundary conditions

    Babuˇ ska, I., Guo, B.Q.: Regularity of the solution of elliptic problems with piece- wise analytic data. Part II: The trace spaces and application to the boundary value problems with nonhomogeneous boundary conditions. SIAM Journal on Mathematical Analysis20(4), 763–781 (1989) https://doi.org/10.1137/0520054

    work page doi:10.1137/0520054 1989

  33. [33]

    Quarteroni, A.: Some results of Bernstein and Jackson type for polynomial approximation inL p-spaces. Japan J. Appl. Math.1, 173–181 (1984) https: //doi.org/10.1007/BF03167866

    work page doi:10.1007/bf03167866 1984

  34. [34]

    SIAM Journal on Numerical Analysis33(2), 789–808 (1996) https://doi.org/10.1137/0733040

    Guo, B.Q., Heuer, N., Stephan, E.P.: Theh-pversion of the boundary element method for transmission problems with piecewise analytic data. SIAM Journal on Numerical Analysis33(2), 789–808 (1996) https://doi.org/10.1137/0733040

    work page doi:10.1137/0733040 1996

  35. [35]

    Springer, Cham, [2021]©2021.https://doi.org/10.1007/978-3-030-56341-7

    Ern, A., Guermond, J.-L.: Finite Elements I: Approximation and Inter- polation. Springer, ??? (2021). https://doi.org/10.1007/978-3-030-56341-7 . https://hal.science/hal-03226049

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

  36. [36]

    PhD thesis, ETH Z¨ urich (2002)

    Wihler, T.P.: Discontinuous galerkin fem for elliptic problems in polygonal domains. PhD thesis, ETH Z¨ urich (2002)

    work page 2002

  37. [37]

    Mathematics of Computation52(1989) https://doi

    Stephan, E.P., Suri, M.: On the convergence of thep-version of the boundary element galerkin method. Mathematics of Computation52(1989) https://doi. org/10.2307/2008651 36

    work page doi:10.2307/2008651 1989