方腔内带Soret效应和Dufour效应的双扩散对流的格子Boltzmann模拟-计算数学专业论文.docxVIP

1 绪 论 引言 流体流动是自然界中的常见现象。流体力学研究流体的运动规律，范围广泛。 经过多年发展，流体力学取得了丰硕的成果，但在很多方面面临着巨大挑战。流体 流动及其复杂，在经典流体力学领域还有不少问题没有解决。同时，流体力学不断 与其他学科交叉，产生了很多传统理论和方法不适用的新问题。 双扩散对流一般指流体中同时存在温度梯度和浓度梯度时，就会有温度扩散和 浓度扩散，并且相互作用，引起流体自然对流。双扩散对流是自然界和工业生产中 的常见现象。如大洋环流、软流层在地壳中的运动、晶体生长等。Soret 效应指温 度场不均匀导致的传质现象，Dufour 效应指浓度场不均匀导致的传热现象。Sr 可 取正值, 负值和零. 正值代表温度梯度的存在产生从高温到低温的传质, 负值代 表温度梯度的存在产生从低温到高温的传质。Dufour 效应类似。 随着计算机技术的蓬勃发展，数值方法逐步成为和理论分析、实验研究同等重 要的科学研究手段。对于流体流动的描述，从不同尺度可以分为：宏观连续模型、 微观分子模型和介观模型。依据所采用的流体模型，流动的数值模拟方法也可以分 为宏观方法、微观方法和介观方法。三种方法在一定条件下是等价的，但由于理论 出发点不同，有各自的适用范围。宏观方法基于连续介质假设，建立流体运动方程， 采用有限元、有限体积、有限差分、谱方法等对微分方程进行离散，得到相应的代 数方程组，再用标准数值方法求解。这是一种发展很成熟和成功的方法，形成了流 体力学的一个分支——计算流体力学（Computational Fluid Dynamics）。但是， 有些情况不能采用连续介质假设，只能采用微观方法和介观方法。微观方法适用数 值方法求解分子运动方程，并确定每个分子在各个时刻的速度和位置。它的假设条 件少，理论上应用范围不受限制，但受现有计算条件的限制时间尺度和空间尺度都 有限。介观方法，把流体离散成一系列的流体粒子。这些粒子比分子级别要大，但 在宏观上无限小。通过构造符合物理规律的演化机制，让流体粒子进行演化计算， 得出符合物理规律的结果。介观方法既具有微观方法假设条件少的优点，又有宏观 方法不关心分子运动细节的优势，从而在处理一些复杂流动问题上有巨大潜力。 格子 Boltzmann 方法（LBM）是近年来最为人们所关注的一种介观方法，在基 础理论和实际应用方面取得了很多成果。格子 Boltzmann 方法已经发展成为流体流 动模拟和物理建模的一种独特的和可信的数值方法。格子 Boltzmann 方法物理意义 清晰、边界处理容易、演化过程清晰、程序易于实现、并行性好等优点，是一种具 有广阔发展前景的数值方法。 格子Boltzmann方法概述 格子Boltzmann方法的发展历史 格子 Boltzmann 方法由格子气自动机（Lattice Gas Automata，LGA）发展而 来。LGA 是更广泛的元胞自动机（Celluar Automata, CA）在流体力学中的应用。 1972 年，法国的 Hardy、Pomeau 和 Pazzis[1]提出了历史上第一个 LGA 模型(命 名为 HPP 模型)。这个模型中，将流体视为由微观粒子所组成，同时时间和空间被 离散到二维正方形网格上。HPP 模型的状态演化分为两个阶段：碰撞过程和迁移过 程。HPP 模型微观上满足质量和动量守恒，但是由于正方形格子缺乏足够的对称性， 其宏观动力学方程不能恢复 Navier-Stokes 方程的非线性项，也不能恢复耗散功。 所以，HPP 模型开始未受到人们的重视，但是其基本思想开创了流体运动建模和模 拟的新思路。 1986 年，法国的 Frisch 和 Pomeau[2]提出了二维正六边形格子的 LGA 模型，即 FHP 模型。FHP 模型网格比 HPP 模型具有更好的对称性，能够恢复到二维不可压 Navier-Stokes 方程。同年，法国的 Humieres、Lallemand 和 Frisch 提出思维面心 一方的格子，能恢复到三维不可压 Navier-Stokes 方程。 LGA 方法有着鲜明的特色: 1、容易从底层刻画流体内部的相互作用；2、粒子 和壁面的相互作用可采用反弹格式，简单直观；3、满足质量、动量、能量守恒；4、 计算稳定性好；5、具有天然的并行性；6、其他。 同时，LGA 方法存在一些缺陷：1、统计噪声：由于碰撞算子含随即因素不可避 免的存在数值噪音。2、碰撞算子的指数复杂性：碰撞算子与离散方向数成指数关 系，对计算和存储有较高要求。3、不满足 Gallieo 不变性：由于对流项前有一个 非单位因子，推导出的动量方程不满足 Gallieo 不变性。 格子气自动机的发展孕育了格子 Boltzmann 方法。正是为了保留 LGA 的优势，