[논문 리뷰] Primary Neutron Spectra in Ion Vlasov-Fokker-Planck Simulations

B. D. Appelbe, W. T. Taitano외 4명|arXiv (Cornell University)|2023. 05. 03.

Laser-induced spectroscopy and plasma인용 수 4

한 줄 요약

이 논문은 이온의 속도 분포 함수나 미분 단면적을 근사하지 않고, 5차원의 완전한 해상도를 갖춘 수치적 방법을 제안하여 이온 Vlasov-Fokker-Planck 시뮬레이션에서 주로 중성자 스펙트럼을 계산한다. 이 방법은 비맥스웰 분포를 띠는 플라즈마에서 중성자 스펙트럼을 정확하게 시뮬레이션할 수 있게 하여, 비평형 이온 분포로 인한 스펙트럼의 상당한 이탈을 드러내며, 특히 충격 구동형 관성형 핵융합 압축에서 중요한 영향을 미친다. 이로 인해 시뮬레이션과 실험 간의 벤치마킹이 향상된다.

ABSTRACT

The energy spectra of unscattered neutrons produced by deuterium-deuterium and deuterium-tritium fusion reactions are an important diagnostic in High Energy Density Physics experiments as the spectra are sensitive to the velocities of reacting ions. Methods exist for calculating these spectra in radiation-hydrodynamic ("hydro") and Particle-in-Cell ("PiC") simulations. The spectra are particularly sensitive to the high energy tail of ion velocity distribution functions since reaction cross sections increase rapidly with the kinetic energy of a reacting ion pair at the energies achieved in laboratory experiments. This means both the hydro and PiC method may not be suitable in certain plasma regimes. The hydro method assumes that the ion velocity distribution is locally Maxwellian, while the PiC method is subject to statistical noise that makes it challenging to accurately simulate finer details of the spectra. In this work, we present a complementary approach: a method for calculating the neutron spectra in ion Vlasov-Fokker-Planck simulations in which the velocity distribution function is fully-resolved. The method is implemented in the spherically-symmetric code iFP which is used to simulate laser-driven Inertial Confinement Fusion experiments. The method is computationally intensive as it requires a five-dimensional numerical integral, but no approximations of the distribution functions or differential cross sections are required. Results show that deviations of the ion distribution functions from Maxwellian can have a noticeable effect on neutron spectra in shock-driven ICF implosions. The method should facilitate more accurate benchmarking of simulations and experiments.

연구 동기 및 목표

• 이온 분포가 맥스웰 분포에서 벗어날 경우 기존의 유체역학적 및 입자-셀 방법이 중성자 스펙트럼을 시뮬레이션하는 데 가지는 한계를 해결하기 위해.
• 이온 속도 분포 함수와 미분 단면적의 완전한 해상도를 유지하는 고정밀도 주로 중성자 스펙트럼 계산 방법을 개발하기 위해.
• 비맥스웰 플라즈마에서 특히 관성형 핵융합에서 관측되는 비맥스웰 이온 분포를 고려한 고에너지밀도 물리 실험에 대한 정확한 합성 진단을 가능하게 하기 위해.
• 중성자 에너지 스펙트럼 계산에서의 근사 없이 시뮬레이션-실험 간의 정밀한 비교를 향상시키기 위해.

제안 방법

• 반응하는 이온의 속도 분포 함수와 중심질량 프레임에서의 미분 단면적에 대한 5차원 수치적 적분을 통해 중성자 스펙트럼을 계산한다.
• 중심질량 프레임에서 산란 각도余弦에 대한 적분을 재귀 관계를 이용해 분석적으로 해결하며, 이는 $ F_j(x) = \int \frac{x^j}{\sqrt{1-x^2}} dx $ 형태의 적분을 기반으로 한다.
• 입자-셀 방법에서 발생하는 통계적 노이즈를 피하고, 복사유체역학 코드에서처럼 局부 맥스웰 분포를 가정하지 않는다.
• 이 방법은 구형 대칭의 iFP 코드에 구현되었으며, 레이저 구동형 관성형 핵융합 실험을 시뮬레이션한다.
• 핵심 방정식은 스펙트럼 커널을 $ \gamma_a, \gamma_b, \gamma_c $, 및 $ p $ 를 포함하는 형태로 변환하며, $ x_{U,L} $ 로 정의된 범위를 통해 각도 적분을 정확하게 평가할 수 있도록 한다.
• 특수한 경우인 $ v_{cm} = 0 $ 또는 $ \sin\alpha = 0 $ 는 델타 함수나 단순화된 표현식을 통해 분석적으로 처리하여 수치적 적분을 대체한다.

실험 결과

연구 질문

• RQ1 충격으로 구동되는 관성형 핵융합 압축에서 비맥스웰 이온 속도 분포 함수는 주로 중성자 스펙트럼의 형태에 어떤 영향을 미치는가?
• RQ2 완전히 해상도를 확보한 Vlasov-Fokker-Planck 방법은 이온 분포나 단면적을 근사하지 않고도 중성자 스펙트럼을 계산할 수 있으며, 이는 유체역학적 또는 입자-셀 방법과 어떻게 비교되는가?
• RQ3 기본적인 수치적 적분에 비해 제안된 5차원 적분 방법의 계산 비용과 수치 정밀도는 어떠한가?
• RQ4 국부적 맥스웰 분포에서의 이탈이 시뮬레이션과 실험 간의 비교에 사용되는 합성 중성자 스펙트럼에 측정 가능한 영향을 미치는 정도는 어느 정도인가?

주요 결과

• 이 방법은 이온 속도 분포 함수나 미분 단면적을 근사하지 않고 주로 중성자 스펙트럼을 계산할 수 있게 하여 스펙트럼 계산의 완전한 정밀도를 보장한다.
• 충격으로 구동되는 ICF 압축에서 이온 분포 함수가 맥스웰 분포에서 벗어날 경우, 특히 고에너지 尾첨에서 중성자 스펙트럼에 눈에 띄는 영향을 미친다.
• 각도 적분에 대한 분석적 해법은 고에너지 반응(예: $ K \sim 10\,\text{MeV} $, 여기서 $ N \sim 10 $)에서 수치적 적분에 비해 계산 비용을 크게 감소시킨다.
• 이 방법은 디테르륨-디테르륨 및 디테르륨-트리티움 핵융합 반응 모두에 적용 가능하며, 임의의 이온 분포 함수를 지원한다.
• 특수한 경우인 $ v_{cm} = 0 $ 및 $ \sin\alpha = 0 $ 를 포함함으로써 특수한 구조에서 효율적이고 정확한 평가가 가능해져 수치적 안정성이 향상된다.
• 더 현실적인 합성 진단을 제공함으로써, 보다 정밀한 시뮬레이션-실험 비교를 가능하게 한다.