CPOTS * 2nd ERASMUS Intensive Program Introduction to Charged
Download
Report
Transcript CPOTS * 2nd ERASMUS Intensive Program Introduction to Charged
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!