[논문 리뷰] Multi-Frequency Implicit Semi-analog Monte-Carlo (ISMC) Radiative Transfer Solver in Two-Dimensions (without Teleportation)

Elad Steinberg, Shay I. Heizler|arXiv (Cornell University)|2021. 08. 05.

Radiative Heat Transfer Studies참고 문헌 31인용 수 16

한 줄 요약

이 논문은 전통적인 암시적 몽테카를로(IMC) 방법에서 유래하는 전이 오류를 제거하기 위해 에너지를 운반하는 물질 입자를 사용하고 암시적 시간 적분을 허용하는 2차원 다주파수 암시적 반사 몽테카를로(ISMC) 방사 전달 해법을 제시한다. ISMC 방법은 1D 및 2D 기하구조에서 회색 및 다주파수 문제 모두에서 정확하고 전이 오류가 없는 해를 도출하며, 광학적으로 두꺼운 영역과 고기울기 투과도 영역에서 IMC보다 빠르게 수렴한다.

ABSTRACT

We study the multi-dimensional radiative transfer phenomena using the ISMC scheme, in both gray and multi-frequency problems. Implicit Monte-Carlo (IMC) schemes have been in use for five decades. The basic algorithm yields teleportation errors, where photons propagate faster than the correct heat front velocity. Recently [Po\"ette and Valentin, J. Comp. Phys., 412, 109405 (2020)], a new implicit scheme based on the semi-analog scheme was presented and tested in several one-dimensional gray problems. In this scheme, the material energy of the cell is carried by material-particles, and the photons are produced only from existing material particles. As a result, the teleportation errors vanish, due to the infinite discrete spatial accuracy of the scheme. We examine the validity of the new scheme in two-dimensional problems, both in Cartesian and Cylindrical geometries. Additionally, we introduce an expansion of the new scheme for multi-frequency problems. We show that the ISMC scheme presents excellent results without teleportation errors in a large number of benchmarks, especially against the slow classic IMC convergence.

연구 동기 및 목표

• 광학적으로 두꺼운 물질에서 물리적으로 불가능한 열파의 전파를 유발하는 전통적인 암시적 몽테카를로(IMC) 방법에서 오랜 기간 지속된 전이 오류 문제를 해결한다.
• 이전에 명시적 방법에 국한되었던 반사 몽테카를로(SMC) 프레임워크를 암시적 형태로 확장하여 정확도를 희생시키지 않고도 큰 시간 간격을 허용한다.
• 다양한 기준 문제를 포함한 2차원 카르테시안 및 원통 기하구조에서 ISMC 방법의 유효성을 검증한다. 이는 모두 회색 및 다주파수(다그룹) 사례를 포함한다.
• ISMC가 고투과도 기울기 또는 광학적으로 두꺼운 물질이 있는 영역에서 고전적 IMC보다 더 빠른 수렴을 달성함을 보여준다.
• IMC가 전이 및 수렴 문제로 실패하는 상황에서, 큰 시간 간격과 낮은 공간 해상도 조건에서도 안정성과 정확도를 확보한다.

제안 방법

• 에너지 운반을 위한 물질 입자를 사용하고, 광자는 존재하는 물질 입자에서만 방출되는 반사 몽테카를로(SMC) 프레임워크를 채택하여 인위적 광자 전파를 제거한다.
• 시간 간격 동안 물질 입자와 광자 간의 에너지 교환을 풀어내는 방식으로 암시적 시간 적분을 구현하며, 시간 간격의 시작 시점에서의 물질 온도를 사용한다.
• 다주파수 문제에 대해 주파수 의존성 투과도 모델을 사용하며, 에너지 의존성 단면적과 온도 의존성 방출 및 흡수 비율을 적용한다.
• 시간 간격당 일정한 입자 수를 사용하여 닫힌 시스템에서 입자 수 조절이 필요 없도록 하며, 방출 및 흡수 추적을 통해 입자 수 보존을 유지한다.
• 이산 공간 메esh를 사용하여 2D에서 방사 전달 방정식을 풀며, 반사 경계 조건과 t=0 시점에 흑체 소스 항을 적용한다.
• 다주파수 기준 문제에서는 Δt = 10⁻¹³ s의 시간 간격을 사용하고 총 입자 수를 2.5×10⁷로 제한하여 수치적 안정성과 수렴성을 확보한다.

실험 결과

연구 질문

• RQ1 ISMC 방법은 광학적으로 두꺼운 물질과 고기울기 투과도 영역이 있는 2차원 방사 전달 문제에서 전이 오류를 제거하는가?
• RQ2 ISMC는 1D 및 2D의 회색 및 다주파수 기준 문제에서 고전적 IMC와 비교해 수렴 속도와 정확도 측면에서 어떻게 성능을 발휘하는가?
• RQ3 IMC가 전이 및 최대 원리 위반으로 실패하는 상황에서, ISMC는 큰 시간 간격과 굵은 공간 해상도 조건에서도 안정성과 정확도를 유지할 수 있는가?
• RQ4 ISMC의 통계적 노이즈는 광학적으로 얇은 영역에서 IMC와 비교해 어떻게 나타나며, 입자 수를 증가시켜 줄일 수 있는가?
• RQ5 ISMC 방법은 불연속적인 물질 특성을 가진 복잡한 2차원 기하구조, 예를 들어 홀라룸 및 격자 구조에서 확장 가능하고 강건한가?

주요 결과

• ISMC 방법은 모든 2차원 기준 문제에서 전이 오류를 성공적으로 제거하였으며, 광학적으로 두꺼운 영역과 고기울기 투과도 영역 모두에서 조차도 굵은 공간 해상도와 작은 시간 간격에서도 효과를 발휘한다.
• 2차원 마샬라크 웨이브 및 홀라룸 문제에서 ISMC는 광학적으로 두꺼운 물질에서 IMC보다 훨씬 더 빠르게 정확한 온도 및 복사 에너지 밀도 프로파일에 수렴한다.
• 알루미늄 블록과 거품이 있는 2차원 다주파수 문제에서 ISMC 결과는 PN 기준 해와 원래 IMC 구현 결과와 뛰어난 일치를 보이며, 관측 가능한 전이 효과가 없었다.
• 이산적 발광 및 흡수로 인해 광학적으로 얇은 영역에서 더 높은 통계적 노이즈가 발생하지만, ISMC는 안정적이고 수렴 가능하며 입자 수를 늘림으로써 노이즈를 줄일 수 있다.
• 동일한 입자 수를 사용할 경우 ISMC는 IMC와 유사한 런타임을 유지하며, 큰 시간 간격에서 최대 원리 위반에 대해 유사한 내성도를 보였다.
• 모든 기준 문제에서 ISMC는 투명하거나 고기울기 영역에서 IMC보다 훨씬 낮은 공간 해상도 요구 조건으로 수렴한 해를 도출하였으며, 이러한 영역에서 뛰어난 효율성을 입증하였다.