[논문 리뷰] Hybrid Solver for the Radiative Transport Equation Using Finite Volume and Discontinuous Galerkin

Vincent Heningburg, Cory D. Hauck|arXiv (Cornell University)|2020. 02. 06.

Computational Fluid Dynamics and Aerodynamics참고 문헌 24인용 수 6

한 줄 요약

이 논문은 복사 전송 방정식(RTE)에 대한 하이브리드 유한체적(FV)/비연속 갈레르킨(DG) 공간 이산화 방법을 제안한다. RTE는 충돌된(산산이 흩트린) 및 충돌하지 않은 성분으로 분리되며, 충돌된 부분은 정확도를 높이기 위해 2차 DG로, 충돌하지 않은 부분은 메모리 및 계산 비용을 줄이기 위해 FV로 해결된다. 이 방법은 전체 DG에 비해 훨씬 낮은 메모리 사용량과 더 빠른 해답 시간을 기록하면서도 유사한 정확도를 확보한다.

ABSTRACT

We propose a hybrid spatial discretization for the radiative transport equation that combines a second-order discontinuous Galerkin (DG) method and a second-order finite volume (FV) method. The strategy relies on a simple operator splitting that has been used previously to combine different angular discretizations. Unlike standard FV methods with upwind fluxes, the hybrid approach is able to accurately simulate problems in scattering dominated regimes. However, it requires less memory and yields a faster computational time than a uniform DG discretization. In addition, the underlying splitting allows naturally for hybridization in both space and angle. Numerical results are given to demonstrate the efficiency of the hybrid approach in the context of discrete ordinate angular discretizations and Cartesian spatial grids.

연구 동기 및 목표

• 산산이 흩트린 영역에서 정확도를 유지하면서 복사 전송 방정식(RTE)에 대한 비연속 갈레르킨(DG) 방법의 메모리 및 계산 비용을 줄이기.
• RTE를 충돌된 및 충돌하지 않은 성분으로 분리함으로써 유한체적(FV) 및 DG 방법을 조합하는 공간 하이브리드화 전략을 개발하기.
• 하이브리드 방법이 점점 더 산산이 흩트리는 영역에서 점점 더 확산되는 근사한 극한을 정확히 포착함으로써 물리적 일관성을 확보하기.
• 기존의 각도 하이브리드 방법을 활용하여 공간 및 각도 이격화를 동시에 적용할 수 있는 자연스러운 확장 가능성을 제공하기.
• 균일한 DG 이격화 대비 메모리 사용량과 해답 시간 측면에서 향상된 계산 효율성 입증하기.

제안 방법

• RTE는 충돌 방정식(산산이 흩트린 소스 포함)과 충돌하지 않은 방정식(산산이 흩트린 없음)으로 연산자 분리되어 별도의 공간 이격화 가능.
• 충돌된 성분은 2차 비연속 갈레르킨(DG) 방법으로 이격화되며, 이는 부족한 메esh에서도 산산이 흩트린을 정확하게 처리할 수 있다.
• 충돌하지 않은 성분은 2차 유한체적(FV) 방법으로 이격화되며, 이는 더 적은 자유도를 요구하고 메모리 사용량을 줄인다.
• FV는 일관된 2차 재구성 기법을 사용하며, 각도 방향에 기반한 업윈드-중심 기반 기울기 제한을 적용하여 표준 FV 대비 정확도 향상.
• 하이브리드 수식은 직각좌표계의 이산 순서 각도 이격화를 사용하며, 기존의 각도 하이브리드 방법과의 원활한 통합 가능.
• 기본 문제와 계산 스케일링 분석을 통해 수치적으로 검증되었으며, 메모리 및 연산 횟수를 명시적으로 추적함.

실험 결과

연구 질문

• RQ1 FV와 DG를 조합한 하이브리드 공간 이격화가 산산이 흩트린 복사 전송 문제에서 정확도를 유지하면서도 메모리 및 계산 비용을 줄일 수 있는가?
• RQ2 하이브리드 FV/DG 방법이 고산산이 흩트린 상태에서의 점점 더 확산되는 극한을 정확히 포착하는가?
• RQ3 기본 문제 전반에 걸쳐 하이브리드 방법의 메모리 사용량과 계산 시간은 균일한 DG 이격화와 비교해 어떻게 되는가?
• RQ4 하이브리드 공간 접근법은 이산 순서 방법과 같은 각도 하이브리드화로 자연스럽게 확장 가능한가?
• RQ5 하이브리드 FV/DG 방법을 사용할 경우 전체 DG 대비 정량적 성능 향상(메모리 및 해답 시간 측면)은 어느 정도인가?

주요 결과

• 하이브리드 FV/DG 방법은 균일한 DG에 비해 메모리 사용량과 계산 시간이 현저히 낮으며, 측정된 메모리 절감은 약 21.7 MB의 오버헤드와 더 낮은 총 벡터 저장소를 기록한다.
• 방법은 정적 상태, 점점 더 산산이 흩트리는 극한에서 2차 정확도를 유지하며, 수치적으로 검증된 바에 따르면 확산 극한을 정확히 포착한다.
• 선원천 기준 문제에서 하이브리드 방법은 자유도 수가 적고 계산 비용이 낮아지면서도 기준 DG 해답과 정확도가 일치한다.
• 계산 스케일링 분석 결과 하이브리드 방법은 메쉬 크기 및 각도 수에 따라 효율적으로 스케일링되며, 연산 횟수와 메모리 사용량은 배열 크기 및 벡터 수로부터 예측 가능하다.
• 방향성 업윈드 기반 재구성은 충돌하지 않은 영역에서 표준 FV 대비 정확도 향상시키며, 더 낮은 메모리로 안정적이고 정확한 해답 가능.
• 하이브리드 접근법은 예를 들어 이전 연구에서 제안된 이산 순서 하이브리드 방법과 같은 각도 하이브리드 방법과 효율적으로 결합 가능하며, 다차원 이격화 전략에서 호환성과 모ularity를 입증한다.