流体热力学的格子Boltzmann方法的理论分析.pdfVIP

( ) 第 15 卷　第 3 期 广西师范大学学报 自然科学版 1. 15　 . 3　 V o N o 　1997 年 9 月 JOU RNAL O F GUAN GX INORM AL UN IV ER S IT Y Sep tem ber　1997　 　 流体热力学的格子Bo ltzm ann 方法的理论分析 刘慕仁　　何　云　　孔令江 (广西师范大学物理与电子科学系, 541001, 桂林; 第一作者 33 岁, 男, 教授) 摘　要　对 F J A lexander 等人提出的流体热力学的格子Bo ltzm ann 方法, 进行了研究, 对模型相 对应的宏观流体热力学方程组给出了理论分析的正确结果和清晰的推导过程 关键词　流体热力学; 格子Bo ltzm ann 方法; 理论分析 分类号　 414 2 O 近几年来的研究表明, 格子 方法成功地模拟了流体力学的 方 Bo ltzm ann N arier Stokes [ 1 ] [ 2 ] 程 但作为沾滞流体, 描述它运动的完备方程组中应含有流体的热力学方程和物态方程 [ 3 ] 人们曾尝试用格子气自动机方法建立这样的方程 , 但由于格子气 自动机方法中, 粒子的分布 函数是 分布函数, 因此不能给出正确的流体力学方程 直到最近, Ferm i D irac F J A lexander [ 4 ] 等人才给出了流体热力学方程 , 为完整地求解流体问题提供了一条新途径 但文献[4 ]在理论分析中存在如下问题: ①导出的方程式中, 多处出现错误; ②推导过程不 明确 鉴于此, 本文重新作了推导, 给出了正确的方程式和清晰的推导过程 1　流体热力学的二维格子Bo ltzm ann 模型 该模型由两个六方格子迭加而成 小格子的边长为 1 个单位, 大格子的边长为 2 个单位, 用单粒子分布函数f (r , ei , t) 表示 r 处, t 时刻具有速度 ei 的粒子出现的几率, = 0, 1, 2 分别 表示大小为 0, 1, 2 三种速度, i= 1, 2, …, 6 表示 6 个可能的速度方向, 粒子的质量为 1 个单位 为方便,f (r , ei , t) 也可简记为f i (r , t) 或者f i , 按统计物理, 它满足如下格子Bo ltzm ann 方 程: 　　 国家自然科学基金及广西自然科学基金资助课题 　　收稿 日期　1997 03 06 © 1994-2008 China Academic Journal Electronic Publishing House. All rights reserved.