ABSTRACT
The immersed boundary method is a general numerical method for modeling elastic boundaries immersed within a viscous, incompressible fluid. It has been applied to biological and engineering systems, including large-scale models of the heart and cochlea. Despite the popularity of the immersed boundary method and the desire to scale the problems to accurately capture the details of the physical systems, parallelization for large-scale distributed memory machine has proven challenging. The primary reason is a classic locality and load balance tradeoff that arises in distributing the immersed boundary data structure across processors. In this paper we describe a parallelized algorithm for the immersed boundary method that is designed for scalability on distributed memory computers. It is implemented using the Titanium language, a Java-based language designed for high-performance scientific computing. Our software package, called IB, takes advantage of the object-oriented features of Titanium to provide a framework for simulating immersed boundaries that separates the generic immersed boundary method code from the specific application features that define the immersed boundary structure and the forces that arise from those structures. We demonstrate the scalability of our design and the feasibility of large-scale immersed boundary computations with the IB package.