[論文レビュー] Fluctuating surface-current formulation of radiative heat transfer: theory and applications

Alejandro W. Rodríguez, M. T. Homer Reid+1 more|DSpace@MIT (Massachusetts Institute of Technology)|Apr 4, 2013

Thermal Radiation and Cooling TechnologiesEngineering被引用数 22

ひとこと要約

本稿では、古典的電磁気学における表面積分方程式（SIE）手法を活用して、任意の形状をした物体間の放射熱伝導を記述する、揺動表面電流（FSC）形式を導入する。熱的揺動を物体の境界面上の仮想的表面電流としてモデル化することで、高対称性幾何における半解析的スペクトル解法と、複雑な構造（有限スラブ、円筒、円錐など）における境界要素法（BEM）による堅牢な数値解法の両方が可能になる。

ABSTRACT

We describe a novel fluctuating-surface current formulation of radiative heat transfer between bodies of arbitrary shape that exploits efficient and sophisticated techniques from the surface-integral-equation formulation of classical electromagnetic scattering. Unlike previous approaches to non-equilibrium fluctuations that involve scattering matrices---relating "incoming" and "outgoing" waves from each body---our approach is formulated in terms of "unknown" surface currents, laying at the surfaces of the bodies, that need not satisfy any wave equation. We show that our formulation can be applied as a spectral method to obtain fast-converging semi-analytical formulas in high-symmetry geometries using specialized spectral bases that conform to the surfaces of the bodies (e.g. Fourier series for planar bodies or spherical harmonics for spherical bodies), and can also be employed as a numerical method by exploiting the generality of surface meshes/grids to obtain results in more complicated geometries (e.g. interleaved bodies as well as bodies with sharp corners). In particular, our formalism allows direct application of the boundary-element method, a robust and powerful numerical implementation of the surface-integral formulation of classical electromagnetism, which we use to obtain results in new geometries, including the heat transfer between finite slabs, cylinders, and cones.

研究の動機と目的

• 任意の形状の物体間の放射熱伝導を扱う一般化されたフレームワークの構築を目的とし、従来の散乱行列法や時領域シミュレーションに起因する制限を克服すること。
• 従来、平衡状態のカシミール力にのみ適用されてきた揺動表面電流（FSC）形式を、異なる温度にある物体間の非平衡放射熱伝導に一般化すること。
• 高対称性物体に対してはスペクトル基底（例：フーリエ級数、球面調和関数）を用いた半解析的解法と、複雑な幾何形状に対しては境界要素法（BEM）による完全な数値解法の両方を可能にすること。
• 非平衡場相関関数と表面積分形式との間の明確な数学的・物理的対応関係を、新規の表面電流表現を通じて確立すること。
• 本手法が、解析的に扱いにくい新しいかつ挑戦的な幾何形状（例：有限スラブ、円筒、円錐）への適用可能性を実証すること。

提案手法

• 物体の境界面上に仮想的表面電流を導入することで、場の波動方程式を解く必要を回避し、放射熱伝導を記述する。
• 古典的電磁気学における表面積分方程式（SIE）フレームワークを用い、グリーン関数作用素を通じて表面電流と電磁界を関連付ける。
• 非平衡系に一般化されたフラクチュエーション・ドレイン理論を適用し、場の相関関数を体積源の積分として表現。その後、SIE形式を用いてこれを表面積分に再定式化する。
• 再帰性と系行列の適切な対称性を保証するため、符号反転行列 $\mathcal{S}$ を用いた複素対称作用素構造を採用する。
• 基底関数（例：BEMにおけるRWG）を用いたガラーキン離散化により、連続的なSIEを線形方程式系 $M\mathbf{x} = \mathbf{b}$ に変換し、ここで $M = W^{-1}$ はSIE行列である。
• 高対称性物体に対してはスペクトル法、任意の幾何形状に対しては境界要素法（BEM）を用い、鋭いコーナーや重ね合わされた構造を含む複雑な形状に対しても対応可能である。

実験結果

リサーチクエスチョン

• RQ1 異なる温度にある任意の形状の物体間の放射熱伝導を、散乱行列や時領域シミュレーションに依存せずにどのように定式化できるか？
• RQ2 従来、平衡状態のカシミール力にのみ適用されてきた揺動表面電流（FSC）形式を、非平衡状態の放射熱伝導に一般化できるか？
• RQ3 非平衡系における体積積分された場の相関関数と、仮想電流による表面積分表現との間の数学的・物理的対応関係は何か？
• RQ4 SIEフレームワークとBEMを用いて、有限スラブ、円筒、円錐といった複雑な幾何形状における熱伝導をどの程度正確に計算できるか？
• RQ5 SIE行列（$M$, $L$, $\Gamma^\star$）の対称性および定性（definiteness）特性が、熱伝導計算における物理的整合性と数値的安定性をどのように保証するか？

主な発見

• FSC形式により、体積離散化を回避し表面未知数に還元することで、任意の幾何形状における放射熱伝導の正確かつ効率的な計算が可能になる。
• スペクトル基底（例：フーリエ級数、球面調和関数）を用いることで、高対称性幾何（例：平面、球面）において収束が速い半解析的結果が得られる。
• 境界要素法（BEM）の実装により、有限スラブ、円筒、円錐といった複雑な構造に対しても、解析的に取り扱いにくい場合でも堅牢な数値解が得られる。
• SIE行列 $M$ および作用素 $L$, $\Gamma^\star$ が再帰性のもとで負半定値かつ複素対称であることが示され、物理的整合性と数値的安定性が保証される。
• 形式は、非平衡系における場相関関数（グリーン関数の積を含む）と表面電流表現との間の明確な数学的リンクを確立し、熱揺動の統一的取り扱いを可能にする。
• 再帰性関係 $\Gamma^{\mathrm{T}}\star = \mathcal{S}(\Gamma\star)\mathcal{S}$ がSIE形式下でも保持され、磁場成分の符号反転を適切に処理することで、源と場の入れ替えが可能になる。