超伝導ケーブルの電磁界シミュレーション

理論と物理

概要 — 超伝導送電の現在地

超伝導ケーブルって実際に電力網で使われてるんですか？ 未来の技術ってイメージが強いんですけど…

いや、もう実用段階に入っている。代表例を挙げると、ドイツ・エッセン市のAmpaCity プロジェクトは10kV・40MVA・約1kmのHTSケーブルを2014年から商用系統に接続して連続運転している。韓国のKEPCOも23kV系統で500m級の実証を完了した。日本でも住友電工が66kV・200m級の実証試験を行っている。

え、もうそんなに進んでるんですか！ でも普通の銅ケーブルと比べて何がそんなに有利なんですか？

一番のポイントは電流密度だ。HTS（高温超伝導）線材は銅の5〜10倍の電流密度で通電できる。これは都市部の地下送電トンネルのスペース節約に直結する。例えば同じ40MVAの送電容量を確保するのに、従来の110kV OFケーブルなら3回線必要なところが、超伝導なら10kV・1回線で済む。変圧器も省略できて設備全体がコンパクトになる。

なるほど！ でも液体窒素で冷やさないといけないんですよね？ その電力はどうなるんですか？

いい質問だ。冷却に必要な電力は典型的にケーブル送電容量の0.5〜1%程度。一方、従来の銅ケーブルの抵抗損失は3〜5%だから、冷却コストを差し引いてもネットで省エネになる。ただしここが設計の核心で、FEM解析ではHTS線材のAC損失と冷却負荷のバランスを精密に評価する必要がある。そのために必要な物理モデルを順に見ていこう。

E-J power law — 超伝導体の構成則

超伝導体って抵抗ゼロなんですよね？ それなのにオームの法則みたいな構成則が要るんですか？

「抵抗ゼロ」というのは理想化された話で、実際のHTS線材は臨界電流密度 $J_c$ 付近で非常に急峻な非線形抵抗を示す。これを記述するのがE-J power lawだ：

ここで $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$ が爆発的に増大する。これは数値的には極端に非線形な構成則を意味する。FEMソルバーにとっては、見かけの抵抗率 $\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$ は2倍以上に増える。逆に3Tの外部磁場下では半分以下に低下する。

じゃあケーブルの設計って、温度と磁場のバランスをどう取るかがキモってことですか？

まさにそう。ケーブル断面内の磁場分布はケーブル自身が作る自己磁場で決まるし、温度分布はAC損失と冷却のバランスで決まる。つまり $J_c$ を決めるには磁場と温度が必要、磁場を決めるには電流（つまり $J_c$ に制約される）が必要、温度を決めるにはAC損失（電流と磁場の関数）が必要——すべてが連成している。だから電磁-熱の連成FEM解析が不可欠なんだ。

AC損失の3成分

「超伝導なのに損失がある」って矛盾してませんか？ 抵抗ゼロなのに発熱するって…

直流なら確かに損失はほぼゼロだ。でも交流の場合、磁束が時間変化するから3種類の損失メカニズムが生じる。これをまとめてAC損失と呼ぶ。超伝導ケーブルの設計において、AC損失の正確な評価は冷却系の設計に直結するから、最も重要な計算項目の一つだ。

3種類って何ですか？ それぞれどのくらいの大きさなんですか？

順番に説明しよう。

1. ヒステリシス損失 $Q_h$ — 超伝導フィラメント内の磁束ピニングに起因する損失で、AC損失の主成分。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倍導電率が高いから、渦電流損失が予想以上に大きくなることがある。特にケーブルのフォーマーが銅製の場合、ここの損失が無視できないケースがある。

全部足し合わせたAC損失が冷却負荷になるわけですね。具体的にはどのくらいの大きさになるんですか？

AmpaCity級の10kV・2.3kA HTSケーブルで、AC損失は1〜3 W/m程度。これを77Kで除去するための冷凍機の所要電力は、カルノー効率を考慮すると損失の10〜15倍（COP ≈ 0.07〜0.1）だから、電気入力換算で10〜45 W/m。1kmのケーブルで冷却系全体が10〜45 kW相当になる。これが銅ケーブルの送電損失（同容量で100〜200 W/m）より十分小さいから経済性が成り立つ。

電磁-熱連成とクエンチ

AC損失で温度が上がると $J_c$ が下がって…って、悪循環になりませんか？

まさにその通り。これがクエンチ（超伝導状態の不可逆的な喪失）の本質だ。電磁-熱連成の支配方程式を書くと：

右辺第2項 $\mathbf{E} \cdot \mathbf{J}$ がAC損失（+常伝導転移時のジュール発熱）、第3項 $q_{cool}$ が液体窒素への冷却。E-J power lawから $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と比較して極めて遅い）。これが問題で、検出が遅れると局所的に数百度に達して線材が焼損する。だからクエンチの検出・伝播解析は安全設計の最重要課題だ。FEM連成解析でクエンチの発生条件と伝播速度を予測することが、保護回路の設計に不可欠になる。

数値解法と実装

FEM定式化 — A-V法とH-formulation

超伝導ケーブルのFEM解析って、普通の電磁界解析と何が違うんですか？ Maxwell方程式を解くだけじゃダメなんですか？

Maxwell方程式を解くこと自体は同じだけど、超伝導体の極端な非線形E-J特性が数値的な難しさを生む。定式化の選択が計算の成否を左右する。主に2つのアプローチがある。

A-V法（磁気ベクトルポテンシャル法）は電磁界解析の標準的手法で、$\mathbf{B} = \nabla \times \mathbf{A}$、$\mathbf{E} = -\partial\mathbf{A}/\partial t - \nabla V$ として磁気ベクトルポテンシャル $\mathbf{A}$ と電気スカラーポテンシャル $V$ を未知数とする：

ただし超伝導体領域では $\sigma$ の代わりにE-J power lawに基づく非線形抵抗率 $\rho_{sc}(J)$ を使う。問題は、この $\rho_{sc}$ が $J$ によって10桁以上変化するため、ヤコビアン行列の条件数が極端に悪化すること。

これに対して、近年超伝導解析で主流になっているのがH-formulationだ：

磁場 $\mathbf{H}$ を直接未知数にするんですか？ なんでそっちのほうがいいんですか？

H-formulationの利点は3つある。第一に、電流密度が $\mathbf{J} = \nabla \times \mathbf{H}$ として自然に得られるため、E-J特性の評価が直接的。第二に、超伝導体の外側（空気領域）に非常に大きな抵抗率（$\rho_{air} \sim 1\,\Omega\cdot\text{m}$）を設定するだけで磁場の侵入を自然に扱える。第三に、2D断面解析ではスカラー変数1つ（$H_z$）になるため自由度が少なく計算効率が良い。

COMSOLのHTSモデリングガイドでもH-formulationが推奨されており、学術論文の大多数がこのアプローチを採用している。

メッシュ戦略 — 薄膜近似とマルチスケール

HTS線材って厚さ1μmのYBCO層があるんですよね？ ケーブル全体は直径100mmとかでしょう？ それを同じメッシュで切れるんですか？

鋭い指摘だ。REBCO線材は典型的に幅4〜12mm、総厚0.1mmで、その中の超伝導YBCO層はわずか1〜2μm。ケーブル外径が100mmだとすると、スケール比は1:100,000。これを直接メッシュで分割するのは非現実的だ。

解決策として薄膜近似（thin shell approximation）を使う。YBCO層を体積要素ではなく面要素として扱い、面内のシート電流密度 $K_s\,[\text{A/m}]$ を自由度とする：

$d_{sc}$ は超伝導層の厚さ。この近似により、1μmの層を厚さ方向に分割する必要がなくなり、メッシュ数を2〜3桁削減できる。COMSOLではこれを「Shell機能」として実装している。ただし超伝導層内の電流密度の厚さ方向分布は捨てることになるから、層厚に対して磁場侵入深さが十分大きい条件（フルペネトレーション）でのみ有効だ。

それでもケーブルに線材が数十本〜数百本巻いてあるんですよね？ 全部モデル化するのは大変そう…

その通り。実際のHTSケーブルは複数層にHTS線材がヘリカルに巻かれ、各層に20〜40本のテープがある。全テープを3D個別にモデル化するのは計算コスト的に非現実的だ。そこで均質化（ホモジナイゼーション）が重要になる。HTS線材層を異方性の等価連続体として扱い、有効な $J_c$、有効抵抗率テンソルを定義する。これについては後のセクションで詳しく説明する。

時間積分と非線形ソルバー

E-J power lawの $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 power lawのヤコビアン（接線抵抗率テンソル）は：

これを数値微分ではなく解析的に与えることで、Newton-Raphsonの収束が大幅に改善する。COMSOLではこれを「Analytical Jacobian contribution」として設定できる。逆にこれを省略すると収束に2〜3倍の反復回数がかかることがある。

均質化モデル — 大規模ケーブルへの適用

さっきの均質化モデルの話を詳しく教えてください。実際にどうやって「等価的な材料」を作るんですか？

ケーブルの1つのHTS線材層を、円筒シェルの異方性連続体として扱う。巻き線方向（ヘリカル方向）の有効臨界電流密度は：

$f_{fill}$ はフィルファクター（その層でHTS線材が占める周方向の充填率、典型的に0.7〜0.9）、$\theta_w$ は巻き角（ケーブル軸に対するヘリカル角、典型的に15°〜35°）。垂直方向（周方向）の電流は流れないと仮定し、異方性抵抗率テンソルとして実装する。

この均質化モデルを使うと、線材1本1本をモデル化するのに比べて自由度が1/100〜1/1000に減り、商用電力ケーブル全長（数百m〜数km）スケールの解析が可能になる。精度は個別テープモデルと比較して電流分布で5%以内、AC損失で10%以内の誤差に収まることが多くの文献で報告されている。

実践ガイド

解析ワークフロー

超伝導ケーブルの電磁界解析って、実際にはどういう手順で進めるんですか？ 最初の一歩から教えてください。

典型的なワークフローはこうだ：

Step 1: ケーブル断面の幾何モデル構築

• クライオスタット外管（ステンレス）→ 断熱真空層 → クライオスタット内管 → LN₂流路 → HTS線材層（複数層）→ フォーマー（中心導体）の同心円構造を構築
• 各HTS層のテープ枚数・巻き角・ギャップを定義
• 2D断面モデルが基本。3Dが必要なのはケーブル端末部（ターミネーション）の解析時

Step 2: 材料特性の定義

• HTS層: E-J power law（$E_0$, $J_c(B,T)$, $n$ 値）を温度・磁場テーブルで定義
• 安定化銅: $\rho_{cu}(T) = 1.7 \times 10^{-8} \times (T/293)$ [Ω·m]（低温での残留抵抗比RRRも考慮）
• LN₂: 熱伝達係数 $h(T)$ を沸騰曲線から取得（核沸騰、膜沸騰の区別）
• 断熱層: $\kappa_{eff} \sim 10^{-4}$ W/(m·K)（多層断熱MLI）

Step 3: メッシュ生成

• HTS層: 薄膜近似を適用する場合はシェル要素、体積モデルなら厚さ方向に最低3層
• 空気領域: ケーブル外径の5〜10倍まで。外周に無限要素（infinite element）を配置
• LN₂流路: 流体解析を連成する場合はプリズム要素のバウンダリーレイヤー

Step 4: ソルバー設定・計算実行

• H-formulation: BDF2、初期時間刻み $10^{-5}$s、最大 $10^{-3}$s
• Newton-Raphson: 相対収束判定 $10^{-4}$〜$10^{-6}$、ダンピング係数0.5〜0.9
• 交流定常状態に達するまで3〜5サイクルは回す（過渡項の減衰待ち）

Step 5: 後処理・AC損失算出

• 定常サイクル（最後の1サイクル）からAC損失を体積積分で算出：$Q = \int_V \int_0^{1/f} \mathbf{E} \cdot \mathbf{J}\,dt\,dV$
• 各層別の損失内訳を分離（HTS層のヒステリシス損失、フォーマーの渦電流損失など）
• 電流分布の均一性（層間電流バランス）を確認

境界条件の設定

H-formulationでの境界条件って、A-V法とは違うんですか？

うん、H-formulationでは以下の境界条件を使う：

電流を印加するのにアンペールの周回積分を使うんですか！ 面白い発想ですね。

そう。H-formulationでは電流密度が直接の未知数ではなく、$\mathbf{J} = \nabla \times \mathbf{H}$ として間接的に決まる。だからケーブルに流す総電流を指定するには、導体を囲む閉路上の $\mathbf{H}$ の線積分を拘束する。COMSOLでは「Coil Group」機能でこれを実現できる。三相ケーブルなら位相差120°の3つの電流を各相に印加する。

検証と妥当性確認

解析結果が正しいかどうか、どうやって確認すればいいですか？ 超伝導ケーブルの実験データってそんなに多くないですよね？

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水準以上でAC損失の収束を確認（誤差 < 2%）
• 時間刻み依存性: 時間刻みを半分にして結果が変わらないことを確認

妥当性確認（Validation）— 計算は物理を正しく再現しているか

• 短尺サンプルのAC損失カロリメトリー測定との比較
• 臨界電流の実測値（自己磁場下、外部磁場下）とモデル予測の比較
• クエンチ電圧の実測波形との比較（あれば）
• HTS-CICC（Cable-In-Conduit Conductor）ではCEA-Cadarache、CERNなどの公開ベンチマークデータが利用可能

よくある失敗パターン

初心者がやりがちなミスってありますか？ 先に知っておきたいです。

実務で頻出する失敗パターンをまとめよう：

ソフトウェア比較

COMSOL — AC/DCモジュール

超伝導ケーブルの解析で一番使われているソフトって何ですか？

学術論文の圧倒的多数がCOMSOL Multiphysicsを使っている。理由は明快で、H-formulationをPDE（偏微分方程式）モジュールで直接定義できるため、超伝導体の非線形構成則を柔軟に実装できる。AC/DCモジュールの標準機能だけでは不十分で、PDEモジュールとの組み合わせが本領を発揮する。

• 強み: マルチフィジクス連成（電磁-熱-流体）がワンツールで完結、PDEの自由な定義、薄膜近似のShell機能
• 弱み: 大規模3D（テープ100本以上）ではメモリ消費が大きい、並列スケーラビリティに限界
• ライセンス: AC/DC + Heat Transfer + CFD + PDE で年間約200万円〜

Ansys Maxwell + Mechanical

Ansysでも超伝導ケーブルの解析はできるんですか？

Ansys MaxwellはA-V法ベースの電磁界解析ソルバーで、非線形B-H曲線の扱いに長けている。ただし、E-J power lawのような電流依存の非線形抵抗率をユーザー定義するのはCOMSOLほど容易ではない。超伝導解析ではカスタムUDFの実装が必要になることが多い。

• 強み: 電磁界解析の精度が高い、渦電流解析の適応的メッシュリファインメント、Workbenchを通じた電磁-熱-構造連成
• 弱み: H-formulationの標準サポートがない、超伝導固有のE-J特性実装に工夫が必要
• 適用例: 超伝導マグネット（MRI、核融合）のクエンチ解析ではMaxwell + Mechanicalの連成が実績多数

オープンソース（GetDP, Elmer, FEniCS）

予算が限られた研究室でも使えるオープンソースのツールはありますか？

超伝導のFEM解析で実績のあるOSSは3つある：

• GetDP（リエージュ大学）: 超伝導研究コミュニティで最も使われているOSS。H-formulation対応、Gmshと統合。CERNの大型マグネット解析にも使用実績あり
• Elmer（CSC, フィンランド）: マルチフィジクスOSSソルバー。超伝導モジュールは限定的だが、カスタム定式化が可能
• FEniCS/FEniCSx: Python FEMフレームワーク。H-formulationを数行で定義でき、研究プロトタイピングに最適。ただし産業規模のメッシング・後処理は別ツールが必要

機能比較マトリクス

結局どれを選べばいいか、一覧で比較できますか？

やっぱりCOMSOLが一強なんですね。でもGetDPも論文で結構見かけますよね？

そう。リエージュ大学のグループ（B. Ramdane, F. Grilli等）がGetDPで精力的に超伝導解析の研究を進めていて、ベンチマーク問題の公開もしている。CERNのFCC（Future Circular Collider）用Nb₃Sn超伝導マグネットの設計でもGetDPが使われている。予算に制約がある研究室なら、まずGetDP + Gmshで始めるのが現実的だろう。

先端技術

機械学習サロゲートモデル

超伝導ケーブルの解析って1回の計算に何時間もかかるんですよね？ もっと速くする方法はないんですか？

確かに、三相HTSケーブルの2D断面の電磁-熱連成解析で1パラメータセットあたり数時間かかることがある。設計最適化で100〜1000通りのパラメータスイープをすると数ヶ月の計算時間になる。

最近注目されているのが機械学習サロゲートモデルだ。FEM計算結果の教師データ（50〜200ケース）からニューラルネットワークやGaussian Processを学習させ、新しいパラメータセットでのAC損失や臨界電流を瞬時に予測する。具体的な成功例としては：

• テープ枚数・巻き角・運転温度を入力 → AC損失を予測するDNN（誤差 < 5%、計算時間1/10,000）
• PINN（Physics-Informed Neural Network）によるH-formulation近似 — 物理法則を損失関数に組み込むことで少ない教師データで汎化

REBCO CORC®ケーブルと大電流化

REBCOケーブルの最新構造って何か新しいものがあるんですか？

最も革新的なのがCORC（Conductor On Round Core）ケーブルだ。従来のHTSケーブルは中空フォーマーに平面的にテープを巻くが、CORCは細い丸棒に複数のREBCOテープをヘリカルに密巻きする。直径5〜10mmの導体で数kA〜数十kAの電流を実現でき、曲げにも強い。

FEM解析の観点からは、CORCは従来ケーブルよりも個別テープの磁場環境が複雑（テープが互いに近接し、磁気的に強く相互作用する）で、均質化モデルの適用が難しい。3Dの個別テープモデルが必要になるケースが多く、計算コストの大幅な増加が課題だ。

直流超伝導送電の展望

さっきAC損失の話で「直流なら損失ゼロ」って言ってましたよね？ じゃあ直流で送電すればいいのでは？

その通りで、直流超伝導送電は究極の低損失送電として注目されている。直流なら時間変動磁場がないからAC損失はゼロ（厳密にはリップル成分による微小損失のみ）。冷却負荷は断熱層からの熱侵入だけで、ケーブル長100kmでも実用的になる可能性がある。

ドイツのSuperLinkプロジェクト（12km、320kV DC）が2030年代の建設を目指している。再生可能エネルギーの大量導入で長距離送電の需要が増しており、洋上風力発電→陸上の送電にDC超伝導ケーブルが検討されている。

FEM解析の観点では、DC超伝導ケーブルの設計課題はクエンチ保護と断熱設計に集中する。AC損失ではなく、接続部（ジョイント）の抵抗発熱や電流リード部からの熱侵入がメインの解析対象になる。

トラブルシューティング

非線形収束の失敗

H-formulationで計算を走らせたら「Newton solver did not converge」って出て止まりました…

超伝導FEMで最頻出のエラーだ。原因と対策を系統的に見ていこう：

原因1: n値が高すぎる初期設定

$n = 30$ でいきなり計算を始めると収束しないことが多い。対策として段階的n値増加法を使う。まず $n = 5$ で計算して安定した初期解を得てから、$n$ を5→10→20→30と段階的に上げていく。COMSOLではパラメトリックスイープでこれを自動化できる。

原因2: 空気領域の抵抗率設定

H-formulationでは空気（非導電領域）にも有限の抵抗率を設定する必要がある。$\rho_{air}$ が小さすぎると空気中に非物理的な渦電流が流れ、大きすぎると行列の条件数が悪化する。推奨値は $\rho_{air} = 1\,\Omega\cdot\text{m}$。

原因3: 時間刻みが粗すぎる

電流がゼロ→ピークに変化する瞬間（$\omega t = 0$ 付近）は磁束侵入が急激に進むため、時間刻みを細かくする必要がある。BDF法の適応的時間刻みを有効にし、最小時間刻みをACサイクルの1/10,000以下に設定する。

AC損失の実験値との乖離

AC損失を計算したんですけど、実験のカロリメトリー値と2倍以上ズレてます。何が間違ってるんでしょう？

AC損失の不一致は超伝導ケーブル解析で最もよく遭遇する問題だ。チェックリストを順に確認してくれ：

• $J_c(B)$ の磁場依存性: 自己磁場下の $J_c$ を全域で使っていないか？ ケーブル断面内の磁場は場所によって0.01〜0.5T程度変化し、$J_c$ が30〜50%変わる
• $J_c$ の異方性: REBCOテープ面に垂直な磁場成分 $B_\perp$ が $J_c$ を大きく下げる。ケーブルの外側層ほど $B_\perp$ が大きい
• 過渡サイクル: 最初の1〜2サイクルは初期磁化の過渡効果が含まれる。必ず3サイクル目以降で評価する
• フォーマー渦電流: 銅フォーマーの渦電流損失を含めているか？ 低温銅は導電率が高いため、この成分が無視できないことがある
• メッシュ解像度: HTS層表面の磁場が急変する領域（侵入フロント）に十分なメッシュ密度があるか

クエンチ伝播速度の過小評価

クエンチ解析をしたら、伝播速度が実験の1/10しか出ないんですけど…

HTSのクエンチ伝播は非常に遅い（典型的に数mm/s〜数cm/s）から、1/10のズレは大きな問題だ。主な原因：

• 電流再分配の無視: クエンチ領域の電流が安定化銅層に転流する効果を正しくモデル化しているか？ この電流再分配による前方のジュール発熱がクエンチ伝播を駆動する
• 接触熱抵抗: REBCO層と基板（ハステロイ）の間の接触熱抵抗を考慮しているか？ これが熱伝導を律速する場合がある
• 2D vs 3D効果: 2D断面モデルではケーブル軸方向の熱伝導（クエンチ伝播方向）が表現できない。少なくとも準3D（2D+軸方向1D）のモデルが必要
• 比熱 $C_p(T)$: デバイモデルに基づく低温比熱を使っているか？ 77K近辺では $C_p$ が急激に変化する

超伝導ケーブルの解析って、本当に多物理の連成が命なんですね。電磁気だけでは全然ダメなんだ…

そう、超伝導ケーブルはCAEの中でも最も連成が強い分野の一つだ。電磁場→発熱→温度変化→$J_c$変化→電流再分配→磁場変化…というループが全て密結合している。だからこそFEMのマルチフィジクス連成解析の真価が問われる、挑戦しがいのある分野だよ。まずはNorris式で検証できる単一テープの2Dモデルから始めて、段階的に複雑化していくのが王道だ。