Pith. sign in

REVIEW 2 major objections 2 minor 31 references

Convergence of the Immersed Interface Method in Linear Elasticity

T0 review · 2 major / 2 minor · reviewed 2026-05-23 · grok-4.3

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The L2 difference between exact and quadratured interface force solutions in linear elasticity equals the quadrature error order.

desk verdict The paper proves the L2 difference between exact-integral and quadrature-forced linear elasticity solutions is order of the quadrature error, using fundamental solutions and singularity removal, for both bounded and unbounded domains. read the letter →

arxiv 2410.10436 v2 pith:ODZNX64T submitted 2024-10-14 math.NA cs.NA

Sabia Asghar , Qiyao Peng , Etelvina Javierre , Fred J. Vermolen This is my paper
classification math.NAcs.NA
keywords immersedinterfacemethodlinearelasticityquadratureerrorconvergenceanalysisfundamentalsolutionssingularityremoval
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper shows that replacing the exact integral force along an interface with a numerical quadrature rule produces a solution whose L2-norm difference from the exact-force solution is controlled directly by the quadrature error. This equivalence is proved for both bounded domains and domains extending to infinity. The argument proceeds by representing solutions via fundamental solutions of linear elasticity, removing the singularities that arise at force points, and applying an extended trace theorem to pass from the interface to the volume. A reader cares because the immersed interface method depends on such quadrature approximations, and the result indicates that solution accuracy tracks quadrature accuracy without additional degradation.

What carries the argument

Fundamental solutions of linear elasticity combined with singularity removal and the extended trace theorem, used to bound the volume L2 difference between exact-integral and quadrature forces.

What would settle it

A concrete calculation on a known interface (for example a circle or sphere) with midpoint quadrature of n points where the L2 solution difference fails to scale as the known quadrature error O(1/n) when n is increased while the domain discretization is held fixed.

Watch

Extended reading notes

Core claim

We consider two linear elasticity problems that differ only in how the distributed normal force on an internal interface is represented: one uses the exact integral over the interface, the other replaces that integral by a quadrature rule. We prove that the L2-norm of the difference between the two solutions is of the same order as the quadrature error. The result holds for both bounded and unbounded domains. The proof rests on the fundamental solutions of linear elasticity together with a singularity-removal principle and the extended trace theorem. Convergence is obtained in the L2-norm on curves and manifolds.

Load-bearing premise

Fundamental solutions of linear elasticity, despite being singular and outside H1, can be combined with singularity removal and the extended trace theorem to equate the solution difference to the quadrature error.

Editorial extensions

If this is right

  • The L2 error bound holds uniformly for both bounded and unbounded computational domains.
  • Convergence of the approximated-force solution to the exact-force solution occurs in the L2-norm on the interface curves and manifolds themselves.
  • The estimated error concerns the exact solutions of the two elasticity problems and does not yet include discretization error from a finite-element scheme.

Reading between the lines

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

  • The same quadrature-to-solution error relation may apply to other elliptic interface problems once analogous fundamental solutions and trace theorems are available.
  • Numerical implementations could therefore focus quadrature refinement on the interface independently of the volume mesh to achieve a target accuracy.
  • The result supplies a theoretical justification for choosing quadrature rules whose error is known a priori when designing immersed-interface codes for elasticity.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 2 minor

Summary. The paper claims to prove that for linear elasticity problems with forces distributed over a closed polyhedral interface (in bounded or unbounded domains), the L²-norm of the difference between the solution with exact integral force and the solution with midpoint-quadrature approximation of that force is of the same order as the quadrature error. The proof strategy invokes the Kelvin fundamental solution (Kelvin tensor), the singularity removal principle, and the Extended Trace Theorem; numerical experiments on exact solutions (not FEM approximations) are shown for both domain types.

Significance. If the central claim holds, the result supplies a quadrature-error analysis for the immersed interface method in linear elasticity that is independent of the underlying discretization (FEM or otherwise). The explicit treatment of both bounded and unbounded domains, together with the use of fundamental solutions to handle the interface singularity, would be a useful addition to the IIM literature.

major comments (2)
  1. [unbounded domain section / abstract] Unbounded-domain analysis (abstract and the section presenting the unbounded case): the claim that ||u − u_q||_{L²(ℝ^d)} is finite and of quadrature order presupposes that the difference decays faster than 1/|x|; this holds only if the quadrature exactly preserves the net force ∫_Γ f dS. The manuscript does not state whether the midpoint rule on a closed polyhedral interface achieves exact net-force preservation, nor does it replace the global L² norm by a weighted or local norm. This is load-bearing for the unbounded-domain half of the theorem.
  2. [proof section] Proof outline (the paragraph invoking fundamental solutions and the Extended Trace Theorem): the argument relies on the representation u(x) = ∫_Γ G(x,y) f(y) dS(y) with G the Kelvin tensor (singular and not in H¹), yet the manuscript provides only a high-level sketch rather than the explicit estimates that convert the quadrature error into an L² bound on u − u_q after singularity removal. A concrete reference to the precise statement of the Extended Trace Theorem used and the resulting constant would be required.
minor comments (2)
  1. [numerical experiments] The abstract states that convergence is demonstrated on curves and manifolds, but the numerical section only reports experiments; a short table or plot caption clarifying the observed rates versus quadrature order would improve readability.
  2. [introduction / problem statement] Notation: the symbol u_q is introduced without an explicit definition in the main text; a displayed equation defining the quadrature force term would remove ambiguity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address each major comment below.

read point-by-point responses

  1. Referee: [unbounded domain section / abstract] Unbounded-domain analysis (abstract and the section presenting the unbounded case): the claim that ||u − u_q||_{L²(ℝ^d)} is finite and of quadrature order presupposes that the difference decays faster than 1/|x|; this holds only if the quadrature exactly preserves the net force ∫_Γ f dS. The manuscript does not state whether the midpoint rule on a closed polyhedral interface achieves exact net-force preservation, nor does it replace the global L² norm by a weighted or local norm. This is load-bearing for the unbounded-domain half of the theorem.

    Authors: We acknowledge the referee's observation. The global L² norm over ℝ^d is finite only when the net force is exactly preserved; otherwise the far-field decay of u − u_q is insufficient for square-integrability. The manuscript does not discuss net-force preservation under the midpoint rule. In revision we will either prove that the midpoint quadrature exactly preserves the net force for the polyhedral closed interface (under the force assumptions of the paper) or replace the global L² norm by a suitable weighted norm. The abstract and unbounded-domain section will be updated accordingly. revision: yes

  2. Referee: [proof section] Proof outline (the paragraph invoking fundamental solutions and the Extended Trace Theorem): the argument relies on the representation u(x) = ∫_Γ G(x,y) f(y) dS(y) with G the Kelvin tensor (singular and not in H¹), yet the manuscript provides only a high-level sketch rather than the explicit estimates that convert the quadrature error into an L² bound on u − u_q after singularity removal. A concrete reference to the precise statement of the Extended Trace Theorem used and the resulting constant would be required.

    Authors: We agree that the proof sketch is high-level. The revised manuscript will expand the relevant section with the explicit estimates that bound ||u − u_q||_{L²} by the quadrature error after singularity removal via the Kelvin representation. We will also cite the precise statement of the Extended Trace Theorem employed and indicate the dependence of the constant on the data. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: derivation invokes standard external fundamental solutions and theorems without self-referential reduction.

full rationale

The paper derives an L2 error bound between exact integral and quadrature-based force representations in linear elasticity by representing solutions via the Kelvin fundamental solution tensor, applying the singularity removal principle, and invoking the Extended Trace Theorem. These are independent, externally established mathematical objects and results (not derived or fitted within the paper). No step reduces the target L2 estimate to a fitted parameter, self-citation chain, or definitional tautology; the quadrature error order is obtained from standard approximation theory applied to the interface integral. The unbounded-domain claim is supported by the same external representation without presupposing the result. This is a self-contained application of known tools rather than a circular construction.

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

Analysis rests on standard domain and interface regularity assumptions plus known singular properties of fundamental solutions; no fitted parameters or new postulated entities appear.

assumptions (2)
  • domain assumption The domain is an open, bounded, simply connected Lipschitz domain in R^d containing a closed polyhedral surface or polygonal contour interface.
    Explicitly stated as the geometric setting for both problems considered.
  • standard math Fundamental solutions of linear elasticity exhibit singular behavior at points of action and do not belong to H1.
    Invoked as a key ingredient of the proof strategy.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Convergence of the Immersed Interface Method in Linear Elasticity." pith.science (2026). https://pith.science/paper/ODZNX64T

@misc{pith2026241010436,
  author       = {Pith},
  title        = {Pith review of: Convergence of the Immersed Interface Method in Linear Elasticity},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ODZNX64T}},
  note         = {Machine review of arXiv:2410.10436}
}
abstract

We consider an open, bounded, simply connected (Lipschitz) domain in $\mathbb{R}^d$, which contains a closed polyhedral surface or polygonal contour, referred to as the interface. From this interface, forces are exerted in the normal direction. The forces are continuously distributed over the interface, resulting in an integral expression. This features an important characteristic of the immersed interface method. Since the integral cannot be resolved exactly, one relies on numerical quadrature rules to approximate the integral. Therefore, we consider two different linear elasticity problems with forces over a curve or surface (interface) that is located within the (open) domain of computation: (1) The force is defined by an integral over the interface; (2) The force is defined by a quadrature approximation of the integral over the interface. We prove that the ${\bf L}^2$-norm of the difference between the solutions from the two elasticity problems is of the same order as the error of quadrature. The results are demonstrated for both bounded and unbounded domains. The proof that we establish relies on the use of: (i) fundamental solutions for linear elasticity, exhibiting singular behaviors (in particular around points of action) and not being in ${\bf H}^1$, and (ii) on the use of singularity removal principle and the Extended Trace Theorem. Convergence is demonstrated in the ${\bf L}^2$-norm on curves and manifolds. We show some numerical experiments on the basis of fundamental solutions with a Midpoint quadrature rule in an unbounded and a bounded domain. We note that the error that we estimate is for the exact solutions and not for finite element solutions. Hence in the numerical finite element-based simulations, the numerical results contain an additional error due to the finite element approach.

Figures

Figures reproduced from arXiv: 2410.10436 by the authors.

Figure 1
Figure 1. a) Geometry and distribution of mesh points in R 2 . The circular cell boundary (red-circular marker) with radius ϵ is divided into finite number of line-segments (mesh points), at the mid-point of each of which point forces are applied (shown with blue-arrows). The outer circular boundary (black circular marker) with radius R over which we evaluate the two-dimensional Green’s function. b) For three dimensions, the … view at source ↗
Figure 2
Figure 2. Displacement field around a circular cell, whose boundary is indicated by the red circle, in two dimensions over a grid of 10 × 10. The displacement vectors are indicated by the blue arrows. The length of the arrows indicate the magnitude of the displacement, but does not reflect the real magnitude of the displacement. The input parameters can be found in [PITH_FULL_IMAGE:figures/full_fig_p017_2.png] view at source ↗
Figure 3
Figure 3. Geometry in R 3 where a spherical cell is centered at the origin i.e. (0, 0, 0). a) A spherical cell and a line segment which is lying outside the cell (non-intersecting). b) A spherical cell and a plane which is lying outside the cell (non- intersecting). c) A line segment with sub-segments (or -divisions) on it (in context of applying the midpoint rule to the line segment). d) A plane segment divided into sub-rect… view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages

  1. [1]

    Adams, J.J

    R.A. Adams, J.J. Fournier, Sobolev spaces, Elsevier, 2003

    work page 2003

  2. [2]

    Bajpai, J

    A. Bajpai, J. Tong, W. Qian, Y. Peng, W. Chen, The interplay between cell-cell and cell-matrix forces regulates cell migration dynamics, Biophys. J. 117 (2019) 1795–1804

    work page 2019

  3. [3]

    Barrow, Portraits of hippocrates, Med

    M.V. Barrow, Portraits of hippocrates, Med. Hist. 16 (1972) 85–88

    work page 1972

  4. [4]

    Benz, Jr, The jeremiah metzger lecture cancer in the twenty-first century: An inside view from an outsider, Trans

    E.J. Benz, Jr, The jeremiah metzger lecture cancer in the twenty-first century: An inside view from an outsider, Trans. Am. Clin. Climatol. Assoc. 128 (2017) 275–297

    work page 2017

  5. [5]

    Bertoluzza, A

    S. Bertoluzza, A. Decoene, L. Lacouture, S. Martin, Local error estimates of the finite element method for an elliptic problem with a dirac source term, Numer. Methods Partial Differ. Equ. 34 (2017) 97–120. 20

    work page 2017

  6. [6]

    Boon, F.J

    W.M. Boon, F.J. Vermolen, Analysis of linearized elasticity models with point sources in weighted sobolev spaces: applications in tissue contraction, ESAIM: Math. Model. Numer. Anal. 57 (2023) 2349–2370

    work page 2023

  7. [7]

    Braess, Finite elements: Theory, fast solvers, and applications in solid mechanics, Cambridge Uni- versity Press, 2001

    D. Braess, Finite elements: Theory, fast solvers, and applications in solid mechanics, Cambridge Uni- versity Press, 2001

    work page 2001

  8. [8]

    Byrne, D

    H. Byrne, D. Drasdo, Individual-based and continuum models of growing cell populations: a comparison, J. Math. Biol. 58 (2008) 657–687

    work page 2008

Show all 31 references
  1. [9]

    J. Chen, D. Weihs, F.J. Vermolen, A model for cell migration in non-isotropic fibrin networks with an application to pancreatic tumor islets, Biomech. Model. Mechanobiol. 17 (2017) 367–386

    2017

  2. [10]

    Correia, The Practical Handbook of Perioperative Metabolic and Nutritional Care, Academic Press, 2019

    M.I.T. Correia, The Practical Handbook of Perioperative Metabolic and Nutritional Care, Academic Press, 2019

    2019

  3. [11]

    Dallon, J.A

    J.C. Dallon, J.A. Sherratt, P.K. Maini, Mathematical modelling of extracellular matrix dynamics using discrete cells: Fiber orientation and tissue regeneration, J. Theor. Biol. 199 (1999) 449–471

    1999

  4. [12]

    Darby, B

    I. Darby, B. Laverdet, F. Bont´ e, D. Alexis, Fibroblasts and myofibroblasts in wound healing, Clin. Cosmet. Investig. Dermatol. 7 (2014) 301–11

    2014

  5. [13]

    D’Angelo, Finite element approximation of elliptic problems with dirac measure terms in weighted spaces: Applications to one- and three-dimensional coupled problems, SIAM J

    C. D’Angelo, Finite element approximation of elliptic problems with dirac measure terms in weighted spaces: Applications to one- and three-dimensional coupled problems, SIAM J. Numer. Anal. 50 (2012) 194–215

    2012

  6. [14]

    Enoch, D

    S. Enoch, D. Leaper, Basic science of wound healing, Surg. (Oxford) 26 (2007) 31–37

    2007

  7. [15]

    Evers, S.C

    J.H. Evers, S.C. Hille, A. Muntean, Modelling with measures: Approximation of a mass-emitting object by a point source, Math. Biosci. Eng. 12 (2015) 357–373

    2015

  8. [16]

    Gjerde, K

    I.G. Gjerde, K. Kumar, J.M. Nordbotten, A singularity removal method for coupled 1d–3d flow models, Comput. Geosci. 24 (2019) 443–457

    2019

  9. [17]

    Janiszewska, M.C

    M. Janiszewska, M.C. Primi, T. Izard, Cell adhesion in cancer: Beyond the migration of single cells, J. Biol. Chem. 295 (2020) 2495–2505

    2020

  10. [18]

    van Kan, A

    J. van Kan, A. Segal, F.J. Vermolen, Numerical methods in scientific computing, VSSD, 2005

    2005

  11. [19]

    Koppenol, Biomedical implications from mathematical models for the simulation of dermal wound healing, Ph.D

    D. Koppenol, Biomedical implications from mathematical models for the simulation of dermal wound healing, Ph.D. thesis, TU Delft, 2017

    2017

  12. [20]

    K¨ oppl, B

    T. K¨ oppl, B. Wohlmuth, Optimal a priori error estimates for an elliptic problem with dirac right-hand side, SIAM J. Numer. Anal. 52 (2014) 1753–1769

    2014

  13. [21]

    Land´ en, D

    N.X. Land´ en, D. Li, M. St˚ ahle1, Transition from inflammation to proliferation: a critical step during wound healing, Cell. Mol. Life. Sci. 73 (2016) 3861–3885. 21

    2016

  14. [22]

    Lions, E

    J.L. Lions, E. Magenes, Non-homogeneous boundary value problems and applications: Vol. 1, volume 181, Springer Science & Business Media, 2012

    2012

  15. [23]

    Peng, S.C

    Q. Peng, S.C. Hille, Quality of approximating a mass-emitting object by a point source in a diffusion model, Comput. Math. Appl. 151 (2023) 491–507

    2023

  16. [24]

    Q. Peng, F. Vermolen, Numerical methods to compute stresses and displacements from cellular forces: Application to the contraction of tissue, J. Comput. Appl. Math. 404 (2022) 113892

    2022

  17. [25]

    Q. Peng, F. Vermolen, Point forces in elasticity equation and their alternatives in multi dimensions, Math. Comput. Simul. 199 (2022) 182–201

    2022

  18. [26]

    Peng, F.J

    Q. Peng, F.J. Vermolen, D. Weihs, Physical confinement and cell proximity increase cell migration rates and invasiveness: A mathematical model of cancer cell invasion through flexible channels, J. Mech. Behav. Biomed. Mater. 142 (2023) 105843

    2023

  19. [27]

    L. Richardson, The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam, Philos. Trans. R. Soc. A 210 (1911) 307–357

    1911

  20. [28]

    Susanne C

    L.R.S. Susanne C. Brenner, The Mathematical Theory of Finite Element Methods, Springer, 2008

    2008

  21. [29]

    Tun¸ c, G.J

    B. Tun¸ c, G.J. Rodin, T.E. Yankeelov, Implementing multiphysics models in fenics: Viscoelastic flows, poroelasticity, and tumor growth, Adv. Biomed. Eng. 5 (2023) 100074

    2023

  22. [30]

    W. Yu, S. Sharma, E. Rao, A.C. Rowat, J.K. Gimzewski, D. Han, J. Rao, Cancer cell mechanobiology: a new frontier for cancer research, J. Natl. Cancer Cent. 2 (2022) 10–17

    2022

  23. [31]

    Yue, Biology of the extracellular matrix: An overview, J

    B. Yue, Biology of the extracellular matrix: An overview, J. Glaucoma 23 (2014) 20–23. 22

    2014

Pith tools

Reviewed May 23, 2026 · model on record in the stance chip above.