UVSim3D의 DOM(Discrete Ordinates Method) 모델이 복사전달방정식을 옳게 푸는지를 복사전열 분야의 표준 문헌 케이스로 검증한 결과입니다.

모델 개요

DOM은 복사전달방정식(RTE)을 유한 개의 이산 방향에 대해 풉니다. 흡수·방출·산란이 있는 매질에서 방향 $\hat{s}$ 를 따르는 복사강도 $I$ 는 다음을 만족합니다.

$$\hat{s}\cdot\nabla I(\mathbf{r},\hat{s}) + (\kappa + \sigma_s)\,I(\mathbf{r},\hat{s}) = \kappa I_b(\mathbf{r}) + \frac{\sigma_s}{4\pi}\int_{4\pi} I(\mathbf{r},\hat{s}')\,d\Omega'$$

여기서 $\kappa$ 는 흡수계수, $\sigma_s$ 는 산란계수, $I_b$ 는 매질의 흑체 복사강도입니다. 관심량인 조사강도(fluence rate) $G$ 는 전방향 적분입니다.

$$G(\mathbf{r}) = \int_{4\pi} I(\mathbf{r},\hat{s})\,d\Omega$$

UV 반응기 해석에서는 물의 산란이 무시할 수준이고 상온에서의 열복사도 UV 파장대에서는 무의미하므로 $\sigma_s = 0$, $\kappa I_b = 0$ 이 되어 식은 순수 감쇠 문제로 축약됩니다. 이때 균질 매질을 지나는 감쇠는 Beer-Lambert 법칙과 일치합니다.

$$I(s) = I_0\,e^{-\kappa s}$$

검증 케이스

Fiveland(1984)가 제시한 세 케이스 중 두 개를 사용했습니다. 하나는 흡수가 지배하는 문제, 다른 하나는 산란만으로 구동되는 문제로, 방정식의 서로 다른 항을 각각 고립시켜 검증합니다.

케이스 A 조사강도 분포 예시 — 네 벽 모두 흡수벽이라 중심이 가장 높고 모서리로 갈수록 낮아지는 대칭 분포

케이스 A — 흡수 매질이 채워진 흑체 밀폐공간

• 정사각 밀폐공간, 네 벽 모두 차갑고 검은 벽(방사율 1, 방출 없음)
• 내부는 순수 흡수 매질이며 단위 방출능 $E_b = 1$ 로 방출, 산란 없음
• 흡수계수 $\kappa$ = 0.1, 1.0, 10.0 m⁻¹ — 광학두께를 100배 범위로 훑습니다
• 비교량은 벽면 복사 열유속의 무차원값 $\hat{Q} = q / E_b$

케이스 B — 순수 산란 매질

• 정사각 밀폐공간, 벽 하나가 뜨거운 벽(단위 방출능), 나머지 셋은 차가운 벽
• 순수 산란 매질 $\sigma_s = 1\ \mathrm{m^{-1}}$, 흡수 없음 ( $\kappa = 0$ )
• 벽 방사율 $\epsilon$ = 1.0, 0.5, 0.1 — 불투명 회색면이므로 반사율은 $\rho = 1 - \epsilon$
• 비교량은 뜨거운 벽에서 빠져나가는 순 열전달률

흡수가 전혀 없으므로 뜨거운 벽이 내보낸 에너지는 전부 다른 벽에 도달해야 합니다. 근사가 아니라 항등식이며, 산란항이 에너지를 만들거나 잃지 않는지를 직접 확인할 수 있습니다.

원 문제는 2차원(무한 길이 덕트)입니다. UVSim3D는 3차원 해석기이므로 한 셀 두께의 격자 양면에 대칭 경계를 주어 무한 덕트를 재현했습니다. 두께 방향으로 균일한 해는 양쪽 대칭면의 경면반사 조건을 자동으로 만족하므로, 이는 근사가 아니라 정확한 재현입니다. 격자는 논문과 같은 한 변 20분할입니다.

결과 1 — 흡수 매질, 벽면 열유속 분포

세 광학두께 모두에서 벽 중앙이 최대, 모서리가 최소인 대칭 분포가 얻어지며, 값의 범위는 논문이 제시한 결과와 일치합니다. 이산 방향 수를 16 → 256으로 늘리면 해가 단조 수렴하며, 144방향과 256방향의 차이는 0.3% 미만입니다.

방향 수가 부족할 때(16방향) 벽 중앙에서 값이 과대평가되는 경향이 보이는데, 이는 DOM의 알려진 특성으로 논문도 최저차 근사는 흡수 매질 문제에 권장하지 않는다고 기술합니다.

결과 2 — 순수 산란 매질

뜨거운 벽에서 빠져나가는 순 열전달률이 세 방사율 모두에서 문헌 결과와 일치합니다.

• $\epsilon = 1.0$ : 0.774 – 0.867 (문헌 약 0.78 – 0.85)
• $\epsilon = 0.5$ : 0.379 – 0.389 (문헌 약 0.38 – 0.40)
• $\epsilon = 0.1$ : 0.0750 – 0.0751 (문헌 약 0.07 – 0.08)

방사율이 낮아질수록 분포가 평평해지는데, 반사가 강한 밀폐공간이 복사장을 균질화하기 때문으로 문헌의 경향과 같습니다.

정량적으로 더 확실한 것은 에너지 쪽입니다. 벽이 내보낸 복사에너지는 계산상 $\epsilon E_b A$ 와 정확히 일치했고 ( $\epsilon$ = 1.0 / 0.5 / 0.1 에 대해 0.0500 / 0.0250 / 0.0050 W), 매질 흡수는 정확히 0, 나머지는 전부 벽에 도달했습니다. 에너지 수지 오차는 세 경우 모두 0.00%입니다. 산란이 복사를 재분배할 뿐 생성하거나 소멸시키지 않는다는 확인입니다.

산란이 있는 경우 수렴에 128회 반복이 필요했습니다(흡수만 있는 경우 14회). 산란이 이산 방향들을 서로 결합시키기 때문이며, 예상되는 거동입니다.

결과 3 — 해석적 극한과의 대조

그래프의 범위 일치만으로는 검증이 약합니다. 아래 두 극한은 참값이 해석적으로 알려져 있어 정량적 대조가 가능합니다.

광학적으로 두꺼운 극한. $\kappa$ 가 충분히 크면 벽은 사실상 흑체를 마주보게 되어 $\hat{Q} \rightarrow 1$ 이 됩니다. $\kappa = 10\ \mathrm{m^{-1}}$ 에서 계산값은 벽 중앙에서 0.9982 로, 이론값과의 차이가 0.18%입니다.

광학적으로 얇은 극한. $\kappa$ 가 작으면 매질이 방출한 복사가 재흡수되지 않고 모두 벽에 도달하므로, 벽 평균 열유속은 $\hat{Q}_\mathrm{avg} = \kappa a$ 로 수렴합니다 ( $a$ 는 한 변 길이 ). $\kappa = 0.1\ \mathrm{m^{-1}}$ 에서 계산값은 0.0933이고 이론값은 0.1인데, 이 차이는 오차가 아닙니다. 같은 계산에서 매질이 자기 방출의 6.7%를 재흡수한 것으로 집계되며, 0.0933은 정확히 0.1의 93.3%입니다. 즉 부족분이 재흡수량과 일치합니다.

에너지 수지. 매질이 방출한 총 복사에너지는 매질에 재흡수되거나 벽에 도달해야 합니다. 모든 계산 조건에서 이 수지가 성립하며(오차 0.00%), 방출 총량도 해석값 $4\kappa E_b V$ 와 정확히 일치합니다. 대칭 경계가 에너지를 만들어내거나 잃지 않는다는 확인이기도 합니다.

결론

• 흡수 매질에서 벽면 열유속 분포가 광학두께 0.1 – 10의 전 범위에서 문헌 결과와 일치하며, 이산 방향 수에 대해 단조 수렴합니다.
• 순수 산란 매질에서도 벽 방사율 0.1 – 1.0 전 범위에서 문헌 결과와 일치합니다.
• 해석해가 존재하는 두 극한에서 오차는 각각 0.18% (두꺼운 극한)와 재흡수량으로 완전히 설명되는 수준(얇은 극한)입니다.
• 에너지 보존이 모든 조건에서 성립합니다(오차 0.00%).

DOM 모델이 복사전달방정식을 옳게 이산화하고 있음을 확인했습니다.

참고문헌

• Fiveland, W. A. (1984). Discrete-Ordinates Solutions of the Radiative Transport Equation for Rectangular Enclosures. Journal of Heat Transfer, 106(4), 699-706.
• Chai, J. C., Lee, H. S., & Patankar, S. V. (1993). Ray Effect and False Scattering in the Discrete Ordinates Method. Numerical Heat Transfer, Part B, 24(4), 373-389.