1 / 15

Study of Plasma Meniscus and Beam Halo in Negative Ion Sources Using 3D3VPIC model

Study of Plasma Meniscus and Beam Halo in Negative Ion Sources Using 3D3VPIC model. Shu Nishioka Faculty of Science and Technology, Keio Univ. 1 st year Master’s degree. B. Simulation model. Geometry of simulation Model. ・ real size of simulation domain:. ・ Scale factor:.

vanna
Télécharger la présentation

Study of Plasma Meniscus and Beam Halo in Negative Ion Sources Using 3D3VPIC model

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Study of Plasma Meniscus and Beam Halo in Negative Ion Sources Using 3D3VPIC model Shu Nishioka Faculty of Science and Technology, Keio Univ. 1st year Master’s degree

  2. B Simulation model Geometry of simulation Model ・real size of simulation domain: ・Scale factor: ・Number of meshes: ・mesh size: ・Initial number of superparticles electron: Ne = 9×106H+: NH+ = 1×107H-: NH- = 1×106 *PG magnetic filter is taken into account and it is parallel to the y-axis. It was based on JT-60U negative ion source.

  3. Simulation model PIC calculation cycle Integration of the motions of the charged particles Solve the Poisson equation First-order particle weighting charged particle

  4. Simulation model Boundary Condition ・ The periodic boundary conditions are used in the y and z directions. ・ The surface produced H- taken to be a parameter and launched at the PG surface by a given number of particles per time step.

  5. Physical quantities Symbol Normalization Physical parameters Symbol Value Time Electron temperature 1 eV First 2D space coordinate Velocity Hydrogen ion temperature 0.25 eV Electrostatic potential Electron density 1018m-3 Magnetic flux density Electron Debye length 7.43×10-6 m Current density Electron plasma frequency 5.64×1010 rad/s density Electron thermal velocity 4.19×105 m/s Simulation model Computatinal resouce CPU : Intel - Core i7-3820 @3.60GHz (4core 2thread): 8 MPI process. GPU : GTX titan : GPU is used only to solve Poisson eq. RAM : 16GB 13s per PIC cycle.

  6. Result

  7. Comparison between 2D3VPIC and 3D3VPIC ~ Magnetic filter field Emittance diagram (at x=17mm) Electron density H- density Magnetic filter field Plasma meniscus θ Fig.1 2D density profile of the 2D3V PIC model Emittance diagram (at x=17mm, z=0mm) H- density Electron density Plasma meniscus Fig. 2 2D density profile of the 3D3V PIC model ・Plasma Meniscus in the 2D model penetrates more deeply into the plasma source region and curvature is larger. ・The beam halo fraction to the total beam current is estimated to be 51.5% in the 2D model while around 6.3% in the 3D model.

  8. Why is Penetration of the Meniscus small in the 3D model Potential Profile during vacuum condition Schematic view of 2Dmodel geometry 2D model 3D model 3Dmodel geometry ・The figure shows potential profile each 2D, 3D model during vacuum condition. In this 2D model, because extraction hole is modeled as a slit, the equipotential surface penetrate more deeply into the plasma source region. Therefore plasma meniscus penetrates deeply than 3D model.

  9. Result ~perpendicular and parallel to PG filter field plane~ Magnetic filter field Emittance diagram (at x=17mm, z=0mm) H- density Electron density Plasma meniscus Fig. 2 2D density profile of the 3D3V PIC model (Z=0) Electron density H- density Emittance diagram (at x=17mm, z=0mm) Plasma meniscus Magnetic filter field Fig. 2 2D density profile of the 3D3V PIC model(Y=0) ・Asymmetry of the electron density profile due to the E×B drift is observed. ・Asymmetry of the plasma meniscus is observed. It induce asymmetry of negative ion current density profile.

  10. Conclusion • The H- beam halo ratio to extraction beam current is dependent on the penetration of plasma meniscus. • 2. The ratio of beam halo is about 6% in the 3D3V-PIC model. This value reasonably agree with the experimental result1,2. • 3. Asymmetricity of the electrons and negative ions due to the E×B drift is observed with 3D3V-PIC model. 1. K. Miyamoto, Y. Fujiwara, T. Inoue, N. Miyamoto, A. Nagase, Y. Ohara, Y. Okumura, and K. Watanabe, AIP conference proceedings, 380, 360 (1996) 2. H.P.L. de Esch, L. Svensson, Fusion Engineering and Design 86, 363 (2011).

  11. Physical parameters Symbol Value Electron temperature Hydrogen ion temperature Electron density Electron Debye length Electron plasma frequency Electron thermal velocity Benchmark problem Geometry 19mm×17mm×17mm Include One PG aperture PG: hole Radius =14mm, thickness=2mm Peak value of the PG filter filed: 50Gauss Physical parameter

  12. Effect of E×B drift Direction and magnitude of the E×B drift Electron density and direction of E×B Magnetic filter field ・The electron transport depends on along the direction of E×B drift close to PG.

  13. ~ J t (Am-2) Temporal Evolution of Extracted current Negative Ion Electron 1750 1500 1250 1000 750 500 250 0 SP start

  14. Definition of Plasma Meniscus The contour surface with defines Plasma Meniscus. This is almost same as the contour along which grad φ=0. The H+ ions cannot go forward beyond this contour and they return back towards the left hand side into the source region. On the other hand, H- ions is accelerated when beyond this contour, then they are extracted.

More Related