Characteristic splitting mixed finite element analysis of Keller-Segel chemotaxis models

详细信息    查看全文

• 作者：Jiansong Zhang ; ^a ; ^b ; ^{jszhang@upc.edu.cn" class="auth_mail" title="E-mail the corresponding author} ; Jiang Zhu^b ; ^{jiang@lncc.br" class="auth_mail" title="E-mail the corresponding author} ; Rongpei Zhang^b ; ^c ; ^{rpzhang@lncc.br" class="auth_mail" title="E-mail the corresponding author}
• 关键词：Splitting mixed system ; Mass conservation ; The method of characteristics ; Keller&ndash ; Segel model
• 刊名：Applied Mathematics and Computation
• 出版年：2016
• 出版时间：31 March 2016
• 年：2016
• 卷：278
• 期：Complete
• 页码：33-44
• 全文大小：1258 K

文摘

Chemotaxis models are described by a system of coupled nonlinear partial differential equations: a convection–diffusion type cell density equation and a parabolic or reaction–diffusion type chemical concentration equation. This paper is focused on development of a new numerical method for these models. In this method, a splitting mixed element technique is used to solve the chemical concentration equation, where the LBB condition is not necessary and the flux equation is separated from the chemical concentration equation. Moreover a mass-conservative characteristic finite element method is used to solve the cell density equation, avoiding nonphysically oscillations and keeping mass balance globally. The coefficient matrices of the discrete system are symmetric and positive definite. The convergence of this method is analyzed and the corresponding error estimates are derived for pre-blow-up time since we assume boundedness of the exact solution. Numerical simulations are performed to verify that our method can obtain accurate, nonnegative, oscillation-free solutions and has the ability to recover blowing-up solutions.