[논문 리뷰] Numerical modeling of neutron transport in SP3 approximation by finite element method

A. V. Avvakumov, V. F. Strizhov외 2명|arXiv (Cornell University)|2019. 03. 27.

Nuclear reactor physics and engineering참고 문헌 19인용 수 4

한 줄 요약

이 논문은 Reactor core 분석을 위한 SP₃ 중성자 이동 근사의 유한요소법(FEM) 구현을 제시하며, 계산 비용을 최소화하면서도 확산 이론에 비해 정확도를 향상시킴을 보여준다. IAEA-2D 및 HWR 기준 문제에서 방법을 검증한 결과, 특정 구성에서 확산 이론과 SP₃ 해 모두 복소수 고유값을 유도할 수 있으며, 특히 HWR의 경우가 그러하다. 그러나 k- 및 α-고유값에 대해서는 기준 운반 해석 결과와 근접한 일致성을 유지함을 확인하였다.

ABSTRACT

The SP3 approximation of the neutron transport equation allows improving the accuracy for both static and transient simulations for reactor core analysis compared with the neutron diffusion theory. Besides, the SP3 calculation costs are much less than higher order transport methods (SN or PN). Another advantage of the SP3 approximation is a similar structure of equations that is used in the diffusion method. Therefore, there is no difficulty to implement the SP3 solution option to the multi-group neutron diffusion codes. In this work, the application of the SP3 methodology based on solution of the λ- and α-spectral problems has been tested for the IAEA-2D and HWR reactor benchmark tests. The FEM is chosen to achieve the 3D geometrical generality, using GMSH as a generic mesh generator. The results calculated with the diffusion and SP3 methods are compared with the reference transport calculation results. It was found for the HWR reactor test that some eigenvalues are complex when calculating using both diffusion and SP3 options.

연구 동기 및 목표

• 3D 기하 구조에서 중성자 이동 방정식의 SP₃ 근사를 해결하기 위한 유한요소법(FEM) 개발 및 검증.
• 전체 핵심 반응로 시뮬레이션에서 SP₃의 정확도 및 계산 효율성과 중성자 확산 이론 간의 비교 평가.
• 반응로 스펙트럼 문제에서 확산 및 SP₃ 해에 복소수 고유값이 발생하는 조건과 그 물리적 의미 탐구.
• 기준 운반 해석 결과와의 비교를 통해 기준 문제에서의 k-효율성 및 α-고유값에 대한 SP₃ 결과 분석.
• 방정식의 구조 유사성 유지를 통해 기존 확산 코드 프레임워크에 SP₃ 통합 가능하도록 설계.

제안 방법

• 구면 조화 함수 전개를 이용해 중성자 이동 방정식을 두 개의 결합된 확산형 방정식으로 감소시킨 SP₃ 근사 적용.
• 3D 기하 구조의 유연성 확보를 위해 GMSH를 사용한 메시 생성을 통한 고차수 라그랑주 유한요소 방법을 이용한 공간 이산화.
• 대규모 고유값 계산을 위한 SLEPc 라이브러리를 활용한 λ- 및 α-스펙트럼 문제 해법.
• 고유값 계산의 정확도 검증 및 수렴 제어를 위한 압축 메시 사용.
• 실제 반응로 동역학 파rameter를 반영한 프롬프트 및 지연 중성자 모델 구현을 통한 α-고유값 계산.
• 정밀한 메시 운반 계산에서 유도된 고밀도 기준 해석 결과와 SP₃ 및 확산 방법의 결과 비교.

실험 결과

연구 질문

• RQ1 전체 핵심 반응로 시뮬레이션에서 SP₃ 근사가 확산 이론에 비해 정확도를 크게 향상시키며, 계산 비용을 최소화할 수 있는가?
• RQ2 IAEA-2D 및 HWR와 같은 표준 기준 문제에서 SP₃ 방법이 기준 운반 결과의 k-효율성 및 α-고유값을 재현할 수 있는가?
• RQ3 확산 및 SP₃ 해에서 복소수 고유값이 어떤 조건에서 발생하며, 그 물리적 의미는 무엇인가?
• RQ4 대칭 및 비대칭 구성에서 SP₃의 스펙트럼 특성(고유값 및 고유함수)은 확산 이론과 어떻게 비교되는가?
• RQ5 구조적 유사성 덕분에 SP₃는 기존 확산 코드 아키텍처에 얼마나 원활하게 통합될 수 있는가?

주요 결과

• SP₃ 방법은 계산 시간을 확산 이론의 약 1.5배로 줄이며, 특히 핀 단위 및 제어봉 반응도 계산에서 정확도 향상을 보였다.
• HWR 기준 문제에서 확산 및 SP₃ 해 모두 작은 허수부를 가진 복소수 고유값을 생성하여 중성자 플럭스 동역학에서 진동하는 행동을 나타냄.
• 반면 IAEA-2D 기준 문제에서는 비대칭 케이스조차도 실수 고유값만 나타내어, 복소수 고유값은 구성 조건에 따라 달라짐을 시사함.
• SP₃ 방법의 기본 α-고유값은 고해상도 메시(예: n=96, p=3)에서 기준값 0.04395 s⁻¹과 매우 유사한 0.04395 s⁻¹을 기록하였으며, 상대 오차는 0.1% 미만.
• 프롬프트 중성자 생성 시간 Λ_pr는 확산 및 SP₃ 모두에서 일관되게 1.917×10⁻⁴ s로 추정되어, 동역학적 행동 모델링의 일致성 확인.
• 소수 고유값(예: α₂ ~ α₅, α₉ ~ α₁₀)은 두 방법 모두에서 작은 허수 성분을 가진 복소수 쌍을 나타내어 중성자 플럭스에서 감쇠 또는 증폭되는 진동 모드 존재를 시사함.