pith. machine review for the scientific record. sign in

arxiv: 2605.08514 · v1 · submitted 2026-05-08 · 🧮 math.NA · cs.NA

Recognition: 2 theorem links

· Lean Theorem

On a PDE-based material parameter identification problem with contact constraints

Ekaterina Sherina, Simon Hubmer, Stefan Kindermann

Pith reviewed 2026-05-12 00:48 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords parameter identificationcontact constraintsinverse PDE problemuniquenessobstacle problemmembrane modelnumerical reconstructionmaterial coefficient
0
0 comments X

The pith

Contact constraints in a PDE parameter identification problem create both unique and non-unique reconstructions of the scalar coefficient.

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

The paper studies the recovery of a scalar material coefficient inside an elliptic PDE whose solution is forced to stay above an obstacle, modeling a membrane pressed against a fixed barrier. It shows that the contact condition can make the coefficient uniquely recoverable from data in some geometries and completely non-unique in others. Reconstruction methods are constructed and run on synthetic examples that illustrate both the identifiable and the ambiguous regimes. The setup is offered as a clean benchmark for testing algorithms on more realistic contact inverse problems.

Core claim

We consider the identification of a scalar coefficient in a PDE-based parameter estimation problem with contact constraints. The considered problem can be used as an idealized model of a membrane under forces, constrained by a barrier or indenter. We discuss both the forward and inverse parameter estimation problems, as well as uniqueness and non-uniqueness issues caused by the contact constraints. Furthermore, we consider the design and implementation of reconstruction approaches which we test on numerical examples, illustrating both uniqueness and non-uniqueness as well as parameter identifiability.

What carries the argument

The obstacle-constrained PDE together with the inverse map from observed data to the unknown scalar coefficient.

If this is right

  • When contact occurs over a sufficiently large region the coefficient becomes non-unique because variations inside the contact set are invisible to the data.
  • In geometries where contact is only partial the coefficient remains uniquely determined by the measurements.
  • Numerical reconstruction schemes can be built that succeed precisely when the identifiability analysis predicts uniqueness.
  • The same framework supplies a controlled test bed for developing regularization strategies that handle the non-unique regime.
  • The model isolates the effect of the contact constraint without additional modeling complications.

Where Pith is reading between the lines

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

  • Varying the shape or position of the barrier could be used to design experiments that restore uniqueness even when full contact would otherwise occur.
  • The same non-uniqueness mechanism is likely to appear in time-dependent or nonlinear contact problems, suggesting a general regularization principle based on contact-region detection.
  • The benchmark could be used to compare different data-acquisition strategies, such as multiple indenter positions, for their ability to guarantee identifiability.

Load-bearing premise

The idealized PDE model with contact constraints captures the essential features of more complex physical contact problems.

What would settle it

A pair of distinct coefficients that produce identical data under the contact constraint, or a case in which the data uniquely determines the coefficient even though contact occurs over a positive-measure set.

Figures

Figures reproduced from arXiv: 2605.08514 by Ekaterina Sherina, Simon Hubmer, Stefan Kindermann.

Figure 2.1
Figure 2.1. Figure 2.1: Illustrative example of a constraint function (= indenter [PITH_FULL_IMAGE:figures/full_fig_p003_2_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: illustrates an instance of [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3.1
Figure 3.1. Figure 3.1: Left: constraint function h. Right: solution of div(a∇u) = 0. Inside the annulus contact region, the coefficient a is not identifiable. u(x) Kc a(x) c1 c2 [PITH_FULL_IMAGE:figures/full_fig_p010_3_1.png] view at source ↗
Figure 3.2
Figure 3.2. Figure 3.2: Illustration of a non-unique parameter a in a 1D contact problem: The solution is piecewise constant over the constraint given by h (blue dashed object). Below is an illustration of the non-uniqueness of a: It is piecewise constant with smooth connection inside the contact area. The constants c1, c2 > 0 can be chosen arbitrarily. is piecewise constant with contact area Kc. The Lagrange multipliers are δ￾… view at source ↗
Figure 6.1
Figure 6.1. Figure 6.1: Log of err∞ (left) and log of errKKT (right) vs. the number of iterations for TestCase 1 and a discretization of n = 40, for the barrier method (blue), Nesterov￾accelerated projected gradient method (red), and projected gradient method (yellow) [PITH_FULL_IMAGE:figures/full_fig_p019_6_1.png] view at source ↗
Figure 6.2
Figure 6.2. Figure 6.2: Indenter h (left) and solution (right) for TestCase 2. function of a circle with center (0.5, 0.5) and radius 0.25, and “cut off” its top by the linear function g(x1, x2) = x1; see [PITH_FULL_IMAGE:figures/full_fig_p019_6_2.png] view at source ↗
Figure 6.3
Figure 6.3. Figure 6.3: Experiment 1: Top: ground truth and 3 different indente [PITH_FULL_IMAGE:figures/full_fig_p022_6_3.png] view at source ↗
Figure 6.4
Figure 6.4. Figure 6.4: Experiment 2: reconstruction from partial observatio [PITH_FULL_IMAGE:figures/full_fig_p024_6_4.png] view at source ↗
Figure 6.5
Figure 6.5. Figure 6.5: Experiment 3: multiple measurements. Left: ground tru [PITH_FULL_IMAGE:figures/full_fig_p024_6_5.png] view at source ↗
Figure 6.6
Figure 6.6. Figure 6.6: Experiment 3: convergence curves for the Nesterov- [PITH_FULL_IMAGE:figures/full_fig_p025_6_6.png] view at source ↗
read the original abstract

We consider the identification of a scalar coefficient in a PDE-based parameter estimation problem with contact constraints. The considered problem can be used as an idealized model of a membrane under forces, constrained by a barrier or indenter. More generally, it serves as a benchmark for the analysis of more complex contact problems and the development of corresponding reconstruction algorithms. In this paper, we discuss both the forward and inverse parameter estimation problems, as well as uniqueness and non-uniqueness issues caused by the contact constraints. Furthermore, we consider the design and implementation of reconstruction approaches which we test on numerical examples, illustrating both uniqueness and non-uniqueness as well as parameter identifyability.

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

0 major / 3 minor

Summary. The paper formulates and analyzes a scalar coefficient identification problem for an elliptic PDE subject to contact constraints, modeled via a variational inequality. It examines the forward problem, derives conditions for uniqueness and non-uniqueness in the inverse problem induced by the contact set, designs reconstruction algorithms, and validates the theory through finite-element numerical experiments that exhibit both identifiable and non-identifiable regimes.

Significance. If the claims hold, the work supplies a clean benchmark problem illustrating how inequality constraints alter identifiability in PDE parameter estimation. The combination of variational-inequality analysis with concrete numerical tests on unique and non-unique cases is useful for algorithm development in contact mechanics and PDE-constrained optimization. The manuscript ships reproducible numerical examples that directly illustrate the theoretical distinctions.

minor comments (3)
  1. [Abstract] Abstract: the phrase 'parameter identifyability' contains a typographical error and should read 'parameter identifiability'.
  2. [Numerical experiments] Numerical section: the description of the finite-element discretization and the choice of solver for the variational inequality would benefit from explicit statements of mesh size, polynomial degree, and convergence tolerance to allow direct reproduction of the reported uniqueness/non-uniqueness transitions.
  3. [Uniqueness and non-uniqueness] Uniqueness analysis: the examples demonstrating non-uniqueness could be strengthened by an explicit statement of the contact set (or active set) on which the coefficient becomes invisible to the data.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation for minor revision. The provided summary accurately captures the formulation, analysis of the forward and inverse problems, uniqueness/non-uniqueness results induced by the contact set, and the numerical validation via finite-element experiments.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained

full rationale

The paper formulates the forward problem as a variational inequality for the contact-constrained PDE, analyzes identifiability via explicit examples of uniqueness and non-uniqueness induced by the constraints, and validates reconstruction methods through finite-element numerical tests that exhibit the predicted regimes. No load-bearing step reduces by construction to a fitted parameter renamed as prediction, a self-citation chain, or an imported uniqueness theorem; the central claims rest on direct theoretical analysis of the variational inequality and independent numerical illustration rather than self-referential definitions or ansatz smuggling.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Review based on abstract only; no specific free parameters, axioms, or invented entities are described in sufficient detail to populate the ledger.

pith-pipeline@v0.9.0 · 5402 in / 993 out tokens · 39754 ms · 2026-05-12T00:48:38.866642+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

Reference graph

Works this paper leans on

64 extracted references · 64 canonical work pages

  1. [1]

    Alessandrini

    G. Alessandrini. On the identification of the leading coefficient of an elliptic equation. Boll. Un. Mat. Ital. C (6) , 4(1):87–111, 1985

    work page 1985

  2. [2]

    Ammari, E

    H. Ammari, E. Bonnetier, Y. Capdeboscq, M. Tanter, and M. Fink . Electrical impedance tomog- raphy by elastic deformation. SIAM Journal on Applied Mathematics , 68(6):1557–1573, 2008

    work page 2008

  3. [3]

    Attouch, G

    H. Attouch, G. Buttazzo, and G. Michaille. Variational analysis in Sobolev and BV spaces , volume 17 of MOS-SIAM Series on Optimization . Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA; Mathematical Optimization Society, Philade lphia, PA, second edition,

    work page

  4. [4]

    Applications to PDEs and optimization

    work page

  5. [5]

    G. Bal, F. Monard, and G. Uhlmann. Reconstruction of a fully aniso tropic elasticity tensor from knowledge of displacement fields. SIAM J. Appl. Math. , 75(5):2214–2231, 2015

    work page 2015

  6. [6]

    G. Bal, W. Naetar, O. Scherzer, and J. Schotland. The Levenbe rg-Marquardt iteration for nu- merical inversion of the power density operator. J. Inv. Ill-Posed Problems , 21(2):265–280, 2013

    work page 2013

  7. [7]

    Bal and G

    G. Bal and G. Uhlmann. Reconstruction of coefficients in scalar se cond-order elliptic equa- tions from knowledge of their solutions. Communications on Pure and Applied Mathematic , 66(10):1629–1652, 2013

    work page 2013

  8. [8]

    H. H. Bauschke and P. L. Combettes. Convex analysis and monotone operator theory in Hilbert spaces. Springer, 2017. 26

    work page 2017

  9. [9]

    A. Beck. First-order methods in optimization , volume 25 of MOS-SIAM Series on Optimiza- tion. Society for Industrial and Applied Mathematics (SIAM), Philadelph ia, PA; Mathematical Optimization Society, Philadelphia, PA, 2017

    work page 2017

  10. [10]

    Bensaada, E.-H

    A. Bensaada, E.-H. Essoufi, and A. Zafrar. Primal-dual formula tion for parameter estimation in elastic contact problem with friction. Appl. Math. Sci. Eng. , 32(1):Paper No. 2367025, 25, 2024

    work page 2024

  11. [11]

    Bonnet and A

    M. Bonnet and A. Constantinescu. Inverse problems in elasticit y. Inverse Problems , 21(2):R1– R50, 2005

    work page 2005

  12. [12]

    A. P. Calder´ on. On an inverse boundary value problem. Computational & applied mathematics , 25(2-3):133–138, 2006

    work page 2006

  13. [13]

    Capdeboscq, J

    Y. Capdeboscq, J. Fehrenbach, F. de Gournay, and O. Kavian . Imaging by modification: Numeri- cal reconstruction of local conductivities from corresponding po wer density measurements. SIAM Journal on Imaging Sciences , 2(4):1003–1030, 2009

    work page 2009

  14. [14]

    Cheney, D

    M. Cheney, D. Isaacson, and J. C. Newell. Electrical impedance tomography. SIAM review , 41(1):85–101, 1999

    work page 1999

  15. [15]

    Chicone and J

    C. Chicone and J. Gerlach. A note on the identifiability of distribut ed parameters in elliptic equations. SIAM J. Math. Anal. , 18(5):1378–1384, 1987

    work page 1987

  16. [16]

    Christof and G

    C. Christof and G. Wachsmuth. Differential sensitivity analysis o f variational inequalities with locally Lipschitz continuous solution operators. Appl. Math. Optim. , 81(1):23–62, 2020

    work page 2020

  17. [17]

    Christof and G

    C. Christof and G. Wachsmuth. Semismoothness for solution op erators of obstacle-type variational inequalities with applications in optimal control. SIAM J. Control Optim. , 61(3):1162–1186, 2023

    work page 2023

  18. [19]

    Clason and V

    C. Clason and V. H. Nhu. Bouligand-Landweber iteration for a no n-smooth ill-posed problem. Numer. Math. , 142(4):789–832, 2019

    work page 2019

  19. [20]

    Constantinescu and N

    A. Constantinescu and N. Tardieu. On the identification of elast oviscoplastic constitutive laws from indentation tests. Inverse Problems in Engineering , 9(1):19–44, 2001

    work page 2001

  20. [21]

    H. W. Engl, M. Hanke, and A. Neubauer. Regularization of inverse problems. Dordrecht: Kluwer Academic Publishers, 1996

    work page 1996

  21. [22]

    H. W. Engl, K. Kunisch, and A. Neubauer. Convergence rates f or tikhonov regularisation of non-linear ill-posed problems. Inverse Problems , 5(4):523, 1989

    work page 1989

  22. [23]

    Fernandez-Zelaia, V

    P. Fernandez-Zelaia, V. J. Roshan, S. R. Kalidindi, and S. N. Melk ote. Estimating mechanical properties from spherical indentation using bayesian approaches . Materials & Design , 147:92–105, 2018

    work page 2018

  23. [24]

    H. Gfrerer. On a globally convergent semismooth* Newton meth od in nonsmooth nonconvex optimzation. Comp. Optim. Appl. , 91:67–124, 2025

    work page 2025

  24. [25]

    Gfrerer, S

    H. Gfrerer, S. Hubmer, and R. Ramlau. On SCD Semismooth ∗ Newton methods for the effi- cient minimization of Tikhonov functionals with non-smooth and non-c onvex penalties. Inverse Problems, 41(7):075002, 2025. Gold OA

    work page 2025

  25. [26]

    Gfrerer, M

    H. Gfrerer, M. Mandlmayr, J. V. Outrata, and J. Valdman. On t he SCD semismooth* Newton method for generalized equations with application to a class of static contact problems with Coulomb friction. Comp. Optim. Appl. , 86:1159–1191, 2022

    work page 2022

  26. [27]

    Gfrerer and J

    H. Gfrerer and J. V. Outrata. On a semismooth* Newton metho d for solving generalized equations. SIAM J. Optim. , 31(1):489–517, 2021. 27

    work page 2021

  27. [28]

    M. S. Gockenbach, B. Jadamba, and A. A. Khan. Equation erro r approach for elliptic inverse problems with an application to the identification of Lam/’e parameter s. Inverse Problems in Science and Engineering , 16(3):349–367, 2008

    work page 2008

  28. [29]

    M. S. Gockenbach and A. A. Khan. An abstract framework for elliptic inverse problems. I. An output least-squares approach. Math. Mech. Solids , 12(3):259–276, 2007

    work page 2007

  29. [30]

    M. S. Gockenbach and A. A. Khan. An abstract framework for elliptic inverse problems. II. An augmented Lagrangian approach. Math. Mech. Solids , 14(6):517–539, 2009

    work page 2009

  30. [31]

    J. Gwinner. An optimization approach to parameter identificatio n in variational inequalities of second kind. Optim. Lett. , 12(5):1141–1154, 2018

    work page 2018

  31. [32]

    J. Gwinner. An optimization approach to parameter identificatio n in variational inequalities of second kind—II. In Deterministic and stochastic optimal control and inverse p roblems, pages 112–130. CRC Press, Boca Raton, FL, 2022

    work page 2022

  32. [33]

    Gwinner, B

    J. Gwinner, B. Jadamba, A. A. Khan, and F. Raciti. Uncertainty quantification in variational inequalities—theory, numerics, and applications . CRC Press, Boca Raton, FL, [2022] ©2022

    work page 2022

  33. [34]

    Gwinner, B

    J. Gwinner, B. Jadamba, A. A. Khan, and M. Sama. Identificatio n in variational and quasi- variational inequalities. J. Convex Anal. , 25(2):545–569, 2018

    work page 2018

  34. [35]

    Hanke, A

    M. Hanke, A. Neubauer, and O. Scherzer. A convergence ana lysis of the Landweber iteration for nonlinear ill-posed problems. Numerische Mathematik , 72(1):21–37, 1995

    work page 1995

  35. [36]

    Harrach and M

    B. Harrach and M. Ullrich. Monotonicity-based shape reconstr uction in electrical impedance tomography. SIAM J. Math. Anal. , 45(6):3382–3403, 2013

    work page 2013

  36. [37]

    Hofmann, B

    B. Hofmann, B. Kaltenbacher, C. P¨ oschl, and O. Scherzer. A convergence rates result for Tikhonov regularization in Banach spaces with non-smooth operators. Inverse Problems , 23(3):987–1010, 2007

    work page 2007

  37. [38]

    Holder and A

    D. Holder and A. Adler. Electrical impedance tomography: methods, history and app lications. Series in medical physics and biomedical engineering. CRC Press, sec ond edition edition, 2022

    work page 2022

  38. [39]

    Hubmer, K

    S. Hubmer, K. Knudsen, C. Li, and E. Sherina. Limited-angle aco usto-electrical tomography. Inverse Problems in Science and Engineering , 27(9):1298–1317, 2019

    work page 2019

  39. [40]

    Hubmer and R

    S. Hubmer and R. Ramlau. Convergence analysis of a two-point g radient method for nonlinear ill-posed problems. Inverse Problems , 33(9):095004, 2017

    work page 2017

  40. [41]

    Hubmer, E

    S. Hubmer, E. Sherina, and R. Ramlau. Characterizations of ad joint Sobolev embedding operators with applications in inverse problems. Electron. Trans. Numer. Anal. , 59:116–144, 2023

    work page 2023

  41. [42]

    V. Isakov. Inverse Problems for Partial Differential Equations. Applied Mathematical Sciences. Springer, New York, NY, 2006

    work page 2006

  42. [43]

    Jadamba, A

    B. Jadamba, A. A. Khan, F. Raciti, and M. Sama. A variational ine quality based stochastic approximation for estimating the flexural rigidity in random fourth- order models. Commun. Nonlinear Sci. Numer. Simul. , 111:Paper No. 106406, 11, 2022

    work page 2022

  43. [44]

    Jin and Y

    Q. Jin and Y. Zhang. Bouligand-Newton type methods for non-s mooth ill-posed problems. Inverse Problems, 40(10):Paper No. 105011, 28, 2024

    work page 2024

  44. [45]

    Kaltenbacher

    B. Kaltenbacher. Regularization based on all-at-once formulat ions for inverse problems. SIAM J. Numer. Anal. , 54(4):2594–2618, 2016

    work page 2016

  45. [46]

    Kaltenbacher

    B. Kaltenbacher. Minimization based formulations of inverse pro blems and their regularization. SIAM J. Optim. , 28(1):620–645, 2018. 28

    work page 2018

  46. [47]

    Kaltenbacher, A

    B. Kaltenbacher, A. Kirchner, and B. Vexler. Goal oriented ad aptivity in the IRGNM for pa- rameter identification in PDEs: II. all-at-once formulations. Inverse Problems , 30(4):045002, 33, 2014

    work page 2014

  47. [48]

    Kaltenbacher, A

    B. Kaltenbacher, A. Neubauer, and O. Scherzer. Iterative regularization methods for nonlinear ill-posed problems. Berlin: de Gruyter, 2008

    work page 2008

  48. [49]

    Kaltenbacher, F

    B. Kaltenbacher, F. Sch¨ opfer, and T. Schuster. Iterative methods for nonlinear ill-posed problems in Banach spaces: convergence and applications to parameter iden tification problems. Inverse Problems, 25(6):065003 (19pp), 2009

    work page 2009

  49. [50]

    A. A. Khan and D. Motreanu. Inverse problems for quasi-varia tional inequalities. J. Global Optim., 70(2):401–411, 2018

    work page 2018

  50. [51]

    Kikuchi and J

    N. Kikuchi and J. T. Oden. Contact problems in elasticity: a study of variational ineq ualities and finite element methods , volume 8 of SIAM Studies in Applied Mathematics . Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 1988

    work page 1988

  51. [52]

    Kindermann

    S. Kindermann. On the tangential cone condition for electrical impedance tomography. Electron. Trans. Numer. Anal. , 57:17–34, 2022

    work page 2022

  52. [53]

    R. V. Kohn and M. Vogelius. Relaxation of a variational method fo r impedance computed to- mography. Comm. Pure Appl. Math. , 40(6):745–777, 1987

    work page 1987

  53. [54]

    Y. Ma, F. Liu, Z. Si, Q. Gerile, K. Yang, H. Li, and Y. Cheng. Param eter identification of contact model for tank saddle ring system based on genetic algorithm and bp neural network. Journal of Physics: Conference Series , 2478(12):122017, jun 2023

    work page 2023

  54. [55]

    Mueller and S

    J. Mueller and S. Siltanen. Linear and Nonlinear Inverse Problems with Practical Appli cations. Society for Industrial and Applied Mathematics, Philadelphia, PA, 20 12

    work page

  55. [56]

    Nesterov

    Y. Nesterov. A method of solving a convex programming problem with convergence rate O(1/k 2). Soviet Mathematics Doklady , 27(2):372–376, 1983

    work page 1983

  56. [57]

    Radov´ a and J

    J. Radov´ a and J. Machalov´ a. Parameter identification in cont act problems for Gao beam. Non- linear Anal. Real World Appl. , 77:Paper No. 104068, 22, 2024

    work page 2024

  57. [58]

    Rauls and G

    A.-T. Rauls and G. Wachsmuth. Generalized derivatives for the s olution operator of the obstacle problem. Set-Valued Var. Anal., 28(2):259–285, 2020

    work page 2020

  58. [59]

    G. R. Richter. An inverse problem for the steady state diffusion equation. SIAM J. Appl. Math. , 41(2):210–221, 1981

    work page 1981

  59. [60]

    A. Rieder. On the regularization of nonlinear ill-posed problems via inexact Newton iterations. Inverse Problems , 15(1):309–327, 1999

    work page 1999

  60. [61]

    Sahin, M

    T. Sahin, M. von Danwitz, and A. Popp. Solving forward and inver se problems of contact mechan- ics using physics-informed neural networks. Advanced Modeling and Simulation in Engineering Sciences, 11:Paper No. 11, 30, 2024

    work page 2024

  61. [62]

    Scherzer, M

    O. Scherzer, M. Grasmair, H. Grossauer, M. Haltmeier, and F. Lenzen. Variational methods in imaging, volume 167 of Applied Mathematical Sciences . Springer, New York, 2009

    work page 2009

  62. [63]

    C. Sun, Q. Lin, B. Wang, J. Huang, and J. Chen. Modeling and inve rse identification method for the characterization of elastic-plastic contact behavior durin g flat punch indentation. Eur. J. Mech. A Solids , 96:Paper No. 104663, 15, 2022

    work page 2022

  63. [64]

    Tardieu and A

    N. Tardieu and A. Constantinescu. On the determination of elas tic coefficients from indentation experiments. Inverse Problems , 16(3):577–588, 2000

    work page 2000

  64. [65]

    Zhang and L

    H. Zhang and L. V. Wang. Acousto-electric tomography. In Proceedings SPIE, Photons Plus Ultrasound: Imaging and Sensing, 5320(9) , page 145, 2004. 29

    work page 2004