An improved immersed boundary-lattice Boltzmann method for simulating three-dimensional incompressible flows

The recently proposed boundary condition-enforced immersed boundary-lattice Boltzmann method (IB-LBM) [14] is improved in this work to simulate three-dimensional incompressible viscous flows. In the conventional IB-LBM, the restoring force is pre-calculated, and the non-slip boundary condition is not enforced as compared to body-fitted solvers. As a result, there is a flow penetration to the solid boundary. This drawback was removed by the new version of IB-LBM [14], in which the restoring force is considered as unknown and is determined in such a way that the non-slip boundary condition is enforced. Since Eulerian points are also defined inside the solid boundary, the computational domain is usually regular and the Cartesian mesh is used. On the other hand, to well capture the boundary layer and in the meantime, to save the computational effort, we often use non-uniform mesh in IB-LBM applications. In our previous two-dimensional simulations [14], the Taylor series expansion and least squares-based lattice Boltzmann method (TLLBM) was used on the non-uniform Cartesian mesh to get the flow field. The final expression of TLLBM is an algebraic formulation with some weighting coefficients. These coefficients could be computed in advance and stored for the following computations. However, this way may become impractical for 3D cases as the memory requirement often exceeds the machine capacity. The other way is to calculate the coefficients at every time step. As a result, extra time is consumed significantly. To overcome this drawback, in this study, we propose a more efficient approach to solve lattice Boltzmann equation on the non-uniform Cartesian mesh. As compared to TLLBM, the proposed approach needs much less computational time and virtual storage. Its good accuracy and efficiency are well demonstrated by its application to simulate the 3D lid-driven cubic cavity flow. To valid the combination of proposed approach with the new version of IBM [14] for 3D flows with curved boundaries, the flows over a sphere and torus are simulated. The obtained numerical results compare very well with available data in the literature.

[1]  Wei Shyy,et al.  On the Finite Difference-Based Lattice Boltzmann Method in Curvilinear Coordinates , 1998 .

[2]  Shiyi Chen,et al.  LATTICE BOLTZMANN METHOD FOR FLUID FLOWS , 2001 .

[3]  Z. Feng,et al.  Proteus: a direct forcing method in the simulations of particulate flows , 2005 .

[4]  P. Colella,et al.  A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains , 1998 .

[5]  M. Thompson,et al.  From spheres to circular cylinders: non-axisymmetric transitions in the flow past rings , 2004, Journal of Fluid Mechanics.

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

[7]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[8]  F. White.  Viscous Fluid Flow , 1974 .

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

[10]  V. C. Patel,et al.  Flow past a sphere up to a Reynolds number of 300 , 1999, Journal of Fluid Mechanics.

[11]  V. Babu,et al.  Numerical solution of the incompressible three-dimensional Navier-Stokes equations , 1994 .

[12]  J. Wu,et al.  Implicit velocity correction-based immersed boundary-lattice Boltzmann method and its applications , 2009 .

[13]  Chang Shu,et al.  Simulation of flows around an impulsively started circular cylinder by Taylor series expansion-and least squares-based lattice Boltzmann method , 2003 .

[14]  Chang Shu,et al.  TAYLOR SERIES EXPANSION AND LEAST SQUARES-BASED LATTICE BOLTZMANN METHOD: THREE-DIMENSIONAL FORMULATION AND ITS APPLICATIONS , 2003 .

[15]  J. Boon.  The Lattice Boltzmann Equation for Fluid Dynamics and Beyond , 2003 .

[16]  H. S. Udaykumar,et al.  A Sharp Interface Cartesian Grid Methodfor Simulating Flows with ComplexMoving Boundaries , 2001 .

[17]  Fotis Sotiropoulos,et al.  A numerical method for solving the 3D unsteady incompressible Navier-Stokes equations in curvilinear domains with complex immersed boundaries , 2007, J. Comput. Phys..

[18]  W. Shyy,et al.  A multi‐block lattice Boltzmann method for viscous fluid flows , 2002 .

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

[20]  F. Sotiropoulos,et al.  A hybrid Cartesian/immersed boundary method for simulating flows with 3D, geometrically complex, moving bodies , 2005 .

[21]  Chi-Wang Shu,et al.  Taylor-series expansion and least-squares-based lattice Boltzmann method: Two-dimensional formulation and its applications. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[22]  M. Thompson,et al.  From spheres to circular cylinders: the stability and flow structures of bluff ring wakes , 2003, Journal of Fluid Mechanics.

[23]  R. Verzicco,et al.  Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations , 2000 .

[24]  Y. Qian,et al.  Lattice BGK Models for Navier-Stokes Equation , 1992 .

[25]  Yan Peng,et al.  A momentum exchange-based immersed boundary-lattice Boltzmann method for simulating incompressible viscous flows , 2006 .

[26]  E. Balaras,et al.  A general reconstruction algorithm for simulating flows with complex 3D immersed boundaries on Cartesian grids , 2003 .

[27]  So-Hsiang Chou,et al.  FINITE VOLUME SCHEME FOR THE LATTICE BOLTZMANN METHOD ON UNSTRUCTURED MESHES , 1999 .

[28]  Mark C. Thompson,et al.  Computations of the drag coefficients for low-Reynolds-number flow past rings , 2005, Journal of Fluid Mechanics.

[29]  Jinhee Jeong,et al.  On the identification of a vortex , 1995, Journal of Fluid Mechanics.

[30]  Y. T. Chew,et al.  Application of multi-block approach in the immersed boundary-lattice Boltzmann method for viscous fluid flows , 2006, J. Comput. Phys..

[31]  Xiaoyi He,et al.  Lattice Boltzmann method on a curvilinear coordinate system: Vortex shedding behind a circular cylinder , 1997 .