Commun. Comput. Phys., 12 (2012), pp. 433-461.


Simulating an Elastic Ring with Bend and Twist by an Adaptive Generalized Immersed Boundary Method

Boyce E. Griffith 1, Sookkyung Lim 2*

1 Leon H. Charney Division of Cardiology, Department of Medicine, New York University School of Medicine, 550 First Avenue, New York, New York 10016, USA.
2 Department of Mathematical Sciences, University of Cincinnati, 839 Old Chemistry Building, Cincinnati, Ohio 45221, USA.

Received 19 February 2011; Accepted (in revised version) 6 August 2011
Available online 20 February 2012
doi:10.4208/cicp.190211.060811s

Abstract

Many problems involving the interaction of an elastic structure and a viscous fluid can be solved by the immersed boundary (IB) method. In the IB approach to such problems, the elastic forces generated by the immersed structure are applied to the surrounding fluid, and the motion of the immersed structure is determined by the local motion of the fluid. Recently, the IB method has been extended to treat more general elasticity models that include both positional and rotational degrees of freedom. For such models, force and torque must both be applied to the fluid. The positional degrees of freedom of the immersed structure move according to the local linear velocity of the fluid, whereas the rotational degrees of freedom move according to the local angular velocity. This paper introduces a spatially adaptive, formally second-order accurate version of this generalized immersed boundary method. We use this adaptive scheme to simulate the dynamics of an elastic ring immersed in fluid. To describe the elasticity of the ring, we use an unconstrained version of Kirchhoff rod theory. We demonstrate empirically that our numerical scheme yields essentially second-order convergence rates when applied to such problems. We also study dynamical instabilities of such fluid-structure systems, and we compare numerical results produced by our method to classical analytic results from elastic rod theory.

AMS subject classifications: 65M06, 65M50, 74S20, 76D05
Key words: Immersed boundary method, Kirchhoff rod theory, adaptive mesh refinement.

*Corresponding author.
Email: boyce.griffith@nyumc.org (B. E. Griffith), limsk@math.uc.edu (S. Lim)
 

The Global Science Journal