[논문 리뷰] Simplified P$_N$ Equations for Nonclassical Transport with Isotropic Scattering

Richard Vasques, Rachel Slaybaugh|arXiv (Cornell University)|2016. 10. 14.

Nuclear reactor physics and engineering참고 문헌 6인용 수 5

한 줄 요약

이 논문은 등방성 산산산란을 가진 비고전적 운반 방정식에 대해 점근 분석을 사용하여 새로운 확산 근사—비고전적 SP_N 방정식—을 유도한다. 이 방법은 자유경로 분포의 모멘트를 통합함으로써 고전적 단순화된 P_N (SP_N) 방정식을 비고전적 매질로 일반화하며, 수치 결과는 비고전적 SP_N 해가 시스템이 확산 한계에 접근할수록 기준 결과로 수렴하고, N과 시스템 크기 M이 증가함에 따라 정확도가 향상됨을 확인한다.

ABSTRACT

An asymptotic analysis is used to derive a set of diffusion approximations to the nonclassical transport equation with isotropic scattering. These approximations are shown to reduce to the simplified P$_N$ equations under the assumption of classical transport, and therefore are labeled nonclassical SP$_N$ equations. In addition, the nonclassical SP$_N$ equations can be manipulated into a classical form with modified parameters, which can be implemented in existing SP$_N$ codes. Numerical results are presented for an one-dimensional random periodic system, validating the theoretical predictions.

연구 동기 및 목표

• 공간적으로 상관된 산산산란 중심을 가진 랜덤 매질에 적용 가능한, 등방성 산산산란을 가진 비고전적 운반 방정식에 대한 고차수 점근 근사법을 개발하는 것.
• 자유경로 분포의 모멘트를 통합하여 고전적 단순화된 P_N (SP_N) 프레임워크를 비고전적 운반에 확장하는 것.
• 비고전적 SP_N 방정식을 고전적 형태로 재구성함으로써 기존 SP_N 코드와의 호환성을 확보하는 것.
• 1차원 랜덤 주기적 시스템에서 앙상블 평균 기준 솔루션을 사용하여 이론적 프레임워크를 수치적으로 검증하는 것.
• 향후 다차원 시스템, 각도 의존 자유경로 분포, 비등방성 산산산란으로의 확장을 위한 기반을 마련하는 것.

제안 방법

• 소수의 매개변수 ε를 사용하여 비고전적 운반 방정식에 점근 전개를 수행하고, 총 단면적, 산산산란 비율, 소스를 스케일링하여 고차수 근사를 유도한다.
• ε² 항까지의 점근 급수 전개에서 항을 일치시켜 비고전적 SP_1 (확산), SP_2, SP_3 방정식의 명시적 형태를 유도한다.
• 비고전적 SP_N 방정식이 N = M일 때 유효하기 위해 자유경로 분포 p(s)의 첫 2M차 순모멘트가 유한해야 한다는 조건이 필요하다.
• 비고전적 SP_N 방정식을 수정된 총 단면적과 흡수 계수를 가진 고전적 SP_N 형태로 재구성하여 표준 Marshak 진공 경계 조건을 적용할 수 있도록 한다.
• 층이 번갈아가며 존재하는 1차원 랜덤 주기적 시스템을 사용하며, 층 두께 ℓ = 0.5는 평균 자유경로 척도를 보장하고, M = 1/ε는 확산 한계를 제어한다.
• 랜덤 시스템의 다수의 물리적 실현을 통해 몬테카를로 방식의 평균을 취하여 기준 솔루션을 생성하며, S_16 이산 순서와 다이아몬드 차분을 사용하여 운반 방정식을 수치적으로 해석한다.

실험 결과

연구 질문

• RQ1 등방성 산산산란을 가진 비고전적 운반 방정식에 대한 고차수 점근 전개가 고전적 SP_N 방정식을 일반화하는 일관된 확산 근사 집합을 도출할 수 있는가?
• RQ2 자유경로 분포가 지수 분포로 감소할 경우 비고전적 SP_N 방정식은 어떻게 행동하는가?
• RQ3 비고전적 SP_N 방정식을 수정된 매개변수를 가진 고전적 SP_N 형태로 재구성할 수 있는가? 이를 통해 기존 SP_N 코드에 구현이 가능한가?
• RQ4 점근 이론이 예측한 바와 같이, 시스템 크기가 증가하고 확산 한계에 접근함에 따라 비고전적 SP_N 근사의 정확도가 향상되는가?
• RQ5 1차원 랜덤 주기적 시스템에서 앙상블 평균 기준 솔루션에 대한 비고전적 SP_N 해의 수렴 거동은 어떠한가?

주요 결과

• 비고전적 SP_N 방정식은 등방성 산산산란을 가진 비고전적 운반 방정식에 대한 고차수 점근 근사로 유도되었으며, 이는 자유경로 분포의 첫 2N차 순모멘트가 유한할 경우에 유효하다.
• 자유경로 분포가 지수 분포일 경우(고전적 운반의 경우), 비고전적 SP_N 방정식은 정확히 고전적 SP_N 방정식으로 감소하여 기존 이론과의 일관성을 확인한다.
• 비고전적 SP_N 방정식은 수정된 총 단면적과 흡수 계수를 가진 고전적 SP_N 형태로 변환될 수 있으며, 이는 표준 Marshak 진공 경계 조건을 적용하고 기존 SP_N 코드에 통합할 수 있음을 의미한다.
• 1차원 랜덤 주기적 시스템에 대한 수치 결과는 x=0에서 스칼라 플럭스 추정치의 상대 오차가 M 증가함에 따라 감소함을 보여주며, 확산 한계로의 수렴을 확인한다.
• 각 N에 대해 오차는 M 증가함에 따라 감소하고, 더 높은 차수의 SP_N 근사(예: SP_3)는 더 낮은 차수의 근사(예: SP_1)보다 낮은 오차를 보이며, 점근 이론이 예측한 정확도 향상 현상이 검증된다.
• 기준 솔루션은 랜덤 시스템의 다수의 실현에 대해 앙상블 평균을 취함으로써 확보되었으며, S_16 이산 순서와 다이아몬드 차분을 사용하여 수치적으로 해석되었으며, 오차 평가의 기준으로 사용된다.