Numerical study of incompressible fluid dynamics with nonuniform density by the immersed boundary method

We apply the immersed boundary method to the dynamics of an incompressible fluid with a nonuniform density. In order to take into account both the inertial and gravitational effects of the variable density, the penalty immersed boundary (pIB) method is used [Y. Kim and C. S. Peskin, Phys. Fluids 19, 053103 (2007)]. Incompressible fluid motion with a nonuniform density has been extensively explored both experimentally and computationally. We show that the pIB method is a robust and efficient numerical tool for the simulation of fluids with variable density by conducting computations of some example problems: The falling of a heavier fluid surrounded by a lighter fluid and the Rayleigh–Taylor instabilities in two dimensions and three dimensions and the dynamic stabilization of the Rayleigh–Taylor instability.

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

[2]  T. Gedeon,et al.  Modeling arthropod filiform hair motion using the penalty immersed boundary method. , 2008, Journal of biomechanics.

[3]  Chi-Wang Shu Total-variation-diminishing time discretizations , 1988 .

[4]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[5]  Gretar Tryggvason,et al.  Computations of three‐dimensional Rayleigh–Taylor instability , 1990 .

[6]  R. Verzicco,et al.  Combined Immersed Boundary/Large-Eddy-Simulations of Incompressible Three Dimensional Complex Flows , 2006 .

[7]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[8]  Z. Feng,et al.  The immersed boundary-lattice Boltzmann method for solving fluid-particles interaction problems , 2004 .

[9]  Yongsam Kim,et al.  Penalty immersed boundary method for an elastic boundary with mass , 2007 .

[10]  G. Wolf DYNAMIC STABILIZATION OF THE INTERCHANGE INSTABILITY OF A LIQUID-GAS INTERFACE. , 1970 .

[11]  M. Michaelis,et al.  An inverted liquid demonstration , 1991 .

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

[13]  Jean-Luc Guermond,et al.  Approximation of variable density incompressible flows by means of finite elements and finite volumes , 2001 .

[14]  Boyce E. Griffith,et al.  An adaptive, formally second order accurate version of the immersed boundary method , 2007, J. Comput. Phys..

[15]  G. Tryggvason Numerical simulations of the Rayleigh-Taylor instability , 1988 .

[16]  Ronald Scherer,et al.  Modeling coupled aerodynamics and vocal fold dynamics using immersed boundary methods. , 2006, The Journal of the Acoustical Society of America.

[17]  Howard L. Schreyer,et al.  Fluid–membrane interaction based on the material point method , 2000 .

[18]  E Weinan,et al.  Small‐scale structures in Boussinesq convection , 1998 .