An immersed boundary energy-based method for incompressible viscoelasticity

Dharshi Devendran, Charles S. Peskin

Research output: Contribution to journalArticlepeer-review


Incompressible viscoelastic materials are prevalent in biological applications. In this paper we present a method for incompressible viscoelasticity in which the elasticity of the material is described in Lagrangian form (i.e. in material coordinates), and Eulerian (spatial) coordinates are used for the equations of motion and to enforce the incompressibility condition. The elastic forces are computed directly from an energy functional without the use of stress tensors, and the immersed boundary method is used to communicate between Lagrangian and Eulerian variables. The method is first applied to a warm-up problem, in which a viscoelastic incompressible material fills a two-dimensional periodic domain. For this problem, we study convergence of the velocity field, the deformation map, and the Eulerian force density. The numerical results indicate that the velocity field and deformation map converge strongly at second order and the Eulerian force density converges weakly at second order. Incompressibility is well maintained, as indicated by area conservation in this 2D problem. Finally, the method is applied to a three-dimensional fluid-structure interaction problem with two different materials: an isotropic neo-Hookean model and an anisotropic fiber-reinforced model.

Original languageEnglish (US)
Pages (from-to)4613-4642
Number of pages30
JournalJournal of Computational Physics
Issue number14
StatePublished - May 20 2012


  • Energy-based method
  • Fluid-structure interaction
  • Immersed boundary method
  • Incompressible viscoelasticity

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • General Physics and Astronomy
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'An immersed boundary energy-based method for incompressible viscoelasticity'. Together they form a unique fingerprint.

Cite this