{"id":"51c7a5bf-fe13-4b72-856e-b92e1daccace","arxiv_id":"2608.10217","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A volumetrically stabilized mixed finite element immersed boundary formulation that enforces Lagrangian incompressibility in fully incompressible hyperelastic solids via a weak solid pressure field.","lead":"The paper introduces a new mixed formulation for the finite element immersed boundary method that adds a solid pressure field to enforce volume conservation in incompressible hyperelastic solids. It is a numerical methods paper that could improve simulations of heart valves, parachutes, and soft materials in fluids.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central 'restores volumetric stability' claim is not established as stated: no inf-sup analysis for the solid (U_h, S_h) pair, and the paper's own radial-ring and inflating-ring tests show J=1 enforcement depends on hs/hx and can leak or lock.","rationale":"The paper's central contribution is a mixed solid pressure intended to make the discrete solution respect J=1. That requires a stable mixed finite element pair, but the paper provides no inf-sup analysis and leaves the element choice vague. More tellingly, the numerical results in Section 4.2 and Section 4.4 are direct evidence that the constraint enforcement is fragile: it depends on hs/hx and can either leak (radially stretched ring at large ratio) or lock (inflating ring at small ratio). Because the authors openly disclose these cases, this is not a charge of concealment; it is an internal inconsistency between the unqualified abstract statement and their own data. The reader's concern about unconditional stability of the time scheme is legitimate but secondary: even a stable time discretization would not rescue the central claim if the discrete volumetric constraint is not robustly satisfied. The method is a credible and useful contribution for the demonstrated parameter window, and the three benchmarks provide real supporting evidence for that narrower claim. Keeping the conditional verdict is appropriate: the paper should be accepted only with the caveat that the broad volumetric-stability claim requires a discrete stability analysis and a precise characterization of the admissible hs/hx regime.","tokens_in":13286,"tokens_out":7323,"duration_ms":79481,"concrete_test":"Repeat the radially stretched ring case C (hs/hx ~ 1/2) and the inflating-ring cases II and the hs/hx ~ 1/2 locking case with a solid pressure space one order lower than the displacement space (e.g., P2 displacement / P1 pressure, or P1/P0), keeping all other parameters identical; in parallel, compute the discrete inf-sup constant of the current (U_h, S_h) pair on the solid mesh. If the sustained pressure difference appears and the inflating-ring locking disappears with the certified stable pair, the central claim was limited by the space choice; if the failures persist, the formulation itself is the limitation and the abstract's claim must be explicitly restricted to the demonstrated regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires the discrete constraint (20d) to enforce J=1 without locking across the intended regime. Section 2.5 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 displacement–pressure 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. The paper's own results indicate the constraint is not robustly satisfied: 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, meaning J=1 is not effectively enforced; only case C (hs/hx ~ 1/2) sustains the pressure difference, and even there 'residual pressure leakage remains.' Section 4.4 reports that the inflating ring at hs/hx ~ 1 retains 'residual volumetric instabilities,' and reducing to hs/hx ~ 1/2 'produced volumetric locking.' Thus the abstract's claim that the mixed formulation 'restores volumetric stability' is supported only in a narrow mesh-ratio window, and the mechanism that would make it general — a stable mixed finite element pair — is not demonstrated. The unconditional-stability assertion inherited from [7] for the enlarged system (19) is also unproven, but the discrete stability of the new (U_h, S_h) pair is the more immediate threat to the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":13584,"tokens_out":5594,"duration_ms":54303,"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":[{"comment":"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.","section":"Section 2.5"},{"comment":"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.","section":"Abstract and Section 2.4"},{"comment":"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.","section":"Sections 4.2 and 4.4"}],"minor_comments":[{"comment":"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.","section":"Section 2.5"},{"comment":"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":"Eq. (17) and (18d)"},{"comment":"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.","section":"Section 4.3"},{"comment":"There is a spacing issue in the title: 'A V olumetrically Stabilized' should read 'A Volumetrically Stabilized'.","section":"Title page"},{"comment":"Reference [16] is missing volume and page numbers; the reference entry should be completed.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope, but the abstract overstates the central claim. The most serious gap is the missing discrete inf-sup analysis for the solid pressure pair and the unverified unconditional-stability assertion; these should be addressed or the claims qualified before publication. The authors' honest reporting of the inflating-ring limitation is a positive feature, but it underscores that the method's robustness is not yet established."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a solid incremental contribution: it grafts a mixed solid-pressure formulation onto the distributed-Lagrange-multiplier FE-IBM, making the volumetric constraint a genuine mixed unknown in a fully variational setting. That combination is new, and the derivation is clean and easy to follow. The implementation in GRINS is reusable, and the three benchmarks are appropriate. The falling-disk result matching an empirical terminal velocity to within 1% is good practical evidence, and I credit the authors for reporting the inflating-ring failure and the mesh-ratio sensitivity instead of hiding them.\n\nThe main problem is the abstract's claim that the method \"restores volumetric stability.\" The paper's own results show this is only true for sufficiently small solid-to-fluid mesh ratio (hs/hx ~ 1/2), and even then residual pressure leakage remains. At coarser ratios the constraint leaks; at finer ratios the inflating ring locks. So the method stabilizes the volumetric response in a limited window, not across the intended regime. The stress-test note is right that no discrete inf-sup analysis is given and the element pair for (U_h, S_h) is never specified. That is a real gap for a method whose whole purpose is to avoid locking. The paper should either provide a stability analysis, run a numerical inf-sup test, or explicitly frame the method as requiring a careful mesh-ratio choice. The unconditional-stability claim inherited from Boffi et al. also deserves at least a remark that the enlarged system is not covered by that proof; it is likely fine in practice but is still an unverified assumption.\n\nMinor but worth fixing: the abstract calls the terminal velocity \"analytical,\" but Eq. (21) is an empirical correlation. That distinction matters.\n\nOverall, this is a useful methods paper for researchers using DLM-FE-IBM for incompressible solids, especially in the GRINS ecosystem. It deserves a serious referee and, after revisions that qualify the central claim and discuss discrete stability, acceptance. I would not cite it in my own immediate work, but I would bring it to a reading group focused on FSI or mixed methods.","headline":"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.","tokens_in":14133,"tokens_out":2551,"would_cite":false,"duration_ms":28608,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M60","65M12","74S05","76M10","74B20"],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["immersed boundary method","finite element method","fluid-structure interaction","distributed Lagrange multiplier","volumetric stabilization","incompressible hyperelasticity","mixed finite element formulation","fictitious domain"],"falsifier":"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.","tokens_in":13069,"feed_emoji":"🧮","tokens_out":9503,"duration_ms":81560,"temperature":0.7,"pith_summary":"The paper is trying to establish that the distributed-Lagrange-multiplier finite element immersed boundary method can be made volumetrically stable for fully incompressible hyperelastic solids by treating a new solid pressure field as a genuine mixed unknown. The failure it addresses is concrete: when only the deviatoric solid stress is discretized, the Lagrangian volume constraint $J=1$ is satisfied only approximately, and immersed rings and other structures collapse or lock. The proposed cure follows the classical mixed treatment of nearly incompressible hyperelasticity: augment the solid stress with a volumetric term derived from a dilatational energy, introduce a solid pressure that weakly enforces $J=1$, and discretize the enlarged system with finite elements. On three benchmarks—an elliptically displaced ring, a radially stretched ring, and a disk falling under gravity—the mixed formulation removes the collapse and reproduces the falling disk's terminal velocity to within about 1%. A fourth, pressure-loaded inflating-ring case is reported as an open challenge because it remains sensitive to the solid-to-fluid mesh-size ratio.","feed_headline":"Solid pressure field stops immersed elastic solids from collapsing","feed_subtitle":"Weakly enforcing J=1 restores volume conservation and matches the falling-disk terminal velocity within 1 percent.","key_machinery":"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.","core_discovery":"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).","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the base distributed-Lagrange-multiplier finite element immersed boundary method formulation that this paper extends, and the semi-implicit scheme whose unconditional stability is inherited.","marker":"[7]"},{"why":"Provides the hyperelastic finite element immersed boundary method variational formulation and stress/density split that the mixed formulation builds on.","marker":"[2]"},{"why":"Introduces the volumetric stabilization idea (dilatational energy enforcing $J=1$) for the immersed finite element/finite difference method, which the paper adapts into a genuine mixed field.","marker":"[13]"},{"why":"Supplies the mixed displacement–pressure theory of nearly incompressible finite elasticity used to justify the solid pressure as a Lagrange multiplier.","marker":"[15]"},{"why":"Provides the inf-sup stable fluid element theory and the bilinear and trilinear forms used in the discretization.","marker":"[17]"},{"why":"Presents the falling-disk immersed finite element method benchmark whose setup the paper follows for the terminal-velocity test.","marker":"[20]"},{"why":"Gives the empirical terminal-velocity formula used as the quantitative reference for the falling disk.","marker":"[22]"},{"why":"Describes divergence-free composite B-spline kernels, the alternative approach the paper identifies for future work on the open mesh-ratio challenge.","marker":"[14]"}],"fun_headline_variants":["Mixed formulation halts spurious collapse in immersed solids","Weak solid pressure enforces incompressibility in FSI","Immersed method gains pressure field to stop volumetric locking","Restoring volume conservation in immersed hyperelastic FSI","Pressure Lagrange multiplier stabilizes immersed elastic solids"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Mixed formulation halts spurious collapse in immersed solids","Weak solid pressure enforces incompressibility in FSI","Immersed method gains pressure field to stop volumetric locking","Restoring volume conservation in immersed hyperelastic FSI","Pressure Lagrange multiplier stabilizes immersed elastic solids"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000439,"raw_usage":{"total_tokens":2335,"prompt_tokens":1161,"completion_tokens":1174,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":777,"completion_tokens_details":{"reasoning_tokens":1098}},"tokens_in":777,"tokens_out":1174,"duration_ms":9499,"temperature":1.0,"reasoning_tokens":1098,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:09:58.830923+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Boffi, N","cited_arxiv_id":null,"evidence_quote":"Supplies the base distributed-Lagrange-multiplier finite element immersed boundary method formulation that this paper extends, and the semi-implicit scheme whose unconditional stability is inherited."},{"cited_title":"Boffi, L","cited_arxiv_id":null,"evidence_quote":"Provides the hyperelastic finite element immersed boundary method variational formulation and stress/density split that the mixed formulation builds on."},{"cited_title":"Vadala-Roth, S","cited_arxiv_id":null,"evidence_quote":"Introduces the volumetric stabilization idea (dilatational energy enforcing $J=1$) for the immersed finite element/finite difference method, which the paper adapts into a genuine mixed field."},{"cited_title":"Brink, E","cited_arxiv_id":null,"evidence_quote":"Supplies the mixed displacement–pressure theory of nearly incompressible finite elasticity used to justify the solid pressure as a Lagrange multiplier."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the inf-sup stable fluid element theory and the bilinear and trilinear forms used in the discretization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Presents the falling-disk immersed finite element method benchmark whose setup the paper follows for the terminal-velocity test."},{"cited_title":"Clift, J","cited_arxiv_id":null,"evidence_quote":"Gives the empirical terminal-velocity formula used as the quantitative reference for the falling disk."},{"cited_title":"Local Divergence-Free Immersed Finite Element-Difference Method Using Composite B-Splines","cited_arxiv_id":"2412.15408","evidence_quote":"Describes divergence-free composite B-spline kernels, the alternative approach the paper identifies for future work on the open mesh-ratio challenge."}],"review_version":1}