[论文解读] 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 Technologies被引用 22

一句话总结

本文提出了一种用于任意形状物体之间辐射传热的波动表面电流（FSC）公式，利用经典电磁学中的表面积分方程（SIE）方法。通过在物体边界上以虚构表面电流建模热涨落，该方法能够在高对称性几何结构中实现半解析谱解，同时在复杂构型（包括有限平板、圆柱体和圆锥体）中通过边界元法（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$）的对称性与半负定性特征在传热计算中如何确保物理解释一致性与数值稳定性？

主要发现

• FSC公式通过将问题简化为表面未知量，避免体积分割，实现了在任意几何形状中对辐射传热的精确高效计算。
• 在高对称性几何形状（如平面、球形）中，通过使用傅里叶级数或球谐函数等谱基，该方法可实现快速收敛的半解析结果。
• 边界元法（BEM）实现可对复杂构型（包括有限平板、圆柱体和圆锥体）提供稳健的数值解，这些构型在解析方法下难以处理。
• SIE矩阵 $M$ 与算子 $L$、$\Gamma^\star$ 在互易性条件下被证明为负半定且复对称，确保了物理解释一致性与数值稳定性。
• 该形式化建立了非平衡场关联（涉及格林函数乘积）与表面电流表示之间的严格联系，实现了对热涨落的统一处理。
• 在SIE公式下，互易关系 $\Gamma^{\mathrm{T}}\star = \mathcal{S}(\Gamma\star)\mathcal{S}$ 得到保持，使得源与场可互换，仅需对磁分量进行适当的符号翻转。