用格子Boltzmann方法模拟非线性热传导方程.PDFVIP

　 ， ３６ １１ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１５年１１月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３６，Ｎｏ．１１，Ｎｏｖ．１５，２０１５ 文章编号：１０００⁃０８８７（２０１５）１１⁃１１５８⁃０９ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 用格子 Ｂｏｌｔｚｍａｎｎ方法模拟 非线性热传导方程∗ １ ２ 刘　 芳 ，　 施卫平 （１． 长春工业大学 基础科学学院，长春 １３００１２； ２． 吉林大学 数学学院，长春 １３００１２） 摘要：　 对具有非线性源项和非线性扩散项的热传导方程建立格子Ｂｏｌｔｚｍａｎｎ求解模型．在演化方 程中增加了两个关于源项分布函数的微分算子，对演化方程实施 Ｃｈａｐｍａｎ⁃Ｅｎｓｋｏｇ展开．通过对演 化方程的进一步改进，恢复出具有高阶截断误差的宏观方程．对不同参数选取下的非线性热传导方 程进行了数值模拟，数值解与精确解吻合得很好．该模型也可以用于同类型的其他偏微分方程的数 值计算中． 关　 键　 词：　 格子Ｂｏｌｔｚｍａｎｎ模型；　 非线性热传导方程；　 ＬＢＧＫ模型；　 非线性源项 中图分类号：　 Ｏ２４１．８　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．３８７９／ ｊ．ｉｓｓｎ．１０００⁃０８８７．２０１５．１１．００４ 引　 　 言 格子Ｂｏｌｔｚｍａｎｎ方法（ｌａｔｔｉｃｅ Ｂｏｌｔｚｍａｎｎ ｍｅｔｈｏｄ，简称ＬＢＭ）是２０世纪８０年代末期提出的 ［１⁃３］ 流体力学的一种新兴的数值计算方法 ．它由格子气自动机演化而来，是从分子运动论层次 描述基本物理规律的动力学方法．该方法的基本思想是构造一个简化的动力学模型（Ｂｏｌｔｚｍａｎｎ 方程），该模型包含了微观或介观过程的物理本质，使其统计平均结果满足所研究的宏观方程 （Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程）．由于ＬＢＭ 具有许多独特的优势，自提出之日起就受到流体力学、物理 学、数学、计算机技术等领域的众多专家学者的关注，并逐渐发展成为模拟复杂流动和多重物 理化学传递过程的强有力的数值工具，并成为国际研究的一大热点．近年来，一些学者根据 ＬＢＭ 的建模机理对原有模型进行改造（改变演化方程和平衡态分布函数的形式；对各阶动量 矩方程灵活限定等），构造出很多类型的偏微分方程的格子Ｂｏｌｔｚｍａｎｎ求解模型．具体工作可见 文献［４⁃１０］等． ［１１］ 本文主要研究用ＬＢＭ求解一般形式的一维非线性热传导方程 ，其形式如下： ２ ｎ ∂ｕ ∂ ｕ ｎ ＝ ＋ － 　 　 ａ ２ ｕ ｕ ，　 　 ｎ ＞ １． （１） ∂ｔ ∂ｘ ＝ 这里ｎ是有理数，ａ是常数，ｕ ｕ（ｘ，ｔ）表示温度，ｘ 和ｔ分别表示空间和时间变量．非线性热传 ∗ 收稿日期：　 ２０１５⁃０７⁃０１；修订日期：　 ２０１５⁃０８⁃２０ 基金项目：　 国家自然科学基金；吉林省教育厅“十二五”科学技术研究计划（２０１４⁃４７８） 作者简介：　 刘芳（１９８２—），女，辽宁铁岭人，讲师，博士（通讯作者． Ｅ⁃ｍａｉｌ：ｆａｉｒｆａｎｇ＠ｓｉｎａ．ｃｏｍ）； 施卫平（１９６３—），男，长春人，教授，博士（Ｅ⁃ｍａｉｌ：ｓｈｉｗｐ＠ｊｌｕ．ｅｄｕ．ｃｎ）．