電池熱暴走シミュレーション
理論と物理
バッテリーの熱暴走って爆発するやつですか？ シミュレーションで予測できるんですか？
ざっくり言うと、リチウムイオン電池は内部温度が約130°Cを超えるとSEI（固体電解質界面）層が分解し始める。するとむき出しになった負極リチウムが電解液と直接反応して発熱し、さらに温度が上がる。200°C以上になると正極からの酸素放出と電解液の発熱反応が暴走して、数秒で800°Cを超えることもある。
え、数秒で800°C!? それが車に積んであるバッテリーで起きたらヤバいですよね…
そのとおり。1セルの熱暴走が隣接セルに熱伝播すると、モジュール全体→パック全体と連鎖的に火災が広がる。これが「延焼（thermal propagation）」で、EV火災の主因だ。だからこそCAEで伝播時間を予測して、断熱バリアの厚さや冷却設計を最適化するのが重要になる。
熱暴走のメカニズム
もう少し詳しく、温度が上がっていく「段階」を教えてもらえますか？
温度上昇は4段階の連鎖反応で進行する。シミュレーションモデルでもこの4段階をそれぞれArrhenius式で表現するのが標準的だよ。
| 段階 | 温度域 | 反応 | 発熱量の目安 |
|---|---|---|---|
| Stage 1 | 80〜120°C | SEI膜の分解 | ~250 J/g |
| Stage 2 | 120〜200°C | 負極-電解液反応（インターカレートLiの露出） | ~1700 J/g |
| Stage 3 | 180〜250°C | 正極分解（NMC, LFPで大きく異なる） | ~800 J/g |
| Stage 4 | 200°C〜 | 電解液分解（有機溶媒の気化・燃焼） | ~400 J/g |
つまり、SEI層の分解が引き金になって次々と反応が連鎖するんですね。LFPはNMCより安全って聞くのは、Stage 3の正極分解が穏やかだから？
そのとおり。LiFePO₄（LFP）はFe-O結合が強くて酸素放出温度が高い。だから正極起因の発熱が少なく、熱暴走に至りにくい。ただし「絶対に暴走しない」わけではなくて、外部短絡や過充電で強制的にStage 2以降に進むケースはある。シミュレーションでは正極材ごとにArrhenius定数を変えてモデル化する。
支配方程式
熱暴走シミュレーションの支配方程式はどうなるんですか？
基本はエネルギー保存式（非定常熱伝導方程式）に、化学反応による発熱項とジュール発熱項を加えたものだ：
ここで $\rho$ は密度、$c_p$ は比熱、$k$ は熱伝導率。右辺第1項が熱伝導、$Q_{\text{rxn}}$ が化学反応発熱、$Q_{\text{ohm}}$ がジュール発熱だ。
Arrhenius反応速度モデル
$Q_{\text{rxn}}$ の中身はどうなっているんですか？ Arrhenius式ってよく聞くけど…
各段階 $i$ の反応速度をArrhenius式で記述して、反応熱 $H_i$ を掛けて合算する。これが $Q_{\text{rxn}}$ だ：
各パラメータの意味はこうだ：
- $A_i$：頻度因子（プリエクスポネンシャル因子）[1/s]
- $E_{a,i}$：活性化エネルギー [J/mol]
- $R$：気体定数 = 8.314 J/(mol・K)
- $c_i$：反応物の無次元濃度（0〜1、消費されると減少）
- $n_i$：反応次数
- $H_i$：反応エンタルピー [J/m³]
$\exp(-E_a/RT)$ の部分がミソなんですね。温度 $T$ が上がると指数的に反応速度が上がって、発熱が増えて、さらに温度が上がる… まさに暴走ループ！
そう、それがSemenov理論で言う「熱爆発」の本質だ。放熱が追いつかなくなった瞬間に温度がポジティブフィードバックで急上昇する。
| 段階 | $E_a$ [kJ/mol] | $A$ [1/s] | $H$ [J/g] |
|---|---|---|---|
| SEI分解 | 135 | 1.67×10¹⁵ | 257 |
| 負極-電解液 | 110〜140 | 2.50×10¹³ | 1714 |
| 正極分解（NMC） | 170〜200 | 1.75×10⁹ | 792 |
| 電解液分解 | 140〜160 | 5.14×10²⁵ | 420 |
これらの値はDSC（示差走査熱量計）やARC（加速速度熱量計）の実験データからフィッティングして求めるのが一般的だ。
電池の発熱モデル
$Q_{\text{ohm}}$ のジュール発熱のほうはどうなるんですか？ 通常運転時の発熱と暴走時の発熱は別物？
通常運転時の発熱モデルはBernardi式と呼ばれる。これも熱暴走シミュレーションの初期条件として重要だ：
- 第1項 $I^2 R_{\text{int}}$：ジュール発熱（抵抗発熱）。常に正で、大電流ほど大きい。内部抵抗 $R_{\text{int}}$ はSOCと温度の関数。
- 第2項 $I \cdot T \cdot dU_{\text{OCV}}/dT$：エントロピー発熱（可逆熱）。充放電で符号が変わる。$dU_{\text{OCV}}/dT$ はOCV（開回路電圧）の温度係数。
例えば急速充電（2C〜3C）だと第1項のジュール発熱が支配的で、セル温度が60〜80°Cまで上昇する。冷却が不十分だとそこからSEI分解の引き金温度に近づいてしまう。
なるほど、通常運転時の発熱モデルと暴走時のArrhenius反応モデルを両方組み合わせるんですね。
セル間熱伝播モデル
1セルが暴走したあと、隣のセルに火が移るのをどうモデル化するんですか？
セル間の伝播モデルでは、3つの伝熱モードをすべて考慮する：
- 熱伝導 $q_{\text{cond}}$：セル間の直接接触やバリア材を介した伝導。最も支配的。バリア材の熱伝導率 $k$ と厚さ $d$ で制御する。
- 対流 $q_{\text{conv}}$：高温ガスのベント（噴出）によるモジュール内の強制対流。ベントガス温度は600〜1000°Cに達する。
- 輻射 $q_{\text{rad}}$：$\sigma T^4$ に比例するステファン・ボルツマン輻射。800°C超で無視できなくなる。
実務では「トリガーセル」（釘刺しや内部短絡）を1つ設定し、そこから隣接セルの温度が暴走開始温度を超えるまでの時間を「伝播時間」として評価する。この伝播時間が乗員の脱出時間（通常5分以上）を確保できるかどうかが設計基準だ。
バリア材の厚さをちょっと変えるだけで伝播時間が大きく変わりそうですね。まさにシミュレーション向きの問題だ。
そう。ある自動車メーカーの事例では、セラミックファイバー系バリアを5mmから8mmに変えるだけで、伝播時間が約40秒から120秒以上に延びた。その80秒の差が乗員の安全を左右する。こういうパラメトリック解析はFEMの得意技だ。
熱暴走の「引き金」——130°Cで何が起きるのか
リチウムイオン電池の熱暴走を理解する鍵は、温度と反応の「連鎖」です。SEI膜が分解し始める80〜120°C、セパレータが溶融して内部短絡が生じる130〜150°C、正極からの酸素放出が起きる180〜200°C——これらの発熱反応が次々とトリガーし合い、温度は数秒で800°Cを超えることもあります。しかし発熱のタイミングは電池の製造ばらつきや劣化状態に強く依存するため、「どの電池でも同じ温度で暴走が始まる」わけではありません。この個体差のモデル化が、熱暴走シミュレーションで最も難しい課題の一つです。
各項の物理的意味
- 蓄熱項 $\rho c_p \partial T / \partial t$：電池材料（正極・負極・電解液・集電体）の熱慣性を表す。ジェリーロール構造では異方性を持つため、面方向と厚さ方向で異なる実効熱容量を設定する場合がある。
- 熱伝導項 $\nabla \cdot (k \nabla T)$：円筒セルの面方向（ $k_r \approx$ 30 W/(m・K)）と径方向（$k_z \approx$ 0.7 W/(m・K)）で熱伝導率が40倍以上異なる強い異方性を持つ。パウチセルでも面方向と厚み方向で大きな差がある。
- 化学反応項 $Q_{\text{rxn}}$：Arrhenius型の4段階反応モデル。温度上昇→反応加速→さらなる温度上昇のポジティブフィードバックが熱暴走の本質。
- 電気的発熱項 $Q_{\text{ohm}}$：通常運転時はBernardi式（$I^2R + I T \, dU/dT$）で記述。内部短絡時は短絡抵抗を通じた局所的な大電流ジュール発熱が支配的。
仮定条件と適用限界
- 均質化仮定：ジェリーロール内部の電極積層構造を均質異方性材料として扱う。セル内部の局所的な反応不均一は無視。
- Arrhenius型反応速度：各段階を単純な1次反応で近似。実際の反応は複数の並行パスを持つが、DSC/ARCデータへのフィッティングで有効なモデルが得られる。
- 反応物濃度の仮定：無次元濃度 $c_i$（0〜1）の消費を追跡。未反応率で表現する場合もある。
- 内部短絡モデル：釘刺し試験の再現では、短絡部の抵抗値と接触面積が結果に大きく影響。実験的キャリブレーションが必須。
- 適用限界：ガス噴出後の流体力学（ベントガス挙動）は別途CFDで解くか、簡易的な対流熱伝達係数として与える必要がある。
次元解析と単位系
| 変数 | SI単位 | 典型値・注意点 |
|---|---|---|
| 活性化エネルギー $E_a$ | J/mol | SEI: 135 kJ/mol、正極: 170〜200 kJ/mol |
| 頻度因子 $A$ | 1/s | 10⁹〜10²⁵のオーダー。段階により大きく異なる |
| 反応エンタルピー $H$ | J/kg または J/m³ | 体積基準と質量基準の変換に注意 |
| 内部抵抗 $R_{\text{int}}$ | Ω | SOC・温度依存。1〜100 mΩ程度 |
| セル間熱伝導率 | W/(m・K) | バリア材: 0.03〜0.5、アルミ筐体: 150〜200 |
数値解法と実装
理論はわかってきました。この方程式をコンピュータで実際にどう解くんですか？
電気化学-熱連成の非定常問題だから、空間離散化にFEM（有限要素法）、時間離散化に陰的スキーム（後退Euler法やCrank-Nicolson法）を使うのが標準だ。Arrhenius項の強い非線形性が数値的に厄介で、ここが熱暴走シミュレーション特有の難しさだよ。
FEM定式化
FEMの定式化はどうなりますか？ 熱伝導のFEMは構造解析と似てるんですか？
形状関数 $N_i$ で温度場を近似して、弱形式に変換するところまでは通常の熱伝導FEMと同じだ。ガラーキン法を適用すると：
- $[C] = \int_\Omega \rho c_p [N]^T [N] \, d\Omega$：熱容量マトリクス
- $[K] = \int_\Omega [B]^T k [B] \, d\Omega$：熱伝導マトリクス
- $\{Q\}$：発熱ベクトル（$Q_{\text{rxn}} + Q_{\text{ohm}}$を積分したもの）
ポイントは右辺の $\{Q\}$ が温度 $\{T\}$ に強く依存すること。Arrhenius項 $\exp(-E_a/RT)$ があるから、発熱量が温度の非線形関数になる。だから反復法（Newton-Raphson法）が必須だ。
時間積分スキーム
暴走が始まると温度が急激に変化しますよね。時間刻みはどう決めるんですか？
これが実務で最もよく悩むポイントだ。熱暴走前は $\Delta t = 1\text{〜}10$ 秒で十分だけど、暴走開始後は温度勾配が急峻になるから $\Delta t = 0.001\text{〜}0.01$ 秒まで細かくしないと追いつかない。
| スキーム | 安定性 | 精度 | 推奨場面 |
|---|---|---|---|
| 後退Euler（1次陰的） | 無条件安定 | $O(\Delta t)$ | ロバストだが拡散が大きい。初期スクリーニング向き |
| Crank-Nicolson（2次陰的） | 無条件安定 | $O(\Delta t^2)$ | 精度重視。暴走前の緩やかな昇温に適する |
| BDF2（2次後退差分） | A-安定 | $O(\Delta t^2)$ | stiff系に強い。暴走遷移期に有効 |
実務では適応的時間刻み制御が必須だ。温度変化量 $|\Delta T|$ がしきい値（例えば5°C/ステップ）を超えたら自動的に $\Delta t$ を半分にし、緩やかなフェーズでは増大させる。COMSOL Multiphysicsなどはこの機能が標準で使える。
非線形連成の取り扱い
Newton-Raphson法で収束させると言いましたが、Arrhenius項が指数関数だから収束しにくいとかありますか？
鋭い質問だ。Arrhenius項の温度に対する感度は極めて高い。ヤコビアン（接線剛性行列）に反応速度の温度微分が入るから：
暴走領域では $\partial Q/\partial T$ が非常に大きくなってヤコビアンの条件数が悪化する。対策としては：
- 荷重ステップの細分化：暴走開始検知で自動的にサブステップを追加
- ライン探索法：Newton更新量にダンピングをかけてオーバーシュートを防ぐ
- 反応物濃度のクランプ：$c_i$ を [0, 1] の範囲に強制拘束して非物理的な値を防ぐ
メッシュ戦略
メッシュの切り方で注意すべきことは？ 電池ってジェリーロール構造で層が何百もありますよね。
実務でジェリーロールの全層を個別にメッシュするのは計算コスト的に不可能だ。だから以下の多段階アプローチを使う：
| スケール | モデル | 典型的な要素数 | メッシュ種類 |
|---|---|---|---|
| 電極層（μm） | 1D電気化学モデル（P2Dモデル） | 100〜1,000 | 1D線分要素 |
| 単セル（mm） | 均質異方性体3D FEM | 10,000〜50,000 | 六面体推奨 |
| モジュール（cm） | 3D FEM（セル+バリア+筐体） | 100,000〜500,000 | 四面体/六面体混合 |
| パック全体（m） | 等価回路モデル + 簡易3D | 500,000〜2,000,000 | 粗い六面体 |
セル間バリアの厚さ方向には最低4層以上の要素を入れること。厚さ1〜3mmの薄いバリアに対して温度勾配が急峻になるから、ここをケチると伝播時間の予測精度が大幅に落ちる。
多段階モデリングの考え方
電池の熱暴走解析は「マトリョーシカ人形」のような入れ子構造になっている。一番小さな人形（電極層スケール）で電気化学反応の発熱量を求め、その結果を次のサイズの人形（単セルスケール）にパラメータとして渡し、さらにモジュール→パックへとスケールアップしていく。各スケールで使うソルバーや要素タイプが異なるため、データの受け渡し精度がシミュレーション全体の信頼性を左右する。
実践ガイド
解析フロー
実際に「電池が暴走しないか」を解析するとき、どんな手順で進めるんですか？
典型的なフローは以下の5段階だ：
- セルレベルキャリブレーション：単セルのARC試験データでArrhenius定数（$A_i, E_{a,i}, H_i$）をフィッティング
- 通常運転時の熱モデル検証：充放電サイクルでのセル温度実測値とBernardi式モデルの一致確認
- トリガーイベントのモデル化：釘刺し（局所短絡抵抗の設定）、過充電（SOCの強制上昇）、外部加熱（ヒーター境界条件）
- モジュールレベル伝播解析：トリガーセルからの隣接セル温度の時系列評価
- パラメトリックスタディ：バリア材厚さ・材質、冷却流量、セル間ギャップ等の設計変数を振って伝播時間を最適化
Step 1のキャリブレーションが一番大変そうですね。実験データがないと始まらない？
そのとおり。ARCかDSCの実験は最低限必要だ。ただし論文データも蓄積があって、例えばNMC111/NMC622/NMC811/LFPの各正極材については代表的なArrhenius定数が文献で公開されている。最初はそれを使って検証し、自社セル固有のデータで微調整するのが現実的だ。
材料パラメータの取得
Arrhenius定数以外に、どんな材料パラメータが必要ですか？
| パラメータ | 取得方法 | 注意点 |
|---|---|---|
| 熱伝導率（面方向/厚み方向） | レーザーフラッシュ法、ホットディスク法 | 異方性が非常に大きい（40倍差） |
| 比熱 $c_p$ | DSC | 温度依存性を考慮（特に電解液の相変化） |
| 内部抵抗 $R_{\text{int}}(T, \text{SOC})$ | EIS（電気化学インピーダンス分光） | SOCと温度の2変数マップが必要 |
| OCV温度係数 $dU/dT$ | 低速充放電でのOCV-T測定 | 正極材の種類で大きく異なる |
| バリア材の熱物性 | メーカースペックシート + 検証試験 | 高温（>300°C）での劣化を考慮 |
バーチャル試験（UL 2580 / GB 38031）
UL 2580とかGB 38031って規格試験ですよね？ シミュレーションで代替できるんですか？
規格認証そのものには実物試験が必須だから、完全な代替はできない。ただしバーチャル試験が非常に有効な場面が3つある：
- 事前スクリーニング：10パターンの断熱設計を全部実物で試したら莫大なコストがかかる。シミュレーションで3候補に絞ってから実試験する。
- 最悪ケースの特定：「どの位置のセルが暴走したら最も延焼が速いか」をFEMで網羅的に探索する。実試験では全位置を試すのは不可能。
- 規格変更への対応：新規格が追加されたとき、設計変更の方向をシミュレーションで即座に検討できる。
| 規格 | 対象 | 主な試験項目 | シミュレーション活用 |
|---|---|---|---|
| UL 2580 | EV用電池パック | 釘刺し、過充電、外部火炎 | 伝播時間予測、ベントガス挙動 |
| GB 38031 | 中国市場向けEVバッテリー | 熱的伝播試験（5分間非延焼） | バリア材パラメトリック設計 |
| UN ECE R100 | 国際統一基準 | 火災安全性 | 外部火炎暴露シミュレーション |
| SAE J2464 | 試験手順ガイドライン | 乱用試験方法の標準化 | 試験条件の設定根拠 |
GB 38031の「5分間非延焼」って、つまり1セルが暴走しても5分以内にパック全体に火が回らなければOKってことですか？
そう。正確には「トリガーセルの熱暴走後、パック外部に火炎や爆発が発生するまでに5分以上の猶予があること」が求められる。乗員が車両から脱出するための最低限の時間という考え方だ。シミュレーションではこの5分を目標値としてバリア設計を最適化する。
断熱バリア設計
断熱バリアの設計って、厚くすれば安全だけど重くなるし、エネルギー密度が下がるジレンマですよね？
まさにそれが実務の核心だ。シミュレーションで最適解を探る典型的なパラメトリック設計を見てみよう：
| バリア材 | 熱伝導率 [W/(m·K)] | 密度 [kg/m³] | 厚さ5mmでの伝播時間 |
|---|---|---|---|
| エアロゲルシート | 0.015 | 150 | 約180秒 |
| セラミックファイバー | 0.03 | 200 | 約120秒 |
| マイカシート | 0.4 | 2800 | 約50秒 |
| バリアなし（空気層のみ） | 0.026（空気） | — | 約30秒 |
エアロゲルは断熱性能が最高だけどコストが高い。セラミックファイバーはバランスが良く多くのOEMが採用している。これにアクティブ冷却（冷却水路）を組み合わせると、さらに伝播時間を延ばせる。
EV火災の「延焼防止」——パック設計者の闘い
熱暴走シミュレーションを実務に活かす最大の目的は「1セルが暴走してもパック全体への延焼を防ぐ」設計です。ある国内自動車メーカーの開発担当者は、シミュレーション結果を元にバリア材の材質をセラミックファイバーからエアロゲル複合材に変更し、厚さも8mmから5mmに薄くすることで、パック重量を3kg削減しながら伝播時間を2倍に延長できたそうです。しかし最終的には1,000回以上のシミュレーションケースと、20回以上の実物釘刺し試験を経て設計を確定しました。「シミュレーションなしでは20回の試験すら打てなかった。100回以上の実物試験が必要だったはず」というのが担当者の実感です。
バリア設計の「熱抵抗」の考え方
断熱バリアの効果は「熱抵抗」$R = d / (k \cdot A)$ で直感的に理解できる。厚さ $d$ を2倍にすれば熱抵抗は2倍、熱伝導率 $k$ を半分にしても熱抵抗は2倍になる。つまり「高性能な断熱材を薄く使う」のと「普通の断熱材を厚く使う」のは等価に思えるが、実際には高温域での材料劣化、機械的強度、コストが絡むため、シミュレーションで複合的に評価する必要がある。
ソフトウェア比較
ツール比較
電池の熱暴走シミュレーションを実際にやるには、どのソフトが使われているんですか？
主要なツールと、それぞれの得意分野を比較するとこうなる：
| ツール | 開発元 | 強み | 電池熱暴走対応 |
|---|---|---|---|
| COMSOL Multiphysics | COMSOL AB | マルチフィジクス連成が容易。Arrhenius反応をGUIで定義可能 | Battery Design Module（電気化学+熱+構造） |
| Ansys Fluent | Ansys Inc. | CFDとの連成。ベントガス挙動の解析 | MSMD（Multi-Scale Multi-Domain）バッテリーモデル |
| Simcenter STAR-CCM+ | Siemens | ポリヘドラルメッシュ。大規模パック解析 | Battery Simulation Module（ESC/NTGK モデル） |
| GT-AutoLion | Gamma Technologies | 電気化学の詳細モデル（P2D） | セルレベルの高精度予測。システムモデルとの連携 |
| LS-DYNA | LSTC / Ansys | 衝突+熱暴走の連成。機械的トリガーの再現 | *BATTERY_ECHEM + 熱ソルバー連成 |
LS-DYNAが入っているのが意外です。衝突解析のソフトですよね？
車両衝突でバッテリーケースが変形→セルが圧壊→内部短絡→熱暴走、という一連のプロセスを1つのモデルで追跡できるのがLS-DYNAの強みだ。機械的なトリガーから熱暴走までをシームレスに扱えるツールは他にあまりない。自動車OEMの衝突安全チームはこれを使うことが多い。
選定ガイド
結局どれを選べばいいですか？ 予算も限られてるんですけど…
用途で選ぶのが一番だ：
- セルレベルの電気化学+熱を詳しくやりたい → COMSOL（直感的なGUI）またはGT-AutoLion（P2Dモデルの精度）
- モジュール/パックの延焼解析が主目的 → Ansys FluentまたはSTAR-CCM+（大規模3D熱解析+CFD）
- 衝突→熱暴走の一気通貫 → LS-DYNA一択
- コストを抑えたい → OpenFOAMのcustomSolverでArrhenius項を自作する手もある。ただし開発工数は覚悟すること
マルチスケール解析のツール使い分け
熱暴走シミュレーションのツール選定で見落とされがちなポイントが「スケールをまたいで計算できるか」です。セル内部の電気化学反応はマイクロメートルスケール、モジュール全体の熱伝導はセンチメートルスケール、パックの冷却流路はミリメートルスケール——これらを1つのモデルで扱うのは計算コスト的に非常に厳しい。そのため実務では「電気化学はGT-AutoLionのP2Dモデル、延焼はFluent/STAR-CCM+の3Dモデル、システム全体はMATLAB/Simulink」というようにスケールごとにツールを使い分け、データを受け渡す手法が一般的です。
先端技術
機械学習サロゲートモデル
最近はAI/機械学習で熱暴走を予測する研究もあるんですか？
はい、大きく2つの方向がある：
- サロゲートモデル：FEMで数千ケースのパラメトリック計算を回してデータセットを作り、ニューラルネットワーク（DNN）やガウス過程回帰（GPR）で近似モデルを構築する。パラメータ最適化の評価関数として使うと、FEMを毎回回す必要がなくなり100〜1000倍速くなる。
- PINN（Physics-Informed Neural Network）：エネルギー保存式を損失関数に組み込んで物理法則を満たす解を学習する。実験データが少ない場合でも物理的に妥当な外挿ができる。ただし暴走域の急激な温度変化をPINNで捉えるのはまだ研究段階だ。
サロゲートモデルでバリア材の最適設計をやったら面白そうですね！
デジタルツイン・BMS連携
実際に車が走っているときにリアルタイムで熱暴走リスクを予測できるんですか？
それがデジタルツインの考え方だ。BMS（Battery Management System）がリアルタイムで取得する電圧・電流・温度データを、事前に構築したサロゲートモデルに入力して、「あと何分で熱暴走に至るか」を推定する。
たとえばセル温度が異常上昇した場合、サロゲートモデルが「現在の条件では120秒後にSEI分解開始」と予測したら、BMSが急速放電停止→冷却強化→乗員への警告を自動実行する。ここまで行けばCAEがリアルタイムの安全制御に直結する。
すごい、シミュレーションが設計だけじゃなく運用フェーズまでつながるんですね！
トラブルシューティング
よくあるエラーと対策
先生、実際に熱暴走シミュレーションを回すと、どんなトラブルに遭遇しますか？
一番多いのは以下の4つだ：
| 症状 | 原因 | 対策 |
|---|---|---|
| 暴走開始付近で収束失敗 （"Time step too small"） | Arrhenius項の急激な温度感度でヤコビアンが悪化 | 適応的時間刻みの下限を $10^{-5}$ s 程度に設定。Newton反復の最大回数を50→200に増加。ライン探索法を有効化 |
| 温度が負の値になる | 反応物濃度 $c_i$ がアンダーフローして発熱が異常 | $c_i$ を [0, 1] にクランプするロジックを追加。COMSOL: "Constraints" で下限設定 |
| 伝播時間が実験と大きくずれる | バリア材の高温劣化（熱伝導率の温度依存性）を無視 | 300°C以上での $k(T)$ カーブを入力。セラミックファイバーは700°C超で焼結し熱伝導率が急増する |
| セル温度上昇が遅すぎる | ジェリーロールの熱伝導異方性の設定ミス（面方向と厚み方向を逆にしている） | 面方向 $k_{\parallel} \approx 30$ W/(m·K)、厚み方向 $k_{\perp} \approx 0.7$ W/(m·K) を確認 |
異方性の設定ミスはありがちですね… 面方向と厚み方向で40倍も違うのか。
そう。銅箔・アルミ箔の集電体が面方向に熱を通すから高い値になる。厚み方向はセパレータや電解液の低い熱伝導率が支配的だ。この異方性を等方性で近似すると、暴走開始時刻が30%以上ずれることがある。
最後に「これだけは覚えておけ」というアドバイスはありますか？
3つだけ。
- まずキャリブレーション：自社セルのARC/DSCデータなしに信頼できる熱暴走シミュレーションは不可能。文献値は初期検証用にしか使えない。
- バリアの温度依存性を忘れるな：室温での熱物性値だけでモデルを組むのは最も危険な失敗。高温域で材料特性が激変する。
- 最悪ケースを探せ：「真ん中のセルが暴走するのが最も危険」とは限らない。冷却構造の近くで暴走すると冷媒が蒸発して冷却喪失になるケースもある。パラメトリック解析で漏れなくチェックすること。
ありがとうございます！ 熱暴走シミュレーションの全体像がよくわかりました。まずはCOMSOLのBattery Moduleでチュートリアルから始めてみます！
いいね！ COMSOLのApplication Galleryに "Lithium-Ion Battery Pack Thermal Runaway Propagation" というチュートリアルがある。まさにこの記事で説明した4段階Arrheniusモデルをそのまま実装しているから、最初の一歩として最適だよ。
関連トピック
なった
詳しく
報告