第39卷第6期2022年11月计算物理CHINESEJOURNALOFCOMPUTATIONALPHYSICSVol.39,No.6Nov.,2022文章编号:1001⁃246X(2022)06⁃0666⁃11收稿日期:2021-12-21;修回日期:2022-02-24基金项目:湖南省自然科学基金(2021JJ30548)、湖南省教育厅科学研究重点项目(21A0329)、吉首大学研究生校级科研基金(Jdy21014)和吉首大学研究生科研创新项目(JGY2022075)资助第一作者:冯舒婷(1997-),女,硕士研究生,研究方向为微分方程数值解,E⁃mail:fst202111@163.com∗通信作者:戴厚平(1979-),男,博士,副教授,研究方向为微分方程数值解及其应用,E⁃mail:daihouping@163.com二维空间分数阶反应扩散方程组的格子Boltzmann方法冯舒婷,戴厚平∗,宋通政(吉首大学数学与统计学院,湖南吉首416000)摘要:提出一种二维空间分数阶反应扩散方程组的格子Boltzmann方法,对方程组的分数阶积分项进行离散化处理。通过Chapman⁃Enskog多尺度技术和Taylor展开技术,从所建立的模型中恢复出二维空间分数阶反应扩散方程组,推导出各个速度方向上的平衡态分布函数。通过数值算例对所构造的LB模型进行检验,数值计算结果与算例的精确解吻合较好。关键词:Riemann⁃Liouville分数阶导数;反应扩散方程组;格子Boltzmann方法;Chapman⁃Enskog多尺度技术中图分类号:O241.82文献标识码:ADOI:10.19596/j.cnki.1001⁃246x.84950引言分数阶微分方程被广泛的应用于物理力学[1-4]、自动控制[5-6]、信号处理[7-9]、生物医学[10-11]等领域。然而,大多数分数阶微分方程的解析解不能显式得到,因此,寻找有效的数值方法求解分数阶微分方程尤其重要。目前,对于分数阶微分方程的数值解法有谱方法[12]、有限差分方法[13-14]、有限元方法[15-16]等。例如Zheng等[17]基于时空谱方法,提出了一种求解多时间分数阶扩散方程的高阶数值方法。Fazio等[18]在非均匀网格上用无条件稳定的隐式有限差分方法研究了含源项的Caputo导数的时间分数阶对流扩散方程的数值解...