The Blob Projection Method for Immersed Boundary Problems

A new finite difference numerical method for modeling the interaction between flexible elastic membranes and an incompressible fluid in a two-dimensional domain is presented. The method differs from existing methods in the way the forces exerted by the membranes on the fluid are modeled. These are described by a collection of regularized point forces, and the velocity field they induce is computed directly on a regular Cartesian grid via a smoothed dipole potential. Comparisons between this method and the immersed boundary method of Peskin and McQueen are presented. The results show that the method proposed here preserves volumes better and has a higher order of convergence.

[1]  Alexandre J. Chorin,et al.  On the Convergence of Discrete Approximations to the Navier-Stokes Equations , 1969 .

[2]  C. Peskin.  Numerical analysis of blood flow in the heart , 1977 .

[3]  A. Fogelson.  A MATHEMATICAL MODEL AND NUMERICAL METHOD FOR STUDYING PLATELET ADHESION AND AGGREGATION DURING BLOOD CLOTTING , 1984 .

[4]  Andrew J. Majda,et al.  High order accurate vortex methods with explicit velocity kernels , 1985 .

[5]  C. R. Anderson.  A method of local corrections for computing the velocity field due to a distribution of vortex blobs , 1986 .

[6]  L. Fauci,et al.  A computational model of aquatic animal locomotion , 1988 .

[7]  C. Peskin,et al.  A three-dimensional computational method for blood flow in the heart. 1. Immersed elastic fibers in a viscous incompressible fluid , 1989 .

[8]  L. Fauci.  Interaction of oscillating filaments: a computational study , 1990 .

[9]  Ole H. Hald,et al.  Convergence of vortex methods , 1991 .

[10]  E Weinan,et al.  A numerical resolution study of high order essentially non-oscillatory schemes applied to incompressible flow , 1992 .

[11]  T. F. Buttke.  Lagrangian Numerical Methods Which Preserve the Hamiltonian Structure of Incompressible Fluid Flow , 1992 .

[12]  Charles S. Peskin,et al.  Improved Volume Conservation in the Computation of Flows with Immersed Elastic Boundaries , 1993 .

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

[14]  Tomas F. Buttke,et al.  Velicity methods: Lagrangian numerical methods which preserve the Hamiltonian structure of incompres , 1993 .

[15]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[16]  T. Hou,et al.  Removing the stiffness from interfacial flows with surface tension , 1994 .

[17]  Charles S. Peskin,et al.  A three-dimensional computer model for fluid flow through a collapsible tube , 1994 .

[18]  L. Fauci,et al.  Sperm motility in the presence of boundaries. , 1995, Bulletin of mathematical biology.

[19]  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.

[20]  R. Cortez.  An Impulse-Based Approximation of Fluid Motion due to Boundary Forces , 1996 .

[21]  B. Wetton,et al.  Analysis and computation of immersed boundaries, with application to pulp fibres , 1997 .

[22]  R. Cortez,et al.  The Dynamics of an Elastic Membrane Using the Impulse Method , 1997 .

[23]  J. Thomas Beale,et al.  A Method for Computing Nearly Singular Integrals , 2000, SIAM J. Numer. Anal..