超导电缆的电磁场仿真
理论与物理
概述 — 超导输电的现状
超导电缆真的已经在电网中应用了吗？印象中总觉得是未来的技术…
不，已经进入实用阶段了。举几个代表性例子：德国埃森市的AmpaCity项目自2014年起就将一条10kV、40MVA、长约1km的高温超导（HTS）电缆接入商用电网并持续运行。韩国的KEPCO也在23kV电网中完成了500米级的实证。日本住友电工也进行了66kV、200米级的实证试验。
诶，已经发展到这个程度了吗！但是和普通的铜电缆相比，到底有什么优势呢？
最关键的一点是电流密度。HTS（高温超导）线材能以铜的5到10倍的电流密度导电。这直接关系到城市地下输电隧道空间的节省。例如，要确保相同的40MVA输电容量，传统的110kV充油电缆需要3回路，而超导电缆只需10kV、1回路即可。变压器也可以省略，使整个设备更紧凑。
原来如此！但是需要用液氮冷却对吧？那部分电力消耗怎么办呢？
问得好。冷却所需的电力通常约为电缆输电容量的0.5%到1%。另一方面，传统铜电缆的电阻损耗是3%到5%，所以即使扣除冷却成本，总体上仍然是节能的。不过，这里正是设计的核心，有限元分析中需要精确评估HTS线材的交流损耗与冷却负荷之间的平衡。为此，我们按顺序来看一下所需的物理模型。
E-J 幂律 — 超导体的本构关系
超导体不是电阻为零吗？那为什么还需要像欧姆定律那样的本构关系呢？
“电阻为零”是理想化的说法，实际的HTS线材在临界电流密度 $J_c$ 附近会表现出非常陡峭的非线性电阻。描述这个现象的就是E-J 幂律：
其中 $E_0 = 1\,\mu\text{V/cm}$（超导判定标准电场），$n$ 是表示陡峭程度的指数。$n$ 值越大，越接近理想的超导体。YBCO系（REBCO）线材的典型值为 $n \approx 20\text{--}40$，Bi-2223带材的典型值为 $n \approx 10\text{--}25$。
$n$ 值变大有什么影响？在解析上有什么意义？
$n$ 值越大，当 $J < J_c$ 时，实际上 $E \approx 0$（超导状态），但当 $J$ 略微超过 $J_c$ 的瞬间，$E$ 会急剧增大。这在数值上意味着极端非线性的本构关系。对于有限元求解器来说，表观电阻率 $\rho = E/J = E_0 J^{n-1}/J_c^n$ 会随电流密度变化好几个数量级，因此牛顿-拉夫森法的收敛会变得非常困难。实际上，当 $n > 30$ 时，如果没有阻尼或自适应时间步长，往往无法收敛。
E-J 幂律的物理背景与等效电阻率
E-J 幂律源于超导体内部的磁通量子（通量子）的钉扎和热激活蠕变。当外部电流产生的洛伦兹力超过钉扎力时，磁通量子开始移动（通量流动），从而产生感应电场。n值反映了钉扎力分布的均匀性。
等效电阻率张量为：
$$ \rho_{sc}(J) = \frac{E_0}{J_c} \left(\frac{|\mathbf{J}|}{J_c}\right)^{n-1} $$该电阻率在 $J \ll J_c$ 时 $\rho_{sc} \to 0$，在 $J = J_c$ 时 $\rho_{sc} = E_0/J_c \approx 10^{-14}\,\Omega\cdot\text{m}$（以REBCO，$J_c \sim 10^{10}\,\text{A/m}^2$ 为例）。相比之下，铜的电阻率为 $1.7 \times 10^{-8}\,\Omega\cdot\text{m}$，因此超导状态下的电阻率比铜低6个数量级以上。
临界电流的温度与磁场依赖性
刚才的公式里出现了 $J_c$，它是一个固定值吗？
不，$J_c$ 强烈依赖于温度 $T$ 和外部磁场 $B$ 两者。这是超导电缆设计最大的约束条件。实际中常用的是Kim-Anderson型模型：
$J_{c0}$ 是基准条件（自场，77K）下的临界电流密度，$T_c$ 是临界温度（YBCO约为92K），$\alpha$ 是温度依赖性的幂指数（典型值为1.5到2.0），$B_0$ 是磁场依赖性的尺度参数。具体数值感来说，REBCO线材在77K、自场下 $J_c \approx 3 \times 10^{10}\,\text{A/m}^2$，但如果冷却到65K，$J_c$ 会增加两倍以上。相反，在3T的外部磁场下会下降到一半以下。
那么电缆设计的关键，就在于如何平衡温度和磁场吗？
正是如此。电缆截面内的磁场分布由电缆自身产生的自场决定，而温度分布则由交流损耗和冷却的平衡决定。也就是说，确定 $J_c$ 需要知道磁场和温度，确定磁场需要知道电流（即受 $J_c$ 约束），确定温度需要知道交流损耗（电流和磁场的函数）——所有这些都是耦合的。因此，电磁-热耦合的有限元分析是必不可少的。
更精确的临界电流模型（包含各向异性）
REBCO线材因其晶体结构而具有强烈的各向异性，$J_c$ 随磁场方向变化很大。区分平行于ab面的磁场（$B_\parallel$）和c轴方向的磁场（$B_\perp$）的各向异性模型：
$$ J_c(B_\parallel, B_\perp, T) = J_{c0}(T) \cdot \frac{B_0}{\sqrt{(k \cdot B_\parallel)^2 + B_\perp^2} + B_0} $$其中 $k$ 是各向异性参数（$k \approx 0.1\text{--}0.3$）。在实际电缆中，由于HTS线材呈螺旋状缠绕，磁场相对于带材平面的入射角会随周向位置变化。这种各向异性效应会导致电流分布不均匀，直接影响设计裕度。
电缆的临界电流 $I_c$ 与裕度设计
电缆整体的临界电流，作为所有HTS线材并联的总和：
$$ I_c = \sum_{i=1}^{N} J_{c,i}(B_i, T_i) \cdot A_{sc,i} $$其中 $N$ 是线材数量，$A_{sc,i}$ 是各线材的超导层截面积。实际应用中，运行电流 $I_{op}$ 通常设定为 $I_c$ 的60%到70%（$I_{op}/I_c \approx 0.6\text{--}0.7$）。这个裕度由交流损耗和稳定性之间的权衡决定。
交流损耗的三种成分
“超导却有损耗”这不矛盾吗？电阻为零却发热…
直流情况下损耗确实几乎为零。但在交流情况下，由于磁通随时间变化，会产生三种损耗机制。这些统称为交流损耗。在超导电缆设计中，准确评估交流损耗直接关系到冷却系统的设计，因此是最重要的计算项目之一。
三种是哪三种？各自大概有多大呢？
按顺序说明一下。
1. 磁滞损耗 $Q_h$ — 源于超导细丝内的磁通钉扎，是交流损耗的主要成分。用Bean模型（临界状态模型）近似：
其中 $d_f$ 是细丝宽度（REBCO线材即导体带材宽度，典型为4到12mm），$\Delta B$ 是外部磁场的波动幅度。关键是 与 $d_f$ 成正比 —— 也就是说，如果将细丝细分（多丝化），就能显著降低磁滞损耗。REBCO线材的多丝化技术研究活跃，正是出于这个原因。
2. 耦合损耗 $Q_c$ — 细丝间基体金属（银包套或哈氏合金基板）中流动的耦合电流引起的损耗：
$\tau_c$ 是耦合时间常数，$\tau_c = \mu_0 l_p^2 / (8\pi^2 \rho_m)$。$l_p$ 是绞距长度，$\rho_m$ 是基体的电阻率。绞距越短，$\tau_c$ 越小，耦合损耗也越小。NbTi线材的绞距为10到20mm，REBCO线材则通过换位（transpose）结构等效地缩短绞距。
3. 涡流损耗 $Q_e$ — HTS线材的稳定化铜层或电缆骨架（中心导体）等常导金属部分感应的涡流损耗：
$d$ 是导体厚度，$\rho_{cu}$ 是铜的电阻率。77K下铜的电导率比室温高约7倍，因此涡流损耗可能比预期更大。特别是当电缆骨架为铜制时，这里的损耗有时不可忽视。
所有损耗加起来就是冷却负荷对吧？具体大概有多大呢？
AmpaCity级别的10kV、2.3kA HTS电缆，交流损耗约为每米1到3瓦。在77K下移除这部分热量所需的制冷机功率，考虑卡诺效率后约为损耗的10到15倍（COP ≈ 0.07到0.1），因此按电力输入换算约为每米10到45瓦。1公里电缆的整个冷却系统相当于10到45千瓦。这远小于铜电缆的输电损耗（同等容量下为每米100到200瓦），因此具有经济性。
电磁-热耦合与失超
交流损耗导致温度上升，$J_c$ 下降… 这不是会形成恶性循环吗？
正是如此。这就是失超（超导状态不可逆丧失）的本质。写出电磁-热耦合的控制方程：
右边第二项 $\mathbf{E} \cdot \mathbf{J}$ 是交流损耗（+转变为常导状态时的焦耳发热），第三项 $q_{cool}$ 是向液氮的冷却。根据E-J幂律，$E \cdot J = E_0 J^{n+1}/J_c^n$，因此当 $J$ 接近 $J_c$ 时，发热急剧增加→温度上升→$J_c$ 下降→发热进一步增加，进入正反馈循环。这就是失超。
好可怕… 电缆失超了会怎么样？
局部发生失超时，该区域会转变为常导状态并产生大量焦耳热。HTS线材热导率差（类似非晶态陶瓷层），因此失超传播速度慢（数毫米/秒到数厘米/秒，与NbTi/Nb₃Sn的数十米/秒相比极其缓慢）。这成了问题，因为检测延迟会导致局部温度达到数百度，烧毁线材。因此，失超的检测与传播分析是安全设计的最重要课题。通过有限元耦合分析预测失超的发生条件和传播速度，对于保护电路的设计至关重要。
AmpaCity — 世界首个商用HTS电缆系统
德国埃森市的AmpaCity项目（2014年至今）在城市地下铺设了1公里长的10kV、40MVA HTS电缆，取代了原有的3回路110kV充油电缆，是世界首个商用案例。值得注意的是，通过从110kV降压至10kV，无需变电站，从而释放了城市中心区的宝贵空间。冷却由两端的2台斯特林制冷机完成，采用冗余设计，即使一台故障也能继续输电。该项目的设计过程中，全面应用了基于COMSOL Multiphysics的HTS电缆电磁-热-流体三物理场耦合分析。可以说，没有有限元分析就无法实现这个项目。
数值解法与实现
有限元公式化 — A-V法与H-formulation
超导电缆的有限元分析，和普通的电磁场分析有什么不同？只解麦克斯韦方程组不行吗？
解麦克斯韦方程组本身是相同的，但超导体的极端非线性E-J特性带来了数值上的困难。公式化的选择左右着计算的成败。主要有两种方法。
A-V法（磁矢势法）是电磁场分析的标准方法，设 $\mathbf{B} = \nabla \times \mathbf{A}$，$\mathbf{E} = -\partial\mathbf{A}/\partial t - \nabla V$，以磁矢势 $\mathbf{A}$ 和电标势 $V$ 为未知数：
但在超导体区域，需使用基于E-J幂律的非线性电阻率 $\rho_{sc}(J)$ 代替 $\sigma$。问题在于，这个 $\rho_{sc}$ 会随 $J$ 变化超过10个数量级，导致雅可比矩阵的条件数极端恶化。
与此相对，近年来超导分析的主流是H-formulation：
直接以磁场 $\mathbf{H}$ 作为未知数吗？为什么这种方法更好？
H-formulation有三个优点。第一，电流密度自然地由 $\mathbf{J} = \nabla \times \mathbf{H}$ 得到，因此E-J特性的评估更直接。第二，只需在超导体外部（空气区域）设置非常大的电阻率（$\rho_{air} \sim 1\,\Omega\cdot\text{m}$），就能自然地处理磁场的侵入。第三，在二维截面分析中，它变为单个标量变量（$H_z$），因此自由度少，计算效率高。
COMSOL的HTS建模指南也推荐H-formulation，大多数学术论文都采用这种方法。
网格策略 — 薄膜近似与多尺度
HTS线材不是有1μm厚的YBCO层吗？电缆整体直径得有100mm左右吧？能用同样的网格划分吗？
指出的很尖锐。REBCO线材典型宽度为4到12mm，总厚度0.1mm，其中的超导YBCO层仅1到2μm。如果电缆外径为100mm，那么尺度…