Pith. sign in

REVIEW 2 major objections 4 minor 30 references

An interaction potential method for passive and active dynamics of hyperelastic materials

T0 review · 2 major / 4 minor · reviewed 2026-08-01 · deepseek-v4-flash

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

Pith's one-line read This paper claims that the discrete energy of an interaction-potential tetrahedron can be made identical to the continuum hyperelastic strain-energy density, by using the six edge strains as exact coordinates for the Green-Lagrange strain t

desk verdict Solid, honest methods paper — the edge-strain construction is exact and the benchmarks are real, but the nodal-volume incompressibility trick is unproven and the cardiac claims lean on it. read the letter →

arxiv 2607.19227 v1 pith:3ZFWAOOJ submitted 2026-07-21 physics.comp-ph

Francesco De Vita , Filippo Caruso Lombardi , Roberto Verzicco , Francesco Viola This is my paper · ORCID
classification physics.comp-ph MSC 74B2074S0565M60
keywords interactionpotentialmethodhyperelasticityactivestraincardiacmechanicstetrahedralmeshesincompressibilityedgeGPUsimulation
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

This paper proposes an interaction-potential method for hyperelastic materials in which the discrete strain energy of each tetrahedron is exactly the continuum strain-energy density U(F), not a spring-network surrogate. The key step is a change of variables: the six edge strains of a tetrahedron map linearly and invertibly to the six independent components of the Green-Lagrange strain tensor, so a constitutive law written in terms of E can be re-expressed exactly in terms of edge strains. The payoff is that nodal forces remain local and edge-wise, ideal for explicit, GPU-style solvers, while preserving the energetic coupling between edges that classical mass-spring models lose. Active contraction is included through the multiplicative decomposition F = F_e F_a, realized as a time-dependent activated reference configuration, and near-incompressibility is handled by a nodal-volume penalty intended to avoid volumetric locking. If correct, the method gives particle-based solvers the fidelity of constant-strain tetrahedral finite elements, validated on standard cardiac-mechanics benchmarks and a realistic ventricle simulation.

What carries the argument

The central object is the edge Green-Lagrange strain ε_jk = 1/2(l^2/l_0^2 - 1). For each tetrahedron, the six edge strains assemble into a vector ε that is a linear function of the six independent components of the continuum Green-Lagrange tensor E, via a 6×6 matrix D built from the reference edge unit vectors. As long as the tetrahedron is non-degenerate, D is invertible, so U(E) can be rewritten exactly as U(D^{-1}ε). This change of variables is what converts a continuum constitutive law into local edge-wise forces while keeping all inter-edge couplings; material-specific derivatives reduce to ∂U/∂ε times a purely geometric factor depending on which node the force acts on.

What would settle it

Construct a non-uniform tetrahedral mesh of a nearly incompressible block, apply large shear or torsion while increasing k_v, and measure the distribution of per-element determinants J_c - 1. If element-level volume errors do not vanish, or locking reappears, even as the nodal control volumes remain near their reference values, then the nodal-volume penalty is not enforcing incompressibility as claimed.

Watch

Extended reading notes

Core claim

The paper's central claim is that the discrete strain energy of a tetrahedron can be made identical to the continuum hyperelastic strain-energy density U(F), not merely an approximation, by changing variables from the six components of the Green-Lagrange strain tensor to the six edge stretches of the tetrahedron. Because a linear tetrahedron's deformation gradient is exactly determined by its edge vectors, the mapping E_v = D^{-1} ε is exact for non-degenerate elements. Consequently, nodal forces computed from edge strains reproduce the internal force of a constant-strain tetrahedral finite element, with the same energetic coupling among edges that mass-spring models lose. Active contraction

Load-bearing premise

The method's near-incompressible tissue behavior depends entirely on the idea that constraining the volume of the lumped region around each node, rather than the local stretch of each tetrahedron, is enough to make the material incompressible without stiff locking; that assumption is supported only by one beam test in the paper.

Editorial extensions

If this is right

  • If the equivalence holds, the method offers a route from arbitrary hyperelastic constitutive laws to particle-style explicit solvers without fitting spring parameters.
  • The active-strain decomposition becomes a per-step update of the reference edge lengths, making active contraction as cheap as a reference-configuration change, with active power entering the energy balance.
  • For near-incompressible tissue, the nodal-volume penalty is claimed to enforce incompressibility without volumetric locking, at first-order convergence in the bulk modulus.
  • Benchmarks on beam, ventricle inflation, cube compression, and active cube place the method's accuracy on par with finite element reference solutions for those test cases.
  • The realistic ventricle simulation produces longitudinal shortening around 20%, wall thickening around 32%, and ejection fraction about 45%, consistent with physiology.

Reading between the lines

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

  • The exactness hinges on linear (P1) tetrahedra; the same change of variables does not carry over to higher-order elements, so the method's route is limited unless a generalized mapping is found.
  • The nodal-volume incompressibility constraint is a lumped constraint; a natural test is whether per-element volume drift remains bounded on highly distorted or non-uniform meshes, where lumped control volumes may misrepresent local Jacobians.
  • Because forces are local and the mass matrix diagonal, the method should extend naturally to explicit coupling with electrophysiology or fluid solvers in a heart model, but such coupling is not demonstrated in this paper.
  • The energy-balance check suggests a diagnostic for any interaction-potential method: whether dE/dt balances active power P_a; the paper shows second-order convergence in that residual.
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 / 4 minor

Summary. The paper proposes an interaction-potential formulation of hyperelasticity for tetrahedral meshes. The key idea is to relate the six edge Green-Lagrange strains of each tetrahedron to the Voigt strain vector through ε = D E_v; since D is invertible for nondegenerate tetrahedra, the continuum strain-energy density U(E) can be evaluated exactly as a function of edge strains. Nodal forces are then expressed as local, edge-wise contributions, making the method suitable for explicit time integration with a diagonal mass matrix. Active contraction is incorporated through the multiplicative decomposition F = F_e F_a, which becomes a time-dependent active reference configuration and requires recomputing D^{-1} each step. To model near-incompressible tissue, the authors add a penalty on lumped nodal control volumes (Eq. 44) rather than per-element Jacobians. The method is validated on beam deflection, ventricle inflation, cube compression, an active cube, an energy-balance test, and a realistic ventricle simulation.

Significance. The central derivation is clean: the edge-strain identity gives an exact energetic equivalence between the discrete model and a constant-strain tetrahedral FEM element with one-point quadrature, and the active-strain treatment is elegant. The open-source GPU implementation is a concrete strength, and the paper makes a useful bridge between continuum constitutive laws and interaction-potential solvers. However, the near-incompressibility strategy is heuristic and, as it stands, not supported by element-level evidence; because the cardiac applications and the 'no volumetric locking' claim depend on it, the significance of those results is currently limited. The active-cube validation is also weakened by the hand-chosen parameter k0.

major comments (2)
  1. [§3.4, Eq. (44)] The nodal-volume penalty is not equivalent to continuum incompressibility. It constrains only Nv lumped volumes, whereas J_c is piecewise constant on ~Nc ≈ 6Nv tetrahedra; the constraint manifold therefore admits displacement fields that leave every V_i unchanged while changing individual J_c (checkerboard volume-transfer modes). Such modes carry zero volumetric energy for any k_v. The validation in §5.1 reports only a global |V−V0|/V0 error, and §5.2 reports total ventricle volume change (~10^-4); no element-wise J statistics are reported anywhere. The claim in §3.4 that the method 'avoids volumetric locking' and the cardiac applications that depend on near-incompressibility are therefore not yet supported. Please report the distribution of J_c over elements for the cube and ventricle cases and test a perturbation that excites checkerboard modes, or replace the constraint with one that
  2. [§5.4] The active-cube validation uses λ_n = k0 λ_f with k0 = 4 chosen by hand, without independent calibration. With one free parameter, matching the wall-thickening time history from [27] is not strong evidence of predictive accuracy. The paper should provide a sensitivity analysis, determine k0 from material parameters, or compare additional output quantities (e.g., deformation gradients or stress measures) against the reference solution.
minor comments (4)
  1. [§1] Typo: 'both for the passive ad for the active dynamics' should be 'both for the passive and for the active dynamics'.
  2. [Eq. (5) and elsewhere] 'V oigt' should be 'Voigt'.
  3. [§3.2, Eq. (32)] Clarify that D_a is built from the active unit edge vectors x_a/l_a, not from F_a itself; this is implicit but easy to misunderstand.
  4. [§5.1] The bottom-right panel of Fig. 2 is labeled 'error of the incompressibility constraint' but the definition of this error is not given. State explicitly whether it is the total-volume error or an element-wise quantity; this is directly relevant to the main concern about Eq. (44).

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the edge-strain reconstruction is an exact change of variables and the validation benchmarks are external.

full rationale

The central derivation chain is self-contained. Equation (8) defines the six edge Green-Lagrange strains as ε = D E_v for a tetrahedron, and Eq. (9) inverts the 6×6 matrix D for a non-degenerate element, giving E_v = D^{-1}ε. Substitution into the continuum strain-energy density, Eq. (10), is an exact coordinate transformation, not a fitted surrogate; the paper explicitly states that the discrete energy is 'the continuum strain-energy density U(F) itself, evaluated on a discrete kinematics defined by the mesh edges.' The active-strain extension via F = F_e F_a likewise yields ε_e = D_a E_{e,v} and U(E_{e,v}) = U(D_a^{-1}ε_e), again by direct substitution rather than by construction of the target output. The mechanical energy balance Eq. (41) is an algebraic identity: P_a is defined as the portion of dU/dt arising from reference-configuration changes, so the numerical energy-balance test checks consistency of the time integrator, not a predicted law. Validation is performed against external benchmarks (Land et al. [19] and Rossi et al. [27]), none of which involve the present authors or are subsumed by the paper's own equations. The stated limitations—linear P1 tetrahedra and explicit time-integration stability—are acknowledged in the conclusion and do not create circularity. Two non-circular caveats remain: the nodal-volume incompressibility penalty (Eq. 44) is a heuristic whose pointwise/elementwise J behavior is not proven, and the active-cube/ventricle parameter k0 is hand-set without derivation; these limit the strength of the incompressibility and cardiac validation, but they do not make the passive or active constitutive derivation reduce to its inputs by construction.

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

The method introduces no invented physical entities. The central claim rests on standard continuum mechanics (hyperelasticity, multiplicative active strain) plus one ad hoc numerical assumption (nodal-volume penalty) and two chosen numerical parameters (k_v, k0).

free parameters (2)
  • k0 (active cross-fiber stretch coefficient) = 4 (cube), 5 (ventricle)
    Prescribed active deformation uses lambda_n = k0 lambda_f; k0 is chosen by hand and affects the active-cube wall-thickening comparison and ventricle metrics.
  • k_v (volumetric penalty stiffness) = 10 MPa (ventricle inflation), 10^5 kPa (cube), varied 10^0-10^4 kPa (beam)
    Penalty parameter for nodal incompressibility; results are shown to converge as k_v increases, so it is a numerical parameter rather than a fitted physical constant.
assumptions (3)
  • domain assumption The tetrahedral mesh remains non-degenerate so D (and D_a) is invertible.
    Equation (9) requires inversion of the 6x6 edge-strain transformation; sliver or degenerate cells would make the change of variables singular.
  • domain assumption Active strain enters multiplicatively at the edge level as lambda = lambda_e lambda_a, equivalent to F = F_e F_a uniform inside each tetrahedron.
    Section 3.2; requires F_a to be element-wise constant so a single activated reference configuration exists per cell.
  • ad hoc to paper The nodal control-volume energy (44) enforces material incompressibility and avoids volumetric locking.
    Section 3.4; no proof of equivalence to J=1 or locking-free convergence is given; only numerical evidence in Section 5.1.

how reviews work

0 comments
Cite this review

Pith. "Pith review of An interaction potential method for passive and active dynamics of hyperelastic materials." pith.science (2026). https://pith.science/paper/3ZFWAOOJ

@misc{pith2026260719227,
  author       = {Pith},
  title        = {Pith review of: An interaction potential method for passive and active dynamics of hyperelastic materials},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3ZFWAOOJ}},
  note         = {Machine review of arXiv:2607.19227}
}
read the original abstract

Simulating active biological tissues, such as the myocardium, requires constitutive models that are both physically faithful and computationally efficient. The most common approach relies on finite element methods that accurately discretize the underlying continuum hyperelastic problem which, in turn, require global nonlinear solve at each time step. On the other hand, fast, interaction-potential methods replace the continuum with a network of independent links approximating the mechanical response. We propose an interaction potential formulation for simulating active biological tissues that bridges this gap. The method recasts continuum hyperelastic constitutive laws in terms of tetrahedral edge strain. Unlike classical mass-spring models, the proposed formulation does not approximate the tissue as independent spring elements but preserves the energetic coupling between adjacent edges. Passive tissue mechanics is described by hyperelastic constitutive laws, while active contraction is incorporated through the active-strain multiplicative decomposition. Within the edge-based formulation, the active strain is incorporated through a time-dependent activated reference configuration. We further introduce a strategy for enforcing the material incompressibility constraint while avoiding volumetric locking. The resulting method can be interpreted as an edge-strain representation of a constant-strain tetrahedral continuum element, providing a bridge between continuum mechanics and discrete interaction potential solvers. Numerical tests demonstrate the capability of the method to effectively simulate different hyperelastic constitutive laws and to properly preserve the energy balance equation. Finally, the proposed method is applied to the active deformation of a realistic ventricle, reproducing longitudinal shortening and wall-thickening values consistent with the literature.

Figures

Figures reproduced from arXiv: 2607.19227 by the authors.

Figure 1
Figure 1. Sketch of the problem: a three-dimensional deformable soft body (such as the human heart) is discretized [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Results for the beam deflection test case: (top panel) shape of the deformed fiber initially at [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Results for the ventricle inflation test case: the left half of the plot report the deformation of the sampled [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: (left column) Cube configuration at maximum contraction with colorbar showing the displacement field in [PITH_FULL_IMAGE:figures/full_fig_p012_4.png]
Figure 5
Figure 5. Figure 5: (upper left) Time history of kinetic (K/E0), potential (U/E0) and mechanical (E/E0) energy, normalized with the initial mechanical energy (E0 = E(t = 0), with E = K +U) for the passive dynamics test case. (upper right) Time history of kinetic (K), potential (U) and mec…
Figure 6
Figure 6. Figure 6: Realistic ventricle test case: (left panel) initial stress-free configuration; (right panel) maximum contraction [PITH_FULL_IMAGE:figures/full_fig_p014_6.png]
Figure 7
Figure 7. Figure 7: (left panel) Time history of the longitudinal shortening (LS) and the wall thickening (WT); (right panel) time [PITH_FULL_IMAGE:figures/full_fig_p014_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

30 extracted references

  1. [27]

    Simone Rossi, Toni Lassila, Ricardo Ruiz-Baier, Adélia Sequeira, and Alfio Quarteroni. Thermodynamically consistent orthotropic activation model capturing ventricular systolic wall thickening in cardiac electromechan- ics.European Journal of Mechanics - A/Solids, 48:129–142, 2014

    2014

  2. [1]

    lifex-fiber: an open tool for myofibers generation in cardiac computational models.BMC Bioinformatics, 24(1):143, 2023

    Pasquale Claudio Africa, Roberto Piersanti, Marco Fedele, Luca Dedè, and Alfio Quarteroni. lifex-fiber: an open tool for myofibers generation in cardiac computational models.BMC Bioinformatics, 24(1):143, 2023

    2023

  3. [2]

    Active stress vs

    D Ambrosi and Simone Pezzuto. Active stress vs. active strain in mechanobiology: constitutive issues.Journal of Elasticity, 107(2):199–212, 2012

    2012

  4. [3]

    Electromechanical coupling in cardiac dynamics: the active strain approach.SIAM Journal on Applied Mathematics, 71(2):605–621, 2011

    Davide Ambrosi, Gianni Arioli, Fabio Nobile, and Alfio Quarteroni. Electromechanical coupling in cardiac dynamics: the active strain approach.SIAM Journal on Applied Mathematics, 71(2):605–621, 2011

    2011

  5. [4]

    John Wiley & Sons, 2014

    Ted Belytschko, Wing Kam Liu, Brian Moran, and Khalil Elkhodary.Nonlinear finite elements for continua and structures. John Wiley & Sons, 2014. 2https://gitlab.com/gssi-fluids/iphem 15 Active strainA PREPRINT

    2014

  6. [5]

    Cardiac excitation–contraction coupling.Nature, 415(6868):198–205, 2002

    Donald M Bers. Cardiac excitation–contraction coupling.Nature, 415(6868):198–205, 2002

    2002

  7. [6]

    Cambridge University Press, 1997

    Javier Bonet and Richard D Wood.Nonlinear continuum mechanics for finite element analysis. Cambridge University Press, 1997

    1997

  8. [7]

    A fast computational model for the electrophysiology of the whole human heart.Journal of Computational Physics, 457:111084, 2022

    Giulio Del Corso, Roberto Verzicco, and Francesco Viola. A fast computational model for the electrophysiology of the whole human heart.Journal of Computational Physics, 457:111084, 2022

    2022

Show all 30 references
  1. [8]

    Systematic coarse-graining of spectrin-level red blood cell models.Computer Methods in Applied Mechanics and Engineering, 199(29-32):1937–1948, 2010

    Dmitry A Fedosov, Bruce Caswell, and George Em Karniadakis. Systematic coarse-graining of spectrin-level red blood cell models.Computer Methods in Applied Mechanics and Engineering, 199(29-32):1937–1948, 2010

    1937

  2. [9]

    Locking-proof tetrahedra.ACM Trans- actions on Graphics, 40(2):1–17, 2021

    Mihai Frâncu, Arni Asgeirsson, Kenny Erleben, and Mads JL Rønnow. Locking-proof tetrahedra.ACM Trans- actions on Graphics, 40(2):1–17, 2021

    2021

  3. [10]

    Springer Science & Business Media, 2013

    Yuan-cheng Fung.Biomechanics: mechanical properties of living tissues. Springer Science & Business Media, 2013

    2013

  4. [11]

    A new computa- tional framework for electro-activation in cardiac mechanics.Computer Methods in Applied Mechanics and Engineering, 348:796–845, 2019

    Emilio Garcia-Blanco, Rogelio Ortigosa, Antonio J Gil, Chun Hean Lee, and Javier Bonet. A new computa- tional framework for electro-activation in cardiac mechanics.Computer Methods in Applied Mechanics and Engineering, 348:796–845, 2019

    2019

  5. [12]

    Approximate simulation of elastic membranes by triangulated spring meshes.Journal of Graphics Tools, 3(2):21–41, 1998

    Allen Van Gelder. Approximate simulation of elastic membranes by triangulated spring meshes.Journal of Graphics Tools, 3(2):21–41, 1998

    1998

  6. [13]

    The generalized Hill model: a kinematic approach towards active muscle contraction.Journal of the Mechanics and Physics of Solids, 72:20–39, 2014

    Serdar Göktepe, Andreas Menzel, and Ellen Kuhl. The generalized Hill model: a kinematic approach towards active muscle contraction.Journal of the Mechanics and Physics of Solids, 72:20–39, 2014

    2014

  7. [14]

    Finite element stress analysis of left ventricular mechanics in the beating dog heart.Journal of Biomechanics, 28(10):1167–1177, 1995

    Julius M Guccione, Kevin D Costa, and Andrew D McCulloch. Finite element stress analysis of left ventricular mechanics in the beating dog heart.Journal of Biomechanics, 28(10):1167–1177, 1995

    1995

  8. [15]

    Mass-spring model for simulation of heart valve tissue mechanical behavior.Annals of Biomedical Engineering, 39(6):1668–1679, 2011

    Peter E Hammer, Michael S Sacks, Pedro J Del Nido, and Robert D Howe. Mass-spring model for simulation of heart valve tissue mechanical behavior.Annals of Biomedical Engineering, 39(6):1668–1679, 2011

    2011

  9. [16]

    Gerhard A Holzapfel and Ray W Ogden. Constitutive modelling of passive myocardium: a structurally based framework for material characterization.Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 367(1902):3445–3475, 2009

    1902

  10. [17]

    V olume conserving finite element simulations of de- formable models.ACM Transactions on Graphics, 26(3):13–es, 2007

    Geoffrey Irving, Craig Schroeder, and Ronald Fedkiw. V olume conserving finite element simulations of de- formable models.ACM Transactions on Graphics, 26(3):13–es, 2007

    2007

  11. [18]

    Computational modelling and analysis of the coupled aero- structural dynamics in bat-inspired wings.Journal of Fluid Mechanics, 1010:A53, 2025

    Sushrut Kumar, Jung-Hee Seo, and Rajat Mittal. Computational modelling and analysis of the coupled aero- structural dynamics in bat-inspired wings.Journal of Fluid Mechanics, 1010:A53, 2025

    2025

  12. [19]

    Sander Land, Viatcheslav Gurev, Sander Arens, Christoph M Augustin, Lukas Baron, Robert Blake, Chris Bradley, Sebastian Castro, Andrew Crozier, Marco Favino, et al. Verification of cardiac mechanics software: benchmark problems and solutions for testing active and passive mate...

    2015

  13. [20]

    A constraint-based formulation of stable Neo-Hookean materials

    Miles Macklin and Matthias Müller. A constraint-based formulation of stable Neo-Hookean materials. InPro- ceedings of the 14th ACM SIGGRAPH conference on motion, interaction and games, pages 1–7, 2021

    2021

  14. [21]

    XPBD: position-based simulation of compliant constrained dynamics

    Miles Macklin, Matthias Müller, and Nuttapong Chentanez. XPBD: position-based simulation of compliant constrained dynamics. InProceedings of the 9th International Conference on Motion in Games, pages 49–54, 2016

    2016

  15. [22]

    Spring-network-based model of a red blood cell for simulating mesoscopic blood flow.International Journal for Numerical Methods in Biomedical Engineering, 29(1):114–128, 2013

    Masanori Nakamura, Sadao Bessho, and Shigeo Wada. Spring-network-based model of a red blood cell for simulating mesoscopic blood flow.International Journal for Numerical Methods in Biomedical Engineering, 29(1):114–128, 2013

    2013

  16. [23]

    Mathematical model of geometry and fibrous structure of the heart.American Journal of Physiology-Heart and Circulatory Physiology, 260(4):H1365–H1378, 1991

    PM Nielsen, IJ Le Grice, BH Smaill, and PJ Hunter. Mathematical model of geometry and fibrous structure of the heart.American Journal of Physiology-Heart and Circulatory Physiology, 260(4):H1365–H1378, 1991

    1991

  17. [24]

    An orthotropic active–strain model for the myocardium mechanics and its numerical approximation.European Journal of Mechanics - A/Solids, 48:83–96, 2014

    S Pezzuto, D Ambrosi, and A Quarteroni. An orthotropic active–strain model for the myocardium mechanics and its numerical approximation.European Journal of Mechanics - A/Solids, 48:83–96, 2014

    2014

  18. [25]

    An orthotropic electro- viscoelastic model for the heart with stress-assisted diffusion.Biomechanics and Modeling in Mechanobiology, 19(2):633–659, 2020

    Adrienne Propp, Alessio Gizzi, Francesc Levrero-Florencio, and Ricardo Ruiz-Baier. An orthotropic electro- viscoelastic model for the heart with stress-assisted diffusion.Biomechanics and Modeling in Mechanobiology, 19(2):633–659, 2020

    2020

  19. [26]

    Implicit finite incompressible elastodynamics with linear finite elements: A stabilized method in rate form.Computer Methods in Applied Mechanics and Engineering, 311:208–249, 2016

    S Rossi, N Abboud, and G Scovazzi. Implicit finite incompressible elastodynamics with linear finite elements: A stabilized method in rate form.Computer Methods in Applied Mechanics and Engineering, 311:208–249, 2016. 16 Active strainA PREPRINT

    2016

  20. [28]

    Springer Science & Business Media, 2012

    Masao Tanaka, Shigeo Wada, and Masanori Nakamura.Computational biomechanics: theoretical background and biological/biomedical problems, volume 3. Springer Science & Business Media, 2012

    2012

  21. [29]

    GPU-accelerated digital twins of the human heart open new routes for cardiovascular research.Scientific Reports, 13(1):8230, 2023

    Francesco Viola, Giulio Del Corso, Ruggero De Paulis, and Roberto Verzicco. GPU-accelerated digital twins of the human heart open new routes for cardiovascular research.Scientific Reports, 13(1):8230, 2023

    2023

  22. [30]

    High-fidelity model of the human heart: an immersed boundary implementation.Physical Review Fluids, 8(10):100502, 2023

    Francesco Viola, Giulio Del Corso, and Roberto Verzicco. High-fidelity model of the human heart: an immersed boundary implementation.Physical Review Fluids, 8(10):100502, 2023. 17

    2023

Pith tools

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