超导电缆的电磁场仿真
超导电缆的电磁场理论基础
概述 — 超导输电的现状
超导电缆真的在电力网中使用了吗？ 我印象中这是未来的技术…
不，已经进入实用阶段了。举几个典型例子吧，德国埃森市的AmpaCity项目从2014年开始，一条约1km长的10kV、40MVA HTS电缆连续运行在商业电网中。韩国KEPCO也完成了23kV系统中500m级的实证。日本住友电工进行了66kV、200m级的实证试验。
哇，已经这么先进了！ 但与普通铜电缆相比优势在哪里呢？
最关键的是电流密度。HTS（高温超导）线材的电流密度比铜高5～10倍。这直接关系到城市地下输电隧道的空间节约。比如同样40MVA的输电容量，传统110kV OF电缆需要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$ 会随电流密度变化几个数量级。对有限元求解器来说，Newton-Raphson法的收敛会变得极其困难。实际上，当 $n > 30$ 时，如果没有阻尼或自适应时间步，收敛往往会失败。
临界电流的温度·磁场依赖性
那个 $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$ 限制）决定；温度由交流损耗（取决于电流和磁场）和冷却决定——全部耦合在一起。因此电磁-热耦合有限元分析是必不可少的。
交流损耗的3个成分
"超导体却有损耗"，这不矛盾吗？ 阻抗为零却发热…
直流时确实损耗接近零。但交流时，磁通的时间变化引起3种损耗机制，统称交流损耗。在超导电缆设计中，准确评估交流损耗对冷却系统设计有直接影响，是最重要的计算项目之一。
那3种损耗分别是什么？ 各占多大的比例？
逐个说明：
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线则用转位（转置）结构等效缩短螺距。
3. 涡流损耗 $Q_e$ — HTS线的稳定化铜层或电缆导体（如中心芯）等常规金属中感应的涡流损耗：
$d$ 是导体厚度，$\rho_{cu}$ 是铜的电阻率。77K的铜导电率是室温的约7倍，所以涡流损耗可能大于预期。特别是电缆中心芯为铜时，这部分损耗不能忽视。
三种加起来的交流损耗大概多大？
AmpaCity级10kV、2.3kA的HTS电缆，交流损耗约1～3 W/m。这点热量要在77K下除去，考虑卡诺效率，冷冻机所需电力大约是损耗的10～15倍（COP ≈ 0.07～0.1），电气输入换算后为10～45 W/m。1km电缆的冷却系统功耗约10～45 kW。相比之下，同容量铜电缆的输电损耗是100～200 W/m，所以经济上是可行的。
电磁-热耦合与淬火
交流损耗温度升高，$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线材热传导差（陶瓷阿莫夫状），淬火传播速度很慢（数mm/s～数cm/s，远低于NbTi/Nb₃Sn的数十m/s）。检测延迟会导致局部温度达数百度，线材烧伤。因此淬火检测与传播分析是安全设计的重中之重。用有限元耦合分析可以预测淬火的发生条件和传播速度，这对保护电路设计至关重要。
AmpaCity — 世界首条商用HTS电缆系统
德国埃森市的AmpaCity项目（2014年起）在市区地下敷设了10kV、40MVA的HTS电缆1km，替代了原来的110kV OF电缆3回路，是全球首例商业运行的HTS输电案例。注意的是，从110kV降压到10kV消除了降压变电所，节约了城市中心宝贵的用地。冷却由两台斯特林冷冻机负责，一台故障时系统仍可运行，设计上有冗余。这个项目的设计过程中，COMSOL Multiphysics的HTS电缆电磁-热-流体三物理耦合分析被全面采用。没有有限元分析，这个项目无法实现。
超导电缆的电磁场数值计算方法
有限元定式 — A-V法与H-格式
超导电缆的有限元分析和普通电磁场分析有什么不同？ 求解Maxwell方程不就行了吗？
Maxwell方程本身没变，但超导体极端的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$ 变化十几个数量级，雅可比矩阵的条件数恶化得极其严重。
近年来超导分析的主流改为H-格式：
磁场 $\mathbf{H}$ 作为直接未知量？ 为什么这样更好呢？
H-格式有三大优点。第一，电流密度自然得到 $\mathbf{J} = \nabla \times \mathbf{H}$，E-J关系直接而明确。第二，超导体外（空气）只需设置极大的电阻率（$\rho_{air} \sim 1\,\Omega\cdot\text{m}$）就能自然处理磁通不侵入。第三，2D截面模型中只有一个标量变量（$H_z$），自由度少、效率高。COMSOL的HTS建模指南也推荐H-格式，学术论文里绝大多数都用这个方法。
网格策略 — 薄膜近似与多尺度
HTS线材里有1微米厚的YBCO层，电缆总直径100毫米吧？ 一个网格怎么网得下来？
精准的观察。REBCO线材通常宽4～12mm、总厚0.1mm，其中YBCO超导层仅1～2微米。电缆外径100mm的话，尺度比是1:100,000。直接网格划分这个是不现实的。
解决办法是薄膜近似（thin shell approximation）。YBCO层不当体积单元，而当面单元处理，引入面内的面电流密度 $K_s\,[\text{A/m}]$ 作为自由度：
$d_{sc}$ 是超导层厚度。这样一来，就不用在厚度方向细分1微米的层，网格数可以削减2～3个数量级。COMSOL中这叫"Shell功能"。代价是放弃了超导层厚度方向的电流密度分布，但在磁通完全侵入条件下（侵入深度大于层厚）这个近似就成立。
可电缆里有几十条甚至几百条线材螺旋绕着…… 每条都要单独模型吗？
如果全部单独建模那也不现实。这时均质化（homogenization）就派上用场了。在后面的章节详细讲。把某一HTS线层看作异方性的连续介质，定义有效的 $J_c$ 和有效阻抗率张量，就能大幅降低自由度，使大规模电缆（数百米～数千米）的分析成为可能。精度上，个别线材模型对比，电流分布误差5%以内，交流损耗误差10%以内。
时间积分与非线性求解器
$n$ 值大的时候求解器收敛困难，具体怎么对付？
有几个技巧：
- BDF（后退差分）时间积分：1～5阶BDF格式为标准。当 $n > 20$ 时，2阶BDF在稳定性和精度间取得较好平衡
- 自适应时间步：电流通过 $J_c$ 附近时自动细分时间步。典型 $\Delta t_{min} \sim 10^{-6}\,\text{s}$，$\Delta t_{max} \sim 10^{-3}\,\text{s}$
- Newton-Raphson阻尼：限制更新步长（$\delta H_{max} < 0.1 H_{peak}$）防止振荡
- 初值策略优化：用斜坡函数缓慢施加外部磁场，避免初期磁通突然侵入
另一个关键是雅可比的解析计算。E-J幂律的雅可比（切向阻抗率张量）为：
用解析式而不是数值微分来计算，Newton-Raphson的收敛速度能显著提升。COMSOL中可设置"Analytical Jacobian contribution"。如果这个被忽略，收敛迭代次数会增加2～3倍。
均质化模型 — 大规模电缆应用
均质化模型的细节教一下。实际上怎样把"一层胶带"变成"等效材料"？
把电缆的一层HTS线材看作圆柱形壳的异方性连续体。沿螺旋绕线方向的有效临界电流密度为：
$f_{fill}$ 是填充因子（该层HTS线在周向占的比例，通常0.7～0.9），$\theta_w$ 是绕线角（相对电缆轴的螺旋角，通常15°～35°）。在垂直（周向）上不流电，作为各向异性的阻抗率张量实现。
这个均质化模型能将自由度降低100～1000倍，使得电力电缆（数百米～数千米长）级的分析成为可能。精度上与单条胶带模型对标，电流分布误差 <5%，交流损耗误差 <10%，大多数文献都有这样的报告。
超导电缆的电磁场实务应用
分析工作流程
超导电缆的电磁场分析，实际流程是什么？ 最开始怎么做？
典型的工作流程如下：
第1步：电缆断面几何建模
- 杜瓦瓶外管（不锈钢）→ 断热真空层 → 杜瓦瓶内管 → 液氮流路 → HTS线层（多层）→ 中心导体的同心圆结构
- 每层HTS胶带数、绕线角、间距定义
- 2D截面模型为主。3D仅在电缆终端（接线端子）分析时必需
第2步：材料参数定义
- HTS层：E-J幂律（$E_0$、$J_c(B,T)$、$n$ 值）用温度-磁场表格定义
- 稳定化铜：$\rho_{cu}(T) = 1.7 \times 10^{-8} \times (T/293)$ [Ω·m]（低温时也需考虑剩余电阻比RRR）
- 液氮：从沸腾曲线读取热交换系数 $h(T)$（核沸腾、膜沸腾的分别）
- 断热层：$\kappa_{eff} \sim 10^{-4}$ W/(m·K)（多层断热MLI）
第3步：网格生成
- HTS层：应用薄膜近似时用壳单元；用体单元模型时在厚度方向至少3层
- 空气域：电缆外径的5～10倍范围。外周用无限单元
- 液氮流路：若连成流动分析，需用棱柱单元的边界层
第4步：求解器设置·计算执行
- H-格式：BDF2、初始时间步 $10^{-5}$s、最大 $10^{-3}$s
- Newton-Raphson：相对收敛判定 $10^{-4}$～$10^{-6}$、阻尼系数0.5～0.9
- 至少运行3～5个周期直到到达交流稳定（过渡项衰减）
第5步：后处理·交流损耗计算
- 用稳定周期（最后一个周期）的数据，以体积积分算交流损耗：$Q = \int_V \int_0^{1/f} \mathbf{E} \cdot \mathbf{J}\,dt\,dV$
- 分离各层的损耗明细（HTS层的滞后损、中心导体的涡流损等）
- 检验电流分布均匀性（层间电流平衡）
边界条件设置
H-格式的边界条件和A-V法是不一样的吗？
是的，H-格式的边界条件如下：
| 边界 | 条件 | 物理意义 |
|---|---|---|
| 计算域外周 | $\mathbf{n} \times \mathbf{H} = 0$（Neumann） | 磁感线垂直射入边界（远场磁场衰减） |
| 对称面 | $\mathbf{n} \times \mathbf{H} = 0$ 或 $\mathbf{n} \cdot \mathbf{H} = 0$ | 对称/反对称条件，计算量减半 |
| HTS导体边界 | 自然边界（连续条件） | $H$ 的切线分量连续（物理自然） |
| 电流施加 | 积分约束 $\oint \mathbf{H} \cdot d\mathbf{l} = I(t)$ | 安培定律强制电流 |
电流用安培定律的环形积分来施加！ 有意思的想法呢。
对。在H-格式中，电流密度不是直接未知量，而是 $\mathbf{J} = \nabla \times \mathbf{H}$ 的结果。所以要指定流过导体的总电流，就要对包围该导体的闭合路径约束 $\mathbf{H}$ 的线积分。COMSOL中用"Coil Group"功能实现。三相电缆的话，对各相分别施加相位相差120°的电流。
验证与合理性确认
分析结果对不对，怎样确认？ 超导电缆的实验数据也不多吧？
V&V（验证与合理性确认）的步骤如下：
验证（Verification）— 代码是否正确求解方程
- 单根胶带的Bean模型解析解对比（Norris式：$Q_h = \frac{\mu_0 I_c^2}{\pi}[f \ln(1-f) + (1-f)\ln(1-f) + f]$，$f = I_{peak}/I_c$）
- 网格收敛性：3个及以上密度级，交流损耗收敛度 < 2%
- 时间步依赖性：时间步减半后结果无变化
合理性确认（Validation）— 计算是否真实再现物理
- 短尺度样品的交流损耗卡路里测量值对比
- 临界电流的实测值（自磁场、外部磁场）与模型预测对比
- 淬火电压的实测波形对比（如果有）
- 对于HTS-CICC（Cable-In-Conduit Conductor），有来自CEA-Cadarache、CERN等公开的基准数据可用
常见失败模式
新手容易犯的错误有哪些？ 事前知道就好了。
实务中常见的失败模式汇总如下：
| 失败模式 | 原因 | 对策 |
|---|---|---|
| 交流损耗是实验值的2～5倍 | 忽视 $J_c(B)$ 的磁场依赖性，全区用自磁场下 $J_c$ | 正确计算各处的 $J_c(B,T)$ |
| 层间电流分配不合理 | 各HTS层的绕线角、胶带数输入有误 | 对照设计值与各层的电感和 $I_c$ 核实 |
| 不收敛 | 过渡计算用作结果 | 至少计算3～5周期，用最后周期 |
| 求解不收敛 | 空气域 $\rho_{air}$ 太小或太大 | $\rho_{air} = 1\,\Omega\cdot\text{m}$（$\rho_{sc}$ 最大值的约100倍） |
| 淬火温度不准 | 用室温 $C_p$ 代替HTS线的低温比热 | 用低温的 $C_p(T)$ 表格（77K约为室温的1/3） |
"过冷问题" — 液氮亚冷产生的意外设计裕度
液氮沸点为77.3K（1 atm），但加压亚冷可降至65～70K时，$J_c$ 增加接近两倍。这能使电缆的电流容量翻番，或同容量下胶带用量减半。AmpaCity二代的设计中曾评估过5气压、67K亚冷运行。只是温度越低，冷冻机的COP越差（65K时与77K比约需1.5倍的能耗），所以有限元分析必须定量评估$J_c$增益与冷却成本增加的权衡。
超导电缆的电磁场软件比较
COMSOL — AC/DC模块
超导电缆分析用得最多的软件是什么？
学术论文的绝大多数采用COMSOL Multiphysics。理由很清楚：AC/DC模块的标准功能不够，但PDE模块能直接定义H-格式，对超导体的非线性E-J特性可以灵活实现。
- 优点：多物理耦合（电磁-热-流体）一站式完成、PDE自由定义、薄膜近似的Shell功能
- 缺点：大规模3D（胶带100+条）内存消耗大、并行扩展有限
- 许可：AC/DC + Heat Transfer + CFD + PDE 年约200万人民币
Ansys Maxwell + Mechanical
Ansys也能做超导电缆分析吗？
Ansys Maxwell是基于A-V法的电磁场求解器，对非线性B-H曲线处理得当。但E-J幂律之类电流相关的非线性阻抗就不如COMSOL容易，通常需自定义UDF才行。
- 优点：电磁场精度高、涡流分析的自适应网格加密、Workbench通过电磁-热-结构耦合
- 缺点：没有标准H-格式支持、E-J特性实现要费工夫
- 实绩：超导磁体（MRI、核融合）的淬火分析有大量Maxwell + Mechanical应用
开源软件（GetDP、Elmer、FEniCS）
经费有限的研究室也能用的开源工具有吗？
超导有限元分析有实绩的开源软件有3种：
- GetDP（比利时列日大学）：超导研究社区中用得最广的开源。H-格式对应，和Gmsh集成。有CERN大型磁体分析的实绩
- Elmer（芬兰CSC）：多物理开源求解器。超导模块有限，但可自定义
- FEniCS/FEniCSx：Python有限元框架。H-格式数行代码搞定，研究原型制作理想。但产业级网格生成、后处理需配套工具
功能对比矩阵
要选哪个，一览表对比一下？
| 功能 | COMSOL | Ansys Maxwell | GetDP | FEniCSx |
|---|---|---|---|---|
| H-格式 | ◎（经PDE） | △（UDF） | ◎（标准） | ◎（直接定义） |
| E-J幂律 | ◎ | △ | ◎ | ◎ |
| 薄膜近似 | ◎（Shell） | △ | ○ | ○（手工） |
| 电磁-热耦合 | ◎ | ○（Workbench） | ○ | ○（手工） |
| 3D大规模 | ○ | ◎ | ○ | △ |
| GUI/学习成本 | ◎（低） | ○ | △（脚本） | △（Python） |
| 费用 | 高 | 高 | 免费 | 免费 |
| 超导文献量 | 非常多 | 中等 | 较多 | 增加中 |
COMSOL一家独大呢。不过GetDP在论文里也常见。
对。列日大学的团队（B. Ramdane、F. Grilli等）用GetDP做超导分析的研究持续推进，还公开了基准问题。CERN的FCC（未来圆形对撞机）用Nb₃Sn超导磁体设计中也用GetDP。经费紧张的研究室，先用GetDP + Gmsh开始是务实的选择。
超导电缆的电磁场前沿研究
机器学习代理模型
超导电缆分析一次要跑好几小时吧？ 快的办法没有吗？
确实，三相HTS电缆的2D截面电磁-热耦合分析，一个参数组合要几小时。设计最优化要扫几百到几千个参数组合，就成了月级以上的计算。
最近注目的是机器学习代理模型。从FEM计算结果的教师数据（50～200个案例）训练神经网络或高斯过程，新参数下的交流损耗、临界电流瞬间预测。具体成功案例：
- 输入胶带数、绕线角、运行温度 → 预测交流损耗的DNN（误差 < 5%、计算时间1/10000）
- PINN（物理信息神经网络）用H-格式近似——物理律变成损失函数，教师数据少也能泛化
REBCO CORC®电缆与大电流化
REBCO电缆有新的结构吗？
最革新的是CORC（Conductor On Round Core）电缆。传统HTS电缆是在中空芯上平面螺绕胶带，CORC则在细的圆棒（直径5～10mm）上密集螺绕多条REBCO胶带。直径小的导体却能做到数千安～数万安的电流，且抗弯。
从有限元角度，CORC的胶带磁场环境比传统复杂（胶带间近距离、磁耦合强），均质化模型适用困难。多数需3D个别胶带模型，计算成本大增。
直流超导输电的展望
说好直流没损耗，为什么不直流输电呢？
完全对。直流超导输电作为终极低损输电备受关注。直流无时间变化磁场，交流损耗零（只有涟波分量的微小损耗）。冷却负荷仅来自断热层的热侵，100km长度的电缆也实用化可能。
德国SuperLink项目（320kV DC、12km）目标是2030年代施工。可再生能源大规模接入驱动长距输电需求增加，洋上风力 → 陆地的直流超导电缆被检讨。
有限元分析的重点从交流损耗转向淬火保护和断热设计。接点（接头）的接触热阻发热、电流引线的热侵入成为主要分析对象。
核融合与超导 — ITER到FCC
核融合反应堆ITER的环向磁场线圈是全球最大的超导磁体（Nb₃Sn、13T、68kA），其设计倾注了几十年的电磁-热-结构耦合有限元分析。下一代FCC（100km周长的圆形对撞机）考虑16T级Nb₃Sn二极管磁石，更先还要用HTS（REBCO）插入线圈。核融合加速器领域培养的超导有限元技术如今反过来应用于输电电缆。特别是淬火保护分析的专用代码（THEA、SUPERMAGNET等）也被移植到输电分析。
超导电缆的电磁场故障排除
非线性收敛失败
H-格式计算跑了一会儿就"Newton solver did not converge"了…
超导有限元最多见的错误。原因和对策体系看：
原因1：n值初设过高
$n = 30$ 一上来就跑会收敛失败。对策是分阶段提升n值。先 $n = 5$ 算出稳定初值，再 5→10→20→30 逐步提升。COMSOL的参数扫描可自动化。
原因2：空气阻抗率设置
H-格式中空气（绝缘域）也需有限阻抗。$\rho_{air}$ 太小则空气产生非物理涡流，太大则矩阵条件数恶化。推荐 $\rho_{air} = 1\,\Omega\cdot\text{m}$。
原因3：时间步粗
电流从零→峰值瞬间（$\omega t = 0$ 附近）磁通侵入急剧，需细的时间步。启用BDF自适应时间步，最小设为交流周期的1/10000以下。
交流损耗与实验值偏差
算出来的交流损耗和实验的卡路里值差2倍以上…… 哪里出问题了？
交流损耗的不符是超导电缆分析最常见的问题。检查清单逐项确认：
- $J_c(B)$ 磁场依赖：全区用自磁场下 $J_c$ 吗？ 断面内磁场变化0.01～0.5T，$J_c$ 变30～50%
- $J_c$ 各向异性：REBCO胶带垂直于面的磁场分量 $B_\perp$ 大幅降低 $J_c$。外层HTS尤其 $B_\perp$ 大
- 过渡周期：前1～2周是初值磁化过渡。必须3周期后评价
- 中心芯涡流：铜芯的涡流损耗包了吗？ 低温铜导电率高，这项不可忽
- 网格分辨率：HTS面表磁场陡变的侵入锋面（penetration front）有足够网格吗
淬火传播速度被低估
淬火分析做了，可传播速度实验的1/10……
HTS淬火很慢（数mm/s～数cm/s），1/10的差是大问题。主要原因：
- 电流再分配忽视：淬火域的电流转流入稳定铜的效果合理模型化了吗？ 这个向前的焦耳热驱动淬火传播
- 接触热阻：REBCO层与镍钴基板间的接触热阻考虑了吗？ 这往往是热传导的瓶颈
- 2D vs 3D：2D截面模型无法表现沿电缆轴向的热传导（淬火传播方向）。至少准3D（2D+轴向1D）才行
- 比热 $C_p(T)$：用的是Debye模型的低温比热吗？ 77K附近 $C_p$ 急变
超导电缆分析，要多物理耦合真的很关键呢。电磁场单独就完全不够……
就是，超导电缆是CAE中最强耦合的领域之一。电磁场 → 发热 → 温度升高 → $J_c$ 降低 → 电流重分配 → 磁场改变…… 这整个环形全部强耦合。所以多物理耦合有限元分析的真价值在这里充分体现。单根胶带的2D模型从Norris式能验证开始，逐步增加复杂度，这是王道。加油！