布朗运动模拟器 返回
统计物理·扩散现象

布朗运动·随机游走模拟器

自由调整粒子数、步长、轨迹显示,实时观察扩散过程。验证爱因斯坦的MSD=2dDt公式,估算扩散系数D。

参数设置
粒子数
步长 σ
每帧步数
显示设置
轨迹长度
颜色模式
边界条件
预设
计算结果
0
MSD (px²)
0
D估计值
0
经过步数
粒子轨迹

模拟自动开始。改变预设或参数来观察扩散过程。

理论·主要公式
$$\langle r^2 \rangle = 2dDt$$

d=维数(2D),D=扩散系数,t=经过步数

爱因斯坦-斯托克斯式:

$$D = \frac{k_B T}{6\pi\eta r}$$

布朗运动·随机游走概述

🙋
布朗运动?屏幕上的粒子只是随意动来动去而已,有什么特别的吗?
🎓
这就是妙处所在啊!这个"随意"的运动实际上证明了分子的存在。1905年,爱因斯坦用数学理论解释了这种现象——看不见的水分子在轰击花粉,导致它随机运动。用这个模拟器,把"粒子数"增加到多个,观察它们平均行为的规律性,就会明白了。
🙋
规律性?随机的动作怎么会有规律呢?屏幕右边的图标记着"MSD",这是什么?
🎓
完全正确!MSD(平均二乘位移)就是粒子离开初始位置有多远的平方值的平均。看起来随机,但MSD值会随时间均匀增长——这就是规律!如果你把"步长 σ"滑块调大,一次"被推的距离"就变大,MSD图表的斜率(代表扩散系数D)就会增陡。试试看!
🙋
哦!那这个"扩散系数D"在真实世界中是怎么由什么决定的呢?
🎓
绝妙的提问!在现实中,它由爱因斯坦-斯托克斯方程决定。比如一种药物粒子在血液中扩散,就用这个概念。试试把"边界条件"改成"反射"——就能模拟粒子被容器限制的情况。在工程中,这些理论用来预测细胞内分子的运动,半导体掺杂扩散,污染物扩散……基础就是这个模型。

常见问题

MSD是平均二乘位移(Mean Squared Displacement)的缩写,表示粒子从初始位置扩散的程度。遵循爱因斯坦关系⟨r²⟩=2dDt,随时间线性增长。从图表的斜率可以推算出扩散系数D。
增加粒子数能提高统计精度,MSD曲线的抖动减少。增大步长会让单次移动距离变大,扩散加快。但步长过大会破坏随机游走的前提,需要在合理范围内调整。
观察MSD图表中线性段的斜率。对于二维(d=2),斜率等于4D,因此斜率÷4就是D。例如斜率为0.8,则D=0.2。粒子数越多推算越准确。
基于理想随机游走模型,定性上能再现爱因斯坦关系。但真实情况受粒子大小、流体粘性、温度等影响,需要用爱因斯坦-斯托克斯方程调整参数才能定量相符。

现实应用

药物开发·靶向输送:纳米级药物粒子在血液中如何扩散、到达患处的速度,都可用扩散理论预测。粒子大小与扩散速度的关联直接影响用药设计。

半导体工艺:硅晶圆中杂质(掺杂元素)的热扩散过程需要精确控制,使用扩散方程来预测杂质分布,其基础正是布朗运动理论。

环境科学:大气微粒(PM2.5)和水体污染物的扩散都遵循这一原理。复杂环境中会叠加湍流影响,但出发点都是随机扩散的概念。

金融工程:股票价格变动常用几何布朗运动建模,期权定价的黑-舒尔茨模型就基于随机游走理论,用来评估金融风险。

常见误解与注意事项

CAE初学者用这个模拟器时容易踩几个坑。第一个就是想当然地认为"粒子越多越精准"。确实统计抖动会减小,但计算负荷激增。实际工程中要权衡精度和成本——概览扩散系数可能只需100个粒子,但要精确得到分布形状可能需要1000+。这种判断能力就是专业性。

第二个常见错误是混淆"步长σ和扩散系数D的因果关系"。在模拟器里σ大则D大,但真实世界中D由物质属性(温度、粘性、粒子大小)决定,是"结果"。要模拟真实的D,需要反推合适的σ和时间步Δt。公式是$$D = \frac{\sigma^2}{2 \Delta t}$$。理解这个"建模"环节很关键,否则就只是玩游戏。

第三个坑是"反射边界的陷阱"。用反射条件能模拟容器,很方便,但现实中分子碰壁不是完全弹性反射。靠近壁面会产生摩擦甚至化学吸附。在模拟中用反射只是"第一近似"。如果要模仿细胞这样复杂的环境,可能还要给边界加特殊规则(比如某个概率被吸收),这样模型才更贴近实际。

使用指南

  1. 设置粒子数在1~100范围,从10个开始较为合适
  2. 调整步长(sigmaVal)在0.5~5.0像素范围。分子运动模拟建议1.0~2.0
  3. 设置采样频率(spfVal),决定MSD计算的时间分辨率。通常每步都采样
  4. 调整轨迹长度(trailVal),显示过去n步的轨迹。为节省内存,建议50~200
  5. 启动模拟后,监测实时MSD值和D估算值

具体计算示例

设粒子数n=20、步长σ=1.5px、总步数1000,初期(步0~100)MSD≈150px²,中期(步200~500)MSD≈850px²递增。用MSD=4Dt公式,t=500时MSD=850,反推D≈0.425px²/step。现实物理中,水中红细胞(直径7μm)标准D≈100μm²/s。若设模拟时间标度为0.1ms/step,就能转换到物理值。

工程实务注意

  1. 初期MSD曲线(0~50步)线性性差,不宜用于D推算。需等待统计数据足够
  2. 粒子少于5个时个别轨迹波动大,统计误差大。建议n≥15提高可信度
  3. 步长和采样频率组合会改变表观D值。与物理系统对应时,需用量纲分析归一化
  4. 无壁面摩擦或漂移力,仅模拟纯随机游走。与实验对比需另评估外场影响