REVIEW 3 major objections 5 minor 23 references
A Volumetrically Stabilized Mixed Formulation of the Finite Element Immersed Boundary Method for Fluid Structure Interaction with Fully Incompressible Hyperelastic Solids
T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Adding a weakly enforced solid pressure as a Lagrange multiplier for $J=1$ removes the volumetric collapse and locking that afflict the distributed-Lagrange-multiplier finite element immersed boundary method for fully incompressible…
desk verdict A credible mixed DLM-FE-IBM extension with honest reporting, but the volumetric-stability claim is only true in a narrow mesh-ratio window and the abstract overstates it. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing mechanism is the solid pressure $p_s$, introduced as a Lagrange multiplier for the Lagrangian incompressibility constraint $J=1$. It enters through a volumetric energy $U(J)=\frac{1}{2}(\ln J)^2$ with $U'(J)=\ln J/J$, enforced weakly as $\int_B U'(J) q_s \, dX = 0$ for all test functions $q_s$ on the solid mesh; the resulting volumetric stress $p_s J F^{-T}$ is added to the isochoric neo-Hookean stress. The time discretization is the semi-implicit scheme inherited from the original method, in which the solid position in the fluid–solid coupling term is frozen at the previous time level while the elastic stress is evaluated at the current configuration, so that the fluid shape functions need not be re-evaluated at every Newton iteration; the paper relies on the cited unconditional stability of that scheme. The spatial discretization uses inf-sup stable fluid elements and continuous piecewise-polynomial spaces for displacement, solid pressure, and multiplier, with extra quadrature used to integrate the non-matching fluid–solid coupling terms accurately.
What would settle it
Run the radially stretched ring at a solid-to-fluid mesh-size ratio of about 1/2 with a much finer mesh than case E and monitor the enclosed area over time; if the sustained pressure difference and the stretched configuration decay toward the reference configuration even while the weak constraint residual is zero, then the discrete solution is not actually enforcing $J=1$. Alternatively, use the falling-disk case with very large time steps and check whether the kinetic energy grows; that would indicate the inherited unconditional stability does not hold for the mixed formulation.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is that volumetric failure in the immersed setting is not an intrinsic limitation of the immersed boundary method but a missing Lagrange multiplier. The paper derives a mixed formulation in which the first Piola–Kirchhoff stress is replaced by $P = p_s J F^{-T} + \mu_s J^{-2/3}(F - I_1/3 F^{-T})$, with $p_s$ the new solid pressure, and enforces the Lagrangian incompressibility constraint through the weak condition $\int_B U'(J) q_s \, dX = 0$ with $U(J)=\frac{1}{2}(\ln J)^2$. Here $p_s$ is a finite element field living on the solid mesh, not a penalty parameter, and in the fully incompressible limit $U'(J)=0$ is imposed exactly in the weak sense. The paper shows by numerical verification that this restores the expected physics: an elliptically displaced thick ring returns to its circular equilibrium at every refinement level, a radially stretched ring sustains a physical pressure difference across its wall when the solid mesh is fine enough relative to the fluid mesh, and a falling disk matches the empirical terminal velocity to within 1% over a range of densities and viscosities. The authors also report the inflating-ring case as a remaining limitation, where the method either relaxes back (for coarse solid meshes) or locks (for very fine solid meshes).
Load-bearing premise
The load-bearing assumption is that the unconditional stability of the semi-implicit time discretization carries over from the original distributed-Lagrange-multiplier immersed boundary method to the enlarged system with the new solid pressure unknown and the nonlinear $J=1$ constraint; the paper cites the earlier stability proof but gives none for the new system.
Editorial extensions
If this is right
- Fully incompressible hyperelastic solids immersed in viscous fluids can be simulated without remeshing and without spurious volume collapse, as long as the solid mesh is sufficiently fine relative to the fluid mesh.
- Physical pressure jumps across solid walls become computable: the radially stretched ring benchmark shows a sustained pressure difference between interior and exterior when the solid-to-fluid mesh-size ratio is about 1/2.
- Quantitative settling problems are within reach: the falling-disk terminal velocity matches the empirical value to within about 1% for multiple densities and viscosities.
- A practical resolution requirement is identified: non-matching coupling integrals need extra quadrature (order 4–8) before results converge, which any implementation must supply.
- The method's remaining open challenge is the pressure-loaded inflating ring, where the correct physical response is obtained only in a narrow range of mesh-size ratios; outside it the method either fails to hold volume or locks.
Reading between the lines
- The paper does not prove stability of the enlarged mixed system; if the inherited stability proof does not carry over, the unconditional-stability claim could fail for large time steps in regimes where the solid pressure couples strongly to the flow. Testing the scheme on the falling disk with increasingly large time steps would settle this.
- The mesh-ratio sensitivity suggests the pair (solid displacement, solid pressure) may need its own discrete inf-sup condition in the immersed setting; one testable extension is to vary the polynomial order of the solid pressure at fixed mesh ratio and see whether the inflating-ring behaviour improves.
- The volumetric-energy route is general: replacing $U(J)$ with other functions would extend the same mixed framework to different incompressible or nearly incompressible constitutive laws.
- A practical design rule implicit in the results is to keep the solid-to-fluid mesh-size ratio at or below about 1/2 for volume-conserving simulations; this could be verified independently on the radially stretched ring benchmark.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a mixed formulation of the distributed-Lagrange-multiplier finite element immersed boundary method (DLM-FE-IBM) for fluid–structure interaction with fully incompressible hyperelastic solids. It introduces an additional solid pressure field that acts as a Lagrange multiplier for the Lagrangian incompressibility constraint J=1, and augments the solid stress with a volumetric contribution derived from a dilatational strain energy. The formulation is discretized with a semi-implicit time scheme and finite elements, implemented in the GRINS multiphysics framework, and verified on three benchmarks: an elliptically displaced thick ring, a radially stretched ring, and a disk falling under gravity. The paper reports that the mixed formulation removes the volumetric failure of the unstabilized method on the first two benchmarks and reproduces the empirical terminal velocity of the falling disk to within about 1%, but it also reports an inflating-ring case that exposes a remaining sensitivity to the solid-to-fluid mesh-size ratio, which the authors describe as an open challenge.
Significance. If the central claim holds, the method would be a useful extension of the DLM-FE-IBM to fully incompressible hyperelastic solids, which are prone to volumetric locking and spurious volume collapse. The paper gives a clean derivation of the mixed formulation, describes a reusable implementation in an established multiphysics framework, and provides quantitative verification against an independent empirical correlation for the falling disk. The authors also deserve credit for transparently reporting the inflating-ring failure. However, the central claim is stronger than the evidence: the numerical tests show that volumetric stability is restored only in certain mesh-ratio regimes, and no discrete inf-sup analysis or stability proof is given for the enlarged mixed system. These gaps are load-bearing because they directly concern the method's claimed generality and reliability.
major comments (3)
- [Section 2.5] The discrete stability of the solid mixed pair (U_h, S_h) is not addressed. The text specifies only that the solid displacement, solid pressure, and multiplier are discretized with continuous piecewise-polynomial elements, without stating the element pair or verifying a discrete inf-sup condition for the mixed system (20c)–(20d). For incompressible mixed elasticity, the displacement and pressure spaces must satisfy an inf-sup condition; arbitrary equal-order pairs are known to lock or oscillate. Since the central claim is that the mixed formulation restores volumetric stability, the absence of this analysis leaves the claim unestablished at the discrete level.
- [Abstract and Section 2.4] The assertion that the semi-implicit time discretization is unconditionally stable is not supported for the enlarged system (19). The cited stability proof in [7] applies to the original DLM-FE-IBM, which does not include the solid pressure unknown p_s or the nonlinear constraint (19d). No stability analysis is given for the new system, so the claim should either be proved or qualified as an inherited property that has not been verified for the mixed formulation.
- [Sections 4.2 and 4.4] The numerical results show that volumetric stability is not restored in general, contradicting the unqualified claim in the abstract. In Section 4.2, cases A and B (hs/hx ~ 2 and ~ 1) allow the radially stretched ring to relax and the pressure difference to decay, indicating that J=1 is not effectively enforced; only case C (hs/hx ~ 1/2) sustains the pressure difference, and even there residual pressure leakage remains. In Section 4.4, the inflating ring at hs/hx ~ 1 retains residual volumetric instabilities, and reducing to hs/hx ~ 1/2 produces volumetric locking. The abstract should qualify the claim that the mixed formulation restores volumetric stability to the tested regimes, and the paper should discuss the conditions under which the stabilization is effective.
minor comments (5)
- [Section 2.5] Please specify the actual finite element spaces used in the computations (e.g., the polynomial degree of the continuous piecewise-polynomial elements for U_h, S_h, and Lambda_h), rather than only stating that they are continuous piecewise-polynomial.
- [Eq. (17) and (18d)] The constraint (18d) is written as the integral of U'(J) against q_s; since U'(J)=ln J / J, its vanishing is equivalent to J=1 only when J is positive. This equivalence should be stated explicitly.
- [Section 4.3] The empirical formula (21) is for a rigid disk; the paper should justify its applicability to an incompressible neo-Hookean disk, for example by reporting the maximum deformation or by showing that the solid modulus makes the disk effectively rigid under the tested conditions.
- [Title page] There is a spacing issue in the title: 'A V olumetrically Stabilized' should read 'A Volumetrically Stabilized'.
- [References] Reference [16] is missing volume and page numbers; the reference entry should be completed.
Circularity Check
No significant circularity: the mixed formulation is self-contained and its verification uses external, parameter-free benchmarks.
full rationale
The paper's central claim is that adding a solid pressure unknown and the weak constraint (18d) enforces the Lagrangian incompressibility constraint J=1 in the DLM-FE-IBM. This is a formulation construction, not a prediction obtained from fitted inputs: the volumetric energy U(J)=1/2(ln J)^2 is adopted from [13] as a modeling choice, and the solid pressure is defined by p_s=κU'(J) and then imposed weakly, which is the standard mixed treatment of incompressible elasticity rather than a circular reuse of the target result. Material parameters (μ_s, ρ_f, μ_f) are prescribed for the benchmarks, not fitted to the outputs. The falling-disk terminal velocity is compared with the external empirical formula (21) from [22], and the ~1% agreement is an independent validation. No step in the derivation reduces to its own output by construction. The only caveat is that the unconditional-stability assertion in Sections 2.4 and 5 is inherited from Boffi et al. [7] for the original DLM-FE-IBM without the new solid-pressure unknown and constraint (19d); that is an unverified assumption and a correctness risk, and the paper itself flags the mesh-ratio sensitivity in the inflating-ring test (Section 4.4), but neither constitutes circularity. Citations to the authors' own GRINS and libMesh work are implementation references and do not carry the mathematical claim, so there is no load-bearing self-citation.
Assumptions & free parameters
assumptions (5)
- domain assumption The Cauchy stress additively splits into fluid and solid parts, sigma = sigma_f in the fluid region and sigma_f + sigma_s in the solid region.
- domain assumption The solid strain energy splits into volumetric and isochoric parts, W = kappa U(J) + Wbar(C).
- domain assumption The volumetric energy U(J) = 1/2 (ln J)^2 is chosen.
- domain assumption The weak enforcement of p_s = kappa U'(J) via integral_B (U'(J) - p_s/kappa) q_s = 0 gives the correct incompressible limit as kappa goes to infinity.
- ad hoc to paper The semi-implicit time discretization is unconditionally stable for the mixed formulation.
Cite this review
Pith. "Pith review of A Volumetrically Stabilized Mixed Formulation of the Finite Element Immersed Boundary Method for Fluid Structure Interaction with Fully Incompressible Hyperelastic Solids." pith.science (2026). https://pith.science/paper/LOEQLXNE
@misc{pith2026260810217,
author = {Pith},
title = {Pith review of: A Volumetrically Stabilized Mixed Formulation of the Finite Element Immersed Boundary Method for Fluid Structure Interaction with Fully Incompressible Hyperelastic Solids},
year = {2026},
howpublished = {\url{https://pith.science/paper/LOEQLXNE}},
note = {Machine review of arXiv:2608.10217}
}
read the original abstract
The finite element immersed boundary method (FE IBM) with a distributed Lagrange multiplier is an attractive framework for fluid structure interaction (FSI) because it couples an Eulerian description of an incompressible fluid to a Lagrangian description of an immersed solid on two independent, nonconforming meshes, avoiding the costly remeshing required by body fitted arbitrary Lagrangian Eulerian methods. When the immersed solid is modelled as a fully incompressible hyperelastic material, however, a direct finite element discretization of the deviatoric solid stress fails to enforce the Lagrangian incompressibility constraint, and the computed structure exhibits spurious volumetric instabilities and locking. In this work we present a mixed formulation of the distributed Lagrange multiplier FEIBM that restores volumetric stability. Following the theory of nearly incompressible hyperelasticity, we augment the solid stress with a volumetric contribution derived from a dilatational strain energy and introduce an additional solid pressure field, enforced weakly, that plays the role of the Lagrange multiplier for the Lagrangian incompressibility constraint J = 1. The resulting formulation is discretized in space by finite elements and in time by an unconditionally stable semi implicit scheme, and is implemented as a reusable Immersed Boundary physics kernel within the GRINS multiphysics framework, built on the libMesh finite element library. The method is verified on three FSI benchmarks, an elliptically displaced thick ring returning to equilibrium, a radially stretched incompressible ring, and a disk falling under gravity in a viscous fluid, for which the mixed formulation removes the volumetric failure observed with the unstabilized formulation and reproduces the analytical terminal velocity of the falling disk to within 1%.
Figures
Figures from the paper (13 more)
Reference graph
Works this paper leans on
- [7]
-
[1]
T. J. R. Hughes, W. K. Liu, T. K. Zimmermann, Lagrangian–Eulerian finite element formulation for incompressible viscous flows, Computer Methods in Applied Mechanics and Engineering 29 (3) (1981) 329–349
work page 1981
- [2]
-
[3]
C. S. Peskin, The immersed boundary method, Acta Numerica 11 (2002) 479–517
2002
-
[4]
Heltai, The finite element immersed boundary method, Ph.D
L. Heltai, The finite element immersed boundary method, Ph.D. thesis, Universit `a di Pavia, Dipartimento di Matematica (2006)
work page 2006
- [5]
- [6]
-
[8]
V . Girault, R. Glowinski, T. W. Pan, A fictitious domain method with distributed Lagrange multipliers for the Stokes problem, Applied Nonlinear Analysis (Kluwer/Plenum) (1999) 159–174
work page 1999
Show all 23 references
-
[9]
Glowinski, Y
R. Glowinski, Y . Kuznetsov, Distributed Lagrange multipliers based on fictitious domain method for second order elliptic problems, Computer Methods in Applied Mechanics and Engineering 196 (8) (2007) 1498–1506
2007
-
[10]
Alshehri, D
N. Alshehri, D. Boffi, F. Credali, L. Gastaldi, Advances on finite element discretization of fluid–structure interaction problems, Arabian Journal of Mathematics 15 (1) (2025) 27–51
2025
-
[11]
P. T. Bauman, R. H. Stogner, GRINS: A multiphysics framework based on the libMesh finite element library, SIAM Journal on Scientific Computing 38 (5) (2016) S78–S100
2016
-
[12]
B. S. Kirk, J. W. Peterson, R. H. Stogner, G. F. Carey, libMesh: A C++library for parallel adaptive mesh refinement/coarsening simulations, Engineering with Computers 22 (3–4) (2006) 237–254
2006
-
[13]
Vadala-Roth, S
B. Vadala-Roth, S. Rossi, S. Prakash, N. A. Barrett, B. E. Griffith, Stabilization approaches for the hyperelastic immersed boundary method for problems of large-deformation incompressible elasticity, Computer Methods in Applied Mechanics and Engineering 365 (2020) 112978
2020
-
[14]
L. Li, C. Gruninger, J. H. Lee, B. E. Griffith, Local divergence-free immersed finite element–difference method using composite B-splines, arXiv preprint arXiv:2412.15408
-
[15]
Brink, E
U. Brink, E. Stein, On some mixed finite element methods for incompressible and nearly incompressible finite elasticity, Computational Mechanics 19 (1) (1996) 105–119
1996
-
[16]
Karabelas, G
E. Karabelas, G. Haase, G. Plank, C. M. Augustin, Versatile stabilized finite element formulations for nearly and fully incompressible solid mechanics, Computational Mechanics 65 (1)
-
[17]
M. D. Gunzburger, Finite Element Methods for Viscous Incompressible Flows: A Guide to Theory, Practice, and Algorithms, Academic Press, 1989
1989
-
[18]
Gonzalez, A
O. Gonzalez, A. M. Stuart, A First Course in Continuum Mechanics, Cambridge University Press, 2008
2008
-
[19]
Heltai, F
L. Heltai, F. Costanzo, Variational implementation of immersed finite element methods, Computer Methods in Applied Mechanics and Engineering 229–232 (2012) 110–127
2012
-
[20]
S. Roy, L. Heltai, F. Costanzo, Benchmarking the immersed finite element method for fluid–structure interaction problems, Computers & Mathematics with Applications 69 (10) (2015) 1167–1188
2015
-
[21]
L. T. Zhang, M. Gay, Immersed finite element method for fluid–structure interactions, Journal of Fluids and Structures 23 (6) (2007) 839–857
2007
-
[22]
Clift, J
R. Clift, J. R. Grace, M. E. Weber, Bubbles, Drops, and Particles, Academic Press, New York, 1978
1978
-
[23]
Turek, J
S. Turek, J. Hron, Proposal for numerical benchmarking of fluid–structure interaction between an elastic object and laminar incompressible flow, in: Fluid–Structure Interaction, V ol. 53 of Lecture Notes in Computational Science and Engineering, Springer, 2006, pp. 371–385. 18
2006
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.