{"id":"1c80c6b2-f13e-4080-81c3-4d9738874f0d","arxiv_id":"2505.21776","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For interface finite element problems, the penalty method converges only at half-order while Nitsche's method, derived as a simple correction, converges optimally.","lead":"This paper compares two numerical methods for coupling separate regions in a simulation, the penalty method and Nitsche's method, and shows the penalty method hits an accuracy wall that Nitsche's method does not. It is a compact, readable guide for engineers who need a reliable method for contact or interface problems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The O(h^{1/2}) penalty ceiling is proved only for the Poisson mortar model; the elastic-contact claim in the title and abstract rests on an unproved transfer of the modeling-error bound (6) to a free-boundary variational inequality.","rationale":"Good-faith reading: the paper is a teaching-style note whose core calculation is correct for the Poisson mortar problem. I verified the key inequality (6) in that setting: testing the penalized error w = u - u^ε with w, using the interface relation ∂w_i/∂n_i + ε^{-1}(w_i - w_j) = -∂u_i/∂n_i, and summing gives ||w||_E^2 = -∫_Γ ∂u1/∂n1 (w1 - w2) ≤ ||∂u1/∂n1||_{0,Γ} ε^{1/2}||w||_E, hence ||w||_E ≲ ε^{1/2}. Thus the model's O(h^{1/2}) ceiling is internally valid, modulo the typo in Eq. (9), which should read h^{s-1} rather than h^s. The genuinely load-bearing assumption is that this estimate survives the passage to elastic contact with an unknown contact zone. The paper says this only adds technicalities and defers to [7], but [7] is cited for Nitsche's method, not for a proof of the penalty method's suboptimality in the contact setting. Since the abstract and title make the elastic-contact claim, this is a real gap, not a manufactured one. My proposed test would settle whether the scaling changes; if it does not, the paper remains a correct model-problem exposition but should be reworded to avoid claiming the full contact result without proof.","tokens_in":7807,"tokens_out":10985,"duration_ms":107992,"concrete_test":"Re-derive the modelling-error estimate (6) for the penalized Signorini/elastic-contact variational inequality, isolating the free-boundary term near the contact-set boundary; determine whether the energy norm of u - u^ε is bounded by C ε^{1/2}||p||_{0,Γ} with p the contact pressure, or whether an O(ε) term from the transition zone appears. If the latter, then with ε≈h and p=1 the penalty total error is O(h), contradicting the abstract. A supplementary numerical run on a manufactured 2D contact solution with ε = 10^{-2}, 10^{-3}, 10^{-4} on a fixed fine mesh would provide corroborating evidence for the ε-rate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2 explicitly narrows the analysis: \"All the questions we are going to address can be answered by considering mortaring of the Poisson equation\" and defers unknown-contact technicalities to [7]. The quantitative statement that the penalty method cannot be optimally convergent -- Eq. (10) -- is built on the modeling-error estimate (6), ||u - u^ε||_E ≲ ε^{1/2}||∂u/∂n||_{0,Γ}, which is asserted but not proved. Within the Poisson mortar problem this estimate is true (testing with u - u^ε and using the interface condition gives the stated ε^{1/2} rate), so the internal argument for the model is sound. The load-bearing gap is the leap to elastic contact: there the constraint is an inequality, the active contact set is unknown, and Section 2 provides no analogue of (6) for that setting. If the free-boundary/complementarity effects change the scaling to O(ε) or alter the optimal ε-h relation, then with ε≈h the penalty method could reach O(h) for linear elements and the central conclusion would not follow. The numerics in Section 4 are also performed on the Poisson mortar problem, not on an elastic contact problem, so they do not close the gap. A secondary typo: Eq. (9) states h^s while Eq. (8) gives h^{s-1} for the interpolation contribution; this does not change the O(h^{1/2}) ceiling but should be corrected.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper discusses the penalty method and Nitsche's method for a problem motivated by elastic contact, but the actual analysis is carried out on a two-domain Poisson 'mortaring' problem. For the penalty method with penalty parameter ε ≈ h, the authors argue that the total error is limited by a modeling error of order O(h^{1/2}), so the method cannot be optimally convergent. They then derive Nitsche's method as a simple consistency correction of the penalty method, prove stability and an optimal a priori error estimate, and present numerical experiments confirming the predicted rates for both methods. The paper positions itself as a short didactic note rather than a comprehensive analysis of contact problems.","tokens_in":8104,"tokens_out":7410,"duration_ms":70968,"significance":"For the Poisson interface model, the paper's analysis is essentially correct and the numerical experiments match the predicted rates. The derivation of Nitsche's method from the penalty formulation is clean and likely useful for teaching, and the explicit constant for linear elements is a practical strength. However, the abstract and title claim results for elastic contact problems, while the analysis and all numerical experiments concern a linear Poisson interface problem with a known interface. The central 'penalty cannot be optimal' conclusion for actual contact problems therefore depends on an unproved transfer from the linear model. This is a significant gap between the presentation and the content, but it is fixable by either narrowing the claims or supplying the missing argument.","major_comments":[{"comment":"The title and abstract promise results for elastic contact problems, but the analysis in Section 2 is explicitly restricted to the Poisson mortar model: 'All the questions we are going to address can be answered by considering \"mortaring\" of the Poisson equation.' The numerical experiments in Section 4 also solve Poisson interface problems, not contact problems. For true contact, the constraint is an inequality and the contact set is unknown, and the paper provides no analogue of estimate (6) in that setting. The load-bearing conclusion that the penalty method cannot be optimally convergent for elastic contact therefore rests on an unproved transfer from the linear model. Please either revise the title and abstract to describe the actual scope, or give a rigorous argument (or a precise theorem from [7]) that the Poisson-mortar estimates imply the same scaling for the contact variational inequality.","section":"Abstract and Section 2"},{"comment":"Equation (6), ||u - u^ε||_E ≲ ε^{1/2} ||∂u/∂n||_{0,Γ}, is the key modeling-error estimate that produces the O(h^{1/2}) ceiling in (10). It is asserted with 'It is not too difficult to show' but no proof or reference is supplied. Since this estimate is load-bearing for the central negative result, please include a short proof (e.g., testing the error equation with u - u^ε and using the interface condition) or provide a precise citation.","section":"Section 2, Eq. (6)"},{"comment":"There is a quantitative error in the interpolation estimate: substituting ε≈h into (8) gives a bound of order h^{s-1}||u^ε||_s, not h^s||u^ε||_s as stated in (9). The same h^s appears in (10). While this correction does not change the O(h^{1/2}) ceiling (for s≥2, h^{s-1} is still smaller than h^{1/2}), the displayed estimates and the phrase 'optimal convergence rate' in (9) should be corrected to h^{s-1}, or the total bound should be stated as O(h^{min(1/2, s-1)}).","section":"Section 2, Eqs. (8)-(10)"}],"minor_comments":[{"comment":"The title contains typographical artifacts ('CONT ACT', 'PENAL TY'); also, in Figures 2, 5, and 8 the convergence-rate labels appear as 'O(N 1/2)' and 'O(N 1)' without minus signs, which is confusing. Please ensure the exponents are displayed correctly.","section":"Title and figures"},{"comment":"The sentence 'the condition number also increases from O(h^{-2})' is vague; please specify the dependence on ε (or on the ratio γ) so that the reader knows the actual growth for the recommended choice ε ≈ h.","section":"Section 2, paragraph after Eq. (10)"},{"comment":"The a posteriori estimator in Eq. (11) appears to contain a missing factor: the first term should presumably be h_E times a residual or jump norm squared. Please check and clarify the notation.","section":"Section 2, last paragraph"},{"comment":"In the concluding remarks, 'it is pity' should be 'it is a pity'.","section":"Section 5"},{"comment":"The deferral to [7] for variational-inequality technicalities would be more useful if a specific theorem or section number from [7] were cited.","section":"Section 2, unknown-contact discussion"}],"recommendation":"major_revision","confidential_remarks":"The core analysis of the Poisson mortar problem appears sound, and the numerical experiments support the conclusions for that model. The main issue is the mismatch between the announced elastic-contact setting and the actual linear-model analysis. If the authors restrict the title and abstract to the interface/mortaring problem and add a proof of (6), the paper would be much stronger; otherwise, the elastic-contact claim is not established. The self-citation to [6] and [8] is appropriate and not a concern."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper is a clear, expository note. The core message is correct: for the Poisson mortar problem, the penalty method with ε ≈ h gives a total error O(h^{1/2}) because the modeling error dominates, while Nitsche's method, obtained by a simple consistency correction, is optimally convergent. The derivation of Nitsche from the penalty form is elegant and explicit, and the stability recipe with the elementwise constant for linear elements is practically useful. The numerics match the predicted rates, including the adaptive experiments that show refinement along the interface for penalty but not for Nitsche.\n\nWhat is actually new is not the mathematical result — the suboptimality of penalty and optimality of Nitsche are well established, including in the authors' own cited work. The contribution is the clean, unified presentation and the explicit stability threshold. That has real value for practitioners and students who want to understand why the penalty parameter cannot be tuned to get optimal convergence.\n\nNow the soft spots. The paper's title and abstract promise elastic contact, but the analysis is entirely for the Poisson mortar problem. Section 2 says so explicitly, and the numerics are on Poisson problems, not contact. The leap to elastic contact is deferred to reference [7], which is a separate published paper, so this is a scope mismatch rather than a technical error. Still, as written, the headline claim is not supported by the analysis in this paper. Second, the modeling error bound (6), which is the load-bearing step for the O(h^{1/2}) ceiling, is asserted without proof. Within the Poisson mortar model it is true, as the stress-test note confirms, but it should be proved or at least clearly cited. Minor typo: equation (9) says h^s while equation (8) gives h^{s-1} for the interpolation contribution; this does not affect the conclusion.\n\nThe paper is for readers who want a concise, reliable account of why the standard penalty method is limited and how Nitsche fixes it. For researchers already in contact mechanics, there is little new. But as an expository piece with correct analysis within its stated model, it deserves a serious referee, provided the authors either extend the analysis to contact or adjust the title and abstract and add the missing proof of (6). I would engage with it, and I would recommend conditional acceptance after revision.","headline":"A clean, well-written note on why the penalty method caps at O(h^{1/2}) for Poisson mortar problems and Nitsche restores optimality, but the title and abstract overstate the elastic-contact scope.","tokens_in":8645,"tokens_out":1505,"would_cite":true,"duration_ms":16709,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N30","65N15","74M15"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that the finite-element penalty method for elastic contact, when tuned with $\\varepsilon$ proportional to the mesh size $h$, cannot do better than an $O(h^{1/2})$ energy-norm error even for smooth solutions, whereas…","keywords":["penalty method","Nitsche method","elastic contact","mortar method","finite element error analysis","interface consistency","adaptive refinement","trace inverse inequality"],"falsifier":"In the paper's own Poisson mortar test, measure the energy-norm error $\\|u-u^\\varepsilon_h\\|_E$ for $\\varepsilon = ch$ over successive uniform refinements of a smooth solution: the bound (10) predicts a slope of $1/2$, so each halving of $h$ should reduce the error by only about 30%. A measured slope closer to 1 would refute the dominating-modelling-error ceiling. For the contact claim itself, run the same experiment on a genuine two-body unilateral contact problem with a known smooth solution and check whether the penalty error slope remains $1/2$.","tokens_in":7597,"feed_emoji":"🧮","tokens_out":6792,"duration_ms":70392,"temperature":0.7,"pith_summary":"The paper argues that, for finite element treatment of elastic contact problems, the classical penalty method cannot deliver optimal convergence no matter how the penalty parameter is chosen, and that Nitsche's method repairs the defect with a small modification. The argument is carried out on the two-domain Poisson 'mortaring' problem, which the authors state is notationally equivalent to elasticity; true unknown-contact technicalities are deferred. With the balancing choice $\\varepsilon \\approx h$, the total error is shown to be dominated by an $O(h^{1/2})$ modelling, or inconsistency, term, so the approximation error's optimal $O(h^s)$ rate is masked. Nitsche's method adds a consistency term and its symmetric counterpart to the penalty form, yielding a stable method with optimal error $O(h^{s-1})$ in the energy norm. Numerical experiments with uniform and adaptive meshes confirm slopes $O(h^{1/2})$ versus $O(h)$.","feed_headline":"Penalty contact methods never converge optimally; Nitsche's fix does","feed_subtitle":"With mesh size h the penalty error stalls at the sqrt of h; Nitsche's correction reaches the full order h^s.","key_machinery":"The engine of the argument is the pair of competing error contributions in the penalty method: the modelling error $\\|u-u^\\varepsilon\\|_E \\lesssim \\varepsilon^{1/2}\\|\\partial u/\\partial n\\|_{0,\\Gamma}$ (equation (6)), which shrinks as $\\varepsilon \\to 0$, and the finite-element interpolation error (8), which contains an $\\varepsilon^{-1/2}h^{s-1/2}$ interface term that grows as $\\varepsilon$ shrinks. Their trade-off fixes $\\varepsilon \\approx h$ and caps the total error at order $h^{1/2}$. Nitsche's machinery is the consistency correction: adding $-\\int_\\Gamma \\frac{\\partial u_1}{\\partial n_1}(v_1-v_2)\\,ds$, valid because the exact solution satisfies $u_1=u_2$ on $\\Gamma$, followed by its symmetric counterpart, makes the discrete form consistent with the exact problem. Stability then comes from the elementwise trace inverse inequality (28), which controls the boundary normal-derivative trace by the element interior gradient scaled by $h_E^{-1/2}$, and this dictates the local parameter choice $\\varepsilon|_E \\lesssim h_E$.","core_discovery":"For the penalized Poisson mortar problem, the exact penalized solution $u^\\varepsilon$ satisfies the modelling-error bound $\\|u-u^\\varepsilon\\|_E \\lesssim \\varepsilon^{1/2}\\|\\partial u/\\partial n\\|_{0,\\Gamma}$ (equation (6)), while the finite-element interpolation bound (8) contains an interface term proportional to $\\varepsilon^{-1/2}h^{s-1/2}$; balancing these two forces the choice $\\varepsilon \\approx h$ and leaves the total error at $h^{1/2}\\|\\partial u/\\partial n\\|_{0,\\Gamma} + h^s\\|u^\\varepsilon\\|_s$ (equation (10)). Hence optimal order is unattainable for the penalty method. Nitsche's method is derived by adding the flux consistency term $-\\int_\\Gamma \\frac{\\partial u_1}{\\partial n_1}(v_1-v_2)\\,ds$ to the penalty form and then symmetrizing with $-\\int_\\Gamma \\frac{\\partial v_1}{\\partial n_1}(u_1-u_2)\\,ds$; it is consistent by construction and, with $\\varepsilon|_E = h_E/(2\\gamma)$ for $\\gamma > C_I$ (or the explicit ratio (36) for linear elements), coercive with respect to the energy norm, yielding the optimal error bound (32).","pith_inferences":["The paper proves the $O(h^{1/2})$ ceiling for the Poisson mortar problem, not for unilateral contact with an unknown contact boundary; the headline claim transfers to elasticity and contact only through the Section 2 assertion, so the transfer should be checked against the referenced contact analysis before relying on it in a production contact code.","The failure mode suggests a repair the authors do not explore: since the modelling error is concentrated on the interface $\\Gamma$, one could keep $\\varepsilon \\approx h$ only there and refine in space, or drive $\\varepsilon$ by the a posteriori estimator; whether this restores optimality is untested here.","The same consistency-correction recipe should apply to any non-conforming penalty form, such as artificial Robin boundary conditions, yielding a family of Nitsche-like methods from a single template.","A practical diagnostic follows from the numerics: an adaptive penalty computation refines the entire interface, while Nitsche refines only singular corners, so the refinement pattern itself reveals which method is in use."],"forward_implications":["For any fixed polynomial degree $p$ and smooth data, the penalty method with $\\varepsilon \\approx h$ has total energy error $O(h^{1/2})$; halving $h$ reduces the error only by a factor of about $\\sqrt{2}$, so the method cannot deliver high accuracy.","Nitsche's method achieves $O(h^{s-1})$ for $1 \\le s \\le p+1$, the same rate as the underlying finite-element interpolation, with the penalty-like parameter chosen locally as $\\varepsilon|_E = h_E/(2\\gamma)$, $\\gamma > C_I$.","For piecewise linear elements the stability constant is explicit: choosing $\\varepsilon|_E = \\alpha (\\int_K dx)/(2\\int_E ds)$ with $\\alpha < 1$ gives uniform stability without solving an eigenvalue problem.","The a posteriori error estimate for the penalty method contains the interface term (11), and that term is efficient---bounded by the true error---only when $\\varepsilon \\approx h$, giving a second, estimator-based justification for the same tuning.","The Nitsche formulation can be written as a one-sided master-slave method or, by symmetrization with average fluxes, in a form that handles differing material parameters in the two subdomains."],"supporting_citations":[{"why":"Supplies the interpolation estimate (8) with the $\\varepsilon^{-1/2}h^{s-1/2}$ interface term that forces the $\\varepsilon \\approx h$ balance.","marker":"[1]"},{"why":"Supplies the a priori error analysis of Nitsche's mortar method and the reliability and efficiency of the interface error estimator.","marker":"[6]"},{"why":"The paper to which the extension to elastic contact and unknown-contact technicalities is deferred, and the source for the convex-flux formulation for different materials.","marker":"[7]"},{"why":"Gives the a posteriori interface term (11) for the penalty method and the result that it is efficient only when $\\varepsilon \\approx h$.","marker":"[8]"},{"why":"The original Nitsche method that the paper rederives as a consistency correction of the penalty method.","marker":"[11]"},{"why":"Provides the explicit constant in the trace inverse inequality (37) used to set the stability parameter for higher-order elements.","marker":"[17]"}],"fun_headline_variants":["Penalty contact error limited to sqrt(h); Nitsche achieves full order","Penalty method can't reach optimal order; Nitsche's correction does","Elastic contact: penalty stalls at half-order, Nitsche converges optimally","Why penalty contact fails to converge optimally; Nitsche's fix works","Optimal convergence impossible with penalty; Nitsche fixes it"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion that penalty cannot be optimal rests on a stated-but-unproved estimate of the modelling error (equation (6)), together with the Section 2 assertion that elasticity and unknown-contact technicalities add only notation to the Poisson interface analysis.","fun_headline_variants_meta":{"raw":{"variants":["Penalty contact error limited to sqrt(h); Nitsche achieves full order","Penalty method can't reach optimal order; Nitsche's correction does","Elastic contact: penalty stalls at half-order, Nitsche converges optimally","Why penalty contact fails to converge optimally; Nitsche's fix works","Optimal convergence impossible with penalty; Nitsche fixes it"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000788,"raw_usage":{"total_tokens":3452,"prompt_tokens":902,"completion_tokens":2550,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":518,"completion_tokens_details":{"reasoning_tokens":2453}},"tokens_in":518,"tokens_out":2550,"duration_ms":19985,"temperature":1.0,"reasoning_tokens":2453,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:23:37.906852+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In the paper's own Poisson mortar test, measure the energy-norm error $\\|u-u^\\varepsilon_h\\|_E$ for $\\varepsilon = ch$ over successive uniform refinements of a smooth solution: the bound (10) predicts a slope of $1/2$, so each halving of $h$ should reduce the error by only about 30%. A measured slope closer to 1 would refute the dominating-modelling-error ceiling. For the contact claim itself, run the same experiment on a genuine two-body unilateral contact problem with a known smooth solution and check whether the penalty error slope remains $1/2$.","supporting_citations":[{"cited_title":"The finite element method with penalty.Math","cited_arxiv_id":null,"evidence_quote":"Supplies the interpolation estimate (8) with the $\\varepsilon^{-1/2}h^{s-1/2}$ interface term that forces the $\\varepsilon \\approx h$ balance."},{"cited_title":"Error analysis of Nitsche’s mor- tar method","cited_arxiv_id":null,"evidence_quote":"Supplies the a priori error analysis of Nitsche's mortar method and the reliability and efficiency of the interface error estimator."},{"cited_title":"On Nitsche’s method for elastic contact problems","cited_arxiv_id":null,"evidence_quote":"The paper to which the extension to elastic contact and unknown-contact technicalities is deferred, and the source for the convex-flux formulation for different materials."}],"review_version":1}