# American Institute of Mathematical Sciences

September  2020, 28(3): 1207-1225. doi: 10.3934/era.2020066

## A multiple-relaxation-time lattice Boltzmann method with Beam-Warming scheme for a coupled chemotaxis-fluid model

 1 Department of Applied Mathematics, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong, China 2 School of Mathematics and Computational Science, Hunan First Normal University, Changsha, China

* Corresponding author: Xuguang Yang

Received  April 2020 Revised  June 2020 Published  July 2020

Fund Project: Z. Qiao's work is partially supported by Hong Kong Research Council GRF grant 15300417. X. Yang's work is supported by the Natural Science Foundation of China (Grant No. 11802090) and was also supported by the Hong Kong Polytechnic University Postdoctoral Fellowships Scheme (No. 1-YW1D)

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. Electronic Research Archive, 2020, 28 (3) : 1207-1225. doi: 10.3934/era.2020066
##### References:

show all references

##### References:
Numerical solutions of bacteria concentrations $n(\mathbf{x}, t)$ of (27)-(28) by the standard LBM (A = 0): (a) $t = 5.0\times10^{-6}$, (b) $t = 1.0\times10^{-5}$, (c) $t = 4.5\times10^{-5}$, (d) $t = 1.15\times10^{-4}$
Numerical solutions of bacteria concentrations $n(\mathbf{x}, t)$ of (27)-(28) by the present LBM with B-W scheme (A = 0.5), (a) $t = 5.0\times10^{-6}$, (b) $t = 1.0\times10^{-5}$, (c) $t = 4.5\times10^{-5}$, (d) $t = 1.15\times10^{-4}$
1D profiles of $n(\mathbf{x}, t)$ along y = 0 at $t = 1.15\times10^{-4}$: (a) by the standard LBM (A = 0); (b) by the present LBM with B-W scheme (A = 0.5)
(a) Vertical profiles of the oxygen $c$ and cell densities $n$ and (b) streamlines at t = 0.22
Vertical profiles of (a) oxygen $c$ and (b) cell densities $n$ at t = 0.22 with $\alpha = 5$, $\delta = 1$ and different $\beta$
Vertical profiles of (a) oxygen $c$ and (b) cell densities $n$ at t = 0.22 with $\beta = 10$, $\delta = 1$ and different $\alpha$
Bioconvection patterns of the cell density $n(x, y)$ at (a) $t = 0.32$, (b) $t = 0.37$, (c) $t = 0.42$, and (d) $t = 0.47$
For $\delta = 5$, cell concentration $n(x, y)$ at (a) $t = 0.27$, (b) $t = 0.34$, and (c) $t = 0.39$
For $\delta = 25$, cell concentration $n(x, y)$ at (a) $t = 2.9$, (b) $t = 3.4$, and (c) $t = 3.9$
For $\delta = 50$, cell concentration $n(x, y)$ at (a) $t = 1.9$, (b) $t = 2.4$, and (c) $t = 2.9$
Convergence test with different lattice spacings, $\delta t = 1.0\times10^{-5}$
 $Meshes$ $E_n$ order $E_c$ order $E_\mathbf{u}$ order $64\times 64$ $1.3934 \times 10^{-4}$ – $1.6327 \times 10^{-4}$ – $2.0029 \times 10^{-4}$ – $128\times 128$ $3.6892 \times 10^{-5}$ 1.9172 $4.3355 \times10^{-5}$ 1.9130 $5.3150 \times 10^{-5}$ 1.9139 $256\times 256$ $9.4955 \times 10^{-6}$ 1.9580 $1.1174\times 10^{-5}$ 1.9560 $1.3699 \times 10^{-5}$ 1.9560 $512\times 512$ $2.4090 \times 10^{-6}$ 1.9788 $2.8380\times 10^{-6}$ 1.9772 $3.4775 \times 10^{-6}$ 1.9779
 $Meshes$ $E_n$ order $E_c$ order $E_\mathbf{u}$ order $64\times 64$ $1.3934 \times 10^{-4}$ – $1.6327 \times 10^{-4}$ – $2.0029 \times 10^{-4}$ – $128\times 128$ $3.6892 \times 10^{-5}$ 1.9172 $4.3355 \times10^{-5}$ 1.9130 $5.3150 \times 10^{-5}$ 1.9139 $256\times 256$ $9.4955 \times 10^{-6}$ 1.9580 $1.1174\times 10^{-5}$ 1.9560 $1.3699 \times 10^{-5}$ 1.9560 $512\times 512$ $2.4090 \times 10^{-6}$ 1.9788 $2.8380\times 10^{-6}$ 1.9772 $3.4775 \times 10^{-6}$ 1.9779
 [1] Alina Chertock, Alexander Kurganov, Mária Lukáčová-Medvi${\rm{\check{d}}}$ová, Șeyma Nur Özcan. An asymptotic preserving scheme for kinetic chemotaxis models in two space dimensions. Kinetic & Related Models, 2019, 12 (1) : 195-216. doi: 10.3934/krm.2019009 [2] Christoforidou Amalia, Christian-Oliver Ewald. A lattice method for option evaluation with regime-switching asset correlation structure. Journal of Industrial & Management Optimization, 2021, 17 (4) : 1729-1752. doi: 10.3934/jimo.2020042 [3] Haibo Cui, Haiyan Yin. Convergence rate of solutions toward stationary solutions to the isentropic micropolar fluid model in a half line. Discrete & Continuous Dynamical Systems - B, 2021, 26 (6) : 2899-2920. doi: 10.3934/dcdsb.2020210 [4] Harumi Hattori, Aesha Lagha. Global existence and decay rates of the solutions for a chemotaxis system with Lotka-Volterra type model for chemoattractant and repellent. Discrete & Continuous Dynamical Systems, 2021  doi: 10.3934/dcds.2021071 [5] Hai-Liang Li, Tong Yang, Mingying Zhong. Diffusion limit of the Vlasov-Poisson-Boltzmann system. Kinetic & Related Models, 2021, 14 (2) : 211-255. doi: 10.3934/krm.2021003 [6] Maolin Cheng, Yun Liu, Jianuo Li, Bin Liu. Nonlinear Grey Bernoulli model NGBM (1, 1)'s parameter optimisation method and model application. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021054 [7] Zhang Chen, Xiliang Li, Bixiang Wang. Invariant measures of stochastic delay lattice systems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (6) : 3235-3269. doi: 10.3934/dcdsb.2020226 [8] G. Deugoué, B. Jidjou Moghomye, T. Tachim Medjo. Approximation of a stochastic two-phase flow model by a splitting-up method. Communications on Pure & Applied Analysis, 2021, 20 (3) : 1135-1170. doi: 10.3934/cpaa.2021010 [9] Dugan Nina, Ademir Fernando Pazoto, Lionel Rosier. Controllability of a 1-D tank containing a fluid modeled by a Boussinesq system. Evolution Equations & Control Theory, 2013, 2 (2) : 379-402. doi: 10.3934/eect.2013.2.379 [10] Rama Ayoub, Aziz Hamdouni, Dina Razafindralandy. A new Hodge operator in discrete exterior calculus. Application to fluid mechanics. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2021062 [11] Xinyuan Liao, Caidi Zhao, Shengfan Zhou. Compact uniform attractors for dissipative non-autonomous lattice dynamical systems. Communications on Pure & Applied Analysis, 2007, 6 (4) : 1087-1111. doi: 10.3934/cpaa.2007.6.1087 [12] Caifang Wang, Tie Zhou. The order of convergence for Landweber Scheme with $\alpha,\beta$-rule. Inverse Problems & Imaging, 2012, 6 (1) : 133-146. doi: 10.3934/ipi.2012.6.133 [13] Roberto Civino, Riccardo Longo. Formal security proof for a scheme on a topological network. Advances in Mathematics of Communications, 2021  doi: 10.3934/amc.2021009 [14] Zaihong Wang, Jin Li, Tiantian Ma. An erratum note on the paper: Positive periodic solution for Brillouin electron beam focusing system. Discrete & Continuous Dynamical Systems - B, 2013, 18 (7) : 1995-1997. doi: 10.3934/dcdsb.2013.18.1995 [15] Bochao Chen, Yixian Gao. Quasi-periodic travelling waves for beam equations with damping on 3-dimensional rectangular tori. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021075 [16] Wolf-Jüergen Beyn, Janosch Rieger. The implicit Euler scheme for one-sided Lipschitz differential inclusions. Discrete & Continuous Dynamical Systems - B, 2010, 14 (2) : 409-428. doi: 10.3934/dcdsb.2010.14.409 [17] Mia Jukić, Hermen Jan Hupkes. Dynamics of curved travelling fronts for the discrete Allen-Cahn equation on a two-dimensional lattice. Discrete & Continuous Dynamical Systems, 2021, 41 (7) : 3163-3209. doi: 10.3934/dcds.2020402 [18] Ying Yang. Global classical solutions to two-dimensional chemotaxis-shallow water system. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2625-2643. doi: 10.3934/dcdsb.2020198 [19] Hong Yi, Chunlai Mu, Guangyu Xu, Pan Dai. A blow-up result for the chemotaxis system with nonlinear signal production and logistic source. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2537-2559. doi: 10.3934/dcdsb.2020194 [20] Xu Pan, Liangchen Wang. Boundedness and asymptotic stability in a quasilinear two-species chemotaxis system with nonlinear signal production. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2021064

Impact Factor: 0.263