PROJECT 06 / E04B / NUMERICAL EXPLORATION

液体进入有限边界后，吸力怎样变化？

完成 18 组主矩阵及2组容差核对，共20组原生 C++ 计算。复用 DRIFT 膜能量，将膜通量与有限边界储液放在同一次隐式更新中；同时测量单膜、共享边界及固定压力对照。

当前是固定界面的简化压力边界算例，不是已解析的膜—边界过渡区，也不是完整泡沫或实验验证。没有表活、气体交换、运动或破裂。E04 的物理交换闭合仍未完成。

2.66e-15主算例最大液量相对误差

-700.0 → -688.9 Pa有限储液边界压力

0.1627%主算例保存暂态场误差 · 对1280参考

三张膜，共享一份实际储液状态

圆盘是实际边界截面积的储液符号；右侧圆点是求解器记录的固定单元，颜色显示 h/h0（红色趋近0，蓝色约1.1）。圆盘及膜厚均放大显示；不表示真实尺寸的泡沫表面。

仅播放 51 个保存状态，不插值，不预设轨迹。物理时长 1 ms。

怎样解读这些结果？

有限边界与固定压力对照具有相同初始吸力，有限边界吸液后压力随截面积更新。这只检验当前压力边界模型的反馈。半单元 Poiseuille 通量是该边界模型的离散化，不是已标定的真实过渡区阻力。半径扫描同时改变初始吸力与容量，不能单独归因为容量。细网格是同方程参照；局部最小值与整体场误差分别保留。

另完成2组1280单元成对计算，将 Newton 容差收紧10倍。排液差异估计的相对变化为 1.201e-11，小于本次0.1%敏感性门槛。因此已完成的20组配置共101,501个物理步。

另有1次收紧100倍的尝试，在1507步 / 0.0941875 ms后线搜索失败，未计为成功。失败记录保留，未自动放宽容差或接受失败步。

算例	单膜单元数	步数	最大液量误差	最小厚度 / μm	耗时 / s
pilot40	40	1000	1.55e-15	0.6508	0.10
main80	80	1000	2.66e-15	0.5891	0.15
space160	160	1000	2.44e-15	0.5706	0.22
reference320	320	4000	5.00e-15	0.5660	1.06
time-coarse	80	500	8.88e-16	0.5887	0.09
time-half	80	2000	2.89e-15	0.5893	0.18
fixed80	80	1000	2.66e-15	0.5834	0.12
radius25	80	1000	2.22e-15	0.3769	0.14
radius100	80	1000	9.99e-16	0.8978	0.13
shared-three	80	1000	3.55e-15	0.6001	0.32
one-step	40	1	0.00e+00	1.9652	0.02
pair160-finite	160	2000	4.88e-15	0.5708	0.30
pair160-fixed	160	2000	1.78e-15	0.5649	0.27
pair320-fixed	320	4000	4.11e-15	0.5601	0.86
pair640-finite	640	8000	6.99e-15	0.5648	3.35
pair640-fixed	640	8000	1.14e-14	0.5589	2.91
pair1280-finite	1280	16000	1.13e-14	0.5645	11.93
pair1280-fixed	1280	16000	1.58e-14	0.5586	11.68
tight10x1280-finite	1280	16000	1.38e-14	0.5645	11.10
tight10x1280-fixed	1280	16000	1.58e-14	0.5586	11.31

参数均为示例输入，未匹配实测材料。新图片保留 PNG 和矢量 PDF。