ITPA St Petersburg 20-24 April 2009 G.Gorini 1
JET results on the determination of thermal/non-thermal fusion yield from
neutron emission spectroscopy
Uppsala University
Istituto di Fisica del Plasma “Piero Caldirola”,CNR & Milano-Bicocca University, Milano, Italy
EFDA JET Collaborators
Giuseppe Gorinion behalf of
ITPA St Petersburg 20-24 April 2009 G.Gorini 2
Issue: Thermal and non-thermal neutron yield (Q value)used nowadays for physics analysis of fusion plasmas (e.g. JET). Not measured.Qnth/Qth desirable measurement on ITER.Can Neutron Emission Spectroscopy determinethe thermal and non-thermal neutron yield and their ratio?
Here: JET results from selected D dischargeswith NBI heating.
ITPA St Petersburg 20-24 April 2009 G.Gorini 3
TOFOR neutron spectromer
neutrons
TOFOR is a new 2.5 MeV Time-Of-Flight neutron spectrometer Optimized for high Rate operation (>200 kHz range)
ITPA St Petersburg 20-24 April 2009 G.Gorini 4
TOFOR neutron spectraMain features in the measured TOF spectrum:
• Broad component (e.g. due to NBI deuterons)
• Low energy tail due to neutron scattering
•Higher energy neutrons: missing in the absence of RF heating
#69242, D plasma, NB heating, single injector tof[ns]
ITPA St Petersburg 20-24 April 2009 G.Gorini 5
TOFOR neutron spectra (linear scale)
#69242, D plasma, NB heating, single injector
ITPA St Petersburg 20-24 April 2009 G.Gorini 6
NBI in JET
• Full model in TRANSP
ITPA St Petersburg 20-24 April 2009 G.Gorini 7
TRANSP• Input experimental data
• Simulation of:– Equilibrium– Trasport– Fast ion dynamics– Radiation emission etc (e.g. for diagnostic comparison)
• NUBEAM
ITPA St Petersburg 20-24 April 2009 G.Gorini 8
NUBEAM
• Motion of particles• Neutral deposition• Ionizzation• Orbits• Collisions and
thermalization
ITPA St Petersburg 20-24 April 2009 G.Gorini 9
Control Room
Control Room is a MonteCarlo code first developed in 1994.
The original system consisted on a C++ library and a few applicationswhose purpose was the calculation of slowing-down ion distributions and the resulting neutron emission in thermonuclear plasmas.
Now it is possible to use the library from within the Python programming language. In fact, one could say the library was designed for being used as a Python module. This gives the user a number of advantages. On theone hand, this allows one to test-drive the library through the interactivePython interpreter and makes it easier to experiment with its features. On the other hand, this enables the user to write simple scripts.
ITPA St Petersburg 20-24 April 2009 G.Gorini 10
Example of spectrum
#69625t = 58 s.
ITPA St Petersburg 20-24 April 2009 G.Gorini 11
TOFOR data
Pulse: 69652
PNBI = 13 MWIp = 2.6 MABT = 2.2 T
T= 1.7 keV
Chi2 = 0.98
ITPA St Petersburg 20-24 April 2009 G.Gorini 12
Pulse
ITPA St Petersburg 20-24 April 2009 G.Gorini 13
Missing Neutrons
ITPA St Petersburg 20-24 April 2009 G.Gorini 14
Missing Neutrons
ITPA St Petersburg 20-24 April 2009 G.Gorini 15
Yth/Ynth
ITPA St Petersburg 20-24 April 2009 G.Gorini 16
Role of sight line• To what extent is the plasma volume seen by TOFOR representative of the TOTAL volume?
Equation: y = a*(x-b) R2 = 0.98364 a = 1.216 ± 0.02391b= 0.00269 ± 0.00369
ITPA St Petersburg 20-24 April 2009 G.Gorini 17
Conclusions
TRANSP-based accurate simulation of neutron emission spectra
with thermal and non-thermal components.
Separation of thermal/non-thermal components possible.
Accuracy depends on shape difference.
On ITER shape will be very different => separation easier.