Modeling elastic shells immersed in fluid

We describe a numerical method to simulate an elastic shell immersed in a viscous incompressible fluid. The method is developed as an extension of the immersed boundary method using shell equations based on the Kirchhoff-Love and the planar stress hypotheses. A detailed derivation of the shell equations used in the numerical method is presented. This derivation, as well as the numerical method, uses techniques of differential geometry. Our main motivation for developing this method is its use in constructing a comprehensive, threedimensional computational model of the cochlea (the inner ear). The central object of study within the cochlea is the basilar membrane, which is immersed in fluid and whose elastic properties rather resemble those of a shell. We apply the method to a specific example, which is a prototype of a piece of the basilar membrane, and study the convergence of the method in this case. Some typical features of cochlear mechanics are already captured in this simple model. In particular, numerical experiments have shown a traveling wave propagating from the base to the apex of the model shell in response to external excitation in the fluid. c � 2004 Wiley Periodicals, Inc.

[1]  Edward Givelberg,et al.  A comprehensive three-dimensional model of the cochlea , 2003 .

[2]  J. Brackbill,et al.  A numerical method for suspension flow , 1991 .

[3]  K Grosh,et al.  Three-dimensional numerical modeling for global cochlear dynamics. , 2000, The Journal of the Acoustical Society of America.

[4]  O. C. Zienkiewicz,et al.  Analysis of thick and thin shell structures by curved finite elements , 1970 .

[5]  D. Chapelle,et al.  The Finite Element Analysis of Shells - Fundamentals , 2003 .

[6]  C S Peskin,et al.  A general method for the computer simulation of biological systems interacting with fluids. , 1995, Symposia of the Society for Experimental Biology.

[7]  C. Peskin,et al.  Heart Simulation by an Immersed Boundary Method with Formal Second-order Accuracy and Reduced Numerical Viscosity , 2001 .

[8]  C. Peskin The immersed boundary method , 2002, Acta Numerica.

[9]  P. Coleman,et al.  Experiments in hearing , 1961 .

[10]  K. Bathe,et al.  On the finite element analysis of shells and their full interaction with Navier-stokes fluid flows , 2002 .

[11]  A. Fogelson,et al.  Truncated newton methods and the modeling of complex immersed elastic structures , 1993 .

[12]  Charles S. Peskin,et al.  Fluid Flow in Collapsible Elastic Tubes: A Three-Dimensional Numerical Model , 2001 .

[13]  Philippe G. Ciarlet,et al.  Introduction to Linear Shell Theory , 1989 .

[14]  Dominique Chapelle,et al.  Some new results and current challenges in the finite element analysis of shells , 2001, Acta Numerica.

[15]  R. Valid The Nonlinear Theory of Shells Through Variational Principles: From Elementary Algebra to Differential Geometry , 1996 .

[16]  G. Békésy,et al.  Experiments in Hearing , 1963 .

[17]  J. P. Beyer A computational model of the cochlea using the immersed boundary method , 1992 .

[18]  L. Taber,et al.  Comparison of WKB calculations and experimental results for three-dimensional cochlear models. , 1979, The Journal of the Acoustical Society of America.

[19]  C. Peskin,et al.  Fluid Dynamics of the Heart and its Valves , 1996 .

[20]  A. Fogelson,et al.  A fast numerical method for solving the three-dimensional Stokes' equations in the presence of suspended particles , 1988 .

[21]  Randall J. LeVeque,et al.  Solution of a two-dimensional Cochlea model with fluid viscosity , 1988 .