[논문 리뷰] Globally Conservative, Hybrid Self-Adjoint Angular Flux and Least-Squares Method Compatible with Void

Vincent Labouré, Ryan G. McClarren외 1명|arXiv (Cornell University)|2016. 05. 17.

Seismic Imaging and Inversion Techniques참고 문헌 8인용 수 5

한 줄 요약

이 논문은 Self-Adjoint Angular Flux (SAAF)와 Conservative Least-Squares (CLS) 수식을 조합하여 전역적으로 보존되고 벽이 없는 영역과도 호환되는 2차 수송 방법을 제안한다. 비공백 영역에서는 SAAF를 사용하고, 공백 영역에서는 수정된 최소제곱(Least-Squares, LS) 수식을 적용하여 전역 보존성을 복원함으로써, 오직 공백 영역에서만 대칭성 손실이 발생하도록 한다. 이 방법은 $P_N$ 및 $S_N$ 이산화에 대해 검증되었으며, 표준 LS 및 SAAF–VT 방법과 비교해 $k$-고유값 문제와 공백이 지배하는 문제에서 뛰어난 정확도와 보존성을 확보하였다.

ABSTRACT

In this paper, we derive a method for the second-order form of the transport equation that is both globally conservative and compatible with voids, using Continuous Finite Element Methods (CFEM). The main idea is to use the Least-Squares (LS) form of the transport equation in the void regions and the Self-Adjoint Angular Flux (SAAF) form elsewhere. While the SAAF formulation is globally conservative, the LS formulation need a correction in void. The price to pay for this fix is the loss of symmetry of the bilinear form. We first derive this Conservative LS (CLS) formulation in void. Second we combine the SAAF and CLS forms and end up with an hybrid SAAF-CLS method, having the desired properties. We show that extending the theory to near-void regions is a minor complication and can be done without affecting the global conservation of the scheme. Being angular discretization agnostic, this method can be applied to both discrete ordinates (SN) and spherical harmonics (PN) methods. However, since a globally conservative and void compatible second-order form already exists for SN (Wang et al. 2014), but is believed to be new for PN, we focus most of our attention on that latter angular discretization. We implement and test our method in Rattlesnake within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework. Results comparing it to other methods are presented.

연구 동기 및 목표

• 기존 SAAF 및 LS 수식의 한계를 해결하기 위해 전역적으로 보존되고 벽이 없는 영역과도 호환되는 2차 수송 방법을 개발하는 것.
• 공백 영역에서 표준 최소제곱(Least-Squares, LS) 수식이 전역 보존성을 상실함으로써 $k$-고유값 문제의 정확도가 떨어지는 문제를 해결하는 것.
• 일정한 교차단면을 가진 근접 공백 영역으로까지 방법을 확장하면서도 전역 보존성을 유지하는 것.
• $P_N$ 및 $S_N$ 각각의 각도 이산화 방식과 모두 호환되도록 하여, 특히 $P_N$ 방법에 대한 빈도를 메우는 것.
• MOOSE 프레임워크 내 Rattlesnake 프레임워크에 구현하고 기존 방법 및 기준 해와의 비교를 통해 하이브리드 SAAF–CLS 방법을 검증하는 것.

제안 방법

• 영역을 공백 영역 ($\mathcal{D}_0$)과 비공백 영역 ($\mathcal{D}_1$)으로 분할하고, $\mathcal{D}_1$에서는 SAAF를, $\mathcal{D}_0$에서는 수정된 LS 수식을 사용한다.
• 공백 영역에서 전역 보존성을 복원하기 위해 비대칭 보정 항을 추가하여 Conservative LS (CLS) 수식을 유도하며, 이는 LS 형태의 이산화 오차를 보정한다.
• 하이브리드 SAAF–CLS 변분 수식은 두 수식을 스케일링 인자로 결합하여 경계 조건 항이 사라지도록 하여 전역 보존성을 유지한다.
• 일정한 교차단면 영역으로 간주하여 근접 공백 영역으로 일반화하며, LS 수식에 수정된 보정 항을 적용한다.
• 각도 이산화에 관계없이 적용 가능한 방법으로, $S_N$ 및 $P_N$ 방법 모두 지원하며, MOOSE 프레임워크 내 Rattlesnake 코드에 구현되었다.
• 비대칭성 손실이 공백 보정으로 인해 발생하더라도, 철저한 변분 수식을 통해 전역 보존성이 유지된다.

실험 결과

연구 질문

• RQ1 기존 방법이 조건이 나쁜 문제로 인해 실패하는 것과는 달리, $P_N$ 각도 이산화를 사용할 때 전역적으로 보존되고 벽이 없는 영역과도 호환되는 2차 수송 방법을 개발할 수 있는가?
• RQ2 표준 LS 수식이 공백 영역에서 비보존적임을 어떻게 보정할 수 있으며, 안정성이나 계산 효율성을 희생하지 않고서도 가능한가?
• RQ3 적절한 스케일링을 통해 SAAF와 CLS 수식을 조합하면 공백 영역과 비공백 영역 간 경계에서 전역 보존성이 유지되는가?
• RQ4 하이브리드 SAAF–CLS 방법은 $k$-고유값 문제와 공백이 지배하는 문제에서 LS–$P_N$, LS–$S_N$, SAAF–VT–$S_N$ 방법과 비교해 성능가능한가?
• RQ5 균일한 교차단면을 가진 근접 공백 영역으로까지 확장할 수 있으며, 전역 보존성과 정확도를 유지할 수 있는가?

주요 결과

• Conservative LS (CLS) 수식은 공백 영역에서 전역 보존성을 성공적으로 복원하였으며, 스칼라 플럭스 오차는 메쉬 및 각도 이산화의 정밀도 향상에 따라 감소하였고, 특히 $P_N$ 및 $S_N$ 방법에서 두드러졌다.
• SAAF–CLS–$P_N$ 방법은 가장 정밀한 메쉬에서 첫 번째 공간 점에서 오차가 $10^{-2}$ cm⁻²s⁻¹ 이하로 떨어졌으며, SAAF–VT–$S_N$ 방법과 유사한 정확도를 확보하였다.
• SAAF–CLS–$S_N$ 방법은 스트리밍 항의 비대칭 결합 항 수가 적어 $P_N$ 방법에 비해 계산 비용이 크게 낮았다.
• 특히 큰 공백 영역이 포함된 $k$-고유값 문제에서 표준 LS–$P_N$ 및 LS–$S_N$ 방법보다 전역 보존성과 정확도 면에서 뛰어난 성능을 보였다.
• 근접 공백 영역에서도 하이브리드 SAAF–CLS 방법은 전역 보존성을 유지하였으며, SAAF–VT–$S_N$ 방법과 유사한 결과를 도출하여 일반화된 보정 항의 타당성을 입증하였다.
• Hypre BoomerAMG 조절자는 SAAF–VT–$S_N$ 방법에 비해 SAAF–CLS–$S_N$ 방법에 대해 덜 효과적이었는데, 이는 메쉬 정밀도 향상 동안 공백 영역에서 비영인 대각항이 유지되기 때문일 것이다.