3차원 임의 형상 근접장 복사를 위한 열 이산 쌍극자 근사법

초록

열 이산 쌍극자 근사법(T‑DDA)은 근접장 열복사를 다루는 새로운 수치 기법으로, 물체를 작은 입방체 서브볼륨으로 분할하고 각 서브볼륨을 전기 점쌍극자로 모델링한다. 열전류의 플럭투에이션‑디스퍼전 정리를 적용해 스토캐스틱 맥스웰 방정식을 풀어 전기장과 쌍극자 모멘트를 구한다. 구형과 입방체 실리카 입자 사이의 전도율을 정확히 재현했으며, 구형에서는 유전율이 클수록 수렴에 필요한 서브볼륨 수가 증가한다는 특징을 보였다.

상세 분석

본 논문은 근접장 열복사 문제를 해결하기 위해 전통적인 이산 쌍극자 근사법(DDA)을 확장한 열 이산 쌍극자 근사법(T‑DDA)을 제안한다. 핵심 아이디어는 열복사에서 발생하는 무작위 전류 밀도 J(r) 를 플럭투에이션‑디스퍼전 정리로 통계적으로 기술하고, 이를 전기 점쌍극자 p 로 변환함으로써 스토캐스틱 맥스웰 방정식의 해를 구하는 것이다.

1. 이론적 기반

   • 맥스웰 방정식에 열 전류 J(r) 를 추가하고, 시간‑조화 형태( e^{-iωt} )로 변환한다.
   • 평균 전류는 0이지만, 전류‑전류 상관함수는 온도 T 와 물질의 복소 유전율 ε(ω) 에 의해 결정된다(식 2.2).
   • 이 상관함수를 이용해 전기장 E(r) 의 두점 상관함수를 구하고, 열전달 플럭스를 평가한다.

2. 볼륨 적분 방정식

   • 전체 전기장을 입사장 E_inc 와 산란장 E_sca 로 분리하고, 자유공간 그린함수 G(r,r′) 를 이용해 적분식(3.7)을 도출한다.
   • 입사장은 열전류에 의해 생성되는 장이며, 산란장은 물체 내부의 유전율 차이에 의해 발생한다.

3. 이산화 및 점쌍극자 모델

   • 물체를 N개의 정육면체 서브볼륨으로 분할하고, 각 서브볼륨을 균일한 전기장과 유전율을 갖는 점쌍극자 p_i 로 가정한다.
   • 클라우시우스‑모소티(Clausius‑Mossotti) 폴라리제이션 α_i 와 복사 폴라리제이션 α_i^{rad} (식 3.23)를 도입해 자기‑상호작용을 포함한다.
   • 서브볼륨 간 상호작용은 자유공간 그린함수 G(r_i,r_j) 를 통해 행렬 A 에 기록되며, 전체 시스템은 3N×3N 복소 선형 방정식(식 3.24)으로 정리된다.

4. 통계적 해석

   • 열전류에 의해 유도된 점쌍극자 p_i 는 평균이 0이지만, 상관함수 ⟨p_i p_j^*⟩ 는 식 3.20에 의해 온도와 물질 손실( Im