M19
轴对称哑铃曲面平均曲率流颈缩奇点的数值测定：Type-I/Type-II 临界参数与收敛性护栏
1 · 研究问题
对单参数哑铃族初值的轴对称平均曲率流（mean curvature flow, MCF），用自建一维径向 PDE 求解器能否复现颈缩（neckpinch）奇点的 Type-I 标度 max|A| ~ (T−t)^(−1/2)，并把"颈缩 vs 整体收缩"的临界参数 λ_c 测到三位有效数字？进一步：λ_c 及标度指数对空间离散方式（参数化 vs 径向图）是否稳健？
2 · 研究背景与空白
MCF 是几何分析核心对象；轴对称情形约化为 1D 抛物 PDE，是笔记本可及的奇点数值实验室。理论：Huisken 的 Type-I 颈缩、Angenent–Velázquez 的退化颈缩 Type-II 精确速率 (T−t)^(−1−1/m) 均已证明。数值：Garfinkle 等《Critical behavior of collapsing surfaces》(arXiv:0903.3369) 在旋转对称哑铃族上测得临界值 λ_c ≈ 0.9076 并观察临界解行为；Isenberg–Wu 合作的系列数值稳定性研究（arXiv:2105.06613, 2206.01890）覆盖非紧 Type-II 情形。
已有工作确立了现象与个例数值。空白在方法可复现性维度：已发表数值工作各用自己的离散化，未见对"λ_c 与标度指数随离散方式/分辨率的稳健性"的系统对照；且 λ_c 仅见个位报道、无跨方法复核（未找到不等于不存在）。把收敛性护栏本身作为交付物，是绕开算力墙的合法课题形态。
3 · 可检验假设
- H1 自建求解器在颈缩区复现 Type-I 指数:−d(log max|A|)/d(log(T−t)) → 0.5，末段拟合值偏差 ≤ 5%（在网格外推意义下）。
- H2 两种独立离散化（径向函数差分 + 弧长参数化）给出的 λ_c 在各自 Richardson 外推误差内一致，且与已发表值 0.9076 的偏差 ≤ 1%——或测得系统性偏差并能归因于初值族定义差异。
4 · 量化验收标准
- 方法学校验（硬门槛）：自建求解器复现两个解析解——球面收缩 R(t)=√(R₀²−4t)（曲面情形系数 2(n−1)）误差 ≤ 0.1%；圆柱颈的自相似线性化增长率与理论值偏差 ≤ 2%。不过关则奇点测量全部无效。
- 收敛性检查（护栏即交付物）:网格 4 档 + 自适应加密策略写明；奇点时刻 T 的估计随分辨率的漂移序列必须呈现收敛趋势并给外推值与余项。
- 守恒/单调量监测:表面积单调下降、avoidance 原理（内含球面下界）逐步校验，违例即判失效。
- 奇点邻域拟合窗口选取规则写死（如 max|A| 上升 2–4 个十倍程区间），窗口敏感性必须报告。
- λ_c 用二分法定位至 Δλ ≤ 5×10⁻⁴，每次判定（颈缩/无颈缩）判据可复现（颈半径低于 10⁻³ 且加密后维持）。
- 代码开源、一键重跑。
5 · 数据与工具
| 用途 | 来源 / 工具 |
|---|---|
| PDE 求解 | 自建:1D 有限差分 + 刚性隐式推进（scipy.integrate.solve_ivp 的 Radau/BDF 内置；自适应网格须自建——这是主要工程量） |
| 解析参照 | 球面/圆柱精确解（手推）；Angenent 颈缩定理表述 |
| 数值对照 | arXiv:0903.3369 的 λ_c ≈ 0.9076 与初值族定义；arXiv:2105.06613 方法细节，仅用于校验与对比，不计入贡献 |
| 符号推导 | sympy 核验轴对称 MCF 方程与线性化 |
| 可视化 | matplotlib 旋转面重建 |
| 算力 | 1D 网格 ≤ 10⁴ 点、单次流演化分钟级；二分 ~50 次 × 2 种离散化 × 收敛档，总 CPU < 50 小时 |
6 · 方法路径
- 推导轴对称 MCF 的径向图方程（sympy 核验），实现差分求解器。
- 球面/圆柱双解析校验（第 1 条验收）。
- 实现自适应加密与奇点邻域诊断，冻结判据。
- 复现单个哑铃颈缩，测 Type-I 指数（H1）。
- 实现第二种离散化（弧长参数化），同一初值互验。
- 二分定位 λ_c，两离散化 + 分辨率外推给合成误差。
- 与发表值对照；差异若超界，逐项排查初值族定义（差异归因写入论文）。
7 · 新颖性边界
本课题不声称发现新奇点类型或证明任何定理；Type-I/Type-II 理论（Huisken、Angenent–Velázquez）与 λ_c ≈ 0.9076 的数值测定（Garfinkle 等）均已发表。本项目的主结论是跨两种离散化的 λ_c 独立复测、误差可控的收敛性护栏、以及一条可操作的奇点数值判据——属"方法层面可复现性贡献"，在计算领域被承认，但必须如此定位，否则会被误读为重复。若复测与发表值偏差超 1%（H2 负向），差异归因（初值族/判据定义）本身构成有效主结论。历届丘奖几何流相关获奖（2020 mention cigar soliton 单调性公式、2022 mention Hamilton 等周比）皆为理论方向，无重合；本题恰好接续该评委口味谱系的数值侧。
8 · 决策门槛（go / no-go）
- 第 4 周末:球面/圆柱校验通过。若隐式推进在颈缩前就失稳，先换 BDF/减小容差；2 周未解决则降级路径 A：放弃逼近奇点，改测"远离奇点段的收敛阶 + 颈半径演化与理论比较"，保留方法学对照框架。
- 第 10 周末:单个颈缩 Type-I 指数可测（拟合窗 ≥ 2 个十倍程）。若自适应网格达不到，降级路径 B：固定环高分辨率网格 + 提前终止（max|A| 达 10³ 即止），指数拟合窗缩短、误差棒放大，H1 保留。
- 第 20 周末:第二种离散化若工程量失控，砍为"单离散化 + 三档分辨率外推"，H2 改述为分辨率稳健性（跨离散化改为跨分辨率），主结论降一档保留。
- 适合对象:对几何分析有真兴趣、PDE 数值有耐心的学生；是几何族中风险最高但评委辨识度最高的一题。