[論文レビュー] Nonclassical Particle Transport in the 1-D Diffusive Limit

Richard Vasques, Rachel Slaybaugh+1 more|arXiv (Cornell University)|Jan 5, 2016

Gas Dynamics and Kinetic Theory参考文献 4被引用数 6

ひとこと要約

本稿では、1次元のランダム周期的拡산系における粒子輸送の非古典的拡散近似を検証し、系が拡散極限に近づくにつれて非古典的輸送方程式の解が非古典的拡散方程式の解に収束することを示している。漸近解析と光学的厚さを変化させた数値シミュレーションを用いて、非指数的パス長分布を示す系、例えばペブルベッド炉のような系において、非古典的拡散モデルが輸送を正確に近似することを確認した。

ABSTRACT

This paper provides numerical results that demonstrate the validity of the nonclassical diffusion approximation to the nonclassical transport equation in certain 1-D diffusive systems. This result provides a more solid foundation in which to improve this theory for relevant nuclear applications.

研究の動機と目的

• 空間的に相関する散乱があるため古典的輸送仮定が成立しない1次元ランダム周期的拡散系における非古典的粒子輸送を調査すること。
• 拡散極限において非古典的輸送方程式が非古典的拡散方程式でよく近似されることを理論的予測するのを検証すること。
• 非古典的ボルツマン方程式の数値的解法を通じて、非古典的輸送の漸近理論を裏付ける最初の計算的証拠を提供すること。
• 光学的厚さが増加するにつれて非古典的スカラーフラックスが非古典的拡散解に収束することを示すこと。
• ペブルベッド炉のような相関したマイクロ構造を持つ核系の粒子輸送のより良いモデル化の基盤を確立すること。

提案手法

• パス長依存の全断面積 $\Sigma_t(s)$ を持つ1次元非古典的線形ボルツマン方程式を定式化し、非指数的衝突統計を考慮する。
• 小さなパラメータ $\varepsilon \ll 1$ を用いた漸近解析を適用し、$1-c = O(\varepsilon^2)$、$Q = O(\varepsilon^2)$、$\partial_x \psi = O(\varepsilon)$ のスケーリングを行い、$\varepsilon$ に関する正則摂動展開を導出する。
• $s$ と $x$ の依存性を分離するために変換されたフラックス $\hat{\psi}^\pm(x,s)$ を導入し、非古典的拡散方程式を一次近似として導出する。
• 擬似時間変数 $s$ と空間変数 $x$ に対して、1次精度のHLL有限体積法を用いて非古典的輸送方程式を数値的に解き、ソース反復法を用い、$s_{\text{max}} = 45$ で打ち切る。
• 光学的厚さが増加するように $M = 10, 15, 20$ の系について、アンサンブル平均スカラーフラックス $\Phi(x) = \int_0^{45} [\psi^+(x,s) + \psi^-(x,s)] ds$ を計算する。
• 拡散係数 $D = \langle s^2 \rangle / (2\langle s \rangle)$ を用いた端点補外境界条件を用いて非古典的拡散方程式を解く。

実験結果

リサーチクエスチョン

• RQ1 非古典的拡散方程式は、1次元拡散極限において非古典的輸送方程式の解を正確に近似するか？
• RQ2 非古典的スカラーフラックスが非古典的拡散解に収束する様子は、光学的厚さにどのように依存するか？
• RQ3 パス長分布 $p(s)$ は非古典的拡散近似の有効性を決定づける役割を果たすか？
• RQ4 非古典的輸送方程式の数値的解法は、相関したマイクロ構造を持つ系において、漸近理論の予測を確認できるか？
• RQ5 非古典的輸送シミュレーションにおいて、散乱比 $c$ がソース反復の収束速度と数値的安定性に与える影響は何か？

主な発見

• $M$ が増加するにつれて非古典的スカラーフラックスが非古典的拡散解に収束する。これは拡散極限における漸近理論の妥当性を確認するものである。
• $M=20$ の場合、非古典的スカラーフラックスと非古典的拡散フラックスは良好な一致を示し、$\langle s \rangle = 2$ および $\langle s^2 \rangle \approx 8.1640$ から、拡散係数 $D \approx 2.041$ が得られた。
• ソース反復の収束には $c = 0.9975$（$M=20$）で936反復が必要であり、理論で予測されたように収縮率が散乱比 $c$ に比例することが確認された。
• 1次元周期的系に対してパス長分布 $p(s)$ が解析的に導出され、長さ2の区間ごとに交互に指数関数的セグメントを持つ非指数的挙動を示した。
• 空間的に相関する散乱子を持つ系、例えばペブルベッド炉において、非古典的拡散方程式が有効な近似であることが示された。古典的輸送は非指数的減衰のため失敗するが、非古典的拡散はその問題を克服する。
• 数値結果により、相関構造を持つランダム媒体において非古典的拡散モデルが正しいフラックス挙動を捉えていることが確認され、このような領域では古典的拡散よりもより正確な代替手段であることが示された。