
Citation: | Sing LEE. Pushing the limits of existing plasma focus towards 1016 fusion neutrons with Q = 0.01[J]. Plasma Science and Technology, 2022, 24(11): 114005. DOI: 10.1088/2058-6272/ac78cc |
Existing conventional megajoule plasma focus machines with 2–3 MA are producing fusion neutron yields of several times 1011 in deuterium operation, the fusion yields predominantly being the beam-gas target. Increasing the current to 10 MA and using 50%–50% D-T mixture will scale the neutron yield towards 1016 D-T fusion neutrons. In this work, we derive the Lawson criterion for plasma focus devices with a beam-target fusion neutron mechanism, so that we may glimpse what future technological advancements are needed for a break-even Q = 1 plasma focus. We perform numerical experiments with a present-day feasible 0.9 MV, 8.1 MJ, 11 MA machine operating in 100 Torr in 50%–50% D-T mixture. The Lee Code simulation gives a detailed description of the plasma focus dynamics through each phase, and provides plasma and yield parameters which show that out of 1.1 × 1019 fast beam ions produced in the plasma focus pinch, only 1.24 × 1014 ions take part in beam-target fusion reactions within the pinch, producing the same number of D-T neutrons. The remnant beam ions, numbering at least 1019, exit the focus pinch at 1.9 MeV, which is far above the 115 keV ion energy necessary for an optimum beam-target cross-section. We propose to regain the lost fusion rates by using a high-pressure D-T-filled drift-tube to attenuate the energy of the remnant beam ions until they reach the energy for the optimum fusion cross-section. Such a fusion enhancement tube would further harvest beam-target fusion reactions by increasing the interaction path length (1 m) at increased interaction density (6 atm). A gain factor of 300 is conservatively estimated, with a final yield of 3.7 × 1016 D-T neutrons carrying kinetic energy of 83.6 kJ, demonstrating Q = 0.01.
A hollow cathode discharge is a unique gas discharge form seen in a cathode cavity. Owing to its high working pressure, high electronic density, and low maintenance voltage [1–4], it was first used as a spectral light source for high-resolution spectral analysis. Later, it was used in metal-vapor ion laser systems, electric propulsion, surface treatment, electrodynamic tether, and other applications [5–9]. For electric propulsion, it is widely used in ion thrusters and Hall-effect thrusters [10–14]. There are usually some discharge instabilities in the discharge process from the discharge chamber of ion thruster, especially under high-current conditions, which significantly affects its industrial application.
When the ion thruster discharge current is a few hundred mA, it enters the arc discharge stage. Domonkos and Williams [15] measured the pulse phenomenon of annular discharge, and found that the pulse of annular discharge is particularly intense under the conditions of high current, low flow, and high voltage.
Subsequently, Goebel [16, 17] reported the presence of high-energy ions in a large current pulse. The potential of these high-energy ions was higher than the voltage applied by the circuit. As the spatial and discharge conditions of the discharge chamber of ion thruster were different, the amplitude and frequency of the pulse were also different. Concurrently, in an experiment designed to calculate the beam current and discharge loss, Goebel found that the ion engine also exhibits discharge instability.
Extensive researches have been conducted to determine the cause of high-current discharge instability, and considerable progresses have been made. For example, Williams and Farnell [18] established a probe system to measure the space potential of the discharge chamber of ion thruster under different anode conditions. They found the presence of a visible discharge pulse, which indicates that the discharge pulse is unrelated to the anode condition. Gabriel [19], through an experimental study, suggested that the discharge instability is mainly related to the filling rate of neutral gas.
Furthermore, Jorns [20] used the average frequency quasi-linear model and experiments to study the acoustic potential pulse oscillations of discharge chamber of ion thruster at a current of 100 A. It was proposed that the ion-acoustic potential pulse oscillations are affected by abnormal electron Landau damping, ion Landau damping, and ion collision frequency, in which the first two terms are more significant. Shunsuke [21] examined the discharge–pulse oscillations by studying the three discharge forms of the hollow cathode (plume, spot, and diffuse). It was identified that the discharge current pulse is affected by the discharge forms of the ion thruster. The effects of the gas flow rate, discharge current, and different types of anode discharge chambers on the ion thruster potential pulse have also been studied [10, 22, 23]. The ion thruster discharge is affected by different types of anode discharge chambers. It was also confirmed that the instability of ionization causes ionization-like high-amplitude potential pulses. Space ions play a decisive role in this kind of space potential pulse. Nelson [24] studied the potential pulse caused by ion-acoustic waves. It was experimentally confirmed that potential pulse oscillation should be universal in a high-current working condition.
The pulse phenomenon of the discharge chamber of ion thruster in the high-current stage is crucial. It exerts a significant influence in practical engineering, but presently, there are no theories that can accurately predict the discharge instability pulse in the arc stage. In this paper, a complete discharge–pulse prediction model is proposed, which combines the equivalent circuit model and the zero-dimensional plasma discharge model, to investigate the discharge–pulse instability in the arc stage. The plasma density and temperature are obtained by using the zero-dimensional model, and then the obtained parameters are brought into the equivalent discharge model to obtain the discharge oscillation curve. The established model is also used to study the effects of different discharge chambers on the pulse characteristics of high-current ion thruster ionization discharge.
The schematic for experimental measurement of the ion thruster pulse is shown in figure 1. A 6 cm Kaufmann ion thruster is used to study, which consists of four parts: the main cathode, the discharge chamber, the optical system and the neutral cathode, and the specific structural parameters are shown in the figure 1. The emission probe is used in this experiment to measure plasma potential. It is a cylindrical tungsten wire with a diameter of 0.127 mm, and the length of the tungsten wire in contact with the plasma is 4 mm. The vacuum chamber has a stainless-steel body with a diameter of 1.8 m and a total length of 4 m. The pumping speed of the diffusion pump is 5000 l s-1, and the ultimate vacuum degree in the chamber is 5×10-4 Pa.
The optical system adopts double grid structure, which is composed of screen grid and acceleration grid. The structural parameters of the discharge chamber and grid are shown in table 1 and table 2.
Magnetic field at the tip of the ring |
1.5 |
Discharge chamber volume |
4.77 |
Discharge chamber length |
1.7 |
Internal surface area of discharge chamber |
4 |
Screen grid area (10-3 m2) | 1.33 |
Acceleration grid area (10-3 m2) | 1.33 |
Grid pitch (10-3 m) | 0.8 |
Grid transparency | 0.75 |
Screen grid plate thickness (10-3 m) | 0.5 |
Acceleration grid plate thickness (10-3 m) | 0.5 |
Screen grid hole radius (10-3 m) | 0.8 |
Acceleration grid hole radius (10-3 m) | 0.8 |
The diagnostic and measuring systems include an emission probe, high-speed camera, current probe (Tektronix TCP0030A), voltage probe (Tektronix TPP1000), the capacitance tester (YTC720, with the range of 2–2000 μF) digital oscilloscope (Tektronix DPO4104B, with a measurement accuracy of ±1%) and automatic capacitance tester. The emission probe is used to measure the plasma potential, a high-speed camera is employed to capture the discharge forms, and no filter is used when capturing the discharge form. The current and voltage probes are used to measure discharge current and discharge voltage, respectively, and the sampling rate is 100 MHz.
Goebel [16, 17] constructed a zero-dimensional plasma discharge model based on the two basic formulae governing neutral gas density and discharge current conservation along with a series of simplifications of ion thruster discharge. The model can be used to solve for plasma density, temperature, primary electron density and neutral gas density.
The plasma density
Ib=12nievaAsTs≈12nee√kTeMAsTs, | (1) |
where
For an unmagnetized plasma, the current
Ii=12nie√kTeMA, | (2) |
where
Ip=n0np⟨σivp⟩V+n0ne⟨σive⟩V, | (3) |
where
The total ion production rate is equal to the total ion absorption rate:
Ip=Ii. | (4) |
Simultaneous equations (1)–(4) can be obtained:
√kTeM⟨σive⟩V+npne⟨σivp⟩V=2n0VA. | (5) |
In the equation (5), 〈σive〉 and 〈σivp〉 are the functions of plasma temperature and density [17]. If
The neutral gas density is calculated as follows:
Qout=Qin-Ibe. | (6) |
In the equation,
Qout=14n0v0AsTsηc, | (7) |
where
ηmd=IbQine. | (8) |
Simultaneous equations (6)–(8) can obtain the expression of neutral gas density as:
n0=4Qin(1-ηmd)v0AsTsηc=4Ibv0eAsTsηc·1-ηmdηmd. | (9) |
The primary electron density in equation (5) can be evaluated from the total primary electron confinement time in the discharge chamber. The emitted current
Ie=npeVτt, | (10) |
where
τp=VvpAp | (11) |
Ap=2rpLc=2B0√2mvpeLc, | (12) |
where
It is assumed that the primary electrons have undergone an inelastic collision with the neutral gas and have lost sufficient energy such that they are then rapidly thermalized with the plasma electrons. The mean time for a collision between the primary and a neutral gas atom to occur is given by:
τc=1n0σvp | (13) |
σ=6.6×10-19(Te4-0.11+Te4)1.6, | (14) |
where
τc=v0eAgTsηcηmd4σvpIb(1-ηmd)=v0AgTsηc4σvpQin(1-ηmd). | (15) |
Finally, primary electrons can also be thermalized by equilibrating with the plasma electrons. The time for primary electrons to slow into a Maxwellian electron population was derived by Spitzer [25] and is given by:
τs=ω2ADl2fG(lfω), | (16) |
where
AD=8πe4nelnΛm2 | (17) |
and
lnΛ=23-0.5ln(10-6neTe3). | (18) |
The function
1τt=1τp+1τc+1τs. | (19) |
The current emitted from the hollow cathode is:
Ie=Id-Is-Ik | (20) |
Is=(1-Ts)2nievaAs=(1-Ts)2nie√kTeMAs, | (21) |
where
np=IeτceV=IeeV[1τp+1τc+1τs]-1=IeeV[vpApV+4σvpQin(1-ηmd)v0AsTsηc+1τs]-1. | (22) |
Assuming that the primary electron loss directly to the anode is negligible, the electron equilibration time is long, and the ion current flowing back to the cathode is small, then the above equation can be written as:
np=Iev0AsTaηc4VσvpIbηmd(1-ηmd)=(Id-Is)v0AsTsηc4VσvpIbηmd(1-ηmd). | (23) |
The plasma density and temperature can be obtained by solving equations (1) and (5).
Aubert et al [26] made significant contributions to the construction of this equivalent-circuit model in the low-current stage. Based on this model, we propose an equivalent-circuit model of the high-current stage. Figure 2 shows an equivalent circuit model of a high-current stage, mainly used to study the ionization-type pulse phenomenon, which is similar to the equivalent model for the low-current stage. The negative electrode of the power supply is at the same potential as ground. The capacitance in the circuit model only depends on the structure of the ion thruster discharge chamber [25], the capacitance tester is used to measure the discharge chamber of ion thruster capacitance, and the measured value is 355 μF.
The basic node current relationship and loop voltage relationship are easy to obtain:
Id=IC+IR | (24) |
V0-VARb=VARp+CdVAdt. | (25) |
According to these two equations, the equivalent circuit model is constructed. The external power supply is represented by a constant voltage source and a current-limiting resistor. The pulse of the discharge current-voltage observed in the experiment is equivalent to the charging and discharging process of the capacitor in the circuit.
Scholars initially used a single sinusoidal function to simulate the discharge voltage pulse. The frequency and amplitude of the sinusoidal function were based on the data obtained in the experiment. Thus, the discharge voltage pulse can be represented as the average of discharge voltage plus this sinusoidal function:
VA=¯VA+A(Te)sin(αt), | (26) |
where
a=2πT, | (27) |
where
The amplitude of the anode voltage is mainly related to the electron temperature, due to its influence on the ionization rate and ionic sound waves. When the electron temperature is high, it promotes the excitation of ionization and the growth of ionic sound waves, leading to the increase of the potential pulse amplitude. When the electron temperature is low, the amplitude of the pulse caused by ionization and ionic sound waves decreases, leading to a decrease in the potential pulse amplitude. In practical experiments, when the discharge is in its plume stage with a low plasma density and high electron temperature, the amplitude of the anode pulse will be relatively large. Conversely, when the discharge form is in its spot stage with a high plasma density and a low electron temperature, the relative pulse amplitude is low. There is an essential relationship between the amplitude of discharge voltage pulse and electron temperature. The results show that there is a linear relationship between the pulse amplitude and the ionization rate [26]. Reference [25] shows that the ionization rate
A(Te)=aTe, | (28) |
where the electron temperature
When the structure of the hollow cathode is fixed, the equivalent charge–discharge resistance of the discharge chamber of ion thruster can be expressed as:
Rp=me(vei+ven)Le2neAs, | (29) |
where,
ven=σn0√8eTeπm | (30) |
vei=2.9×10-12nelnΛTe3/2. | (31) |
The plasma density is related to the gas flow rate, the discharge current, and the discharge voltage working condition parameters through the zero-dimensional plasma discharge model. Therefore, the equivalent resistance value can be obtained from the working condition parameters. Since the plasma density changes with the changes in discharge voltage, the equivalent plasma resistance can be divided into an equivalent plasma charge resistance and an equivalent plasma discharge resistance. In the specific solution process, for the equivalent charging resistance, the discharge current of the input of the model is taken as the value of the maximum. Contrarily, for the equivalent discharge resistance, we take the value of the minimum. For both, the remaining conditions should be consistent. Thus, we can obtain the equivalent charging and discharging resistances, respectively.
These parameters are closely related to microscopic parameters. However, microscopic parameters are coupled with each other, making it impossible to calculate the microscopic parameters corresponding to working conditions directly. The zero-dimensional plasma discharge model adopts two formulae related to conservation of neutral gas density and discharge current. This can be used to solve for neutral and plasma parameters, corresponding to each specific working condition parameter. The corresponding microscopic parameters are then determined. Thus, the function of obtaining microscopic parameters from working condition data is attained.
The equivalent-circuit model can estimate the waveforms of the discharge pulse of the discharge chamber of ion thruster in the high-current stage, and the zero-dimensional plasma discharge model provides a method for calculating the microparameters, such as the plasma density and the temperature, from the working condition parameters. It is, therefore, possible to estimate the discharge–pulse waveforms of the high-current stage, when the input of the equivalent circuit model is obtained from the microscopic parameters. The simulated process is shown in figure 3.
The following set of experimental data is used: the discharge current changes in the range of 5 A, the average of it is 5 A, the beam current is 0.3 A, the supply voltage is 25 V, the value of the current-limiting resistance in the experiment is 1.14 Ω, the gas flow rate of the cathode is 10 sccm, and the discharge supply gas is 2 sccm.
Figure 4 demonstrates the comparison between the experimental measurement waveforms of the discharge voltage pulse and the sinusoidal simulation waveforms. The charging and discharging processes are the same for the sinusoidal simulation, however, this is not aligned with the actual simulation. The increase and decrease in the voltage correspond to the decrease and increase in the internal average potential, respectively. The corresponding physical mechanism is a process of slow accumulation and rapid dissipation of positive ions. Therefore, the times when the anode voltage rises and falls are not equal. However, the simulation of this sinusoidal pulse is highly dependent on the experimental data, and thus, does not have much significance in practical engineering.
The zero-dimensional model and the equivalent circuit model introduced in section 2 are simulated, as a result, the electron temperature is 3 eV, the
In figure 5, the waveforms of the corresponding experimental measurement are compared with those of theoretical simulation. The experimental simulation results fit well with the actual situation, thus addressing the two limitations in the sinusoidal simulation model.
As shown in figure 6, after the application of the fast Fourier transform, the relationship between the discharge current and pulse frequency is obtained at various discharge currents. The five key points are extracted from figure 6 for later simulation, (13 A, 60 kHz), (15 A, 75 kHz), (18 A, 85 kHz), (20 A, 110 kHz), and (25 A, 120 kHz).
When the discharge voltage pulse amplitude, discharge current, current-limiting resistance, and hollow cathode structure parameters are known, the controlled variable method is used by changing only the charging resistance, discharge resistance, and equivalent capacitance. The waveforms are consistent with the frequency measured in the experiment, thereby, revealing the change in the charging resistance, discharging resistance, and frequency of discharge. The simulated results are shown in figures 7–9.
In three sets of equivalent resistance simulation experiments, figures 7 and 8 illustrate the corresponding frequency by adjusting the charging resistance and the discharging resistance, respectively, while figure 9 shows the same frequency by adjusting the capacitance value. The change in charging resistance and discharging resistance is shown in figure 10. The change in the equivalent capacitance is shown in figure 11.
From figures 10 and 11, it can be observed that, with an increase in the discharge current, the plasma equivalent resistance decreases. The plasma resistance changes slowly with discharge current from 13 to 15 A, and 20 to 25 A, and rapidly from 15 to 20 A. However, the frequency of the current pulse decreases with the increase in capacitance, and this change is uniform. The recorded discharge forms are shown in figure 12.
When the discharge current is low, it is in the form of a purple plume. With the further increase in discharge current, the plasma plume grows, reaching a maximum at approximately 15 A. Upon further increasing the discharge current, the plume shrinks to a point near the cathode exit, which is caused by the radial contraction of plasma. When the plasma plume contracts, as seen in figure 12, a bright region is formed between the anode and the cathode. Subsequently, the plasma becomes more compact because the excited ionization region becomes more compact, and when the discharge current increases, the electron temperature decreases rapidly.
The equivalent resistance is inversely proportional to the plasma density. Therefore, when the discharge current increases, the plasma density increases, and thus, the resistance decreases. When the current is between 15 and 18 A, the plasma plume shrinks radially, the plasma density increases rapidly, and thus the resistance changes rapidly. However, the equivalent capacitance depends only on the hollow cathode structure. The change in equivalent capacitance will not change the discharge current and will not cause any change in potential difference. Therefore, changing only the equivalent capacitance will not change the form of the discharge. The capacitance and discharge–pulse frequency will change evenly.
Some experimental waveforms obtained using the experimental system are shown in figure 13. The spatial potential pulse waveforms are synchronized with the oscillation variation of the cathode current waveforms. This indicates that the discharge oscillation is really related to the change in spatial plasma (ion/electron) density.
When other conditions are known, the volume of the discharge chamber is reduced from 47.7 to 10 cm2, the discharge–pulse curve is calculated using the prediction model, and the change of the discharge–pulse curve is shown in figure 14.
When the volume of the discharge chamber decreases from 4770 to 1000 cm2, the discharge voltage pulse frequency increases from 238 to 463 kHz when other conditions remain unchanged. This is due to the increase in the plasma density when the volume of the discharge chamber decreases. Moreover, when other conditions remain unchanged, the decrease in discharge chamber volume will lead to an increase in the electron temperature and thus an increase in the pulse amplitude of discharge anode voltage.
When other conditions are satisfied, the transmittance of the grid plate in the discharge chamber increases from 0.7 to 0.8. The discharge–pulse curve is calculated using the predictive model. The changes in the discharge–pulse curve are illustrated in figure 15.
When the other conditions remain unchanged, and the transmittance of the gate plate increases from 0.75 to 0.85, the pulse frequency of anode voltage decreases from 238 to 206 kHz. This is due to the decrease in the plasma density. Furthermore, the grid plate transmittance has little effect on the electron temperature, thus, the pulse amplitude of the discharge voltage changes slightly.
A theoretical model for predicting the discharge–pulse waveforms is proposed, which aims at addressing the phenomenon of the ionization discharge pulse in the discharge chamber of ion thruster in the arc phase. The effects of plasma density and temperature on the amplitude and frequency of the pulse waveforms are theoretically analyzed. Parameters, such as discharge current, voltage, and gas flow rate, can be used to estimate the amplitude and frequency of the pulse waveforms.
The simulation is performed using the established predictive model, and the waveforms obtained from the simulation are consistent with those obtained from the experiment. It was experimentally found that, when the discharge current is between 15 and 18 A, the plume shrinks radially near the cathode outlet, and the plasma density will increase rapidly in this range. The equivalent resistance in the simulation model showed that the turning point of the resistance change was similar to the current corresponding to the rapid change in plasma density.
The established model is also used to study the effects of different discharge chambers on the pulse characteristics of high-current hollow cathode ionization discharge. When the volume of the discharge chamber reduces, while other parameters are unchanged, the plasma density increases, and the temperature rises. Thus, the discharge–pulse frequency and amplitude tend to increase. When the grid plate transmittance increases and other parameters are constant, the plasma density reduces. Thus, the discharge–pulse frequency tends to reduce.
[1] |
Gribkov V et al 2007 J. Phys. D: Appl. Phys. 40 3592 doi: 10.1088/0022-3727/40/12/008
|
[2] |
Auluck S et al 2021 Plasma 4 450 doi: 10.3390/plasma4030033
|
[3] |
Schmidt A et al 2021 IEEE Trans. Plasma Sci. 49 3299 doi: 10.1109/TPS.2021.3106313
|
[4] |
Lee S 2008 Plasma Phys. Control. Fusion 50 10500
|
[5] |
Lee S and Saw S H 2017 The plasma focus—numerical experiments insights and applications ed R S Rawat Plasma Science and Technology for Emerging Economies Singapore (Berlin: Springer)
|
[6] |
Jackson S L et al 2020 Effect of initial conditions on charged particle acceleration in a dense plasma focus 2020 IEEE International Conference on Plasma Science (ICOPS) (Singapore: IEEE) 376
|
[7] |
Beresnyak A et al 2018 IEEE Trans. Plasma Sci. 46 3881 doi: 10.1109/TPS.2018.2865188
|
[8] |
Decker G et al 1980 Plasma Phys. 22 245 doi: 10.1088/0032-1028/22/3/003
|
[9] |
Struve K W and Freeman B L Dense plasma focus neutron generator for active interrogation SAND2012-0451C (Albuquerque: Sandia National Lab)
|
[10] |
Nukulin V Y and Polukhin S N 2007 Plasma Phys. Rep. 33 271 doi: 10.1134/S1063780X07040022
|
[11] |
Lee S 2009 Appl. Phys. Lett. 95 151503 doi: 10.1063/1.3246159
|
[12] |
Kies W 1986 Plasma Phys. Control. Fusion 28 1645 doi: 10.1088/0741-3335/28/11/003
|
[13] |
Lawson J D 1955 Some criteria for a power producing thermonuclear reactor PhD Thesis Atomic Energy Research Establishment Harwell, Berkshire
|
[14] |
Lee S 2014 J. Fusion Energ. 33 319 doi: 10.1007/s10894-014-9683-8
|
[15] |
Glasstone S and Lovberg R H 1960 Controlled Thermonuclear Reactions (Princeton: Van Nostrand)
|
[16] |
Huba J D 2006 NRL Plasma Formulary (Washington: Naval Research Laboratory)
|
[17] |
Haines M G et al 2006 Phys. Rev. Lett. 96 075003 doi: 10.1103/PhysRevLett.96.075003
|
[18] |
Lee S 1988 Laser Part. Beams 6 597 doi: 10.1017/S026303460000553X
|
[19] |
Lee S 2022 Radiative dense plasma focus computation package: RADPF (www.plasmafocus.net/IPFS/modelpackage/File1RADPF.htm)
|
[20] |
Lee S and Saw S H 2012 Phys. Plasmas 19 2703
|
[21] |
Lee S and Saw S H 2008 Appl. Phys. Lett. 92 021503 doi: 10.1063/1.2827579
|
[22] |
Lee S and Saw S H 2008 J. Fusion Energ. 27 292 doi: 10.1007/s10894-008-9132-7
|
[23] |
Lee S et al 2016 IEEE Trans. Plasma Sci. 44 165 doi: 10.1109/TPS.2015.2497269
|
[24] |
Lee S and Serban A 1996 IEEE Trans. Plasma Sci. 24 1101 doi: 10.1109/27.533118
|
[25] |
Lee S and Saw S H 2010 A course on plasma focus numerical experiments Joint ICTP-IAEA Workshop on Dense Magnetized Plasmas and Plasma Diagnostics (Trieste: ICTP)
|
[26] |
Koh J M et al 2005 Plasma Sources Sci. Technol. 14 12 doi: 10.1088/0963-0252/14/1/002
|
[27] |
Lerner E J, Murali S K and Haboub A 2011 J. Fusion Energ. 30 367 doi: 10.1007/s10894-011-9385-4
|
[28] |
Miller S M et al 2020 Rev. Sci. Instrum. 91 063507 doi: 10.1063/1.5139663
|
[29] |
Stygar W A et al 2015 Phys. Rev. ST Accel. Beams 18 110401 doi: 10.1103/PhysRevSTAB.18.110401
|
[30] |
Lee S 1984 J. Phys. D: Appl. Phys. 17 733 doi: 10.1088/0022-3727/17/4/012
|
[31] |
Lee S and Saw S H 2012 J. Fusion Energ. 31 603 doi: 10.1007/s10894-012-9506-8
|
[32] |
LeChien K et al 2009 A 1-MV, 1-MA, 0.1-Hz linear transformer driver utilizing an internal water transmission line 2009 IEEE Pulsed Power Conference (Washington, DC: IEEE)
|
[33] |
Xu Q et al 2020 Matter Radiat. Extrem. 5 014401 doi: 10.1063/1.5120256
|
[34] |
Qiu H et al 2020 Appl. Sci. 10 8301 doi: 10.3390/app10228301
|
[35] |
Lim L H et al 2015 Phys. Plasmas 22 092702 doi: 10.1063/1.4929856
|
[36] |
Damideh V et al 2019 Vacuum 169 108916 doi: 10.1016/j.vacuum.2019.108916
|
[37] |
Akel M et al 2022 Plasma 5 184 doi: 10.3390/plasma5020014
|
[38] |
Akel M et al 2016 J. Fusion Energ. 35 694 doi: 10.1007/s10894-016-0092-z
|
[39] |
Marciniak Ł et al 2022 Appl. Radiat. Isot. 182 110118 doi: 10.1016/j.apradiso.2022.110118
|
[40] |
Lee S 2022 Plasma focus systems and methods for producing neutrons U.S. Provisional Patent Application No 63/323, 684
|
[41] |
Lee S 2022 Plasma Phys. Control. Fusion 64 082001 doi: 10.1088/1361-6587/ac7b49
|
[42] |
Lee S, Damideh V and Btaiche J C 2022 Plasma focus systems and methods for aneutronic fusion U.S. Provisional Patent Application No 63/323, 697
|
[43] |
Lee S, Damideh V and Btaiche J C 2022 Plasma focus systems and methods with enhanced neutron yield U.S. Provisional Patent Application No. 63/363, 448
|
[1] | Xiaoxi DUAN (段晓溪), Benqiong LIU (刘本琼), Huige ZHANG (张惠鸽), Ben LI (李犇), Jiting OUYANG (欧阳吉庭). Various patterns in dielectric barrier glow discharges simulated by a dynamic model[J]. Plasma Science and Technology, 2019, 21(8): 85401-085401. DOI: 10.1088/2058-6272/ab0d51 |
[2] | Simin ZHOU (周思敏), Xiutao HUANG (黄修涛), Minghai LIU (刘明海). Electrical model and experimental analysis of a double spiral structure surface dielectric barrier discharge[J]. Plasma Science and Technology, 2019, 21(6): 65401-065401. DOI: 10.1088/2058-6272/ab0814 |
[3] | Mohamed MOSTAFAOUI, Djilali BENYOUCEF. Electrical model parameters identification of radiofrequency discharge in argon through 1D3V/PIC-MC model[J]. Plasma Science and Technology, 2018, 20(9): 95401-095401. DOI: 10.1088/2058-6272/aac3cf |
[4] | Chenchen WU (吴辰宸), Xinfeng SUN (孙新锋), Zuo GU (顾左), Yanhui JIA (贾艳辉). Numerical research of a 2D axial symmetry hybrid model for the radio-frequency ion thruster[J]. Plasma Science and Technology, 2018, 20(4): 45502-045502. DOI: 10.1088/2058-6272/aaa8d9 |
[5] | Saravanan ARUMUGAM, Prince ALEX, Suraj Kumar SINHA. Feedback model of secondary electron emission in DC gas discharge plasmas[J]. Plasma Science and Technology, 2018, 20(2): 25404-025404. DOI: 10.1088/2058-6272/aa8e3f |
[6] | Xiang HE (何湘), Chong LIU (刘冲), Yachun ZHANG (张亚春), Jianping CHEN (陈建平), Yudong CHEN (陈玉东), Xiaojun ZENG (曾小军), Bingyan CHEN (陈秉岩), Jiaxin PANG (庞佳鑫), Yibing WANG (王一兵). Diagnostic of capacitively coupled radio frequency plasma from electrical discharge characteristics: comparison with optical emission spectroscopy and fluid model simulation[J]. Plasma Science and Technology, 2018, 20(2): 24005-024005. DOI: 10.1088/2058-6272/aa9a31 |
[7] | ZHANG Junmin (张俊民 ), CHI Chengbin (迟程缤), GUAN Yonggang (关永刚), LIU Weidong (刘卫东), WU Junhui (吴军辉). Simulation of Arc Rotation and Its Effects on Pressure of Expansion Volume in an Auto-Expansion SF6 Circuit Breaker[J]. Plasma Science and Technology, 2016, 18(3): 287-291. DOI: 10.1088/1009-0630/18/3/12 |
[8] | LIN Xin (林莘), WANG Feiming (王飞鸣), XU Jianyuan (徐建源), XIA Yalong (夏亚龙), LIU Weidong (刘卫东). Study on the Mathematical Model of Dielectric Recovery Characteristics in High Voltage SF6 Circuit Breaker[J]. Plasma Science and Technology, 2016, 18(3): 223-229. DOI: 10.1088/1009-0630/18/3/02 |
[9] | GE Lei(葛蕾), ZHANG Yuantao(张远涛). A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges[J]. Plasma Science and Technology, 2014, 16(10): 924-929. DOI: 10.1088/1009-0630/16/10/05 |
[10] | WANG Yan(王燕), LIU Xiang-Mei(刘相梅), SONG Yuan-Hong(宋远红), WANG You-Nian(王友年). e-dimensional fluid model of pulse modulated radio-frequency SiH4/N2/O2 discharge[J]. Plasma Science and Technology, 2012, 14(2): 107-110. DOI: 10.1088/1009-0630/14/2/05 |
1. | Tian, F., Miao, L., Liang, F. et al. Research on Two-dimensional Simulation Model of Discharge Instability of Hollow Cathode Plume Plasma | [空心阴极羽流等离子体放电不稳定性二维仿真模型]. Binggong Xuebao/Acta Armamentarii, 2024, 45(4): 1208-1218. DOI:10.12382/bgxb.2022.0946 | |
2. | Tian, F., Miao, L., Xia, Q. et al. Analysis of the effect of keeper working conditions on hollow cathode performance. Acta Astronautica, 2023. DOI:10.1016/j.actaastro.2023.05.044 |
Magnetic field at the tip of the ring |
1.5 |
Discharge chamber volume |
4.77 |
Discharge chamber length |
1.7 |
Internal surface area of discharge chamber |
4 |
Screen grid area (10-3 m2) | 1.33 |
Acceleration grid area (10-3 m2) | 1.33 |
Grid pitch (10-3 m) | 0.8 |
Grid transparency | 0.75 |
Screen grid plate thickness (10-3 m) | 0.5 |
Acceleration grid plate thickness (10-3 m) | 0.5 |
Screen grid hole radius (10-3 m) | 0.8 |
Acceleration grid hole radius (10-3 m) | 0.8 |