P23
颗粒物料孔口流量的 Beverloo 标度:5/2 指数与截止参数对颗粒形状的稳健性检验
来源说明：本则出自独立撰写的第二批方案。它与前 20 则同样遵循八块结构与硬门槛要求，但撰写时未做英文文献检索，新颖性边界依据的是历届获奖图谱与既有知识，而非当轮查新。因此其「需核实」条目更多，启动前须自行补一轮英文检索。
1 · 研究问题
平底容器孔口直径 D 变化时，颗粒质量流量 W 是否服从 Beverloo 律 W = C ρ√g (D − kd)^(5/2)？指数 5/2 与截止参数 k 在从近球形(小米)到强非球形(大米、红扁豆)的颗粒形状序列上是否稳健？
2 · 研究背景与空白
与液体不同，颗粒从孔口流出的流量几乎不依赖料位高度——孔口上方形成自由落体拱，特征速度 ~√(gD)，由此量纲分析给出 W ∝ D^(5/2)。Beverloo(1961)引入有效孔径修正 (D − kd)(d 为颗粒直径,k≈1.4 对球形颗粒，需核实精确值：搜索关键词 "Beverloo correlation k value spheres")，该律是粮仓、制药下料设计的基石。
已有文献对球形颗粒(玻璃珠)验证充分；非球形颗粒的系统数据较少且分散,k 与指数是否漂移在教学与工程文献间说法不一。
空白在于：用同一装置、同一拟合口径横跨形状序列的对照数据很少，而这恰是低成本、高可重复、正负结论均成立的一年期课题：装置是一块钻孔板加电子秤，难点全在实验纪律与统计。
3 · 可检验假设
- H1：近球形小米的拟合指数与 2.5 的偏差 ≤ 0.1(95% 自助法置信区间覆盖 2.5)。
- H2：非球形颗粒(长径比 ≥ 1.5 的大米)的有效 k 显著大于球形颗粒(置信区间不重叠)；若指数本身漂移超过 0.15，则"Beverloo 形式失效"即为主结论。
4 · 量化验收标准
- 方法学校验(硬门槛)：计量链校准——电子秤用标准砝码(或整包定量食盐)校准误差 ≤ 1%，计时用手机视频帧号(30 fps,±1 帧)校准；随后用小米在单一孔径下重复 10 次，流量变异系数 ≤ 3%，证明装置可重复。此条不过关，后续全部结论无效。
- 孔径序列 ≥ 7 个点，覆盖 D/d = 5–30(下限避堵塞、上限避容器壁效应，容器直径 ≥ 5D)。
- 统计口径预先写死：主拟合为非线性最小二乘同时估计 (C, k, 指数)；另做固定指数 2.5 只估 (C, k) 的对照；两套都报，参数置信区间用 ≥1000 次自助法(以"次"为重抽单元)。
- 每个孔径重复 ≥ 5 次；流量取秤读数-时间曲线线性段斜率，线性段判据(R² ≥ 0.999 的最长窗)写死。
- 颗粒表征：每种颗粒随机取 ≥ 50 粒，游标卡尺测三轴尺寸，报等效径与长径比分布。
- 可复现性：数据表(每次放流的原始秤曲线)、拟合脚本开源。
5 · 数据与工具
| 用途 | 来源 / 工具 |
|---|---|
| 装置 | 大口塑料瓶/亚克力筒 + 可换钻孔底板(孔径 4–25 mm，五金店加工)+ 支架,<80 元 |
| 计量 | 厨房电子秤(0.1 g 分度,<40 元)+ 手机视频计时 |
| 颗粒 | 小米、大米、红扁豆、绿豆(超市),<20 元 |
| 分析 | Python:scipy.optimize.curve_fit + 自写自助法；纯 CPU，秒级 |
| 理论对照 | Beverloo 1961 关联式与球形颗粒 k 值,仅用于校验与对比，不计入本项目的数据贡献 |
6 · 方法路径
- 加工孔板序列，完成秤/计时校准与小米重复性验证(验收 1)。
- 表征全部颗粒的尺寸分布与堆积密度。
- 正式测量:4 颗粒 × ≥7 孔径 × ≥5 次，随机化测量顺序以摊薄湿度漂移。
- 提取流量(判据写死)，拟合两套模型 + 自助法置信区间。
- 核实边界效应：最大孔径下换更大容器复测 2 点，验证壁效应 < 2%。
- 独立交叉校验：用视频颗粒面计数法对小米复测 2 个孔径的流量，与秤法偏差 ≤ 5%。
- 汇总形状序列上的 (指数, k) 图，写双语术语表(orifice flow, free-fall arch)。
7 · 新颖性边界
- 本课题不声称发现 Beverloo 律(1961 年已发表)，不声称建立新的颗粒流理论。
- 已有工作：球形玻璃珠体系的指数与 k 已被多次验证；非球形颗粒有零散数据但口径不一。
- 本项目贡献(主结论)：同装置同口径下跨形状序列的 (指数, k) 稳健性数据 + 公开原始数据。属"参数维度补全 + 可复现性贡献"，定位须讲清，否则会被误读为验证性实验。
- 价值:"k 随长径比上升而指数不动"或"Beverloo 形式失效"都是可发表级信息；负结果同样成立，但须报告参数置信区间证明分辨力(能分辨 k 的 20% 差异)。
- 核对获奖清单：颗粒/沉降类相邻获奖有 "Boycott Effect"(2022 优胜，悬浮液沉降)与 "Falsification of Sand Painting on Induction Cooker"(2023 优胜，振动颗粒图案)——均不涉及孔口流；本赛道尚无 Beverloo 类获奖题。
8 · 决策门槛(go / no-go)
⚠️ 第 4 周末：装置重复性达标(变异系数 ≤ 3%)且最小孔径堵塞率 ≤ 30%。若小孔频繁堵塞，收窄孔径范围至 D/d ≥ 7 并增加大孔端点；若非球形颗粒出现强间歇流(秤曲线台阶化),降级：把间歇性本身升级为研究对象(流量涨落的功率谱与堵塞概率 vs D/d)，保留"形状如何改变孔口流"的主结论框架。 ⚠️ 第 8 周末：核实湿度影响——同批颗粒烘干前后流量差若 > 3%，全部实验加干燥剂储存流程并记录湿度。 ⚠️ 预算裁剪：颗粒种类可从 4 减到 2(球形 + 一种非球形)，主结论仍成立。