模拟组合KdV方程孤立波的LBGK模型.pdfVIP

第 24卷 第 6期 福建师范大学学报 (自然科学版) Vo1．24 No．6 2008年 11月 JournalofFujianNormalUniversity (NaturalScienceEdition) NOV．2008 文章编号：1000—5277(2008)06—0021—04 模拟组合KdV方程孤立波的LBGK模型 黄宝华，马昌凤 (福建师范大学数学与计算机科学学院，福建 福州 350007) 摘要：构造了一个求解组合KdV方程 +6uu 一 “+ “ 一0孤立波解的BGK型格子Bohzmann 模型．用组合KdV方程中的系数，确定了局部平衡分布函数chapman—Enskog展开中的可调参数，并对该方法 进行 了讨论．数值模拟与理论结果相吻合． 关键词：组合KdV方程；LBGK模型；数值模拟 中图分类号：O411．3；O414．2 文献标识码：A A LBGK M odelforSimulatingSolitaryW aves oftheCombinedKdV Equation HUANG Bao—hua，M A Chang—feng (SchoolofMathematicsandComputerScience，FufianNo1malUniversity，Fuzhou350007，China) Abstract：A LBGK modelforsimulatingsolitarywavesofthecombinedKdV equation， + 6uu ～ Ou “+ “一 = 0isestablished．ThetunableparametersinChapman—Enskogex— pansionofthelocal—rquilibrium distributionfunctionaredeterminedbythecoefficientofthe combinedKdV equation． Simulatingresultsfitcloseinwiththetheoreticalresults． Keywords：combinedKdV equation；LBGK model；numericalsimulation 格子Boltzmann方法是近年来发展起来的用于模拟各种物理问题的数值方法．本文考虑如下组合 KdV方程的LBGK模拟： Ut+ 6uu 一 Ou。 + Ux 一 r． r． 0， (1) 其中00是某个常数．最近 3O年，格子Boltzmann方法(LBM)已经成功应用于流 体力学的各个领 域，如流体动力学，多相流，多组分流等[1]．近年来，这一方法也在模拟各种非线性发展方程方面取 得了很大的成功，如Burgers方程[4]，MKDV方程L6 等．组合KdV方程是一类重要的非线性物理方 程，本文建立模拟这类方程的一个带修正项的LBGK模型，数值模拟表明，该模型是十分有效的． 格子 Boltzmann方法 用小参数e来表示空间步长(e可以看作平均 自由程与特征长度的比值)．采用D1Q4模型使用的离 散速度方向 e1，e2，e3，e4]一 [1，一 1，2，一 2]． 设粒子的分布函数 (z，f)(=1，…，4)，定义其演化规律，即格子Boltzmann方程为 (z+ccei~t+e。’)一 (，￡)一一÷[(x,t)一。’(z，￡)]+ehi(，f)． (2) 收稿 日期：2007—1l—O9 基金项 目：国家自然科学基金资助项 目 (2004036133)；福建省科技计划资助项 目 (2008F5019) 作者简介：黄宝华 (1986一 )，女，福建泉州人，硕士研究生． *通讯作者：mac[@fj