「热巧克力效应」中气泡液体声速的时变规律：手机频谱测量与 Wood 方程检验

来源说明：本则出自独立撰写的第二批方案。它与前 20 则同样遵循八块结构与硬门槛要求，但撰写时未做英文文献检索，新颖性边界依据的是历届获奖图谱与既有知识，而非当轮查新。因此其「需核实」条目更多，启动前须自行补一轮英文检索。

1 · 研究问题

搅拌后含气泡的热液体中，敲击杯底音调随时间恢复的频率曲线 f(t)，能否由 Wood 方程(气泡液体声速公式)加上单一气泡消除时间常数的模型定量预测？液体温度与气泡引入方式如何改变该时间常数？

2 · 研究背景与空白

含微小气泡的液体声速远低于纯水：气泡提供压缩性、液体提供惯性,Wood 方程给出体积分数 φ~10⁻³ 量级的气泡即可把声速压低数倍。杯中驻波频率正比于声速，因此敲击音调是气泡含量的实时探针——这就是冲热巧克力/速溶咖啡时音调在几秒内爬升的原因(hot chocolate effect)。

Crawford(1982, Am. J. Phys.)对该效应做了经典的定性与半定量解释，报告音调爬升可达约 3 个八度，并指出 Wood 方程可解释量级。此后该现象常作为教学演示出现。

空白在于：该效应缺乏系统的参数扫描(温度 × 气泡引入方式 × 容器几何)与配套的公开数据集和不确定度预算。手机频谱 app 使逐秒追踪 f(t) 变得零成本，一年期学生课题恰好能补上这块——且"模型拟合好/不好"两种结果都是有效结论。

3 · 可检验假设

• H1:f(t) 恢复曲线可由"Wood 方程 + 气泡体积分数单指数衰减 φ(t)=φ₀e^(−t/τ)"模型拟合，决定系数 R² ≥ 0.95；时间常数 τ 随水温升高(50→90 °C,≥3 档)单调下降。
• H2：由频率反演的初始气泡体积分数 φ₀ 与独立几何测量(搅拌前后液面体积差)在同一数量级(比值落在 0.3–3 之间)；若系统偏离，偏离方向可由气泡尺寸分布假设解释(该解释须标注为推测)。

4 · 量化验收标准

1. 方法学校验(硬门槛)：静置纯水杯的敲击基频与一维驻波模型(f=c/4h，水深 h 实测)预测偏差 ≤ 5%，且用信号发生器 app 播放 440/1000 Hz 标准音校验手机频谱管线的频率读数误差 ≤ 1%。此条不过关，后续全部结论无效。
2. 每组实验条件重复 ≥ 10 次;f(t) 曲线提取判据(功率谱主峰，窗长 0.1 s，峰高阈值)预先写死并写入脚本，不允许逐条人工挑峰。
3. 参数扫描 ≥ 3 个温度 × ≥ 3 种气泡引入方式(搅拌/倾倒/吸管注气)，合计 ≥ 90 条有效曲线;τ 的拟合给自助法(bootstrap,≥1000 次重抽)95% 置信区间。
4. H1 的"单调下降"用 Spearman 秩相关 + 置换检验报 p 值，不用肉眼判断。
5. 物理自洽检查：反演 φ₀ ≤ 5×10⁻²(超出则说明模型失效区，需报告)。
6. 可复现性：全部原始音频、脚本与处理流程开源(GitHub/Gitee)，第三方可一键重跑出全部图表。

5 · 数据与工具

用途	来源 / 工具
声学采集	手机 + Phyphox(开源,phyphox.org)或任意录音 app，采样 44.1 kHz；纯手机即可
频谱分析	Python:numpy / scipy.signal(STFT)；全程纯 CPU，单条曲线处理 < 1 s
温度测量	学校实验室温度计或 <30 元电子测温计；冰水/沸水两点校准
理论对照	Wood 方程(标准两相介质声学公式，教材级);Crawford 1982 结果仅用于校验与对比，不计入本项目的数据贡献
器材	玻璃杯/量筒、吸管、可可粉或洗洁精微泡液，合计 <100 元

6 · 方法路径

1. 搭建测量链：手机固定支架 + 标准敲击器(定高落锤小球)，用标准音完成频率校准(验收 1)。
2. 跑通纯水基线：验证驻波模型、确定敲击可重复性(基频标准差/均值 ≤ 2%)。
3. 编写自动峰提取脚本，判据写死；用 3 条人工标注曲线验证脚本一致性。
4. 执行 3×3 参数扫描，每格 ≥10 次；当日备份原始音频。
5. 拟合 Wood + 指数衰减模型，输出 τ(T) 与自助法置信区间。
6. 独立交叉校验：用液面体积差反演 φ₀，与声学反演对照(H2)。
7. 整理开放数据包 + 双语术语表(bubbly liquid, Wood equation, standing wave)。

7 · 新颖性边界

• 本课题不声称发现新物理：效应机制与 Wood 方程解释已由 Crawford(1982)发表，本项目不得声称为自己的发现。
• 已有工作:Crawford 给出机制与量级估计；教学文献中多为演示性复现。
• 本项目贡献(主结论)：首个系统的温度 × 气泡引入方式参数扫描 + 时间常数 τ 的定量模型检验 + 公开可复现数据集。属"评价维度的转换 + 参数维度补全"。
• 价值：把一个演示级现象升级为带误差棒的定量检验；若模型在某参数区失效，失效边界本身就是结论。
• 核对获奖清单：声学类相邻获奖有 "Singing Tube"(2020 优胜，波纹管气流激励)、"An Investigation into the Sound Amplification in a Music Box"(2024 优胜)、"The Whistling Metal Plate"(2025 铜奖，结构振声)——均为发声机制研究；本题对象是声速作为气泡含量探针，测量量与物理问题不同。
• "τ 与温度无显著关系"同样是有效结论，但须给出 τ 的置信区间证明有能力分辨 2 倍以内的差异。

8 · 决策门槛(go / no-go)

⚠️ 第 4 周末：完成验收 1 校准 + 10 条热液体预实验曲线。可测阈值：敲击频移动态范围 ≥ 1 个八度、单条曲线可追踪时长 ≥ 5 s。若热液体气泡消除过快(<2 s)导致追踪失败,立即降级：改用室温水 + 吸管定量注气体系(气泡量可控、消除更慢),Wood 方程检验的主结论框架完全保留，仅损失"温度维度"。 ⚠️ 第 8 周末：核实项——手机麦克风在 200–4000 Hz 的响应是否足够平直(双 app 交叉验证)；若否，改报相对频移而非绝对频率，分析框架不变。 ⚠️ 预算裁剪顺序：容器几何维度可砍(不影响主结论)；温度与气泡方式两维度砍一保一仍可成文。

---