# 布朗运动与扩散模拟 - 随机行走、均方位移、斯托克斯–爱因斯坦、菲克定律

> 在线布朗运动与扩散模拟：成千上万个随机行走者、朗之万方程、均方位移斜率拟合、列维飞行与亚扩散，斯托克斯–爱因斯坦算扩散系数，菲克定律演示。

- 网址: https://qilingtools.com/zh/physics/random-walk
- 分类: 物理工具 (https://qilingtools.com/zh/physics)
- 标签: [物理](https://qilingtools.com/zh/tags/physics)
- 最后更新: 2026-10-29
- 价格: 免费，无需注册
- 隐私: 在浏览器本地处理，文件不上传服务器

同时模拟上千个从原点出发的行走者：格点随机行走、带阻尼与热噪声的朗之万布朗运动、列维飞行和连续时间随机行走，实时画出轨迹、双对数均方位移曲线（含斜率拟合与理论曲线）和位移分布直方图。斯托克斯–爱因斯坦计算器根据半径、黏度和温度给出花粉、细菌、蛋白质、糖分子乃至气体分子的扩散系数和典型位移；一维扩散方程页把数值解、行走者直方图和菲克通量放在一起对比。

## 使用步骤

1. 在「行走者与均方位移」页选择运动模型，调整阻尼 γ、温度 k_BT 或尾指数 α，点「运行」
2. 观察 log–log 均方位移曲线：正常扩散斜率为 1，朗之万短时为 2，列维飞行大于 1，亚扩散小于 1
3. 在「斯托克斯–爱因斯坦」页选择粒子与介质、设定温度，读出 D、1 秒/1 小时的位移和扩散到 1 mm 所需时间
4. 在「扩散方程与菲克定律」页选择初始分布和边界条件，拖动观测线 x₀ 查看通量 J 与累计穿过量

## 功能特点

- 格点行走、朗之万布朗运动、列维飞行、CTRW 四种模型，最多 5000 个行走者
- 双对数均方位移、最近十倍时间的斜率拟合，以及弹道–扩散过渡的精确解
- 位移直方图与同方差高斯对比，并给出超额峰度
- 斯托克斯–爱因斯坦（水、甘油、空气含坎宁安修正）与气体动理论扩散系数
- 一维扩散方程数值解、行走者直方图与点源解析解同图，反射/吸收边界
- 菲克第一定律通量曲线与可拖动观测点的累计净通量

## 常见问题

### 均方位移为什么和时间成正比？

每一步的位移彼此独立、平均为零，所以总位移的方差等于各步方差之和，随步数线性增长：⟨r²⟩ = 2dDt，d 是维数。这就是爱因斯坦 1905 年对布朗运动的解释，也是扩散的标志：距离只按 √t 增长，扩散 10 倍远需要 100 倍时间。所以在双对数坐标上正常扩散是一条斜率为 1 的直线。

### 朗之万模型里短时间为什么斜率是 2？

真实粒子有质量和速度，速度要经过动量弛豫时间 τ = m/γ 才会被周围分子的碰撞打乱。t ≪ τ 时粒子几乎沿直线运动，⟨r²⟩ ≈ (k_BT/m)·d·t²，斜率为 2（弹道运动）；t ≫ τ 时才变成扩散 ⟨r²⟩ = 2dDt。完整的结果是 ⟨r²⟩ = 2dD[t − τ(1 − e^{−t/τ})]，工具里的理论虚线就是它。对水中半径 1 µm 的颗粒 τ 只有约 0.2 µs，所以显微镜下看到的都是扩散段。

### 斯托克斯–爱因斯坦公式适用于哪些情况？

D = k_BT/(6πηr) 把粒子当作在连续黏性流体中缓慢运动的刚性小球，对胶体、细胞、病毒、蛋白质都很准，对蔗糖这样的小分子也大致可用；但当粒子与溶剂分子差不多大（如水分子自扩散）时会明显低估。空气中的亚微米颗粒要乘坎宁安滑移修正，因为颗粒尺寸接近空气分子平均自由程（约 65 nm）。气体分子之间的扩散则完全不适用斯托克斯定律，工具改用气体动理论 D = 3/(8nd²)·√(k_BT/πm)。

### 什么是反常扩散？

如果 ⟨r²⟩ ∝ t^β 且 β ≠ 1 就叫反常扩散。列维飞行的步长服从重尾分布 P(ℓ) ∝ ℓ^{−1−α}（α < 2），偶尔出现极长的跳跃，典型位移按 t^{1/α} 增长，属于超扩散，常用来描述动物觅食和湍流输运。连续时间随机行走里粒子每跳一次要等待重尾分布的时间（α < 1），越来越容易被长时间「困住」，⟨r²⟩ ∝ t^α，属于亚扩散，见于细胞质中的大分子和多孔介质。

### 菲克定律和随机行走是什么关系？

大量独立行走者的密度 c(x, t) 满足扩散方程 ∂c/∂t = D ∂²c/∂x²，这就是菲克第二定律；对应的粒子流 J = −D ∂c/∂x 是菲克第一定律：物质从浓度高处流向低处，流量与梯度成正比。工具同时用差分法解方程、用高斯步长移动几万个行走者，两者的直方图应当吻合，点源时还与解析解 c = M/√(4πDt)·e^{−x²/4Dt} 一致。想从概率角度理解随机过程，可以看马尔可夫链工具。

## 相关工具

- [单位换算](https://qilingtools.com/zh/math/unit): 免费在线单位换算：长度、面积、体积、重量、温度、时间、速度、数据存储、压强、能量、角度共 11 类，一次列出所有单位的换算结果，含市斤、亩、尺等中国传统单位。
- [向量计算](https://qilingtools.com/zh/math/vector): 免费在线向量计算器：二维/三维向量的加减与数乘、点积与叉积、模长、夹角、单位向量、投影长度与投影向量，并判断两向量是否垂直或平行。
- [有效数字计算器](https://qilingtools.com/zh/math/significant-figures): 免费在线有效数字计算器：数出有效数字位数并逐位标出，按 N 位有效数字或小数修约，支持四舍六入五成双等 6 种规则，给出科学计数法并按规则做加减乘除。
- [频率发生器](https://qilingtools.com/zh/audio/tone-generator): 免费在线音频频率发生器：1～20000 Hz 正弦、方波、锯齿、三角波，对数滑块与 0.1 Hz 微调，按音名设频率，线性/对数扫频，三种噪声，双耳节拍，可下载 WAV。
- [音频频谱图](https://qilingtools.com/zh/audio/spectrogram): 免费在线音频语谱图：上传音频或视频，用 Hann 窗短时傅里叶变换画出时间-频率-强度图，FFT 窗口与配色可调，点击读取时间、频率、音名和分贝，导出高清 PNG。
- [分贝测试仪](https://qilingtools.com/zh/audio/decibel): 免费在线分贝仪：用麦克风实时测量环境噪音，支持 A/C/Z 频率计权、快/慢时间计权、最小/最大/平均/Leq 统计、60 秒曲线和校准偏移，附常见声音分贝对照表。
- [物理公式计算器](https://qilingtools.com/zh/physics/physics): 在线物理公式求解：匀变速运动、抛体、牛顿定律、斜面、功和能、欧姆定律、串并联、RC、压强浮力、波、透镜成像、万有引力，任选已知量求未知量，带单位换算和代入步骤。
- [傅里叶级数与 FFT 可视化](https://qilingtools.com/zh/math/fourier): 免费在线傅里叶工具：方波、锯齿波、三角波、脉冲或自定义函数的 aₙ、bₙ 系数与部分和动画（吉布斯现象），手绘图形生成本轮圆，FFT 幅度谱与相位谱，支持加窗。
