技術要約：車載用パワーエレクトロニクスの過渡熱モデリングに向けたスケーラブルなアプローチ

問題提起

インバータ、整流器、電源装置などの車載用パワーエレクトロニクスにおいて、効率的な熱管理は信頼性と性能の観点から極めて重要である。過度な熱は性能を低下させ、コンポーネントの寿命を縮め、システム故障を引き起こす可能性がある。数値流体力学（CFD）は正確な熱シミュレーションを提供するが、複数の熱源や混合熱伝達モード（伝導、自然対流、強制対流）を含む複雑で大規模なシステムにおいては、計算コストが高く、非常に時間がかかる。

既存の次数低減モデル（ROM）は、精度と効率のバランスを取ろうとしているが、以下のような重大な限界に直面している：

• Krylov部分空間およびバランス・トランクレーション法： 線形時不変（LTI）システムには効果的であるが、非線形性、熱・電気の結合挙動、および時変境界条件への対応が困難である。また、展開点の試行錯誤が必要であり、部分空間が大きくなるにつれて計算負荷が増大する。
• 固有直交分解（POD）： 主要な熱モードを持つシステムには効果的であるが、自動車用途のように熱エネルギーが多くの空間的・時間的スケールに分散している場合、得られる利点は限定的となる。
• 機械学習： リアルタイム予測を進展させているものの、これらデータ駆動型のアプローチは、データの品質への依存、過学習、解釈性の欠如、および安定性の保証という課題を抱えており、クリティカルなシステムにおける安全上の懸念を生じさせる。
• コンパクト熱ネットワーク（CTN）： パラメータ抽出のために広範な測定に依存することが多く、チップ間の熱結合を省略したり、既知の温度依存性を無視して線形特性を仮定したりする場合がある。

したがって、フルCFDシミュレーションの膨大なコストをかけることなく、高速な過渡応答、変動する流速、および非線形な熱伝達メカニズムを扱うことができる、スケーラブルで解釈可能かつ計算効率の高いモデリング手法が求められている。

手法

本論文では、**集中定数線形重ね合わせ（LPLSP）**法を提案する。このハイブリッドアプローチは、集中定数モデリングと線形重ね合わせの原理を組み合わせることで、過渡的な熱挙動を効率的に捉えるものである。この手法は、一般化熱方程式に基づいているが、CFDデータに適合させた経験的な定式化を利用している。

ケーススタディの構成

本手法は、3相システム用のハイサイドおよびローサイドのペアとして配置された6つのMOSFETを含むインバータモジュールを用いて検証されている。このモジュールは、フィン付きヒートシンクに取り付けられたプリント基板実装（PCBA）上に搭載されている。シミュレーションは、Ansys® Icepak™を用いて以下の3つの条件下で実施された：

1. 自然対流
2. 一定速度の強制対流
3. 変動する流速の強制対流

モデルの定式化

LPLSPモデルは以下のステップで導出される：

1. 断熱系： 温度拡散が瞬時であると仮定し、時間の経過とともに温度が線形に上昇すると仮定する（$T(t) = T_0 + (P/C)t$）。
2. 伝導： 空間勾配が顕著な複数の物体へと拡張する。温度は、線形項と、RC回路のような過渡応答（$1 - e^{-\theta t}$）に従う偏差項（$T_D$）としてモデル化される。複数の熱源を考慮するために重ね合わせが適用される。
3. 対流： 流速の影響を組み込む。変動する流速に対しては、パラメトリックスタディを用いて、熱抵抗（$R_{ij}$）、時定数（$\theta_i$）、および流速（$U_f$）の間の関係を確立する。これらのパラメータは、流速のべき乗関数（$R_{ij} = m U_f^\gamma + c$）として表現される。
4. 一般化： 混合モードおよび可変入力を持つシステムにおいて、任意のコンポーネントの温度は、各熱源の寄与を、流速の関数であるそれぞれの熱抵抗および時定数で重み付けして合計することで算出される。

パラメータ推定

モデル定数は、以下のパラメトリックCFDシミュレーションプロセスを通じて決定される：

• 他のMOSFETがオフの状態のとき、一度に1つのMOSFETに対して単位電力（1W）を印加する。
• 過渡応答を捉えるために、短時間（約10秒間）のシロレーションを実行する。
• 熱抵抗と時定数は、CFDデータとモデルによる予測された線形温度曲線からの偏差の二乗誤差の和を最小化するように最適化される。
• 可変速度のシナリオについては、離散的な速度（例：1, 5, 10, 20 m/s）でパラメトリックシミュレーションを実行し、中間値には線形補間を使用する。

主な貢献

• スケーラブルな過渡モデリング： LPLSP法は、複数の熱源や混合熱伝達モード（変動する流速を含む）を持つシステムにおいて、過渡的な熱挙動を正確にモデル化するフレームワークを提供する。
• 非線形性の処理： 従来のLTI-ROMとは異なり、LPLSPは、局所的な線形化が成立する短い時間窓内で応答を特徴付けることにより、非線形領域（例：流動と温度が結合している自然対流）においても高い精度を維持する。
• 計算効率： パラメトリックシミュレーション（一度限りのコスト）が完了すれば、モデルは任意の入力変化（電力消費、流速）を1秒未満で評価でき、これはフルCFDよりも数桁高速である。
• 解釈性： モデルは、物理的パラメータ（形状、材料特性、流動）と熱応答との間の明確な関数関係を提供し、機械学習の「ブラックボックス」的な性質を回避する。
• 最適化戦略： 本論文は、幾何学的形状や材料特性が類似している、あるいは空間的に近接しているコンポーネントは、同様のモデルパラメータを持つことを特定している。コンポーネントをクラスタリングし、時定数が無視できる程度の相互作用を排除することで、必要なパラメトリックシミュレーションの回数をさらに削減できることを示唆している。

結果

LPLSPモデル、LTI-ROM、およびフルCFDシミュレーションの間で、3つのシナリオにわたる比較が行われた：

1. 自然対流：
   • LPLSP： CFDと比較して、平均誤差率は**4.74%**を達成した。
   • LTI-ROM： 自然対流特有の非線形性と時変性により、LTIの仮定が崩れるため、**25.11%**という大幅な誤差を出した。
2. 強制対流（一定速度）：
   • LPLSPおよびLTI-ROMのいずれも、定常流条件下ではシステムが線形に振る舞うため、5%未満の誤差（それぞれ2.48%および4.97%）を達成した。
3. 強制対流（変動速度）：
   • LPLSP： 速度依存パラメータに対して区分線形近似を用いることで、5%未満の誤差（4.22%）を維持した。
   • LTI-ROM： 大きな速度変化は固定速度のインパルス関数を不正確にするため、評価対象外とした。

計算時間：

• フルCFDシミュレーションは、1ケースあたり730～900秒を要した。
• LPLSPのパラメトリック開発には（速度ステップ数に応じて）1440～4800秒を要したが、その後のMATLABでの評価時間は1秒未満であった。
• LTI-ROMの開発には（定常収束が必要なため）4140秒を要し、評価時間は約6秒であった。

意義と主張

本論文は、LPLSP法が自動車用パワーエレクトロニクスの設計サイクルを短縮するための実用的な経路を提供すると主張している。その主な意義は以下の通りである：

• 高速な過渡現象や変動する周囲条件を伴うシステムにおける、迅速な熱性能評価を可能にする。
• 従来のLTIベースのROMが失敗するシナリオ（自然対流、変動流など）においても、高い精度（誤差5%未満）を提供する。
• フルCFDが極めて困難な、初期段階の設計イテレーションや長時間のミッションプロファイルシミュレーションのための堅牢なツールとして機能する。

著者らは限界についても認めており、モデルは固定された機械的レイアウトを前提としており、複雑な多相流（例：ヒートパイプ）や、大きく異なる時定数を持つ多段階の指数関数的上昇を示すシステムには苦慮する可能性があるとしている。今後の課題として、実験データに対するモデルの検証や、相変化ダイナミクスを取り入れたハイブリッド・サブモデルの導入が提案されている。しかし、固定レイアウトの自動車用電子機器という範囲内において、本手法は従来のシミュレーション技術に代わる、スケーラブルで正確かつ計算効率の高い手法として提示されている。