In this work, a multiple-relaxation-time (MRT) lattice Boltzmann method (LBM) is proposed to solve a coupled chemotaxis-fluid model. In the evolution equation of the proposed LBM, Beam-Warming (B-W) scheme is used to enhance the numerical stability. In numerical experiments, at first, the comparison between the classical LBM and the present LBM with B-W scheme is carried out by simulating blow up phenomenon of the Keller-Segel (K-S) model. Numerical results show that the stability of the present LBM with B-W scheme is better than the classical one. Then, the second order convergence rate of the proposed B-W scheme is verified in the numerical study of the coupled Navier-Stokes (N-S) K-S model. Finally, through solving the coupled chemotaxis-fluid model, the formation of falling bacterial plumes is numerically investigated. Numerical results agree well with existing ones in the literature.
Citation: Zhonghua Qiao, Xuguang Yang. A multiple-relaxation-time lattice Boltzmann method with Beam-Warming scheme for a coupled chemotaxis-fluid model[J]. Electronic Research Archive, 2020, 28(3): 1207-1225. doi: 10.3934/era.2020066
In this work, a multiple-relaxation-time (MRT) lattice Boltzmann method (LBM) is proposed to solve a coupled chemotaxis-fluid model. In the evolution equation of the proposed LBM, Beam-Warming (B-W) scheme is used to enhance the numerical stability. In numerical experiments, at first, the comparison between the classical LBM and the present LBM with B-W scheme is carried out by simulating blow up phenomenon of the Keller-Segel (K-S) model. Numerical results show that the stability of the present LBM with B-W scheme is better than the classical one. Then, the second order convergence rate of the proposed B-W scheme is verified in the numerical study of the coupled Navier-Stokes (N-S) K-S model. Finally, through solving the coupled chemotaxis-fluid model, the formation of falling bacterial plumes is numerically investigated. Numerical results agree well with existing ones in the literature.
[1] | The Keller-Segel model for chemotaxis with prevention of overcrowding: linear vs. nonlinear diffusion. SIAM J. Math. Anal. (2006) 38: 1288-1315. |
[2] | Z. Chai, H. Liang, R. Du and B. Shi, A lattice Boltzmann model for two-phase flow in porous media, SIAM J. Sci. Comput., 41 (2019), B746–B772. doi: 10.1137/18M1166742 |
[3] | Z. H. Chai and T. S. Zhao, Lattice Boltzmann model for the convection-diffusion equation, Phy. Rev. E, 87 (2013), 063309. |
[4] | Lattice Boltzmann method for fluid flows. Annu. Rev. Fluid Mech. (1998) 30: 329-364. |
[5] | Sinking, merging and stationary plumes in a coupled chemotaxis-fluid model: A high-resolution numerical approach. J. Fluid Mech. (2012) 694: 155-190. |
[6] | Numerical study of plume patterns in a chemotaxis-diffusion-convection coupling system. Comput. Fluids (2016) 126: 58-70. |
[7] | The Keller-Segel model with logistic sensitivity function and small diffusivity. SIAM J. Appl. Math. (2005) 66: 286-308. |
[8] | C. Dombrowski, L. Cisneros, S. Chatkaew, R. E. Goldstein and J. O. Kessler, Self-concentration and large-scale coherence in bacterial dynamics, Phys. Rev. Lett., 93 (2004), 098103. doi: 10.1103/PhysRevLett.93.098103 |
[9] | (2004) Chemotaxis. London: Imperial College Press. |
[10] | Upwind-difference potentials method for Patlak-Keller-Segel chemotaxis model. J. Sci. Comput. (2012) 53: 689-713. |
[11] | Discontinuous Galerkin methods for the chemotaxis and haptotaxisn models. J. Comput. Appl. Math. (2009) 224: 168-181. |
[12] | New interior penalyt discontinous Galerkin methods for the Keller-Segel chemotaxis model. SIAM J. Numer. Anal. (2008/2009) 47: 386-408. |
[13] | A finit volume scheme for the Patlak-Keller-Segel chemotaxis model. Numer. Math. (2006) 104: 457-488. |
[14] | Bioconvection. Fluid Dynam. Res. (2005) 37: 1-20. |
[15] | Bioconvection in suspensions of oxytactic bacteria: Linear theory. J. Fluid Mech. (1996) 324: 223-259. |
[16] | Lattice-Boltzmann model for bacterial chemotaxis. J. Math. Biol. (2005) 51: 302-332. |
[17] | Numerical investigation of falling bacterial plumes caused by bioconvection in a three-dimensional chamber. Eur. J. Mech. B Fluid. (2015) 52: 120-130. |
[18] | Recent advances in lattice Boltzmann computing.. Annu. Rev. Comput. Phys. (1995) 3: 195-242. |
[19] | Thermodynamic-consistent multiple-relaxation-time lattice Boltzmann equation model for two-phase hydrocarbon fluids with Peng-Robinson equation of state. Int. J. Heat Mass Tran. (2019) 141: 1216-1226. |
[20] | Stability of operator splitting methods for systems with indefinite poerators: advection-diffusion-reaction systems. J. Comput. Phys. (2009) 228: 3508-3516. |
[21] | Conservative upwind finite-element method for a simplified Keller-Segel system modeling chemotaxis. IMA J. Numer. Anal. (2007) 27: 332-365. |
[22] | Numerical investigation of chemotaxic phenomenon in incompressible viscous fluid flow. Comput. Fluids (2014) 103: 290-306. |
[23] | B. C. Shi and Z. L. Guo, Lattice Boltzmann model for nonlinear conveciton-diffusion equations, Phy. Rev. E, 79 (2009), 016701. |
[24] | Bacterial swimming and oxygen transport near contact lines. Proc. Natl. Acad. Sci. (2005) 102: 2277-2282. |
[25] | Fractional step methods applied to a chemotaxis model. J. Math. Biol. (2000) 41: 455-475. |
[26] | A lattice Boltzmann analysis of the conjugate natural convection in square enclosure with a circular cylinder. Appl. Math. Model. (2019) 71: 31-44. |
[27] | A multiple-relaxation-time Lattice-Boltzmann model for bacterial chemotaxis: Effects of initial concentration, diffusion, and hydrodynamic dispersion on traveling bacterial bands. Bull Math. Biol. (2014) 76: 2449-2475. |
[28] | A coupled lattice Boltzmann method to solve Nernst-Planck model for simulating electro-osmotic flows. J. Sci. Comput. (2014) 61: 222-238. |
[29] | T. Zhang, B. C. Shi, Z. L. Guo, Z. H. Chai and J. H. Lu, Genearl bounce-back scheme for concentration boundary condition in the lattice Boltzmann method, Phys. Rev. E, 85 (2012), 016701. |