从离散速度模型到矩方法.pdf

从离散速度模型到矩方法 蔡振宁，樊玉伟，李若 要 为了求解动理学方程，我们通过研究在一维情形下对离散速度模 型的离散速度点的进行自适应的技术，发现可以自然地得到 Grad 矩 方程组。作为一个统一的认识，矩方程组可以看作是对离散速度点自 适应的离散速度模型，而离散速度模型可以看作是取特别形式的 “矩” 的矩方程组。这使得我们可以在一致的框架下来理解离散速度模型和 矩方法，而不是将它们对立起来。为了建立这样的一致框架，最近在 [3] 中发展的正则化理论是根本性的。 关键词：玻尔兹曼方程，离散速度模型，矩方法，双曲性，正则化 1 本文所要考虑的问题是所谓的动理学理论，具体地来说，是要给出对 下面的玻尔兹曼方程进行近似求解的方案。玻尔兹曼方程从介观层次上 来描述同一种类的大量微观粒子的运动，其基本变量是所谓的分布函数 ，其中 ，表示的是位于空间坐标、时 刻 具有速度 的微观粒子的数密度。所有的粒子只有两种基本的运动方 式：它要么遵循牛顿第一定律，以自身的速度 在空间中漂移，要么因为 受到力的作用而改变其运动的速度。粒子的受力具有很多的可能性，比如 最典型的情况就是和另外的粒子发生了相互碰撞，而本文将不太关心粒子 本文是应中国工程物理研究院第九届学术年会力学与数学学科分会的邀请，为了做年 会上的院外特邀报告而作。本文是三位作者的共同工作，最后由李若执笔完成。特别感谢 于翠影女士联系使得我可以使用TEX 而不是Word 来完成这个稿件，否则我会失去动工 的勇气。 北京大学数学科学学院，email: harecat@ 。 北京大学数学科学学院，email: ywfan@ 。 北京大学数学科学学院、应用物理和技术中心，email: rli@ 。 1 具体受力的因素和结果。在这样的情况下，玻尔兹曼方程描述了分布函数 随时间演化的机制 (1) 其中 就是碰撞项，碰撞项会导致粒子改变自己的速度，从而使得取 不同的 的 之间相互交换信息。 麦克斯韦通过干净优美的假设推导出下面的结论：如果让一个封闭的系 统中粒子充分地相互作用，直到状态不再发生改变，此时 ，这些 粒子的系统将会达到所谓热力学平衡态，而处于平衡态的分布函数将会形 为 ( ) exp ( ) (2) 1 其中的变量、 和 是宏观的流体变量，分别为密度、速度和温度 ，满 足下面的关系： ∫ d ∫ d (3) ∫ d 由于碰撞的作用会导致分布函数趋于平衡态，玻尔兹曼碰撞项有时候会使 用下面的BGK 碰撞项 来进行近似