A continuous energy-based immersed boundary method for elastic shells

TL;DR

This paper introduces a novel continuous energy-based immersed boundary method for elastic shells, improving force computation accuracy in fluid-structure interaction simulations involving biological membranes.

Contribution

It presents a new continuous force computation method for elastic shells, enhancing accuracy over previous discretized approaches in immersed boundary simulations.

Findings

More accurate force calculations than previous methods

Improved geometric information such as mean curvature

Successful application to biological cell models

Abstract

The immersed boundary method is a mathematical formulation and numerical method for solving fluid-structure interaction problems. For many biological problems, such as models that include the cell membrane, the immersed structure is a two-dimensional infinitely thin elastic shell immersed in an incompressible viscous fluid. When the shell is modeled as a hyperelastic material, forces can be computed by taking the variational derivative of an energy density functional. A new method for computing a continuous force function on the entire surface of the shell is presented here. The new method is compared to a previous formulation where the surface and energy functional are discretized before forces are computed. For the case of Stokes flow, a method for computing quadrature weights is provided to ensure the integral of the elastic spread force density remains zero throughout a dynamic simulation. Tests on the method are conducted and show that it yields more accurate force computations than previous formulations as well as more accurate geometric information such as mean curvature. The method is then applied to a model of a red blood cell in capillary flow and a 3D model of cellular blebbing.