从排液与粗化到
泡沫的级联崩塌
From Drainage and Coarsening to Cascading Foam Collapse
Bubble-Cascade · 原始⑥ FROTH 方向的展开规划 · SIGGRAPH 2027 目标，按证据决定投稿
研究主线：追踪多泡系统中液膜变薄、Plateau 边界储液、气体扩散粗化，以及有限速破裂后的液体与动量传递；研究这些过程如何共同决定一次破裂会传播，还是被邻近泡沫吸收而终止。
保留原始 HTML 的完整目标：多泡 → 膜变薄与边界液体 → 逐级破裂 → 气体扩散粗化，并形成相互反馈。网络排液是基础模块，不再作为本项目的终点；粗化进入核心闭环，不移至一个未定义的远期阶段。表面活性剂是贯穿这条链的物理变量。
状态：仅有截面热力学、液量账本及工程测试。下文的新方法、机理假设、场景与验收指标均待实现和验证；尚无 FROTH 独立耦合泡沫结果。
01｜忠于原始规划：研究泡沫为何接连崩塌
总体规划 bubble_siggraph_plan.html 对⑥的原文为：“多泡、膜变薄、Plateau 边界里的液体、逐级破裂，以及气体扩散造成的粗化。”本版以此为边界，将此前过度收窄的“液膜—边界守恒耦合”恢复为完整研究问题。
| 原始组成 | 在新规划中的职责 | 必须看见的反馈 |
|---|---|---|
| 多泡 | 共享液膜、边界、顶点以及气体腔室的相互作用，不是若干独立气泡并排渲染 | 一个局部事件改变邻域压力、几何、流动与后续事件 |
| 膜变薄 | 排液、拉伸和界面输运产生空间非均匀厚度；需要时使用经验证的分离压 | 局部薄区、边缘夹缩及界面状态决定对扰动的响应 |
| Plateau 边界液体 | 有限储液、沿边界流动、膜—边界交换；破裂后吸收液体和冲量 | 边界既可能承受并传递扰动，也可能耗散或重新供液；结果不能预设 |
| 逐级破裂 | 区分首孔的形成与已有破裂的传播；回缩、冲击、邻膜响应及拓扑更新 | 邻膜事件不能由预先指定顺序或随机翻转概率制造 |
| 气体扩散粗化 | 缓慢改变泡体积、膜面积和邻接关系；传气受膜厚及液体结构影响 | 粗化改变后续排液与易损性，排液和事件又改变粗化通路 |
目标画面：同一泡簇先排液、形成非均匀厚度并缓慢粗化；受控的一次破裂随后引起回缩和边界响应，邻膜破裂或恢复，留下不同的崩塌范围。给出膜厚、边界含液量、气体量、通量和事件时间，而不只展示彩色泡沫。
完整项目不以赶完一个固定日期为科学完成标准。若核心物理链未成立，则延后独立论文；不能把“静态网络已能运行”重新包装成原始崩塌目标已经实现。
02｜更明确的创新：让慢演化与快破裂真正相互作用
“排液 + 破裂 + 粗化”均有先例。创新应来自新的离散耦合能力、对简化模型失效区间的预测，以及可复验的物理结论，而不是增加机制数量或使用新的缩写。以下三项是候选贡献，不是已证实的首次。
C1．有限速破裂与边界网络的守恒事件耦合
建立液膜—液环—边界—邻膜之间的状态转移接口，使破裂释放的液体、表面活性剂及动量具有明确去向。有限速回缩时允许边界响应、邻膜流动和几何重排同时发生。和“整张膜瞬时删除并均摊液体”的近似比较，测量后续事件时间、空间传播范围及吸收/终止情形。[R02][R03][R13]
要证明的难点：拓扑改变前后的自由度不同、多个输运路径并存，离散转移仍不重复计液体、不凭空产生冲量，并能承受多次事件。真实液环与飞滴由 CODA 提供；⑥负责泡沫网络、接口集成和反馈验证。没有该能力时不能把 A4 的二维压缩边缘带称为三维液环。
C2．排液—粗化留下的空间分布，能否改变级联的传播与终止？
设计具有相同邻接拓扑、近似相同几何、总含液量和平均膜厚，但膜厚—边界储液—表面活性剂空间关联不同的泡簇；由可记录的排液/粗化历史生成，再在同一位置施加同样的受控初始破裂。比较哪些局部液体通路造成传播，哪些造成抑制或恢复。
不额外发明“历史记忆变量”。历史只通过当前解析的物理场体现；完整状态相同就不应产生不同预测。不同初值导致不同结果本身不是发现，必须证明平均状态模型或既有近似遗漏了可量化的关键变量，并在未用于模型选择的工况中提高预测。[R26][R27]
C3．带误差控制的慢—快耦合，不让大步长跳过物理事件
平缓阶段以隐式步推进网络排液与气体交换；接近快速失稳或已出现开孔时缩小步长，用事件窗口推进回缩、边界响应和邻膜动力学，再返回慢阶段。慢阶段不能默认完全冻结几何；事件窗口也不能暂停全部排液与界面输运。
状态交接、步长切换及拓扑转移使用同一收支定义；用时间步细化和逐阶段误差比较选择切换条件。候选优势是相比全程统一小步长，在给定误差下节省成本，又相比完全尺度分离保留决定后续事件的瞬态反馈。变分、隐式或多速率这些词本身不是新意。[R02][R20]
- H1：至少存在一个实测材料和可实现几何区间，有限速事件与边界响应会改变后续事件，而不仅改变播放速度。若差异随加密消失，不能支持此假设。
- H2：在匹配粗粒度描述后，解析的空间液体/界面分布仍是预测传播所必需的信息。若平均状态模型已足够准确，应报告它成立的范围。
- H3：守恒慢—快算法能在相同预测精度下，比全程小步长更经济；若只有极大步长下的数值差异，则没有物理贡献。
推荐论文题目暂用 From Drainage and Coarsening to Cascading Foam Collapse。若最终证据只支持某一种机理，应在成稿时准确限定副标题；项目主线仍保留完整原始方向，不提前把结果写进标题。
03｜与先例逐项比较，不靠“首次”措辞
| 先例 | 已经具备 | 本项目争取的额外证据 |
|---|---|---|
| Saye–Sethian 2013/2016 [R01][R02] | 膜—边界—节点排液、动态边界截面积、运动中守恒搬运、瞬时破裂和重排 | 可动界面下的主动交换；事件窗口中有限速回缩与边界/邻膜反馈；和相同方程简化版作误差比较 |
| MELP 2022、thin-film SPH 2021 [R03][R05] | 粒子薄膜、形变、表面流、活性剂、多泡及拓扑相关能力 | 独立有限液量边界网络与事件传递；不能把粒子表示、黑斑图像和多泡当作独有贡献 |
| Ishida 2020、Tao 2024 [R04][R06] | 变形膜厚与非流形输运，或形变薄膜中的涡流和活性剂动力学 | 液环/边界具有液量和动量的多次事件闭环，以及经过验证的慢—快交接 |
| 局部膜—边界模型与实验 [R09][R10][R11][R12] | 过渡区通量、可动界面 CFD、膜供液和弯月面反馈 | 从局部闭合到运动泡簇传播；几何匹配验证，不把固体板弯月面直接当自由干边界 |
| Yanagisawa–Kurita 2019 [R13] | 回缩膜冲击、飞滴穿膜及含液率相关的集体崩塌 | 分辨并预测相应机制或明确只覆盖其中之一；提出实测可检验的终止条件 |
| 排液—粗化及剪切泡沫研究 [R17][R18][R26][R27] | 排液与粗化相互作用；重排/排液时标竞争维持非平衡厚膜 | 空间分辨地把缓慢准备阶段接到快速失效，而不是声称发现时标竞争 |
| 隐式 ALE 曲面方法 [R20] | 自演化曲面的单体隐式耦合 | 本混合维度网络和变拓扑事件特有的守恒、正性、稳定性与误差优势 |
不作为独立贡献：SI 单位、拉格朗日表示、动态边界含液量、守恒搬运、120° 平衡角、加入粗化、出现级联视频。它们分别是前人能力、必要验证或目标现象。边界是否“既能传递又能缓冲扰动”是待检验的机制解释，不能提前宣布为新定律。
文献依据来自本地已整理的 28 篇文献、24 份 PDF 及 1 篇全文 XML；这里只列与规划直接相关的条目。R10 仍只有出版社摘要；对其算法细节不作确定断言。未检得完全等价方法不构成全球首次证明。
04｜工程现状、模块边界和真实缺口
| 内容 | 当前证据 | 本计划的处理 |
|---|---|---|
| 截面热力学、液量账本 | foam/ 已有；规定转移和导数检查通过 | 复用。1000 次规定转移不是 1000 个物理排液步 |
| 多泡几何与表面动力学 | VIPER 双泡/三泡组件；尚无对应厚度—边界液量耦合 | 消费现有接口；新增通用泡簇构造和泡沫事件集成放 foam/ |
| 单膜、SI 材料与活性剂 | DRIFT 已有；排液实验误差、有限速率交换等仍开放 | 固定版本对照后复用；所需基础修正在③实现再合入，不复制私有热力学补丁 |
| 真实液环、碎裂、飞滴 | CODA 所属能力，当前⑥没有完整实现 | 必须设接口及交付关口；不能默认已可用 |
| 气体交换与破裂后气体账本 | ⑥未实现 | 本轮核心小泡簇必须覆盖；闭合并泡与向外界排气分别处理 |
| 空气惯性与阻尼 | GUST 所属；当前继承组件不含完整空气耦合 | 筛选不需此效应的验证区间；若实测快速响应依赖空气，则列为必要依赖，不能拟合阻尼替代 |
已有 125 项通过属于工程/组件验证；没有独立泡沫排液、真实级联或粗化闭环结果。本版仅调整研究规划和交接文字，不改变历史结果状态。
05｜模型：统一状态与收支，分别验证慢过程和快事件
A．变量与所有权
膜粒子带位置、速度、液体体积 V、全厚度 h = V/a、界面活性剂量及需要时的体相溶质量；边界结点带体积 VPB、截面 A = VPB/ℓ 和相应流动/界面变量；液环带液量、动量和活性剂；每个气泡带气体摩尔数 n。顶点先用无独立储液体积的压力结点；增加顶点液池时再补齐体积、能量与耗散。
S总 = S界面 + S体相（包括各区域）
n总 = Σn气泡；开放系统另记流入/流出与逸散气体
P液 = Σm u；气体/固体/约束作用通过冲量交换或边界收支记录
幽灵粒子不拥有液体。同一份质量不能既留在旧膜上，又进入液环或边界。开放系统的流出、气体逸散、支撑约束冲量和外力做功必须单独记录；有重力或固定框架时，不要求液体总动量恒定，而要求含外部冲量的平衡闭合。
B．能量、耗散与界面闭合
先推导界面自由能、气体自由能、分离势、重力和动能，避免在已有界面自由能之外再次叠加一份张力面积能。张力和化学势来自同一热力学约定；移动界面、吸附/解吸与膜—边界转移分别检查功与溶质收支。全膜厚/半膜厚、单界面 γ/双面膜张力 2γ 全程写清。
E相对/ℓ = γ(π − 2√3)r；∂(E/ℓ)/∂A = −γ/r
以上截面只适用于相应干泡沫、等张力近似。湿边界、显著不等张力或固体接触应换用对应截面/界面计算；不能用同一公式硬套所有场景。导通率的数值系数、界面流变和膜—边界过渡区闭合需有文献/局部求解依据，不能只用 A²/μ 的量纲表达。[R09][R10][R21][R23][R24]
冻结几何的耗散步可单独验能量下降；完整力学系统允许势能转为动能，检查总能量与外功/耗散的平衡。VIPER 现有分裂势存在能量误差；沿用时记录该误差，只称分步耦合。只有推导与验证成立后才称全局最小化或无条件稳定。
C．慢演化：排液、移动界面与粗化
膜内、沿边界及膜—边界交换通量写入同一个守恒更新。几何变化只更新面积和边长，液量由输运决定；重分布、插入和合并使用保守映射。界面活性剂先做不溶性、有限弹性基线，再在 DRIFT 接口验证后增加可溶动力学；不能将恒张力近似作为完整创新模型。
气体从高化学势向低化学势传递，闭合系统成对交换摩尔数；等温理想气体采用 pV = nRT，传气通量需与气体化学势/扩散阻力一致。膜内扩散、Henry 溶解和界面阻力的输入来源逐项列出；湿边界/顶点传气何时不可忽略，按 R18 评估。von Neumann–Mullins、MacPherson–Srolovitz 等只作相应极限下的诊断，不编码为驱动面积/体积的目标律。[R17][R18][R26]
D．首孔与传播分开，避免把阈值制造成机理
第一阶段用受控首孔。指定孔的位置、初始尺寸、注入/移除的液体和触发功，所有对照一致。它只回答已有破裂如何传播，不能预测自然寿命。后续若研究自发成核，再引入由材料、缺陷或涨落支持的模型，并独立报告敏感性；黑膜达到平衡不等于必然开孔。
液环以局部膜厚、张力和动量平衡自然演化，不直接指定 Taylor–Culick 速度。液环碰到边界时，通过有限时间交换/接触求解把液量、活性剂与冲量传入网络；局部碰撞/吸收闭合须有独立验证。不能把整膜液量瞬时均摊、将冲量清零后仍声称解析了破裂动力学。
区分三条路径：重排引起邻膜拉伸/再排液；回缩膜/液环对边界及邻膜的冲击；飞滴穿膜。前两者属于核心事件反馈；飞滴穿膜只有在 CODA 与接触模型可用时纳入。若不含飞滴，则限定研究覆盖的机理，不能拿 R13 的完整快速集体崩塌作已复现实验。
E．拓扑事件与慢—快交接
并泡时合并气体摩尔数，不按新体积重置压力或气量；向外界破裂时跟踪逸散通量。移除失效膜、合并/重连边界和 T1 时，建立旧自由度到新自由度的显式映射与局部收支。回缩期间膜仍有未回缩部分，不能先删除全部几何再“补画”液环。
慢—快切换依据相对状态变化、非线性求解误差和事件接近程度，并以逐步细化确定容差；不预设通用切换阈值。每次交接记录事件前后状态、能量/冲量、最小厚度和通量残差。允许渐近回到平衡后排液，但不能依赖该假设跳过决定级联的瞬态。
06｜开发路线：每一步都服务于完整闭环
| 阶段 | 实现与问题 | 交付物与退出标准 |
|---|---|---|
| M0 · 约定与依赖 | 复用已通过的截面/账本；记录 VIPER/DRIFT 基点、SI 与接点定义；明确 CODA/GUST 接口 | 变量表、导数核对、源版本清单；不重复搭环境，不重建分支 |
| M1 · 真实网络排液 | 单边界+邻膜 → 双泡/三泡；物理导通率、正性、开放/封闭收支 | 1000 个物理步、细网格参考和压力/通量历史；这是数值基线，非最终贡献 |
| M2 · 移动与界面状态 | SI 运动泡沫、主动交换、有限界面弹性；保守重分布；四泡拓扑构造 | 液量/活性剂/冲量闭合；分离被动搬运和主动排液；步长/分辨率比较 |
| M3 · 粗化准备阶段 | 在相同网络上推进气体交换和几何更新；二泡诊断、三/四泡联动 | 至少一个粗化改变排液/薄区分布的小泡簇；对流动、气体和时间尺度的账本独立闭合 |
| M4 · 受控首孔与有限速事件 | CODA 液环接入、边界吸收、邻膜响应、并泡及接点更新；T1 按场景需求实现 | 单事件完整状态转移；瞬时删除作消融；解析到的传播与终止机制明确 |
| M5 · 最小完整闭环 | 四泡/七泡排液与粗化 → 受控首孔 → 后续事件或恢复 → 剩余泡簇继续演化 | 至少一次后续事件与一次终止/恢复对照；无预设事件顺序；必要时如实报告假设不成立 |
| M6 · 验证与规模 | 局部实验、关键消融、稳健性和成本；27 泡作规模扩展 | 能复跑的证据包；27 泡视频不替代四泡/七泡的定量验证 |
粗化、界面输运和事件耦合均属于完整方案。排期上可在 M2 后先做 M3 的小算例，不必等大型泡簇完成。刚性/恒张力版本只用于调试和对照；若最终只有该版本，必须重新评估论文差异。
若 CODA 未能提供守恒的有限速液环模块，不把它标成可选而继续承诺真实级联；应当暂停这一核心结论，协调其所属分支或调整投稿时间。⑥不在自身目录复制另一份未经归属审查的液环/碎裂求解器。
07｜最小实验矩阵：证明反馈，而不是堆出机制
三类计算对照与物理消融分开
- 全程小步长参考：对同一模型、相同初值与边界逐步细化；检验慢—快切换的数值误差。
- 尺度分离对照：重排时被动守恒搬运，平衡后排液；保持可共同保持的材料与通量一致。若刚性界面是额外假设，另列为物理消融，不能混成算法误差。
- 瞬时事件对照：用明确的保守转移规则瞬时删除；和有限速事件比较，避免用不守恒的弱基线抬高优势。
在同一主场景预先指定“同步/分离”“有限速/瞬时”“动态粗化/固定气量”三因素的 2×2×2 组合；报告主效应与交互项。刚性/可动界面另做配对实验以控制工作量。固定液池版本只用于局部诊断，不充当 Saye–Sethian 的完整替代。
四组必要场景
| 场景 | 比较内容 | 观测量 |
|---|---|---|
| E1 单边界与邻膜 | 刚性/可动界面；有无邻膜供液；自由边界与实验证实的固体接触装置分别处理 | 速度剖面、截面、膜厚、压差和通量；验证导通率与交换闭合 [R09][R11][R12][R23][R24] |
| E2 双泡/四泡慢演化 | 移动中排液；有无气体交换；同方程时空收敛 | 膜厚/气体量历史、局部收支、界面状态、粗化影响下的薄区位置 [R18][R26][R27] |
| E3 同位点触发的四泡/七泡 | 保守瞬时删除与有限速回缩；解析或关闭某条物理传递路径 | 事件间隔、首个受影响邻膜、冲量与供液响应、最大受影响范围、传播或终止 [R13][R15] |
| E4 历史与空间分布对照 | 匹配粗粒度状态、改变实际空间关联；提前确定材料/几何扫描和保留工况 | 局部响应与简化模型误差、控制变量匹配误差、结果对分辨率/初扰动的敏感性 |
小泡簇先做机制和确定性响应，不用四/七泡的少数事件拟合雪崩幂律。只有在足够规模、可解释的初始条件/材料不确定性样本和统计收敛后，才讨论事件分布或概率；禁止临时添加随机破裂规则来凑分布。
验收目标：数值检查、物理验证和新结论分开
| 层级 | 目标/规则（均尚待验证） | 不能推导出的结论 |
|---|---|---|
| N0 热力学 | 解析导数/完整离散梯度与 Hessian 检查；初始截面 FD 相对误差目标 1e-8，做差分步长扫描 | 不代表可动/湿边界截面已经验证 |
| N1 固定拓扑液量 | 1000 物理步相对收支误差目标 ≤1e-13；打印绝对残差与总输运量，体积保持正值 | 不能靠负值裁剪或全局归一化达标 |
| N2 单边界 | 三档空间分辨率和至少三档步长；与同方程细网格场误差目标 ≤1% | 同方程参考不是实验验证 |
| N3 运动/事件/气体 | 每次重映射、事件和切换分别记液量、活性剂、气体及含外力冲量的残差；初始相对目标 ≤1e-10，绝对容差按 SI 尺度预先定义，报告累积误差 | 不能把气体逸散当成数值漏量，也不能把固定框架反力忽略 |
| N4 预测收敛 | 步长与粒距细化后膜厚历史目标变化 <2%；快速事件时刻、回缩速度和传播范围单独检验；离散事件顺序改变时扩大工况/集合分析 | 仅均值收敛不能支持事件顺序稳定 |
| P1 局部实验 | 先匹配几何、材料、界面流变、含液率和时间尺度；冻结参数后报告偏差及实验不确定性，不预设“都在5%以内” | 不同固体板/自由边界几何不能互换 |
| P2 快速事件 | 单孔回缩用相应极限的 Taylor–Culick 作诊断；邻膜冲击/吸收或飞滴机制与对应实验逐项核对 | 回缩速度正确不代表级联统计或自然寿命正确 |
| C 核心贡献 | 在保留工况和给定数值误差下，候选机制解释有效且相对合理基线改善预测；或明确给出简化模型成立范围 | 仅不同图像、极大步长差异、后验挑案例不支持新方法 |
拟合若确为必要，只允许独立局部校准，并保留未参与校准的材料或几何；将参数标为拟合值，不再使用“无拟合预测”的表述。临界厚度、界面黏性、边界吸收、传气率、空气阻力等均需来源和不确定性。
08｜论文图与可复现交付
- 主视觉：同一泡簇的排液、粗化、一次受控触发、传播或终止；同时显示真实膜厚与边界含液量，标注物理时间及触发方式。
- 方法：二维膜、一维边界与液环、顶点、气体腔室的自由度和事件转移；用收支式说明质量、溶质、气体及冲量的去向。
- 基础与成本：单边界参考、守恒/正性、误差随步长/分辨率、时间分项与内存；用误差—成本曲线比较多速率和统一小步长。
- 关键反馈：瞬时/有限速、同步/分离、动态粗化/固定气量的交互；不把多种物理同时更改后的一张差异图当消融。
- 空间历史：匹配均值但空间分布不同的泡簇、匹配残差与响应；给出简化模型是否失效及保留工况。
- 实验与适用范围：至少一项局部输运、至少一项有限速事件；列出未覆盖的空气、湿边界或飞滴区间，展示四泡/七泡定量案例，规模允许再加27泡。
新场图按粒子圆盘着色，不用网格插值冒充粒子解析。学术配色以 Nature 风格红蓝为主；新图片保存 PNG，矢量图同时导出 PDF。局部双击可用的 HTML 查看器嵌入数据或提供文件选择器。
每次研究运行记录提交、完整参数、单位、初始条件、边界、随机种子（若有）、日志/帧位置、耗时、迭代数与失败步；大于5MB的帧不进 Git。若训练任务超过2小时，设置每30分钟监视；本规划不启动训练或长计算。训练/拟合不是实现核心物理的前提。
09｜排期与决策：保留目标，按证据决定本轮能否成稿
2026-11-11 再核对 SIGGRAPH 2027 官方页面：表单与冲突信息 2027-01-21 22:00 UTC，论文截止 01-28 22:00 UTC，上传截止 01-29 22:00 UTC。细则尚待发布；内部仍以 01-10 冻结主要结果、01-25 定稿倒排。[S1]
| 时间 | 推进内容 | 必须交付/决策 |
|---|---|---|
| 11/11–11/22 | M0/M1：真实通量、单边界、双泡网络；并行做 CODA/GUST 接口可行性审查 | 能复跑的物理网络；是否具有有限速液环交付路径。后者缺失时不能承诺本轮真实级联 |
| 11/23–12/06 | M2/M3：运动泡膜、界面变量、四泡构造；二泡传气与小泡簇粗化 | 12/06：运动中的排液、气量交换和守恒重映射成立；至少一个可检验的事件反馈场景设计；否则重排期限 |
| 12/07–12/20 | M4：受控首孔、有限速事件、边界响应、并泡账本；局部实验匹配 | 12/20：至少一个单事件验证、传播/终止对照与可信基线；没有有限速反馈就不声称完整原始方向已完成 |
| 12/21–01/03 | M5：排液—粗化—事件闭环，四/七泡机制消融、空间历史工况 | 至少一个含所有核心环节的最小闭环；冻结主假设及保留工况，失败假设如实保留 |
| 01/04–01/10 | M6：收敛、独立复跑、局部实验、误差—成本和论文图 | 冻结可支持的结论；若闭环或物理验证未成立，延后完整论文，不以静态网络代替 |
| 01/11–01/20 | 共同作者审阅、论文与补充材料、版权与复现检查 | 01/18 内部登记材料；明确③/④/⑥贡献与署名/结果归属 |
| 01/21–01/28 | 登记、终稿和上传检查；不再加新机制 | 01/25 内部定稿；按官方正文及上传截止完成 |
资源假设：此时间表是高风险目标，不是工期承诺。有限速液环、泡沫拓扑、粗化和验证需要可并行推进的负责人；若仍由同一人同时承担③与⑥，不默认“30%投入足够完成两篇论文”。11/22 根据接口与成本测量明确人力安排；没有外部交付时延后完整闭环，保持研究目标不变。
三个研究关口
- 继续独立成稿：核心闭环可靠、候选贡献至少一项有清楚证据、局部实验匹配、对先例差异可量化。
- 假设未成立但模型可靠：报告哪些简化成立，继续寻找由预先定义扫描支持的问题；负结果不自动等于 SIGGRAPH 贡献。
- 仅网络或受控删膜可用：作为阶段成果保留并继续研究；不能称为“自然级联崩塌已复现”，也不自动改变③的研究范围。
10｜分支纪律、接口风险和第一阶段行动
继续在 06_FROTH/project 的 p6-codex-development 工作，只推送该分支，不切换③工作树。所有⑥新增求解、事件网络集成和泡簇构造放在 foam/。VIPER、DRIFT、CODA、GUST 的基础修改先在所属分支实现，审查后合入；共享 core/ 仅增量扩展。
- 先写膜/边界/液环/气体及开放外界的所有权和单位表，列出每次通量/事件的供体、受体和功/冲量。不要重做已完成的静态截面与规定转移测试。
- 从文献和局部截面求解核对刚性/可动边界导通率、膜—边界交换及不等张力适用范围；优先用可验证最小几何。
- 实现单边界加邻膜的真实时间推进、正性与收支；同时记录构建/求解成本并检查接口依赖。
- 写 CODA 事件接口：未回缩膜、液环状态、与边界接触、释放/接收冲量、是否支持飞滴。确认实际能力，不把计划接口当成已交付。
- 把二泡传气、四泡构造和匹配粗粒度状态的实验设计写成可复跑配置；不等大型渲染完成后再确定科学问题。
实现阶段按 Working rules 分步提交并推送，显式暂存路径；结果改变时同一提交同步 README、RESULTS 与结果页。每次运行前检查当前分支、工作区和远端。物理参数、模型选择和局限随输出保存。
11｜依据、文献与修订说明
原始范围来自 SPH/bubble_siggraph_plan.html 的⑥ FROTH 卡片；本次不修改总体 HTML 或其他分支。工程基线为⑥提交 2097fc6，状态见 开发验证记录。本版覆盖旧版过窄的论文范围与贡献定位，但不改写历史数值结果。
文献编号沿用本地创新性审查库，便于追溯；分支内仅保存公开引用及规划，不包含第三方 PDF 或依赖机器路径的外部文件。
- [S1] SIGGRAPH 2027 Technical Papers，官方日期于2026-11-11核对。
- [R01] Saye & Sethian, Multiscale Modeling of Membrane Rearrangement, Drainage, and Rupture in Evolving Foams, Science (2013)。
- [R02] Saye & Sethian, Multiscale modelling of evolving foams, JCP (2016)，重点§2.2–2.5。
- [R03] Deng et al., A Moving Eulerian-Lagrangian Particle Method for Thin Film and Foam Simulation, TOG (2022)。
- [R04] Ishida et al., A Model for Soap Film Dynamics with Evolving Thickness, TOG (2020)。
- [R05] Wang et al., Thin-Film Smoothed Particle Hydrodynamics Fluid, TOG (2021)。
- [R06] Tao et al., A Vortex Particle-on-Mesh Method for Soap Film Simulation, TOG (2024)。
- [R09] Breward & Howell, The drainage of a foam lamella, JFM (2002)。
- [R10] Anazadehsayed & Naser, A combined CFD simulation of Plateau borders including films and transitional areas of liquid foams (2017)；本次仅核对出版社摘要。
- [R11] Vigna-Brummer et al., Flowing menisci: coupled dynamics and liquid exchange with soap films, PRL (2025)。
- [R12] Vigna-Brummer et al., Quantifying the liquid flow between a soap film and a vertical meniscus, 2026 预印本，非已确认期刊最终稿。
- [R13] Yanagisawa & Kurita, In-situ observation of collective bubble collapse dynamics in a quasi-two-dimensional foam, Scientific Reports (2019)。
- [R15] Ritacco, Kiefer & Langevin, Lifetime of Bubble Rafts: Cooperativity and Avalanches, PRL (2007)；本次仅核对摘要。
- [R17] Morgan & Cox, Mean-field model for the bubble size distribution in coarsening wet foams, Soft Matter (2026)。
- [R18] Schimming & Durian, Border-Crossing Model for the Diffusive Coarsening of Two-Dimensional and Quasi-Two-Dimensional Wet Foams (2017)。
- [R20] Sauer, A curvilinear surface ALE formulation for self-evolving Navier–Stokes manifolds (2025)。
- [R21] Koehler, Hilgenfeldt & Stone, Foam drainage on the microscale I (2004)；本次仅核对摘要与机构书目。
- [R23] Koehler et al., Drainage of single Plateau borders: Direct observation of rigid and mobile interfaces, PRE (2002)。
- [R24] Koehler et al., Foam drainage on the microscale II: Imaging flow through single Plateau borders (2004)。
- [R26] Hilgenfeldt, Koehler & Stone, Dynamics of Coarsening Foams: Accelerated and Self-Limiting Drainage, PRL (2001)。
- [R27] Saint-Jalmes & Trégouët, Foam coarsening under a steady shear: interplay between bubble rearrangement and film thinning dynamics, Soft Matter (2023)，重点§4.1与Fig.8。
本次修订：恢复原始完整目标；将基础网络与候选创新分开；加入有限速事件、空间分布历史和慢—快误差控制；把粗化纳入最小闭环；增加机制消融、终止对照、依赖关口和失败假设的处理；不声称新结果已经出现。