NET ACCUMULATION / JOINT HENRY–MARANGONI / CONTROLLED MOVING GEOMETRY

液体流入后，
边界体积真正增长。

10组正向算例共270步，另有1步故意漏掉网格运动的负对照。边界尺寸由实际收到的液量决定；浓度、Marangoni流动和三张有限膜的库存进入同一次求解。13项预定检查及5项独立保存数据审计通过。

真实保存状态：20秒联合供液

节点位置随实际边界尺度变化，圆点颜色表示c/c*。坐标用同一尺度绘制，帧间不改变视场；没有插值。0.925%的尺度变化本来就很小，数值读数与原始坐标可直接核对。

矢量PDF · 11组矩阵与13项检查 · 5项独立审计与空间误差 · 材料、模型输入和源码版本 · 实际保存帧

这次解决了什么

过去固定边界要求三口液量代数和为零。这次允许三张膜同时供液：ΔA=−ΔtΣQ，Q是相对运动入口的外向液体流量。边界随净来液增大；绝对不可压速度场同时包含气液界面的法向运动，不再把净来液塞进固定体积。

体相控制体容量随s²变化，Henry界面容量随s变化。膜的远端框架固定，近端随边界移动，因此膜的投影长度、液量、表活量与入口浓度都更新。体相、界面和三个接触口分别上风输运，两个方向的通量不会先抵消。所有新浓度与流函数同时求解，没有摩尔数归一化或负值裁剪。

新实现退回平衡通流、固定几何时，与上一阶段保存结果的浓度扰动范数差为1.84×10⁻¹²。体相均匀对照在半径增长4.54%的应力测试中仍保持均匀至约10⁻¹³；故意漏掉网格通量的算例，总摩尔数仍在舍入误差内守恒，但局部浓度已经错误。这两个检查不能互相替代。

能量：保留几何和约束做功

容量变化与相对流动压缩项可以有正有负。另从原生速度算子检查黏性／SIP耗散、Marangoni功和运动边界反力功；两式结合时，几何与接触约束的剩余功单独列出。这些量直接由容量增量和实际通量计算，并未用能量残差反推修正。

这里尚未得到自由法向形状运动、洁净表面张力能及接触区机械功的完整闭合。因此可以验证本离散参考的收支一致，不能称作完整自由泡沫的能量单调定律。能量增加也不能自动解释成数值不稳定。

实际执行的对照

GCL两组只验证体相几何守恒，不代表完整肥皂膜物理。negative-no-mesh是故意错误的消融，预期未通过局部几何检查；它不混入正向算例。主算例采用q0=8×10⁻¹³m²/s；增长诊断用10q0运行20s。所有数值是单位轴向长度的二维截面量，液体为m²、表活量为mol/m。

独立复核、误差与适用范围

从保存的实际三角形坐标重算边界面积，与三张膜失去的液量核对；从V、L和c重建每张膜的表活量，再逐口对上累计传输。原生BDM速度还独立积分了相对运动壁面／入口的法向流量，以及−∇c·u的Marangoni功。入口误差约10⁻¹⁴，梯度积分与通量伴随功的相对差约10⁻¹³。

最后两档时间步的浓度扰动场差约0.345%；两个网格的共同内部3491点重构场差6.42%，尚未空间收敛。主／细网格SIP数值面项占比约29.6%／24.6%；20s粗网格增长诊断约48%，不能全部称为物理黏性。

入口口径随s缩放，膜厚由V/L得到，两者仍有差别：主算例最大约0.00291%，长时联合诊断约2.88%；体相大幅运动测试约13.8%。长时／大幅变化结果用于验证收支和接口，不用于量化真实泡沫寿命。稀Henry材料也未与具体实验标定。

启动批次在体相应力算例的Newton残差1.54×10⁻¹¹处停滞，未达到更严格的内部尝试值10⁻¹¹。保留失败输出和源码ZIP；现在仅在已经满足预定10⁻¹⁰验收阈值且线搜索不再改善时，允许以roundoff_limited_stop标记停止。验收阈值和物理参数未放宽。成功矩阵约191秒，CPU求解，没有机器学习训练。

依赖复查：DRIFT的边缘记忆／压力端口原型已包含同总量、不同内部状态的对照；它是固定端点、不可动界面的线性毛细条带，没有可溶表活、曲面边界或破裂。FROTH不重复提出通用记忆／降阶，也不把该不可动流动律接进当前可动界面求解器；后续按所属分支审查复用。2025与2026膜—弯月面实验可作为下一阶段的独立参考，不能把实测通量律写进模型后再声称预测成功。

下一步：接到共享膜网络

本阶段不继续追加截面均匀加密。下一步处理真实膜—边界接触匹配和唯一库存网络，逐步移除指定流量与受约束形状，再连接有限速CODA液环吸收和动量传递。创新假设仍是：同总量、不同空间状态是否会在有限速事件后改变邻膜响应；本页尚未验证该假设。四泡／七泡连接及文献驱动的计划。

全部原始速度CSV（无损ZIP）；研究仓库保留各原始文件。