[论文解读] PIC simulation of particle behavior in an Inertial Electrostatic Confinement Fusion device using XOOPIC code

D. Bhattacharjee, S. Adhikari+2 more|arXiv (Cornell University)|Feb 14, 2020

Fusion and Plasma Physics Studies被引用 4

一句话总结

本研究使用开源XOOPIC代码对圆柱形惯性电静电约束聚变（IECF）装置中的离子动力学进行粒子-网格（PIC）模拟，揭示了核心区域形成多个势阱，以及离子密度高达10¹⁶ m⁻³的现象。模拟得到的离子能量分布函数与电位分布与实验数据吻合，验证了模型的准确性，并为优化IECF放电参数提供了指导。

ABSTRACT

The kinetic analyses are quite important when it comes to understand the particle behavior in any device as they start to deviate from continuum nature. In the present study, kinetic simulations are performed using Particle-in-Cell (PIC) method to analyze the behavior of ions inside a cylindrical Inertial Electrostatic Confinement Fusion (IECF) device which is being developed as a tabletop neutron source. Here, the lighter ions, like deuterium are accelerated by applying an electrostatic field between the chamber wall (anode) and the cathode (cylindrical gridded wire), placed at the center of the device. These ions recirculate across the cathode grid, which in turn capable of producing fusion reaction at the central region of the device. An open source PIC code (XOOPIC) is used in our study to simulate the ion dynamics at different experimental conditions. The plasma potential profiles obtained from the simulated results indicate the formation of multiple potential well structures inside the cathode grid depending upon the applied cathode potential (from $-1$ to $-5~kV$). The ion density at the core region of the device is found to be of the order of $10^{16}~m^{-3}$, which closely resembles the experimental observations. Spatial variation of Ion Energy Distribution Function (IEDF) has been measured in order to observe the characteristics of ions at different cathode voltages. Finally, the simulated results are compared and found to be in good agreement with the experimental profiles. The present analysis can serve as a reference guide to optimize the technological parameters of the discharge process in IECF devices.

研究动机与目标

• 理解在不同静电场作用下，惯性电静电约束聚变（IECF）装置中离子的动能行为。
• 采用粒子-网格方法，对圆柱形IECF装置中的离子循环与约束动力学进行建模。
• 将模拟得到的电位与离子密度分布与实验观测结果进行对比，以实现装置优化。
• 分析不同阴极电压下离子能量分布函数（IEDF）的空间变化特性。
• 为通过动能模拟调节IECF装置技术参数提供参考框架。

提出的方法

• 使用开源XOOPIC代码，模拟中心阴极网与腔体壁作为阳极的圆柱形IECF装置中的离子动力学。
• 在阴极上施加-1 kV至-5 kV的静电场，以加速氘离子并诱导其循环运动。
• 追踪离子轨迹与电荷累积，以计算等离子体电位分布与离子密度分布。
• 在不同径向位置计算离子能量分布函数（IEDF），以评估离子能量特性。
• 采用PIC方法求解自洽的泊松方程与麦克斯韦方程，以描述等离子体的动能行为。
• 将模拟结果与实验数据进行对比，以验证模型在不同阴极电压下的准确性。

实验结果

研究问题

• RQ1 在IECF装置中，不同阴极电位下，阴极网上方如何形成多个势阱结构？
• RQ2 在不同施加电压下，核心区域的离子密度空间分布如何？
• RQ3 在阴极电压从-1 kV到-5 kV变化时，离子能量分布函数（IEDF）在径向上如何变化？
• RQ4 模拟得到的等离子体电位与离子密度分布与实验测量结果的吻合程度如何？
• RQ5 使用XOOPIC的动能模拟能否准确预测IECF中子源优化的关键性能指标？

主要发现

• 由于施加的阴极电位，阴极网内部形成多个势阱结构，且势阱深度随电压从-1 kV增至-5 kV而增加。
• 核心区域的离子密度达到约10¹⁶ m⁻³，与实验观测结果高度一致。
• 离子能量分布函数（IEDF）表现出显著的空间变化，表明离子在装置半径方向上呈现能量依赖性行为。
• 模拟得到的等离子体电位分布与实验数据在测试电压范围内具有良好的一致性。
• 模拟的离子动力学与能量分布与实验测量结果高度吻合，验证了基于XOOPIC的模型可靠性。
• 本研究为优化IECF装置放电参数以提升中子产额，提供了经过验证的模拟框架。