颗粒孔口流的 Beverloo 标度与堵塞统计：以连续称重判定 5/2 指数与雪崩分布形态

标签：优先级:高 · 桌面实验/颗粒物质族 · 预算档位:≤600元 · Arduino称重传感器 · 技能取向:自动化采集+统计

1 · 研究问题

自制竖直料仓中颗粒经圆孔的稳态流量 Q 是否满足 Beverloo 定律（Beverloo law）Q = Cρ√g (D-kd)^5/2，拟合指数的 95% CI 是否含 2.5？进一步：当孔径比 D/d 降至 6 以下进入堵塞区时，两次堵塞之间的"雪崩规模"（avalanche size）分布是否为指数型，其特征规模随 D/d 如何增长？

2 · 研究背景与空白

颗粒流量与孔径的 5/2 幂关系（Beverloo 1961）是颗粒物质（granular matter）的基石经验律；沙漏计时与料仓设计都依赖"流量与料高无关"这一反直觉事实。Mankoc 等（2007）给出跨大范围孔径比的修正表达式；Zuriguel 等系统测量了堵塞概率与雪崩分布，发现雪崩规模呈指数分布、且"是否存在临界孔径比"仍有争议（To 等 vs Zuriguel 等）；2024–2025 年仍有 DEM 论文在验证 Beverloo 适用界（arXiv:2512.03698 指出 D/d≈10 时流量还随料高变化）。即：定律本身与堵塞统计均为已发表结果。

空白在：以上实验多用精密球形颗粒与实验室天平；用一套 ≤600 元装置在同一实验中同时覆盖自由流区（标度拟合）与堵塞区（雪崩统计），并检验非球形厨房颗粒（小米/绿豆/玻璃珠对照）对指数与分布形态的影响，未见系统报道。这是"换样本 + 参数维度"型空白，正负结论均可发表级陈述。

3 · 可检验假设

• H1：对近球形玻璃珠，Q–(D-kd) 双对数拟合指数为 2.50 ± 0.15（95% CI），且固定孔径下流量随料高（3 倍变化）漂移 < 3%。
• H2：D/d < 6 区雪崩规模 s 的生存函数在半对数图上为直线（指数分布，K–S 检验 p > 0.05），非球形颗粒（绿豆）的特征规模在同一 D/d 下显著小于球形颗粒。

4 · 量化验收标准

1. 方法学校验（硬门槛）：HX711 + 称重传感器用标准砝码（借学校物理组）五点标定，线性偏差 ≤ 0.5%；再复现"流量与料高无关"（D/d≈8，料高 3 档，流量差 ≤ 3%）。不过关则全部数据无效。
2. 采样率 ≥ 20 Hz 连续记录 m(t)，稳态流量取中段线性拟合；每个孔径重复 ≥ 10 次装填，报均值 ± 标准误、RSD。
3. 指数拟合：ODR + Theil–Sen 交叉验证 + bootstrap 95% CI；k 与指数联合拟合并报告参数相关性；孔径 ≥ 7 档、跨度 ≥ 1 个数量级的 (D-kd)。
4. 雪崩统计：每个 D/d 采集 ≥ 200 次堵塞事件；分布形态判定预先写死为"指数 vs 幂律"的对数似然比 + K–S 检验（α = 0.05）。
5. 误差源量化：环境湿度（颗粒吸潮改变摩擦，RH 记录且 > 70% 天次作废）、装填方式（漏斗定高倾倒，写成 SOP）、传感器温漂（每轮前后零点核查）。
6. 脚本 + 原始 m(t) 数据全部开源。

5 · 数据与工具（含器材清单与估价）

用途	来源 / 工具	估价
透明亚克力管（内径 80 mm，长 400 mm）+ 底板支架	淘宝定制	≈150 元
孔板 ≥ 8 片（3D 打印或亚克力激光切孔 3–30 mm）	学校创客室/淘宝代打	≈100 元
玻璃珠（d≈2 mm 与 4 mm 各 2 kg）、小米、绿豆	淘宝/超市	≈120 元
HX711 + 5 kg 称重传感器 + Arduino Nano	淘宝	≈50 元
DHT22 湿度记录、振动排堵用小电机（堵塞后重启流动）	淘宝	≈40 元
手机 240 fps 慢动作（孔口流态定性记录；局限：无法解析单颗粒轨迹，仅作辅助）	自有	0 元
对照基准	Beverloo 1961、Mankoc 2007 修正式、Zuriguel 雪崩统计（仅校验对比，不计入贡献）	—

合计约 450–600 元。安全：无。

6 · 方法路径

1. 搭料仓 + 完成称重标定与"料高无关性"复现（第 4.1 条）。
2. 写 Arduino/Python 自动采集与堵塞自动判定（流量为零持续 > 5 s 判堵塞，判据写死）。
3. 玻璃珠扫孔径 7 档以上，ODR 拟合 Beverloo 指数与 k。
4. 进入堵塞区，逐 D/d 采 ≥ 200 次雪崩，做分布形态检验。
5. 换绿豆/小米重复第 3–4 步，对照颗粒形状的影响。
6. 独立交叉校验：用秒表 + 定量杯粗测两档流量，与传感器结果比对（系统误差核查）。
7. 视频：装置全景、孔口 240 fps 流态、堵塞拱结构特写、数据屏同框；成片 ≤ 3 min。

7 · 新颖性边界

本课题不声称发现 Beverloo 定律或雪崩指数分布——均已发表（Beverloo 1961；Mankoc 2007；Zuriguel 等）。丘奖历届有颗粒相关获奖（2023 优胜 "Falsification of Sand Painting on Induction Cooker" 为振动颗粒成图；2024 金奖为听声计球），均不涉及孔口流标度。本项目主结论是：(i) 低成本连续称重装置对 5/2 指数的独立复测与不确定度预算；(ii) 厨房级非球形颗粒对指数、截距 k 与雪崩特征规模的形状效应定量对照。指数不漂移则普适性获得独立支持，漂移则量化形状依赖，两种结果都有效。"临界孔径比是否存在"的学术争议只作背景，不声称裁决。

8 · 决策门槛（go / no-go）

• 第 4 周末：称重链标定或料高无关性不过 → 检查支架共振与传感器量程；2 周不过改用"接液瓶置于电子秤 + 手机拍秤面读数"的备胎记录法（时间分辨率降至 1 Hz，只影响短时涨落分析）。
• 第 9 周末：单孔径流量 RSD > 10% → 降级路径 A：孔径档数减至 5、每档重复加倍，只做 H1，放弃雪崩统计；主结论框架（标度独立复测）不变。
• 第 18 周末：堵塞事件采集速率 < 50 次/日（人工重启太慢）→ 启用振动电机自动排堵；仍不足则雪崩统计只做 2 个 D/d 档，分布形态结论降为"单点检验"。
• 预算裁剪顺序：先砍第二种玻璃珠粒径，再砍振动电机；砍后 H2 的形状对照缩为"玻璃珠 vs 绿豆"一组，框架不变。

---