第三讲：随机行走到连续扩散.pdfVIP

第三讲：随机行走到连续扩散 Martin Z. Bazant Department of Mathematics, MIT 2006 年2 月13 日 概述 在前面的讲座中（由Yip教授演讲），我们已经讨论过单个粒子如何复杂的 移动，在液体或气体中粒子的可视随机运动被称为“扩散”。 在空气尘埃粒子 的所谓布朗运动中观察到同样由于碰撞的产生宏观尺寸上的无序运动。本讲座中 展示这个现象的模拟（java 程序），并可以在 /classes/109N/more_stuff/Applets/brownian/brownian.html 中找 到。这个例子说明，可以通过分子动力学在微观尺寸上模拟扩散现象，但是也可 以更简单、更粗糙描述粒子水平上模拟扩散，这就是本次讲座的主题。最简单的 “Monte Carlo ”（或概率）扩散模型是随机行走，这里假设粒子遵循一条独立随 机位移的轨迹。我们将证明，经过很长时间后，随机行走粒子集合的平均特性可 由扩散方程很好的描述出来，这个方程描述连续水平下的扩散，并且为下次讲座 的基础（由Rosales 教授演讲）。 历史 1827年Robert Brown第一次在显微镜下观察到花粉颗粒在水中、似乎一直持 续的无序运动。直到1877年Desaulx假定“布朗运动”的由流体分子的热碰撞引 起的运动。当然，不能应用分子动力学仿真去测试这些假设，但是一个基于随机 行走概念的数学理论即将出现。 “随即行走”最初由Karl Pearson在1905年提出1 。在给Nature 的一封信中， 1 为了引用和讨论有趣的历史，可见B. Hughes, Random Walks and Random Environmental, Vol. 22.00j/18.361j/... 建模与仿真概论 —— 第三讲 2 提出一个描述森林中蚊子横行的简单模型。在每一个时间步，一个蚊子移动一个 固定长度为a，角度任意选择。Pearson想要知道经过多个已采取的步骤后蚊子的 分布。 图1 对于不同的N值Pearson随机行走中P r N  的Rayleigh渐近逼近 这封信由Lord Rayleigh 回复，他早在1880年在不同材料中的声波背景下以更 普遍的形式解决了这个问题。模拟声波通过材料的传播可被认为是对一系列随机 波矢k （振幅相等但相位随机）的求和：材料中声波波长大致相等，但是它们在 材料中的散射点被改变。 我们希望找到声波经过多个已采取的步骤后的概率密度函数。设P R dR   表 N R dR 示N步中移动距离在R 和 之间的概率。对于单位长度的步骤，Lord Rayleigh 表示当N  时，  ~ 2RR N2 / （1） P R e N N 对于几个不同的N值，这个函数如图1所示。我们看到预计的传播距离与步数的 R N2 ~ ，这是典型的“扩散”现象。 平方根成比例，即 Ⅰ, Sec. 2.1 (Oxford, 1995)。 2 22.00j/18.361j/... 建模与仿真概论 —— 第三讲 3 大约在同一时间，随机行走理论也由Louis Bachelier在他真正非凡的博士论 文——1900年出版的La Théor