弹性波动方程的格子Boltzmann理论.pdfVIP

摘 要 格子 Boltzmann 方法 （LBM）源于二十世纪七十年代提出的格子气自动 机 （Lattice Gas Automata，LGA）方法，因此本文首先介绍了格子气自动 机方法的理论基础，LGA 研究的是一个完全离散的动力学系统：流场被离 散成为规则的格子，将流体考虑成离散的粒子，时间按固定的步长离散成 为时间序列，在这个系统中，流体粒子只可以沿网格线运动，粒子在一个 时间步内移动时，只能移动到最相邻的格点，而且粒子的速度是一个离散 的有限集合。 在格子 Boltzmann 模型中，用粒子的分布函数代替 LGA 中粒子本身进 行演化，其演化方程直接采用格子 Boltzmann 方程，并根据分布函数直接 计算流体的密度和速度。通过时-空离散，可以将复杂的 Boltzmann 微分、 积分方程表示成简单离散形式的格子 Boltzmann 方程。 在本文的第一章中详尽的介绍了有关于格子 Boltzmann 方法的基本思 想、基本模型和发展现状。 本文的第二章中首先介绍了弹性波动方程，然后将格子 Boltzmann 方 程应用于弹性波动方程。首先对格子 Boltzmann 方程进行 Taylor 展开，其    次 对 分 布 函 数 f (xt, ) 、 附 加 项 分 布 函 数  (xt, ) 以 及 时 间 t 做   Chapman-Enskog 展开和多尺度分析，然后根据宏观定义下的守恒条件，得 到一系列不同尺度下关于平衡态分布函数的偏微分方程。 在本文的第三章中根据这些平衡态分布函数的偏微分方程做动力学方 程的恢复，根据弹性波动方程的特点，并且在源项为二阶的假设下，通过 适当的定义高阶矩的表达式，便可得到修正了的弹性波动方程。进而通过 修正的波动方程就可以得到方程的截断误差以及误差的形式，并判断误差  eq 的精度，并且根据守恒条件可求得平衡态分布函数f (xt, ) 。  I 关键词： LGA， 格子 Boltzmann 方法，Chapman-Enskog 展开，弹性波动方程， II