IYPT 2026 No. 8 磁力加速器

理论现状与数值解计划 | 更新: 2026-08-11

一句话结论：这是一个有限行程去钉扎问题。发散几何产生有限、空间波纹状的磁势能景观；只有当独立签署的端点落差和每个下游势垒/损耗裕度都允许时，宽→窄设置才会加速。运动将有限落差转化为平动和车轮转动，而涡流阻尼、滚动损耗、接触、抬起、滑移和横向/偏航不稳定性限制了可实现性。

三条承载方程及其物理意义

以下三条方程是评委需要看到的全部公式。每条都有明确的物理意义，它们构成一个自洽循环——不是三条孤立的模型。

方程 1 · 磁体相互作用能（"预算"）

$$U_{12}=\frac{\mu_0}{4\pi}\iint\frac{\sigma_1(\mathbf r)\sigma_2(\mathbf r')}{|\mathbf r-\mathbf r'|}\,dS\,dS'$$
这是什么：两块刚性、均匀磁化磁体之间的静磁相互作用能。$\sigma_i = \mathbf M_i \cdot \hat{\mathbf n}_i$ 是磁体 $i$ 表面的等效磁荷密度（只有极面携带磁荷；内部 $\nabla\cdot\mathbf M=0$）。双重积分遍历两个表面。
物理意义：这是能量预算——磁场中储存了多少势能。它回答："如果我把这两块磁体放在这里，场有多少能量？"
它去哪：对所有车轮-固定磁体对求和，得到总势能景观 $U(s) = \sum_{i,j} U_{12}^{(i,j)}(s) + E_{\rm plate}(s)$。这个 $U(s)$ 是方程 2 的原材料。

方程 2 · 运动方程（"花钱"）

$$M_{\rm eff}\dot v=-\frac{dU}{ds}-\beta(s)v-\xi,\quad \xi\in F_{rr}\,\partial|v|,\qquad \dot s=v$$
这是什么：车的牛顿第二定律，含三个力项。
每项的物理意义：
$-\dfrac{dU}{ds}$ —— 磁力：
这是方程 1 给出的势能景观的梯度。如果 $U$ 随 $s$ 增加而减小（宽→窄），则 $-dU/ds > 0$，车被向前推。它回答："场想往哪推车，用多大力？"
$-\beta(s)v$ —— 涡流阻尼：
车运动时，金属板中变化的磁通感应出电流，以焦耳热形式耗散能量。$\beta(s)$ 是位置依赖的阻尼系数。它回答："板从运动中偷走了多少能量？"
$-\xi$ —— 滚动阻力：
阻碍运动的广义摩擦力。在 $v=0$ 时成为静阈值：只有当 $|dU/ds| > F_{rr}$ 时车才开始动。$\xi \in F_{rr}\,\partial|v|$ 中的 $\partial|v|$ 是 $|v|$ 的次微分：$v \neq 0$ 时取 $\mathrm{sgn}(v)$，$v=0$ 时取 $[-1,1]$ 区间。
它从哪来：方程 2 中的 $U(s)$ 正是方程 1 对所有磁体对的求和，再加上方程 3 的修正（板自身的能量）。
它去哪：解这个常微分方程给出 $v(s)$ 和分类（exit / trapped / reversed）。出口速度 $v_f$ 告诉你预算是否花得明智。

方程 3 · 铁板响应（"反馈"）

$$\mathbf M(\mathbf r)=\chi_p\,\mathbf H_{\rm local}(\mathbf r),\qquad \mathbf H_{\rm local}=\mathbf H_{\rm ext}+\mathbf H_{\rm demag}(\mathbf M)$$
这是什么：铁磁板的构成定律，与退磁场耦合。
物理意义：
它从哪来：外场 $\mathbf H_{\rm ext}$ 来自硬磁体，由方程 1 中出现的相同表面磁荷 $\sigma$ 计算。
它去哪：解出的 $\mathbf M$ 通过板能量 $E_{\rm plate} = \frac{\mu_0}{2}\int \mathbf M \cdot \mathbf H_{\rm demag}\,dV$ 修改方程 2 中的总能量 $U(s)$。它也修改涡流阻尼 $\beta(s)$，因为变化的 $\mathbf B$ 场包含板的响应。

自洽循环

这三条方程不是清单——它们是一个循环：

方程 1（预算） 方程 2（花钱） 方程 3（反馈） ┌─────────────┐ ┌─────────────┐ ┌─────────────┐ │ U₁₂ 成对能 │──求和──→│ U(s) 势能 │ │ M = χH_loc │ │ 景观形状 │ │ -dU/ds 力 │←β(s)──│ H_loc = │ │ ΔU 端点落差 │ │ v(s) 速度 │ │ H_ext+H_demag│ └─────────────┘ │ exit/trapped│ └─────────────┘ └─────────────┘ │ ↑ │ └──── M 修改 U ←───────┘
人话三步：
1
方程 1 问："银行里有多少钱？"（能量预算 $\Delta U$）
2
方程 2 问："我们能多快花掉它？"（力 $-dU/ds$ 推车；阻力和摩擦分走一杯羹）
3
方程 3 问："板在我们花钱时会改变银行余额吗？"（板磁化，改变场，改变能量，改变力）

循环闭合是因为板的磁化（方程 3）依赖于磁体的场（方程 1），而车的运动（方程 2）改变相对位置，从而改变场和板的响应。

附录 A：工程实现细节

数值计算时，耦合本构律 $\mathbf M=\chi_p\mathbf H_{\rm local}$ 被离散化到网格上。设 $f_p(x,y)$ 为冻结梯形内的精确单元面积分数，使用变换未知量 $\mathbf Q=\sqrt{f_p}\,\mathbf M$。离散系统是对称正定（SPD）线性系统：

$$\left[I+\chi_p\sqrt{f_p}\,\mathcal N^{(n_z)}\sqrt{f_p}\right]\mathbf Q=\chi_p\sqrt{f_p}\,\mathbf H_{\rm ext}$$

其中 $\mathcal N^{(n_z)}$ 是 $n_z$ 层退磁算子（由第 6.1 节 Fourier 表示计算的 Hermitian 半正定矩阵），$\sqrt{f_p}$ 是梯形面积分数掩码。

工程形式性质：

本附录是唯一出现矩阵形式的地方。正文使用物理本构律（方程 3）及其耦合场解释。

当前状态