VARIDRAINAGE / DRIFT / 2026-11-11

局部边缘与界面响应：
第一轮探索计算

冻结协议 → 原生回归 → 64组界面响应 → 12组初检与4组延时边缘诊断。白底图示均来自本批实际保存输出。

NUMERICAL DIAGNOSTICS / C1原型前置检查 / C2尚未获得实测动力学支持

135 / 135原生检查通过

64 + 16响应工况 + 边缘工况

52.5%候选局部网格单元减少，仍需运动粒子验证

这批边缘计算使用固定有限体积单元、静止界面和固定负压储库；没有表面流、解析弯月面或边缘再生。面积响应采用未标定的有效交换率，尚未耦合排液。不能将这些检查表述为Monier实验复现、H2机制证实、黑膜形成或完整C1/C2贡献。

本轮推进了什么

E1 · 局部加密与守恒接口账本

200×80 μm、2 μm初始膜，侧向1%扰动；对照均匀细网格与30/60 μm边缘加密区域。每条跨区面只有一个通量，在两侧液量账本中取相反号。计算到1 ms。

候选局部网格2432个单元，对照5120；单元减少52.5%。积分流出误差0.0003868%，厚度场L2误差0.0022%，变薄量L2误差0.056%。

E2 · 独立的均匀面积响应

继承SDS平衡热力学，使用两侧表面量与体相总量的闭合账本。面积改变±0.1%/±1%，检查2/6 μm、四个未标定响应时间和时间步减半。

全部64组通过可行性、表活收支与固定面积耗散门槛；最大表活相对收支1.47e-16。最大粗/细时间步张力差为初始响应的0.465%。

局部边缘：厚度与通量误差

边缘剖面及通量误差对照 — 同一简化方程和相同材料输入。参考解dx=1.25 μm、dy=2.5 μm、dt=2.5 μs；通量10%是内部诊断门槛。矢量PDF

均匀细网格与局部网格的单元中心圆盘场 — 圆盘是有限体积单元中心，不是运动粒子；虚线为60 μm接口。没有网格插值渲染。矢量PDF

实际保存状态回放

仅播放保存状态，不做时间插值；界面是简化计算域边界，非真实弯月面。颜色：蓝1.6 μm → 浅色 → 红2.05 μm。

完整边缘工况表 / 12组

Case	Cells	min h / μm	液量收支	积分流出误差	变薄量L2误差	计时 / s	守恒/耗散
crosscheck_1d	40	1.67135	3.9e-16	0.000%	0.000%	0.19	PASS
uniform10	320	1.71996	1.3e-14	8.688%	11.026%	0.47	PASS
uniform5	640	1.66571	1.0e-14	2.065%	2.741%	0.72	PASS
uniform2p5_y5	1280	1.65064	2.4e-14	0.382%	0.549%	1.16	PASS
uniform2p5	2560	1.65061	1.4e-13	0.382%	0.548%	4.08	PASS
reference_dt5	5120	1.64674	3.3e-13	0.043%	0.033%	8.49	PASS
reference_dt2p5	5120	1.64677	9.3e-13	0.000%	0.000%	17.88	PASS
local5_bulk10	416	1.66571	9.9e-15	2.066%	2.745%	0.53	PASS
local2p5_bulk10	1216	1.65063	1.1e-13	0.384%	0.585%	1.86	PASS
local1p25_bulk5	2432	1.64674	4.7e-14	0.042%	0.066%	3.77	PASS
local1p25_bulk5_dt2p5	2432	1.64678	5.0e-13	0.000%	0.056%	8.17	PASS
local1p25_patch30	1856	1.64698	2.5e-13	0.006%	0.387%	5.58	PASS

一维回归行不作二维误差比较。计时包含建图、求解和账本，不含状态归档；单次进程内运行。单元数减少不等同于已验证的完整求解器加速。

时间延长十倍：保留检查

保持材料、几何、初始场和边界不变，将终止时间从1 ms延长到10 ms。四组工况在执行前冻结：均匀细参考dt=10/5 μs，以及60/30 μm局部区域dt=5 μs；没有依据结果重调参数。

原候选60 μm区域相对较小时间步参考的积分流出差0.0039%，变薄量L2差0.082%。本检验仍限于同一简化方程，不代表运动弯月面或实验验证。

10ms保留时间检查 — 更长时间下的厚度与变薄量。矢量PDF · 结果与源提交 · 冻结协议

Case	Cells	积分流出差	变薄量L2差	液量收支	内部精度门槛
holdout_ref_dt10	5120	0.0249%	0.017%	3.9e-13	PASS
holdout_ref_dt5	5120	0.0000%	0.000%	2.7e-12	PASS
holdout_local60	2432	0.0039%	0.082%	8.8e-13	PASS
holdout_local30	1856	0.0021%	0.302%	8.9e-13	PASS

界面恢复：可能的效应与尚缺的证据

不同交换时间下的张力恢复 — 实线为有限速率；虚线是不交换极限；零线为相同总量的闭合平衡。固定面积后的自由能下降通过检查，面积阶跃的外功另计。矢量PDF

有效正反交换遵循继承热力学的详细平衡；瞬时体相胶束平衡保持不变。τ bath是固定体相的线性时间尺度，并非已测SDS速率，也未必等于闭合膜的精确弛豫时间。本轮只证明数值响应可区分；是否影响实际边缘交换，仍须独立动力学数据和运动/输运耦合。

h0 / μm	面积比	τ bath / ms	初始Δγ / mN·m⁻¹	1ms剩余Δγ	表活收支
2	1.001	0.1	0.3326	0.0000	1.5e-16
2	1.001	1	0.3326	0.1252	1.5e-16
2	1.001	10	0.3326	0.3015	1.5e-16
2	1.001	100	0.3326	0.3293	1.5e-16
6	1.001	0.1	0.3327	0.0000	0.0e+00
6	1.001	1	0.3327	0.1253	0.0e+00
6	1.001	10	0.3327	0.3016	0.0e+00
6	1.001	100	0.3327	0.3294	0.0e+00

表格选取+0.1%面积、较小时间步；全部64组及时间步差见JSON。未根据实验指数筛选交换率。

决策与下一步

以局部网格的守恒和精度结果作为参考，实现运动表面/粒子局部边缘及其输运接口；仍需解析表面流和真实夹缩。
保持C2为条件性扩展：先获得独立速率或跨膜扩散证据，再决定是否接入整体膜。
在局部通量和运动接口收敛前，不启动16×5 mm主装置长跑，不报告条纹指数吻合或黑膜实验验证。

复现与证据

关键结果 12组边缘结果 64组响应与时间步比较运行协议边缘参数运行状态冻结实验计划总计划书

计算源提交 6dd72427c7cda1a96f460e308540fe618eadd438。完整CSV、NPZ与运行命令位置见本地RESULTS；原生回归与边缘文本日志纳入Git。复现：构建test_surface_exchange、tf_response_bench及原有测试后运行 thin_film/scripts/phase2/run.py，再运行report.py。

局部边缘与界面响应：第一轮探索计算