超导体的交流损失
超导体的交流损失理论基础
概要
老师！今天就是讨论超导体的交流损失，是什么呢？
超导导线·电缆在交流磁场下的损失计算。磁滞损、耦合损、涡流损的分离评估。丝径和绞距的最优化。
哦～，超导导线·电缆的话题，超级有趣！ 请多多讲解。
支配方程
我明白了...超导体的交流损失看似简单，但其实内涵很深。
离散化手法
这个方程在计算机上如何实际求解？
使用有限元法（FEM）进行空间离散化。构造单元刚度矩阵，建立全局刚度方程。
矩阵求解算法
矩阵求解算法具体是什么意思？
通过直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程组。对于大规模问题，带预处理的迭代法较为有效。
| 求解法 | 分类 | 内存使用 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小～中规模 |
| Cholesky分解 | 直接法（对称正定） | O(n²) | 小～中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说在有限元法这个地方如果缩水，后面会吃苦头呀。我要铭记在心！
商用工具中的实现
那么做超导体的交流损失要用什么样的软件呢？
| 工具名 | 开发方/现属 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Maxwell | ANSYS Inc. | .aedt, .maxwell |
| JMAG-Designer | JSOL Corporation | .jmag, .jproj |
供应商谱系和产品整合历史
各软件的成立经过很有故事性吧？
COMSOL Multiphysics
请介绍一下"COMSOL Multiphysics"！
1986年在瑞典成立。以MATLAB联动的FEMLAB开始，后来改名为COMSOL。多物理场领域有优势。
现属公司：COMSOL AB
Ansys Maxwell
请介绍一下"Ansys Maxwell"！
Ansoft Maxwell。低频电磁场解析。2008年整合入Ansys。
现属公司：ANSYS Inc.
JMAG-Designer
具体什么是JMAG？
由日本JSOL Corporation开发。电气设备设计专用的电磁场解析工具。
现属公司：JSOL Corporation
啊，原来是这样！ 年在瑞典成立的意思，我终于明白了。
文件格式和互操作性
不同软件之间交换数据时有什么注意点吗？
| 格式 | 扩展名 | 种类 | 概要 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | ISO 10303符合的3D CAD数据交换格式。形状+PMI对应。 |
| IGES | .igs/.iges | 中立CAD | 初期的CAD数据交换规范。曲面数据的互操作性有课题。迁移到STEP正在推进。 |
| VTK | .vtk/.vtu | 可视化 | Visualization Toolkit格式。ParaView等中使用。 |
不同求解器之间转换模型时，要注意单元类型的对应关系、材料模型的兼容性、荷载·边界条件的表现差异。特别是高阶单元或特殊单元（粘聚力单元、用户定义单元等）在求解器间往往无法直接转换。
我明白了…格式看似简单，但其实内涵很深。
实务注意事项
教科书上没有的"现场智慧"有吗？我想事先学到。
网格收敛性的确认、边界条件的妥当性验证、材料参数的敏感性分析特别重要。
超导体的交流损失的全局有所掌握！ 明天尝试在实务中应用。
嗯，很好！ 实际动手才是最好的学习。有疑问随时问我。
超导体的交流损失——"本应零电阻却会发热"的悖论
超导体直流电阻为零，但在交流（时变）磁场·电流下会产生"交流损失"。这是因为磁通量子（通量量子：2.07×10⁻¹⁵ Wb）侵入超导体内部并运动时，克服钉扎力所消耗的能量导致的散热。损失类型有①磁通流动损失（超过临界电流时）、②磁滞损失（即使在临界电流以下也会发生）。Norris模型（1970）给出了圆丝·带材交流损失的解析表达式，至今仍被用作CAE验证的基准式。在HTS线圈中，交流损失成为冷冻负荷的主要因素。
超导体的交流损失数值计算方法
数值方法详解
具体用什么算法求解超导体的交流损失？
我明白了。如果超导体的交流损失够好，基本上没问题？
离散化的定式化
用形状函数 $N_i$ 对未知量进行近似：
用公式表示就是这样。
基础方程的离散形式
用公式表示就是这样。
仅看公式有点模糊...代表什么？
连续体的支配方程离散化后，得到如下代数方程组：
这里$[K]$是全局刚度矩阵（或等效系统矩阵），$\{u\}$是未知节点变量向量，$\{F\}$是外力向量。
啊，原来如此！ 连续体的支配方程离散化后，就是这个工作原理。
单元技术
听说过"单元技术"，但没真正理解...
| 单元类型 | 次数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体2阶 | 二次 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中～高 | 中 |
积分方案
积分方案具体是什么意思？
听到现在，终于理解为什么单元类型这么重要了！
收敛性和稳定性
不收敛了怎么办，首先查什么？
收敛速度：二阶单元以 $O(h^2)$ 的阶数误差减小（光滑解情况）
我明白了…网格细分看似简单，但其实内涵很深。
求解器设置建议
具体用什么算法求解超导体的交流损失？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理手法 | ILU(0) or AMG | 根据问题规模而定 |
| 最大迭代次数 | 1000 | 不收敛时需要调整设置 |
| 内存模式 | In-core | 尽可能采用 |
棱单元（Nedelec单元）
专为电磁场解析的单元。自动保证切向分量连续性，消除虚假模式。3D高频解析的标准。
节点单元
用于标量位势定式化。在静磁场的标量位势法或静电场解析中有效。
FEM vs BEM（边界元法）
FEM：适应非线性材料·非均质媒质。BEM：可自然处理无限区域（开域问题）。混合型FEM-BEM也很有效。
非线性收敛（磁饱和）
B-H曲线的非线性用牛顿-拉夫逊法处理。残差基准：$||R||/||R_0|| < 10^{-4}$是一般做法。
频率域解析
在时间谐波假设下化简为定常问题。需要复数运算，但宽带特性由时间域解析获得。
时间域的时间步长
最高频率分量的1/20以下的时间步长为必要。隐式时间积分允许更大步长，但要注意精度。
超导体的交流损失实务应用
实践指南
老师，请告诉我"实践指南"！
讲解超导体的交流损失实务解析流程和注意事项。
前辈说的"超导体的交流损失一定要好好做"的意思，我现在明白了。
解析流程
从第一步开始教我！ 要从哪里开始？
1. 预处理 (Pre-processing)
- CAD数据导入和形状简化
- 材料特性定义
- 网格生成（单元类型·尺寸的决定）
- 边界条件和荷载条件的设置
2. 求解 (Solving)
- 求解器设置（求解法、收敛基准、输出控制）
- 作业投入和计算实行
- 收敛监控
3. 后处理 (Post-processing)
- 结果可视化（位移、应力、其他物理量）
- 结果验证和妥当性确认
- 报告制作
网格生成最佳实践
网格的好坏怎么判断？
单元品质指标
请介绍"单元品质指标"！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 纵横比 | 1.0 | < 5.0 | 精度下降 |
| 雅可比比 | 1.0 | > 0.3 | 单元退化 |
| 翘曲 | 0° | < 15° | 精度下降 |
| 斜度 | 0° | < 45° | 收敛性恶化 |
| 锥形比 | 0 | < 0.5 | 精度下降 |
网格密度的决定
网格密度的决定具体是什么意思？
边界条件设置指南
我听说边界条件弄错了，全部都废…
啊，原来如此！ 过度约束要注意，就是这么回事。
按商用工具的实现步骤
有好多种软件吧？ 各种的特性请教！
| 工具名 | 开发方/现属 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Maxwell | ANSYS Inc. | .aedt, .maxwell |
| JMAG-Designer | JSOL Corporation | .jmag, .jproj |
COMSOL Multiphysics
请介绍一下"COMSOL Multiphysics"！
1986年在瑞典成立。以MATLAB联动的FEMLAB开始，后来改名为COMSOL。多物理场领域有优势。
现属公司：COMSOL AB
Ansys Maxwell
请介绍一下"Ansys Maxwell"！
Ansoft Maxwell。低频电磁场解析。2008年整合入Ansys。
现属公司：ANSYS Inc.
老师的说明很易懂！ 工具名的雾气散去了。
常见失败和对策
初学者常范哪种错？ 想事先知道！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、不妥当的边界条件 | 网格改善、约束条件见直 |
| 应力异常大 | 应力奇异点、网格依存 | 奇异点回避、局部网格细分 |
| 位移不现实 | 材料常数误差、单位系混乱 | 输入数据确认 |
| 计算时间过长 | 不必要的细分、非效率求解 | 网格最优化、并行计算 |
质量保证检查清单
教科书上没有的"现场智慧"有吗？
超导体的交流损失的全景我掌握了！ 明天试试在实务中应用。
嗯，很好！ 实际动手才是最好的学习。有疑问随时问我。
"MRI用HTS线圈的交流损失"——磁场脉冲左右设计
MRI装置的超导磁铁以永久电流模式运行，直流损失为零，但梯度线圈动作时的急速磁场变化会在HTS中诱导交流损失。梯度线圈的脉冲磁场（dB/dt=100 T/s以上）会在超导线圈中产生数W的交流损失，如果超过冷却机容量，液体氦的蒸发会引发失超（Quench）风险。在CAE中，需要输入实际MRI摄像序列模拟的磁场波形，进行HTS的E-J曲线和热传导的"磁场-热耦合解析"，计算交流损失的时间序列。这个解析已成为高性能MRI开发的必须工程。
超导体的交流损失软件比较
商用工具比较
有好多种软件吧？ 各种的特性请教！
介绍超导体的交流损失对应的主要商用CAE工具的功能比较，以及各产品的历史背景。
我明白了。如果超导体的交流损失够好，基本上没问题？
支持工具列表
那么做超导体的交流损失要用什么样的软件呢？
| 工具名 | 开发方/现属 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Maxwell | ANSYS Inc. | .aedt, .maxwell |
| JMAG-Designer | JSOL Corporation | .jmag, .jproj |
COMSOL Multiphysics
请介绍一下"COMSOL Multiphysics"！
1986年在瑞典成立。以MATLAB联动的FEMLAB开始，后来改名为COMSOL。多物理场领域有优势。
现属公司：COMSOL AB
Ansys Maxwell
请介绍一下"Ansys Maxwell"！
Ansoft Maxwell。低频电磁场解析。2008年整合入Ansys。
现属公司：ANSYS Inc.
JMAG-Designer
具体什么是JMAG？
由日本JSOL Corporation开发。电气设备设计专用的电磁场解析工具。
现属公司：JSOL Corporation
我明白了…年在瑞典成立看似简单，但其实内涵很深。
功能比较矩阵
预算和时间都有限，哪个性价比最高？
| 功能 | COMSOL | Maxwell | JMAG |
|---|---|---|---|
| 基本功能 | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ |
| GPU对应 | △ | △ | ○ |
转换时的风险
转换时的风险具体是什么意思？
啊，原来如此！ 不同工具间的转换就是这么回事。
许可证形式
听说过"许可证形式"，但没真正理解...
| 工具 | 许可证 | 特征 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高价但有官方支持 |
| OpenFOAM | GPL | 免费但支持收费 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的开源求解器 |
选择指南
最后，怎么选，请给判断标准。
在超导体的交流损失工具选择中，考虑以下因素：
超导体的交流损失的全景我掌握了！ 明天试试在实务中应用。
嗯，很好！ 实际动手才是最好的学习。有疑问随时问我。
HTS交流损失解析工具——COMSOL vs 专用Toolbox
HTS交流损失解析的主要工具是COMSOL（AC/DC模块+HTS模型）和开源专用代码。COMSOL可以将H定式化、A定式化、T-Ω定式化作为自定义PDE实现，GUI操作也方便。研究社区可以免费使用CoilCAD（Cambridge大）、MagDyn（CERN）等专用HTS模拟器。商用方面，Infolotica（现Mentor）的Flux HTS具有交流损失计算的专门功能。最近，采用FEniCS和PETSc组合的开源HTS解析框架正在精力充沛地开发中，在大型线圈解析的应用中正在拓展。
超导体的交流损失先进研究
前沿话题和研究趋势
超导体的交流损失领域今后会怎样发展？
来看看超导体的交流损失领域的最新研究动向和先进方法。
我明白了。如果超导体的交流损失够好，基本上没问题？
最新的数值方法
接下来是最新的数值方法。什么内容？
仅看公式有点模糊...代表什么？
高性能计算 (HPC) 的应对
| 并行化手法 | 概要 | 适用求解器 |
|---|---|---|
| MPI (区域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU活用。特别在显式法中有效 | LS-DYNA, Fluent等 |
| 混合型 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
超导体的交流损失故障排除
故障排除
我明白了。如果超导体的交流损失够好，基本上没问题？
常见错误和对策
老师也为超导体的交流损失通宵调试过吗？（笑）
1. 收敛失败
收敛失败具体是什么意思？
症状：求解器在指定迭代次数内未收敛而异常终止
可能原因：
- 网格品质不足（过度扭曲的单元）
- 材料参数设置不妥当
- 初始条件不妥当
- 非线性性太强（荷载步数不足）
对策：
- 进行网格品质检查（纵横比、雅可比比）
- 确认材料参数的单位系
- 荷载分多步骤（增加子步数）
- 放松收敛判定基准（但要注意精度）
也就是说在收敛失败这地方如果缩水，后面会吃苦头。肝心要记！
2. 非物理的结果
接下来是非物理的结果的话题。什么内容？
症状：应力/位移/温度等数值不现实
可能原因：
- 边界条件误设定
- 单位系混用（SI单位与工程单位混淆）
- 单元类型选择不妥当
- 应力奇异点存在
对策：
- 确认反力合计（力的平衡）
- 确认单位系一贯性
- 重新考虑单元类型的适切性
- 特异点消除或局部细分模型
前辈说的"收敛失败一定要好好做"的意思，我现在明白了。
3. 计算时间超过
计算时间超过具体是什么意思？
症状：计算耗时数倍于预想
对策：
- 网格粗密分布的最优化
- 对称性的活用（1/2、1/4模型）
- 求解器设置的最优化（迭代法、预处理的选择）
- 并行计算的活用
4. 内存不足
请介绍一下"内存不足"！
症状：Out of Memory 错误
前辈说的"收敛失败一定要好好做"的意思，我现在明白了。