[논문 리뷰] On the Accuracy of the Non-Classical Transport Equation in 1-D Random Periodic Media

Richard Vasques, Kai Krycki|arXiv (Cornell University)|2014. 12. 09.

Advanced Mathematical Modeling in Engineering참고 문헌 10인용 수 3

한 줄 요약

이 논문은 입자 충돌 거리가 공간 상관성으로 인해 지수 분포가 아닌 분포를 따르는 1차원 무작위 주기적 미디어에서 비고전적 운반 방정식의 수치적 검증을 처음으로 제시한다. 분석적으로 경로 길이 분포 함수 $ p(s) $ 를 유도하고, 이에 대응하는 $ \Sigma_t(s) $ 를 계산하며, 비고전적 운반 방정식을 해결함으로써 저산란 영역 및 소스 근처에서 원자 혼합 모델에 비해 더 높은 정확도를 보임을 입증한다.

ABSTRACT

We present a first numerical investigation of the accuracy of the recently proposed {\em non-classical transport equation}. This equation contains an extra independent variable (the path-length $s$), and models particle transport taking place in random media in which a particle's distance-to-collision is {\em not} exponentially distributed. To solve the non-classical equation, one needs to know the $s$-dependent ensemble-averaged total cross section $Σ_t(s)$, or its corresponding path-length distribution function $p(s)$. We consider a 1-D spatially periodic system consisting of alternating solid and void layers, randomly placed in the infinite line. In this preliminary work, we assume transport in rod geometry: particles can move only in the directions $μ=\pm 1$. We obtain an analytical expression for $p(s)$, and use this result to compute the corresponding $Σ_t(s)$. Then, we proceed to solve the non-classical equation for different test problems. To assess the accuracy of these solutions, we produce "benchmark" results obtained by (i) generating a large number of physical realizations of the system, (ii) numerically solving the transport equation in each realization, and (iii) ensemble-averaging the solutions over all physical realizations. We show that the results obtained with the non-classical equation accurately model the ensemble-averaged scalar flux in this 1-D random system, generally outperforming the widely-used atomic mix model. We conclude by discussing plans to extend the present work to slab geometry, as well as to more general random mixtures.

연구 동기 및 목표

• 공간적으로 상관된 산란 중심이 존재하는 1차원 무작위 주기적 미디어를 통과하는 입자 운반을 모델링하는 데 있어 비고전적 운반 방정식의 정확도를 평가하는 것.
• 교대되는 고체 및 진공 층으로 이루어진 1차원 시스템에서 경로 길이 분포 함수 $ p(s) $ 의 해석적 표현을 유도하는 것.
• 유도된 $ \Sigma_t(s) $ 를 사용하여 비고전적 운반 방정식을 수치적으로 해결하고, 다수의 물리적 실현을 통해 얻은 평균 결과와 비교하는 것.
• 비고전적 모델이 스칼라 플럭스 예측에서 널리 사용되는 원자 혼합 모델에 비해 어떻게 성능을 발휘하는지 평가하는 것.
• 향후 슬래브 기하구조 및 더 일반적인 무작위 미디어로의 확장에 기초를 마련하는 것.

제안 방법

• 무한선 상에 무작위로 위치한 교대되는 고체 및 진공 층으로 이루어진 1차원 공간 주기적 시스템을 모델링하는 것.
• 입자가 방향 $ \mu = \pm 1 $ 에서만 운동하는 막대 기하구조를 가정하여 각도 의존성을 단순화하는 것.
• 재생 이론과 재생 보상 정리를 사용하여 경로 길이 분포 함수 $ p(s) $ 를 해석적으로 유도하는 것.
• 항등식 $ \Sigma_t(s) = p(s) / \left(1 - \int_0^s p(s') ds' \right) $ 를 사용하여 $ p(s) $ 에서 $ s $-의존성 총 단면적 $ \Sigma_t(s) $ 를 계산하는 것.
• 유한 체적 방법과 HLL 리만 해법을 사용하여 비고전적 운반 방정식을 수치적으로 해결하는 것.
• 대규모 물리적 실현의 앙상블을 시뮬레이션하여 각각의 경우에 고전적 운반 방정식을 풀고 스칼라 플럭스를 평균내어 기준 솔루션을 생성하는 것.

실험 결과

연구 질문

• RQ1 공간적으로 상관된 산란이 존재하는 1차원 무작위 주기적 미디어에서 비고전적 운반 방정식은 앙상블 평균 스칼라 플럭스를 정확하게 예측할 수 있는가?
• RQ2 산란 비율과 소스 구성이 다양할 경우 비고전적 모델의 정확도는 원자 혼합 모델에 비해 어떻게 다를까?
• RQ3 비고전적 모델은 원자 혼합 모델에 비해 소스 근처에서 스칼라 플럭스의 진동 구조를 더 잘 유지하는가?
• RQ4 수치적 확산은 비고전적 모델과 원자 혼합 모델 간 상대적 성능에 어떤 영향을 미치는가?
• RQ5 비주기적 무작위 미디어로의 $ p(s) $ 의 해석적 유도는 확장 가능할까?

주요 결과

• 비고전적 운반 방정식은 1차원 무작위 주기적 미디어에서 앙상블 평균 스칼라 플럭스를 정확하게 모델링하며, 앙상블 평균 물리적 실현의 기준 결과와 매우 유사한 해를 제공한다.
• 저산란 비율(예: $ c = 0.1 $)에서 비고전적 모델은 원자 혼합 모델을 능가하며, 특히 소스 근처에서 최대 플럭스를 잘 포착하고 해의 진동형태를 잘 유지한다.
• 문제 세트 A에서 $ c = 0.1 $ 인 경우, 그림 4에 나타나 있듯이 비고전적 모델은 원자 혼합 모델 대비 절대값 기준 상대 오차를 최대 10% 감소시킨다.
• 원자 혼합 모델은 산란 비율 $ c $ 가 증가함에 따라 정확도가 향상되지만, 비고전적 모델의 정확도는 동일 조건에서 약간 감소하는 경향이 있다. 이는 HLL 스킴에서의 수치적 확산으로 인한 것으로 보인다.
• 소스에서 떨어진 영역에서는 특히 저산란 조건에서 비고전적 모델이 기준 결과와 더 나은 정성적 일치를 유지한다.
• 결과는 비고전적 운반의 이론적 프레임워크를 검증하며, 고유값과 고유함수 추정에서 더 높은 정확도가 요구되는 임계성 계산에 그 적용을 지지한다.