Pith. sign in

REVIEW 2 major objections 6 minor 31 references

Elastic Interaction of Pressurized Cavities in Hyperelastic Media: Attraction and Repulsion

T0 review · 2 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash

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

Pith's one-line read Two pressurized cavities in a hyperelastic medium switch from attraction to repulsion once the pressure-to-shear-modulus ratio exceeds a critical value.

desk verdict The paper's new repulsive regime is worth taking seriously, but the finite outer boundary at 50R makes the finding unproven until a nonlinear domain-size convergence check is done. read the letter →

arxiv 2501.07664 v1 pith:LF6QJKB3 submitted 2025-01-13 cond-mat.soft

Ali Saeedi , Mrityunjay Kothari This is my paper · ORCID
classification cond-mat.soft
keywords hyperelasticitypressurizedcavitiescavityinteractionenergydrivingforceneo-HookeanmaterialMooney-RivlinArruda-Boycestrainstiffening
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 establishes that two equal pressurized cylindrical cavities in a 2D hyperelastic medium can repel each other, something linear elasticity says never happens. Below a pressure-to-shear-modulus ratio of about one, the interaction is always attractive; above that ratio, the energy landscape becomes non-monotonic, so cavities closer than a critical separation attract while wider-separated cavities repel. The same transition appears in neo-Hookean, Mooney-Rivlin, and Arruda-Boyce materials, with strain stiffening in the Arruda-Boyce model moving the critical separation. If correct, this gives a quantitative lever for arranging or stabilizing cavities in soft materials and underground structures by choosing pressure and material parameters rather than only geometry.

What carries the argument

The central object is the configurational driving force $F = -d\Pi_{\rm eqm}/d\eta$, the negative derivative of the equilibrium total potential energy with respect to the separation $\eta$ between cavity centers. The authors compute $\Pi_{\rm eqm}$ from finite-element solutions as the strain energy minus the pressure-volume work of the two cavities, fit a smooth curve through the discrete energy values, and read the nature of the interaction from the sign of the derivative. The nonlinear constitutive models are what let the energy landscape acquire a maximum at high pressure, while the linear-elastic landscape stays monotone; strain stiffening, especially in the Arruda-Boyce model, shifts the location of that maximum.

What would settle it

Recompute the neo-Hookean two-cavity energy at $P/\mu=1.5$ with the separation sampled every $\Delta(\eta/R)=0.01$ around the reported $(\eta/R)_{\rm critical}$, and take the driving force directly from finite differences of the raw, unsmoothed finite-element energies; if the derivative never changes sign, the repulsion is an artifact of the smoothing.

Watch

Extended reading notes

Core claim

For two equal pressurized circular cavities in an infinite, incompressible hyperelastic medium under plane strain, the total potential energy as a function of center-to-center separation is monotonically decreasing at low pressures ($P/\mu \lesssim 1$), so the configurational driving force $F = -d\Pi_{\rm eqm}/d\eta$ is always negative and the cavities attract. At higher pressures ($P/\mu \gtrsim 1$), the energy curve develops a maximum at a critical separation $(\eta/R)_{\rm critical}$; the driving force is negative for closer spacings and positive for wider ones, meaning sufficiently separated cavities repel. The critical separation shifts with pressure and with strain-stiffening parameters, and the transition appears in neo-Hookean, Mooney-Rivlin, and Arruda-Boyce models, with Arruda-Boyce showing the strongest dependence on its limiting-stretch parameter.

Load-bearing premise

The sign change in the driving force is read from the derivative of a smooth curve fitted through discrete finite-element energy values, so the entire repulsive regime depends on that fitted curve faithfully representing the true energy landscape.

Editorial extensions

If this is right

  • Below $P/\mu \simeq 1$ the two-cavity system behaves like the linear-elastic case: attraction at every separation, so no equilibrium spacing exists.
  • Above the threshold there is an unstable equilibrium at $(\eta/R)_{\rm critical}$; cavities perturbed inward attract and coalesce, while cavities perturbed outward separate further.
  • The critical separation depends on pressure: higher pressure moves the transition, and in Arruda-Boyce materials a smaller limiting stretch $\lambda_m$ lowers the critical separation.
  • Because all three constitutive models show the transition, repulsion at large separation is a generic consequence of nonlinear elasticity rather than a quirk of one material law.
  • The energy-based driving-force formulation gives a direct design target: match the critical separation to a desired tunnel or bubble spacing by choosing the pressure ratio.

Reading between the lines

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

  • If the transition survives direct differentiation of the raw energies, the energy maximum implies an unstable separatrix; in multi-cavity arrays this could organize cavities into repulsive lattices that resist coalescence without walls.
  • A natural extension beyond the paper is to spherical cavities in 3D, where the same energy-maximum mechanism may regulate cavitation damage and phase-separated droplet spacing in soft solids.
  • The finite outer radius used to approximate infinity should be tested against larger domains at high $P/\mu$; if the energy peak shifts or disappears with domain size, the infinite-domain interpretation needs revisiting.
  • The exact value of the reported critical separation is the softest quantitative output because it comes from a smoothed derivative; direct finite differences of the raw finite-element energies would settle how much of the phase diagram is robust.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 6 minor

Summary. The paper studies two equal pressurized cylindrical cavities in a 2D infinite hyperelastic medium under plane strain, using ABAQUS finite elements for neo-Hookean, Mooney-Rivlin, and Arruda-Boyce materials. The interaction is quantified by the driving force F = -dΠ/dη, where Π is the total potential energy and η is the center-to-center separation. The authors validate their computational pipeline against a corrected bipolar-coordinate linear-elastic series solution and against an analytical single-cavity neo-Hookean solution, and they report mesh-convergence studies and raw energy data. Their central claim is that at low pressure-shear modulus ratios the interaction is always attractive, while at higher ratios the energy becomes non-monotonic: below a critical separation the force is attractive and above it repulsive, a regime absent in linear elasticity. The paper also reports that Mooney-Rivlin results are nearly independent of the stiffening parameter α, while Arruda-Boyce results depend on λ_m, with smaller λ_m giving a smaller critical separation.

Significance. If the repulsive regime is genuine, the result is counterintuitive and relevant for soft-material design, phase-separation problems, and geomechanics. The paper has clear strengths: it validates the numerical method against an external linear-elastic series solution and a single-cavity analytical benchmark, reports raw potential-energy data in Appendix E, and scans material parameters rather than fitting them to obtain the claimed behavior. The central limitation is that the main observable is obtained by differentiating an unspecified smooth-curve fit to discrete energy values, and the infinite-domain claim is never checked by a nonlinear outer-boundary convergence study. Because both issues bear directly on whether the repulsive regime is a property of the infinite hyperelastic medium or an artifact of the numerical procedure, the result is defensible only after those points are addressed.

major comments (2)
  1. [Section 2.4, Figures 4–5] The driving force F = -dΠ/dη is computed from a smooth curve fitted to discrete finite-element energy values, but the paper does not state the fitting function, the smoothing procedure, or the sensitivity of the derivative to the fit. Since the existence of the repulsive regime and the value of (η/R)_critical are defined by a sign change of this numerical derivative, the reader cannot tell whether the transition is a feature of the underlying FE solution or an artifact of the chosen interpolant. Please report the fitting form, fit residuals for each P/μ, and demonstrate that the sign change and critical separation are insensitive to the fitting procedure.
  2. [Section 2.4, Appendix C] The infinite-medium claim is not checked for outer-boundary effects in the nonlinear regime. All simulations fix R_o/R = 50, and the mesh-convergence study in Appendix C varies mesh density at this fixed outer radius. At P/μ = 1.75, Eq. (B.43) gives λ_a ≈ 3.6, so at η/R = 20 the cavity wall is only about 10 deformed radii from the traction-free outer boundary; moreover, as η/R increases from 2.2 to 20, the cavities move from ±1.1R to ±10R, approaching that boundary by roughly 18% of its initial standoff. A traction-free boundary is softer than the infinite medium and attracts a pressurized cavity, producing a contribution to F that is positive—the same direction as the claimed repulsion. A domain-size study with R_o/R = 100 and 200 at high P/μ and large η/R is required to rule out a finite-domain artifact; the linear-elastic validation cannot detect this effect because it is O((R/R_o)^2) at low pressure and small deformation.
minor comments (6)
  1. [Abstract and Section 4] The sentence 'Effect of strain stiffening on these interactions are also analyzed' has a subject-verb agreement error.
  2. [Figure 4] The caption does not identify which curve corresponds to which P/μ value, making it difficult to connect the plotted potential-energy and driving-force curves to the pressure levels discussed in the text.
  3. [Appendix A] The text says the corrected solution and its MATLAB implementation are included, but no MATLAB code appears in the appendix; only the analytical formulas are given.
  4. [Equation (26)] The notation mixing summation and integration symbols (χ0X, πX) is confusing; please clarify that this is a numerical quadrature rule.
  5. [Figure 12] The axes are unlabeled; the plotted 'Area of deformed cavity' should state the normalization and whether it is the deformed area or the area change relative to the undeformed cavity.
  6. [Appendix C] The statement that 'mesh independence was found to hold across pressures and models' is not supported by the displayed results, which show only a neo-Hookean model at P/μ = 1.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the repulsion claim is extracted from independently validated FE energy data, not assumed by construction.

full rationale

The paper's central claim (attraction at low pressure, non-monotonic energy and repulsion at high pressure) is not forced by its inputs. The potential energy is computed from finite-element simulations for a scanned range of P/mu and eta/R, and no parameter is fitted to produce repulsion: the material parameters and geometry are prescribed, not calibrated to the target behavior. The smooth curve fitted to the energy data is used only to compute the derivative F = -dPi/deta, and the raw, non-fitted energy data are reported in Appendix E, so the sign change is visible in the input data rather than manufactured by the fit. The computational framework is validated against two external benchmarks: the bipolar-coordinate series solution for linear elasticity (Appendix A) and the analytical single-cavity neo-Hookean solution (Appendix B). The only self-citations (Refs. [1,2]) appear in the introduction as application context for phase separation and are not load-bearing for the interaction result. The finite outer-boundary and curve-fitting concerns raised in the reader's take are correctness or robustness risks, not circularity: they concern whether the FE data faithfully represent the infinite domain, but the energy values are still computed from first principles rather than assumed. No equation, fitted parameter, or self-citation chain is equivalent by construction to the claimed attraction-to-repulsion transition, so the paper is self-contained with respect to circularity.

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

No constants are fitted to the target result; P/mu, alpha, and lambda_m are scanned inputs. The central claim rests on four stated modeling and computational assumptions, the most fragile being the undisclosed energy-data fit used to compute the driving force.

assumptions (4)
  • domain assumption Plane-strain, quasi-static, incompressible hyperelasticity with constant internal pressure is an adequate model for the cavity pair.
    Used throughout; footnote in Sec 2.1 notes constant pressure excludes gas-filled cavities, so the scope is explicitly limited.
  • standard math The total potential energy Pi = integral(psi dOmega) - 2P Delta A and the driving force F = -dPi/deta correctly quantify configurational interaction.
    This variational definition is standard in mechanics and is benchmarked against the linear-elastic series solution in Sec 3.1.
  • domain assumption The traction-free external boundary at radius 5 with cavity radius 0.1 approximates an infinite domain.
    Finite element truncation; validated only for the linear case in Fig 3, not for the high-pressure nonlinear cases.
  • ad hoc to paper The discrete potential-energy data can be accurately represented by a smooth curve whose derivative yields the driving force.
    Sec 2.4 fits a smooth curve but gives no fitting form, no error analysis, and no derivative validation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Elastic Interaction of Pressurized Cavities in Hyperelastic Media: Attraction and Repulsion." pith.science (2026). https://pith.science/paper/LF6QJKB3

@misc{pith2026250107664,
  author       = {Pith},
  title        = {Pith review of: Elastic Interaction of Pressurized Cavities in Hyperelastic Media: Attraction and Repulsion},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LF6QJKB3}},
  note         = {Machine review of arXiv:2501.07664}
}
read the original abstract

This study computationally investigates the elastic interaction of two pressurized cylindrical cavities in a 2D hyperelastic medium. Unlike linear elasticity, where interactions are exclusively attractive, nonlinear material models (neo-Hookean, Mooney-Rivlin, Arruda-Boyce) exhibit both attraction and repulsion between the cavities. A critical pressure-shear modulus ratio governs the transition, offering a pathway to manipulate cavity configurations through material and loading parameters. At low ratios, the interactions are always attractive; at higher ratios, both attractive and repulsive regimes exist depending on the separation between the cavities. Effect of strain stiffening on these interactions are also analyzed. These insights bridge theoretical and applied mechanics, with implications for soft material design and subsurface engineering.

Figures

Figures reproduced from arXiv: 2501.07664 by the authors.

Figure 1
Figure 1. Schematic of the two-cavity system in a hyperelastic medium under constant internal pressure P. [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. A representative generated mesh: (a) entire domain and (b) locally refined region for [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. (a) Potential energy and (b) driving force for two pressurized cavities in a linear elastic medium. The results are independent of P/µ. 3.2. Nonlinear Elastic Medium To compare and contrast the effects of nonlinear elasticity and strain-stiffening on the interaction between the cavities, the results are presented in a non-dimensionalized form and the three materials models (NH, MR and AB) are taken to have the same … view at source ↗
Figures from the paper (10 more)
Figure 4
Figure 4. Figure 4: Summary of potential energy and driving force for all the models ( [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: Phase diagram of nonlinear elastic material models ( [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: Comparison of (a) potential energy and (b) driving force for nonlinear elastic material models ( [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 7
Figure 7. Figure 7: Driving force variation for different strain-stiffening parameters ( [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]
Figure 8
Figure 8. Figure 8: Bipolar coordinate system and the geometry used in linear elastic problem ( [PITH_FULL_IMAGE:figures/full_fig_p008_8.png]
Figure 9
Figure 9. Figure 9: Single cavity with reference (blue) and deformed (green) coordinates. [PITH_FULL_IMAGE:figures/full_fig_p011_9.png]
Figure 10
Figure 10. Figure 10: Comparison between analytical and computational solutions for single cavity growth in a neo-Hookean medium. [PITH_FULL_IMAGE:figures/full_fig_p012_10.png]
Figure 11
Figure 11. Figure 11: Refinement levels (η/R = 15) and solution convergence (P/µ = 1). D. Deformed Area [PITH_FULL_IMAGE:figures/full_fig_p013_11.png]
Figure 12
Figure 12. Figure 12: shows the variation of the deformed area as a function of the separation between the cavities for P/µ = 1.5. (a) (b) 5 10 15 20 -1.12 -1.1 -1.08 -1.06 -1.04 -1.02 -1 -0.98 η/R=3 η/R=4 η/R=5 η/R=6 5 10 15 20 0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 η/R=3 η/R=4 η/R=5…
Figure 13
Figure 13. Figure 13: Original potential energy data prior to fitting. References [1] M. Kothari, T. Cohen, Effect of elasticity on phase separation in heterogeneous systems, Journal of the Mechanics and Physics of Solids 145 (2020) 104153. [2] M. Kothari, T. Cohen, The crucial role of ela…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

31 extracted references · 30 canonical work pages

  1. [1]

    Kothari, T

    M. Kothari, T. Cohen, Effect of elasticity on phase separation in heterogeneous systems, Journal of the Mechanics and Physics of Solids 145 (2020) 104153

    work page 2020

  2. [2]

    Kothari, T

    M. Kothari, T. Cohen, The crucial role of elasticity in regulating liquid–liquid phase separation in cells, Biomechanics and Modeling in Mechanobiology 22 (2) (2023) 645–654

    work page 2023

  3. [3]

    C. Kooi, A. Verruijt, Interaction of circular holes in an infinite elastic medium, Tunnelling and Underground Space Technology 16 (2001) 59–62

    work page 2001

  4. [4]

    C. W. Barney, C. E. Dougan, K. R. McLeod, A. Kazemi-Moridani, Y. Zheng, Z. Ye, S. Tiwari, I. Sacligil, R. A. Riggleman, S. Cai, et al., Cavitation in soft matter, Proceedings of the National Academy of Sciences 117 (17) (2020) 9157–9165

    work page 2020

  5. [5]

    Jeffries, M

    J. Jeffries, M. Wall, K. Moore, A. Schwartz, He bubble coarsening by migration and coalescence in annealed pu–ga alloys, Journal of nuclear materials 410 (1-3) (2011) 84–88

    work page 2011

  6. [6]

    Zhang, H

    W. Zhang, H. Han, J. Dai, C. Ren, C. Wang, L. Yan, H. Huang, Z. Zhu, Simulation of migration and coalescence of helium bubbles in nickel, Journal of Nuclear Materials 518 (2019) 48–53

    work page 2019

  7. [7]

    K. Xu, Y. Mehmani, L. Shang, Q. Xiong, Gravity-induced bubble ripening in porous media and its impact on capillary trapping stability, Geophysical Research Letters 46 (23) (2019) 13804–13813

    work page 2019

  8. [8]

    G. B. Jeffery, Ix. plane stress and plane strain in bipolar co-ordinates, Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character 221 (582-593) (1921) 265–293

    work page 1921

Show all 31 references
  1. [9]

    Howland, R

    R. Howland, R. Knight, Stress functions for a plate containing groups of circular holes, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences 238 (793) (1939) 357–392

    1939

  2. [10]

    A. E. Green, General bi-harmonic analysis for a plate containing circular holes, Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences 176 (1940) 121–139

    1940

  3. [11]

    C. B. Ling, On the stresses in a plate containing two circular holes, Journal of Applied Physics 19 (1948) 77–82

    1948

  4. [12]

    Toshihiro, M

    I. Toshihiro, M. Kazyu, Stress concentrations in a plate with two unequal circular holes, International Journal of Engi- neering Science 18 (8) (1980) 1077–1090

    1980

  5. [13]

    S. K. Hoang, Y. N. Abousleiman, Extended green’s solution for the stresses in an infinite plate with two equal or unequal circular holes, Journal of Applied Mechanics 75 (2008) 031016

    2008

  6. [14]

    Davanas, Analysis of elastic interactions between holes, Journal of Materials Science 27 (1992) 1589–1598

    K. Davanas, Analysis of elastic interactions between holes, Journal of Materials Science 27 (1992) 1589–1598

    1992

  7. [15]

    R. S. Rivlin, Large elastic deformations of isotropic materials iv. further developments of the general theory, Philosophical transactions of the royal society of London. Series A, Mathematical and physical sciences 241 (835) (1948) 379–397

    1948

  8. [16]

    Mooney, A theory of large elastic deformation, Journal of applied physics 11 (9) (1940) 582–592

    M. Mooney, A theory of large elastic deformation, Journal of applied physics 11 (9) (1940) 582–592

    1940

  9. [17]

    Arruda, M

    E. Arruda, M. Boyce, A three-dimensional constitutive model for the large stretch behavior of rubber elastic materials, J. Mech. Phys. Solids. 41 (1993) 389–412

    1993

  10. [18]

    D. Seol, S. Hu, Y. Li, J. Shen, K. Oh, L. Chen, Computer simulation of spinodal decomposition in constrained films, Acta materialia 51 (17) (2003) 5173–5185

    2003

  11. [19]

    Binder, S

    K. Binder, S. Puri, S. K. Das, J. Horbach, Phase separation in confined geometries, Journal of Statistical Physics 138 (2010) 51–84

    2010

  12. [20]

    Bohinc, V

    K. Bohinc, V. Kralj-Igli vc, S. May, Interaction between two cylindrical inclusions in a symmetric lipid bilayer, The Journal of Chemical Physics 119 (2003) 7435–7444

    2003

  13. [21]

    Kim, Interaction behaviours between parallel tunnels in soft ground, Tunnelling and Underground Space Technology 19 (4-5) (2004) 448

    S.-H. Kim, Interaction behaviours between parallel tunnels in soft ground, Tunnelling and Underground Space Technology 19 (4-5) (2004) 448

    2004

  14. [22]

    J. Fu, J. Yang, L. Yan, S. M. Abbas, An analytical solution for deforming twin-parallel tunnels in an elastic half plane, International Journal for Numerical and Analytical Methods in Geomechanics 39 (5) (2015) 524–538

    2015

  15. [23]

    G. Zeng, H. Wang, M. Jiang, Analytical stress and displacement of twin noncircular tunnels in elastic semi-infinite ground, Computers and Geotechnics 160 (2023) 105520

    2023

  16. [24]

    L. M. Kanner, C. O. Horgan, Inhomogeneous shearing of strain-stiffening rubber-like hollow circular cylinders, International Journal of Solids and Structures 45 (2008) 5464–5482. 14

    2008

  17. [25]

    Syst` emes, Abaqus 2022, simulia, Dassault Syst` emes, Providence, RI, USA (2022)

    D. Syst` emes, Abaqus 2022, simulia, Dassault Syst` emes, Providence, RI, USA (2022). URL https://www.3ds.com/products-services/simulia/products/abaqus/

    2022

  18. [26]

    C. O. Horgan, G. Saccomandi, Phenomenological hyperelastic strain-stiffening constitutive models for rubber, Rubber chemistry and technology 79 (1) (2006) 152–169

    2006

  19. [27]

    Storm, J

    C. Storm, J. J. Pastore, F. C. MacKintosh, T. C. Lubensky, P. A. Janmey, Nonlinear elasticity in biological gels, Nature 435 (7039) (2005) 191–194

    2005

  20. [28]

    K. A. Erk, K. J. Henderson, K. R. Shull, Strain stiffening in synthetic and biopolymer networks, Biomacromolecules 11 (5) (2010) 1358–1363

    2010

  21. [29]

    M. C. Boyce, E. M. Arruda, Constitutive models of rubber elasticity: a review, Rubber chemistry and technology 73 (3) (2000) 504–523

    2000

  22. [30]

    Suchocki, S

    C. Suchocki, S. Jemio lo, Polyconvex hyperelastic modeling of rubberlike materials, Journal of the Brazilian Society of Mechanical Sciences and Engineering 43 (7) (2021) 352

    2021

  23. [31]

    Lucht, Bipolar coordinates and the two-cylinder capacitor, Rimrock Digital Technology 84103 (2015) 61

    P. Lucht, Bipolar coordinates and the two-cylinder capacitor, Rimrock Digital Technology 84103 (2015) 61. 15

    2015

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.