超热电子在稠密等离子体中输运的混合粒子模拟方法

相对论激光等离子体相互作用以及所产生的带电粒子束在高密度等离子体中的输运行为非常重要.该物理问题的数值模拟研究仍面临技术挑战.本文介绍一种粒子/流体混合模拟方法.该方法中超热电子采用动力学方法描述,背景冷的稠密等离子体采用简化的流体方程描述,适合于超热电子密度远小于背景电子密度,超热电子能量远大于背景电子温度.我们的三维并行混合模拟程序HEETS的模拟结果表明：背景材料的电离和电阻率模型至关重要,将严重影响高能电子输运过程的模拟.     

���ӻ��������ȵ��������еĹ��ʳ���

用有限拉莫尔半径(FLR)修正的相对论色散关系进行了数值模拟，研究了电子回旋波在高温、高密等离子体的功率沉积，并将计算结果与应用弱相对论Fokker-Planck方程得到的结果进行了比较。结果表明，在高温、高密等离子体中，波功率的吸收非常集中；平行折射率、极向发射位置和发射波频率的变化都会影响波功率沉积的大小和分布；平行折射率变大后，FLR效应会使波的阻尼减少。     

电子束在相对论返波管振荡器收集极上的能量沉积

在圆周对称的磁场作用下,环形电子束以一定角度轰击在圆柱面的相对论返波管振荡器(RBWO)收集极上并将能量沉积其中,采用蒙特卡罗程序FLUKA,建立了电子的能量沉积分布计算模型,研究了电子能量沉积分布规律;建立了背散射电子的运动模型,模拟了磁场作用下背散射电子的运动轨迹;研究了圆周对称径向磁场的近似方法。研究结果表明:随着磁场强度的增大,最大能量沉积密度增大,背散射电子在更靠近电子束入射区域的位置再次入射并沉积能量,且可能形成一个新的能量沉积峰值。在磁场强度较大时,采用单向的径向磁场即可较好地计算圆周对称径向磁场下背散射电子的能量沉积分布。     

Efective Parameters on Energy Losses of Fast Electron in Fusion Mediums

The stopping and scattering of fast electrons in a dense plasma relevant to inertial confinement fusion(ICF)are investigated numerically with the latest improved cross section equations.Binary and collective efects are considered to determine beam transport parameters such as range,penetration depth,spreading processes as straggling and blooming versus electron energy and plasma parameters.Blooming and straggling efects,which act as consequences of scattering with statistical assumption in collisions,lead to a non-uniform,extended region of energy deposition.Finally the mean angle of deflections is calculated for diferent plasma energies.     

Effective Parameters on Energy Losses of Fast Electron in Fusion Mediums

The stopping and scattering of fast electrons in a dense plasma relevant to inertial confinement fusion (ICF) are investigated numerically with the latest improved cross section equations. Binary and collective effects are considered to determine beam transport parameters such as range, penetration depth, spreading processes as straggling and blooming versus electron energy and plasma parameters. Blooming and straggling effects, which act as consequences of scattering with statistical assumption in collisions, lead to a non-uniform, extended region of energy deposition. Finally the mean angle of deflections is calculated for different plasma energies.     

Heating mechanisms in short-pulse laser-driven cone targets

The fast ignitor is a modern approach to laser fusion that uses a short-pulse laser to initiate thermonuclear burn. In its simplest form the laser launches relativistic electrons that carry its energy to a precompressed fusion target. Cones have been used to give the light access to the dense target core through the low-density ablative cloud surrounding it. Here the ANTHEM implicit hybrid simulation model shows that the peak ion temperatures measured in recent cone target experiments arose chiefly from return current joule heating, mildly supplemented by relativistic electron drag. Magnetic fields augment this heating only slightly, but capture hot electrons near the cone surface and force the hot electron stream into filaments.     

Angular distribution of electrons in the field of a short laser pulse of relativistic intensity

Based on the dynamics of the electron in the field of a laser pulse of relativistic intensity, the electron emission angles with respect to the laser pulse propagation axis are calculated. The angular distribution of accelerated electrons is analyzed together with their energy spectrum. It is shown that fast electrons forming the high-energy spectral region are emitted into a fixed angular cone.     

电子束能量沉积的相对论Fokker-Planck方程的计算方法

针对相对论电子束在高密度等离子体中的能量沉积过程,建立三维动量空间中快电子能量沉积的相对论Fokker-Planck方程的可计算物理模型,构造数值算法并研制数值模拟程序.通过与解析模型和蒙特卡罗模拟相比较,验证数值方法和程序的可靠性.在二维动量空间的模拟基础上,通过计算能量区间为0.5 MeV~3.5 MeV的快电子在背景密度为300 g·cm^-3的氘氚等离子体中的能量沉积过程,发现由于碰撞效应使平均散射角趋近平衡,三维动量空间计算快电子连续射程和穿透深度与二维结果基本一致.     

Collective deceleration of relativistic electrons precisely in the core of an inertial-fusion target

The energy deposition of a relativistic electron beam in a plasma can be managed through turning on or off fast beam-plasma instabilities in desirable regions. This management may enable new ways of realizing the fast-igniter scenario of inertial fusion. Collisional effects alone can decelerate electrons of at most a few MeV within the core of an inertial-fusion target. Beam-excited Langmuir turbulence, however, can decelerate even ultrarelativistic electrons in the core.     

爆推快点火模型分析

 使用简单的爆推模型估算爆推快点火过程及其结果。 首先由ns级主驱动激光直接驱动,形成中心低密度高温热斑,周围为高密度低温主燃料区, 两区压力平衡(等压模型);然后用ps级超短超强激光打入,产生超热电子,其能量在低温主燃料区沉积,主燃料区发生爆炸,一部分向外飞散,一部分向内压缩中心热斑。在这个爆推模型下,热斑体积压缩比为64,中子产率将有极大的提高,相应的中子产额和能量增益得到提高。离子温度因为主燃料区质量过大,提高不大。提高超热电子能量,或者减小低温主燃料区质量,离子温度将会显著提高。不同的初始离子温度对结果有很大影响,较低的初始温度下更容易得到较高的中子产率和产额。     

The impact Of impurity ion in deuterium-tritium fuel on the energy deposition pattern Of The Proton Ignitor Beam

The ignition stage of deuterium-tritium fuel in inertial confinement fusion is a challenging task affected by many undesirable processes especially material mixing processes in the hot-spot region. In this research, an alternative proposal of the enhanced energy deposition in the proton fast ignition has been suggested. It consists of two primary assumptions of the beam-plasma system. In the first place, we have adopted the proton beam generated by TNSA or RPA mechanisms, each described by a Maxwellian or Gaussian energy distributions. Next, a realistic, non-uniform fuel plasma was adopted. Then, the cumulative stopping power of a proton beam of 10 kJ energy, penetrating the low content metal-contaminated deuterium-tritium fuel has been examined. It has been shown that in the case of the very low impurity fractions, irregular spatial fluctuations in the cumulative stopping power relative to pure fuel plasma emerges. However, at the higher concentrations, a systematic pattern becomes visible such that the contribution of the deep layers in the stopping power reduces. We observe the enhanced energy deposition close to the corona/dense core interface. It has been shown that the corona/dense core energy deposition ratio differs by up to 2.5% between pure and contaminated DT plasma. In the contaminated fuel plasma, energy deposition in the TNSA regime will effectively heat the plasma corona. While in the RPA counterpart, at a similar level of contamination, most of the incident beam energy remains inside the core fuel region.     

快点火中相对论电子束能量沉积的动理学研究

针对相对论快电子束在高密度压缩芯区等离子体中的能量沉积过程开展物理建模、程序研制和数值模拟研究。从等离子体粒子碰撞的基本物理出发,综合考虑了高能电子与背景等离子体之间的短程两体碰撞过程和长程集体效应,建立了相对论Fokker-Planck动理学模型,通过采用球谐展开的方法,推导得到了适于数值求解的方程形式并根据方程特点开展相应的数值算法研究及程序研制并完成了物理考核,对快点火能量沉积的典型物理算例进行了模拟研究,并针对即将在神光Ⅱ升级装置上开展的快点火物理实验进行了初步的物理分析。     

High-frequency surface waves in quantum plasmas with electrons relativistic degenerate and exchange-correlation effects

The propagation characteristics of high-frequency surface waves are studied in spin-1/2 quantum plasmas by considering the electron relativistic degenerate and exchange-correlation effects. Using the quantum fluid equations of magnetoplasmas in the presence of the quantum Bohm potential, spin magnetization energy, relativistic degenerate pressure, and exchange-correlation effects, a generalized dispersion relation is derived. The analytical and numerical results show that the relativistic degenerate and exchange-correlation effects significantly modify the propagation properties of high-frequency surface waves. It is found that under the influence of exchange-correlation effects, the frequency spectrum of high-frequency surface waves will be down-shifted. It is also indicated that the dispersion curve shifts up with the increase of relativistic gamma factor. Furthermore, the phase speed of the high-frequency surface waves increases with increasing electron number density. The current research is helpful to understand the propagation of the high-frequency surface waves in quantum plasmas, such as those in dense astrophysical environment.     

Monte Carlo simulation of electron beam air plasma characteristics

A high-energy electron beam generator is used to generate a plasma in atmosphere. Based on a Monte Carlo toolkit named GEANT4, a model including complete physics processes is established to simulate the passage of the electron beam in air. Based on the model, the characteristics of the electron beam air plasma are calculated. The energy distribution of beam electrons (BEs) indicates that high-energy electrons almost reside in the centre region of the beam, but low-energy electrons always live in the fringe area. The energy deposition is calculated in two cases, i.e., with and without secondary electrons (SEs). Analysis indicates that the energy deposition of SEs accounts for a large part of the total energy deposition. The results of the energy spectrum show that the electrons in the inlet layer of the low-pressure chamber (LPC) are monoenergetic, but the energy spectrum of the electrons in the outlet layer is not pure. The SEs are largely generated at the outlet of the LPC. Moreover, both the energy distribution of BEs and the magnitude of the density of SEs are closely related to the pressure of LPC. Thus, a conclusion is drawn that a low magnitude of LPC pressure is helpful for reducing the energy loss in the LPC and also useful for greatly increasing the secondary electron density in dense air.     

Fast electron energy deposition in aluminium foils: Resistive vs. drag heating

The high current electron beam losses have been studied experimentally with 0.7 J, 40 fs, 6 10¹⁹ Wcm^-2 laser pulses interacting with Al foils of thicknesses 10-200 μm. The fast electron beam characteristics and the foil temperature were measured by recording the intensity of the electromagnetic emission from the foils rear side at two different wavelengths in the optical domain, ≈407 nm (the second harmonic of the laser light) and ≈500 nm. The experimentally observed fast electron distribution contains two components: one relativistic tail made of very energetic (T _h ^tail ≈ 10 MeV) and highly collimated (7° ± 3°) electrons, carrying a small amount of energy (less than 1% of the laser energy), and another, the bulk of the accelerated electrons, containing lower-energy (T _h ^bulk=500 ± 100 keV) more divergent electrons (35 ± 5°), which transports about 35% of the laser energy. The relativistic component manifests itself by the coherent 2ω₀ emission due to the modulation of the electron density in the interaction zone. The bulk component induces a strong target heating producing measurable yields of thermal emission from the foils rear side. Our data and modeling demonstrate two mechanisms of fast electron energy deposition: resistive heating due to the neutralizing return current and collisions of fast electrons with plasma electrons. The resistive mechanism is more important at shallow target depths, representing an heating rate of 100 eV per Joule of laser energy at 15 μm. Beyond that depth, because of the beam divergence, the incident current goes under 10¹² Acm^-2 and the collisional heating becomes more important than the resistive heating. The heating rate is of only 1.5 eV per Joule at 50 μm depth.     

Particle-in-cell simulations with charge-conserving current deposition on graphic processing units

We present an implementation of a 2D fully relativistic, electromagnetic particle-in-cell code, with charge-conserving current deposition, on parallel graphics processors (GPU) with CUDA. The GPU implementation achieved a one particle-step process time of 2.52 ns for cold plasma runs and 9.15 ns for extremely relativistic plasma runs, which are respectively 81 and 27 times faster than a single threaded state-of-art CPU code. A particle-based computation thread assignment was used in the current deposition scheme and write conflicts among the threads were resolved by a thread racing technique. A parallel particle sorting scheme was also developed and used. The implementation took advantage of fast on-chip shared memory, and can in principle be extended to 3D.     

Experiments on two-step heating of a dense plasma in the GOL-3 facility

This paper presents the results of experiments on two-stage heating of a dense plasma by a relativistic electron beam in the GOL-3 facility. A dense plasma with a length of about a meter and a hydrogen density up to 10¹⁷ cm⁻³ was created in the main plasma, whose density was 10¹⁵ cm⁻³. In the process of interacting with the plasma, the electron beam (1 MeV, 40 kA, 4 μs) imparts its energy to the electrons of the main plasma through collective effects. The heated electrons, as they disperse along the magnetic field lines, in turn reach the region of dense plasma and impart their energy to it by pairwise collisions. Estimates based on experimental data are given for the parameters of the flux of hot plasma electrons, the energy released in the dense plasma, and the energy balance of the beam-plasma system. The paper discusses the dynamics of the plasma, which is inhomogeneous in density and temperature, including the appearance of pressure waves. Zh. éksp. Teor. Fiz. 113, 897–917 (March 1998)     

激光驱动强流电子束产生和控制

在惯性约束聚变(ICF)电子束快点火物理方案中,需要超强拍瓦激光脉冲驱动MeV能量的强流电子束,并沉积数十kJ能量到压缩氘氚芯区。强流电子束的束流品质是影响点火成功的关键因素之一,为深入了解强流电子束产生物理过程,研制成了三维高性能、适应上万CPU核规模的并行粒子模拟程序,并开展了大规模数值模拟研究,探索了强流电子束的产生机制和输运规律。回顾了近几年来快点火研究团队围绕强流电子束产生和控制开展的研究,介绍了导致束流品质差的两大物理原因:预等离子体效应和束流不稳定性磁场的随机散射。针对这两个物理原因,提出了四种提高强流电子束品质的方法:(1)双层金锥靶减弱预等离子体的负面效应;(2)输运丝产生环向磁场准直强流电子束;(3)外加磁场导引强流电子束提高耦合效率;(4)抑制束流不稳定性以降低随机磁场对电子束流的散射。     

Retarded dielectric theory of electron energy loss in multilayers and superlattices

This work uses computer algebra to develop the dielectric theory of energy loss for fast electrons travelling in multilayer systems and superlattices. The interaction of a relativistic electron beam is analysed in detail for beams travelling normal or parallel to the interfaces.     

Resistive collimation of electron beams in laser-produced plasmas

Intense relativistic electron beams, produced by high-intensity short-pulse laser irradiation of a solid target, have many potential applications including fusion by fast ignition. Using a unique Fokker-Planck code, supported by analytic calculations, we show that fast electrons can be collimated into a beam even when the fast electron source is not strongly anisotropic, and we derive a condition for collimation to occur.