CPOTS * 2nd ERASMUS Intensive Program Introduction to Charged

Report
CPOTS – 2nd ERASMUS Intensive Program
Introduction to Charged Particle
Optics: Theory and Simulation
http://cpots2012.physics.uoc.gr
Dept. of Physics, University of Crete
Aug 19 – Sept 2, 2012
Heraklion, Crete, GREECE
Group Project
Building the International Thermonuclear
Experimental Reactor (ITER)
Ioannis
Jordi
Stefan
Dr. Jason Greenwood
Dept. of Physics and Astronomy
Queen’s University Belfast
E-mail: [email protected]
http://cpots2012.physics.uoc.gr
2
What is ITER?
ITER



Nuclear
fusion research and
engineering project
Expected to produce
ten times the amount
of energy input
Largest Tokamak
reactor
Tokamak


Confine a plasma in
the shape of a torus
Helicoidal field
Fusion


Two or more atomic
nuclei join together to
form a single nucleus
Deuterium, tritium
Parameters of ITER


Central radius R0=6,2m
Inner radius r0=2,0m
Parameters of ITER




18 toroidal coils carrying
up to 80kA giving a
maximum field of 11,8T
(5,3T at outer radius)
Deuterium, tritium
and α-particles
Plasma current 15MA
Plasma temperature 8keV
Building ITER in SIMION
Step 1
 using the solenoid-example
 studying the .lua-file
 Fly´m
Building ITER in SIMION
Step 1
 using the solenoid-example
 studying the .lua-file
 Fly´m
Building ITER in SIMION
Step 2a
 adapting the toroid-example
 studying the .lua-file
 Fly´m
Building ITER in SIMION
Step 2a
 using the toroid-example
 studying the .lua-file
 Fly´m
Building ITER in SIMION
Step 2b
 trying different particle conditions (mass,
charge, angle, energy) and magnetic field
strength
 Trapping is not possible!
Building ITER in SIMION
Step 2b
 trying different particle conditions (mass,
charge, angle, energy) and magnetic field
strength
 Trapping is not possible!
Building ITER in SIMION
Step 2b
 explanation for drift:
 gradient in the
magnetic field
 Non homogeneous-field
 Stronger near
center
Building ITER in SIMION
Step 3a
 adding a current loop in LUA
Building ITER in SIMION
Step 3b
 new parts in the LUA-code:
local poloidal_current = 150000
local field2 = MField.hoop {
current = poloidal_current,
center = MField.vector(0,0,0),
normal = MField.vector(0,0,1),
radius = 6210 }
local field = MField.combined_field{field1,field2}
Building ITER in SIMION
Step 4a
define particles:
Particle
α-Particle
Deuterium
Tritium
Proton
Electron
Symbol
4
2He
1D
3 T
1
2
p
e
Charge Mass Energy
Colour
[e] [amu] [keV]
2+
4
4.000 red
1+
2
100 green
1+
3
100 blue
1+
1
100 pink
1-
5*10-4
100
cyan
Building ITER in SIMION
Step 4b
 Fly´m
α-particle
deuterium
tritium
proton
electron
Results
α-particle
deuterium
tritium
proton
electron
Problems


even small errors in calculation may cause
„untrapping“ (Trajectory quality)
Processing time!
Conclusions



Centered particles are trapped
(theoretically forever) independent from
starting direction!
Particles at the edge aren´t trapped stable
„trapping area“: +0,8m to -0,3m from the
center-radius
Thank you for your
attention!

similar documents