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

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

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

参数设置
300 K
3.0
1.0
80
显示设置
边界条件
预设
实时数值
0.0
扩散系数 D [px²/s]
0
⟨r²⟩ [px²]
0.0
时刻 t [s]
80
粒子数
正常扩散
状态
布朗运动动画
大粒子(示踪粒子) 溶剂分子 轨迹(残影)

“随机游走”=分子碰撞使粒子做不规则运动。增加粒子数可以看到初始团簇扩展开来,也就是扩散。

均方位移 MSD ⟨r²⟩ ∝ t
实测 ⟨r²⟩ 理论 4Dt
计算结果
0
⟨r²⟩ (px²)
0
D 估计值 (px²/s)
0
D 理论值 (px²/s)
理论与主要公式
$$\langle r^2 \rangle = 2dDt = 4Dt\ \ (2\text{D})$$

d=维数(2D),D=扩散系数,t=经过时间。由斜率得 D=⟨r²⟩/(4t)。

爱因斯坦-斯托克斯公式(越大/越冷/越粘→越慢):

$$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. 设置“温度 T”“粒子半径 r”“粘度 η”,扩散系数由爱因斯坦-斯托克斯公式 D=kT/(6πηr) 自动计算
  2. 改变“示踪粒子数”,观察大量粒子的扩散(初始团簇的扩展)
  3. 打开“溶剂分子”显示,可以看到撞击大粒子的高速分子
  4. 在 MSD 图中比较实测 ⟨r²⟩ 与理论线 4Dt,确认实时数值中的 D 估计值是否接近理论值
  5. 用预设(布朗运动/单粒子/扩散/低温·高粘度)可立即再现典型场景

具体计算示例

温度 T=300K、粘度 η=1.0、粒子半径 r=3.0 时,理论扩散系数与 D=kT/(6πηr) 成正比,在本工具中约为 60px²/s。二维随机游走中 ⟨r²⟩=4Dt,因此 t=5s 时期望 MSD≈1200px²。粒子数越多,由 MSD 估计的 D=⟨r²⟩/(4t) 越接近理论值 60px²/s。若要对应到真实物理系统,请另行定义像素长度和 1 秒所对应的物理量后再换算。

实务注意事项

  1. MSD 曲线尚未呈线性的初期阶段(过渡区)不宜用于估计扩散系数,必须等到积累足够的统计数据
  2. 粒子数少时单条轨迹的波动大,统计误差增大;粒子多时可靠性提高
  3. 显示使用归一化单位(px·s)。与物理系统对应时需通过量纲分析进行归一化
  4. 反射边界是理想的弹性反射。与实验对比时需另外评估壁面吸附、摩擦等外部效应