热巧克力效应的定量检验：以同步测量气泡空隙率与共振频率判定 Wood 方程的适用边界

标签：优先级:中 · 桌面实验/声学-两相流族 · 预算档位:≤350元 · phyphox频谱 · 技能取向:信号处理+物理建模

1 · 研究问题

搅动含气液体后敲击杯底所闻音调随时间爬升（热巧克力效应，hot chocolate effect），其瞬时共振频率 f(t) 能否用独立测得的气泡体积分数（void fraction）φ(t) 代入 Wood 方程（Wood's equation，气液混合物声速 c(φ)）定量预测——预测误差随 φ 如何变化，方程在多大 φ、多大气泡尺寸下失效？

2 · 研究背景与空白

液柱敲击音的基频约为四分之一波长驻波 f=c_eff/4h（h 为液深）。少量气泡使混合物压缩性剧增而密度几乎不变，Wood (1930) 方程给出 c 随 φ 非线性骤降——φ∼1% 即可把声速从 1480 m/s 拉到 ~100 m/s，音调因此低几个八度，气泡消退后回升。Crawford（Am. J. Phys. 50, 1982）命名并解释了该效应；J. Sound Vib.（2012）用谱图可视化了频率爬升；多篇教学文献（UT Austin 演示论文等）给出定性对照。即：现象、机理与谱图演示均已发表。

空白在：已有工作几乎都只测 f(t)，把 φ 留作解释性参数；同步、独立地测量 φ(t)（由液面抬升高度换算）并与 f(t) 逐时刻对照 Wood 预测的定量闭环，在低成本装置上未见系统报道。这是"评价维度转换"型空白（从演示到定量检验），正负结论均有效。

3 · 可检验假设

• H1：用碳酸水稀释液配制 φ₀ = 0.2%–2% 的初始空隙率，敲击基频的实测值与"液面抬升法 φ + Wood 方程 + 四分之一波长"链条的预测值在 φ≤1% 区偏差 ≤ 15%。
• H2：预测偏差随气泡半径增大（静置粗化后期 vs 初期）系统增大——当气泡共振频率接近驻波频率时 Wood 低频极限失效，偏差与 φ 无关而与泡径相关（方向性假设，偏差–时间斜率显著为正）。

4 · 量化验收标准

1. 方法学校验（硬门槛）：无泡清水中同一杯的敲击基频换算声速，与 20 ℃ 纯水 1482 m/s 偏差 ≤ 5%（含杯壁弹性修正的讨论；若系统偏低，改用"频率比"口径：一切结论用 f_bubbly/f_clear 表述，消去几何与杯壁因子）；空气中管共振测声速与 343 m/s 偏差 ≤ 1%（phyphox 链路自检）。不过关则全部数据无效。
2. f(t) 提取：等间隔电磁敲击器（Arduino 驱动，1 次/秒，力度固定）+ phyphox/audacity 谱图峰追踪，峰判据（最大谱峰且带宽 < 50 Hz）预先写死；每条件重复 ≥ 8 杯次。
3. φ(t) 独立测量：微距镜头拍液面弯月面刻度，抬升分辨率 ≤ 0.2 mm（对应 φ∼0.2% 量级需液深 ≥ 10 cm 的窄口容器；此为装置设计约束，预先写死）；同步时间戳。
4. 拟合与统计：f–φ 关系与 Wood 曲线的偏差用 bootstrap 95% CI；H2 的斜率检验用 Theil–Sen + bootstrap（对非高斯残差稳健）；温度记录（声速温漂 ~3 m/s/℃）。
5. 误差源量化：敲击力度（电磁敲击器一致性 RSD < 5% 预校）、泡沿壁聚集（换杯重复）、溶解气持续析出（用放置 24 h 的脱气水做零对照）。"偏差与 Wood 预测不可分辨"即 H1 成立，须证明误差棒 ≤ 15% 判据的一半。
6. 音频原始文件与脚本全部开源。

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

用途	来源 / 工具	估价
高窄玻璃量筒/直筒杯（内径 ~40 mm，高 ~250 mm）×3	淘宝	≈80 元
电磁敲击器（电磁铁 + 继电器 + Arduino Nano）	淘宝	≈70 元
苏打水/碳酸水批量 + 量筒移液配比	超市	≈50 元
手机微距贴镜（液面读数）	淘宝	≈40 元
phyphox / Audacity 谱图、Python 峰追踪	免费	0 元
温度计、遮光背景板	自有	≈20 元
对照基准	Crawford 1982、Wood 方程、JSV 2012 谱图（仅校验对比，不计入贡献）	—

合计约 260–350 元。安全：无。

6 · 方法路径

1. 完成声学链路自检与清水基准（第 4.1 条），确定是否采用"频率比"口径。
2. 校准电磁敲击器一致性与谱峰提取脚本。
3. 标定液面法 φ 测量：注入已知体积空气（注射器微量）验证换算，偏差 ≤ 10%。
4. 主实验：4 档初始碳酸水配比 × 每档 ≥ 8 杯次，同步采 f(t) 与 φ(t)。
5. 与 Wood 预测逐时刻对照，量化偏差随 φ、随时间（泡径代理）的结构（H1、H2）。
6. 独立交叉校验：糖粉/可可粉成核方案重复 2 档（更接近原始"热巧克力"场景），检验结论不依赖气源方式。
7. 视频：谱图实时爬升与液面同框录制、装置讲解、敲击慢动作；成片 ≤ 3 min。

7 · 新颖性边界

本课题不声称发现热巧克力效应或 Wood 方程——Crawford (1982) 已命名并解释，谱图演示已发表（JSV 2012）。丘奖 2024 优胜 "Ups and Downs of Objects in Supersaturated Fluid"（重庆育才）研究过饱和液中气泡–固体耦合浮沉，与本题共享"析出气泡"要素但观测量（力学 vs 声学）不同；2024 优胜音乐盒扩音课题属结构声学，均须点名。本项目主结论是：低成本装置上 f–φ 闭环对 Wood 方程的定量检验及其失效边界。若预测精确，给出该演示实验可作"声学空隙率计"的精度评估；若系统偏离，偏差–泡径关系即指向共振修正，两种结果都有效。

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

• 第 4 周末：清水基准偏差 > 5% 且频率比口径仍不稳定（同杯 RSD > 3%）→ 换厚壁玻璃容器并加装杯底硬垫；2 周不过则整题止损转 P15（器材共通）。
• 第 10 周末：液面法 φ 分辨率达不到 0.2%（弯月面抖动）→ 降级路径 A：改用精密称重法（总质量不变、体积变化由溢流管读出）或把最低 φ 档删去，只在 0.5%–2% 区检验，主结论框架不变。
• 第 20 周末：H2 的泡径代理（静置时间）被证明与 φ 强耦合无法分离 → 降级路径 B：放弃 H2，加密 H1 的 φ 档数（4→7 档），把交付物聚焦为"Wood 方程偏差–φ 曲线"，框架不变。
• 预算裁剪顺序：先砍微距镜（用普通微距模式，分辨率降级入误差预算），再砍第三只量筒。

---