Synchrotron emission in Alcator C-Mod: Spectra at three ... · • Relativistic “runaway”...

1
no synchrotron p confined [15] Synchrotron polarization (poloidal projection). B 0 = 3 T, R 0 = 1.75 m, a = 0.4 m, q 0 = 1, r b = 0.15 m, = 50, θ = 0.1 Synchrotron emission in Alcator C-Mod: Spectra at three magnetic fields, visible camera images, and polarization data* RA Tinguely 1 , RS Granetz 1 , M Hoppe 2 , O Embréus 2 , A Stahl 2 , and T Fülöp 2 1 MIT Plasma Science and Fusion Center, Cambridge, MA, USA 2 Chalmers University of Technology, Gothenburg, Sweden *Supported by USDoE award DE-FC02-99ER54512 Abstract: Alcator C-Mod's high magnetic field allows runaway electron synchrotron emission to be observed in the visible wavelength range. Visible spectrometers were used to measure synchrotron spectra at three magnetic fields: 2.7, 5.4, and 7.8 T. Both a test- particle model [1] and kinetic solver CODE (Collisional Distribution of Electrons) [2,3] explore the energy evolution of the runaway population and the impact of magnetic-field- dependent synchrotron radiation as a power loss mechanism. Additionally, distortion- corrected visible camera images capture the spatial distribution and evolution of synchrotron emission in C-Mod. Initial results show good agreement between experiment and the new synthetic diagnostic SOFT (Synchrotron-detecting Orbit-Following Toolkit) [4]. Finally, a first look at synchrotron polarization data is presented. Introduction: In a plasma, a particle’s collision frequency with the background plasma varies as n/v 3 Given a large enough [7,8] electric field, electrons can “run away” to relativistic speeds In a magnetic field, relativistic REs emit synchrotron radiation mc dp dt =− ( , , ,…) Motivation: Relativistic “runaway” electrons (REs) can have energies of ~10s of MeV Full conversion of plasma current to RE current has been seen in TCV, JET [5,6] 15 MA for ITER! Unconfined REs can seriously damage vacuum vessel RE radiation can provide important information about the evolution of REs in momentum and real space 1160824028 Experimental setup: Visible camera (spatial distribution) Motional Stark Effect (polarization, spatial) Absolutely-calibrated spectrometers (energy) HXR Camera (energy, spatial) Synchrotron Orbit-Following Toolkit (SOFT): Synthetic diagnostic of camera AND spectrometer Implements: Camera viewing geometry Magnetic geometry (toroidally symmetric) Momentum space and radial distribution MSE Polarization data: MSE views tangency radii of r/a ~ -0.4 – 0.8 Polarization angle and fraction of polarized light can give information about spatial distribution of runaways Soon to be implemented into SOFT = Summary: RE synchrotron spectra were measured in plasmas with similar driving electric fields (E/Ec) on axis, but three different magnetic fields: 2.7, 5.4, and 7.8 T . From absolute brightness calibrations, the normalized brightness was computed to compare the relative change in spectral shape over time. A test particle model [1,8,11] and kinetic solver CODE [2,3] were used to simulate RE evolution in momentum space, predicting RE energy, pitch, and density. Using appropriate brightness formulas [8,11,12] and synthetic diagnostic SOFT [4], the normalized brightness was calculated and compared to experiment. For B = 2.7 T , the TPM (also with SOFT) reproduces experimental data, indicating that ~35-40 MeV REs are plausible near the plasma core. For B = 5.4 T , the experimental data are best reproduced by the TPM if REs at the q = 3/2 and 2 surfaces contribute significantly to the spectra. For B = 7.8 T , the predicted normalized brightness from the TPM does not match experiment, but the trend in data is captured. First results using the SOFT synthetic camera are able to match experimental synchrotron images well. More work is to be done implementing spatial distributions and spatial evolution of REs. MSE polarization data shows qualitative similarities with theory, but more work is needed for quantitative comparisons. This could provide important spatial information about REs. SOFT CODE Synchrotron emission from REs measured at three magnetic fields: 2.7 T, 5.4 T, 7.8 T Similar RE parameters (E 0 /E c ) [8] on axis For a fixed RE energy and pitch, increasing the magnetic field will shift the synchrotron spectra toward shorter wavelengths B = 2.7 T REs strike limiter in C-Mod + The difficulties with synchrotron spectral measurements: RE density affects amplitude and is hard to reproduce No unique solution, i.e. many combinations of energy and pitch produce same brightness spectra Normalize at one wavelength (675 nm) to observe relative change of spectral shape magnetic axis q = 1, 3/2, 2 Methodology: For each flux surface (magnetic axis; q = 1, 3/2, 2) 1. Obtain flux-surface-averaged plasma parameters (n, T, E, B) 2. Run kinetic solver like CODE [2,3] OR 3. Initiate a RE at many times within the discharge 4. Evolve the test particle model (TPM) for momentum [1] and density [8,11] 5. Check for particle confinement (orbit loss) [12] 6. Determine the maximum allowable RE momentum (expected to dominate synchrotron spectra) 7. Compute the normalized synchrotron brightness from [10,13,14] and using SOFT [4], and compare to experiment B = 5.4 T B = 7.8 T The normalized brightness as calculated from the test particle model using [10,13,14] and normalized power from SOFT are in good agreement with experiment. As the TPM predicts, this suggests that most high energy REs are generated in the plasma core, near the magnetic axis and q = 1 surface. There is little contribution from REs at the outer edge of the plasma. Data is consistent with RE energies of ~35-40 MeV. The normalized brightness predicted from the CODE momentum-space distribution does not match as well with experiment. The reason for this is unclear. The normalized brightness and power calculated from the TPM using [10,13,14] and SOFT match experiment best for the outer q surfaces of q = 2 and 3/2, respectively. REs off axis are contributing significantly to the observed spectra. Data is consistent with RE energies of ~20-25 MeV. If only REs on axis were considered, the predicted spectra would only match the measured spectra for Z eff = 1. Neither the normalized brightness calculated from [10,13,14] nor that from SOFT can reproduce the relative change of spectral shape observed in experiment, even for the highest energy REs at the q = 3/2 and q = 2 surfaces. This could imply higher energy REs than predicted. In fact, the q = 2 surface is just barely observed by the spectrometer. A slight shift in radial position (perhaps due to drifts) could improve the match with experiment. However, the general trend of spectra shifting toward longer wavelengths in time is consistent with experiment. [10] [10] [9] References: [1] JR Martín-Solís, et al., PoP 5 (1998) [2] M Landreman, et al., CPC 185 (2014) [3] A Stahl, et al., NF 56 (2016) [4] M Hoppe, et al., Submitted to NF (2017) [5] G Papp, et al., EX/9-4, 26 th IAEA (2016) [6] C Reux, et al., 30 th MHD ITPA (2017) [7] H Dreicer, Phys Rev 115 (1959) [8] J Connor and R Hastie, NF 15 (1975) [9] J Decker, et al., PPCF (2014) [10] IM Pankratov, Plasma Phys. Reports 25 (1999) [11] M Rosenbluth and S Putvinski, NF 37 (1997) [12] H Knoepfel and DA Spong, NF 19 (1979) [13] JH Yu, et al., PoP 20 (2013) [14] A Stahl, et al., PoP 20 (2013) [15] YM Sobolev, ISSN 1562-6016, BAHT (2013) 0 2.5 5.0 7.5 10 E (MeV) ~ p/2 0 1 2 t (s) B = 2.7 T 2.7 T 5.4 T 7.8 T

Transcript of Synchrotron emission in Alcator C-Mod: Spectra at three ... · • Relativistic “runaway”...

Page 1: Synchrotron emission in Alcator C-Mod: Spectra at three ... · • Relativistic “runaway” electrons (REs) can have energies of ~10s of MeV • Full conversion of plasma current

no synchrotronpconfined

[15]

Synchrotron polarization (poloidal projection).B0 = 3 T, R0 = 1.75 m, a = 0.4 m, q0 = 1,

rb = 0.15 m, 𝛾 = 50, θ = 0.1

Synchrotron emission in Alcator C-Mod: Spectra at three magnetic fields, visible camera images, and polarization data*

RA Tinguely1, RS Granetz1, M Hoppe2, O Embréus2, A Stahl2, and T Fülöp2

1MIT Plasma Science and Fusion Center, Cambridge, MA, USA 2Chalmers University of Technology, Gothenburg, Sweden*Supported by USDoE award DE-FC02-99ER54512

Abstract: Alcator C-Mod's high magnetic field allows runaway electron synchrotron

emission to be observed in the visible wavelength range. Visible spectrometers were usedto measure synchrotron spectra at three magnetic fields: 2.7, 5.4, and 7.8 T. Both a test-particle model [1] and kinetic solver CODE (Collisional Distribution of Electrons) [2,3]explore the energy evolution of the runaway population and the impact of magnetic-field-dependent synchrotron radiation as a power loss mechanism. Additionally, distortion-corrected visible camera images capture the spatial distribution and evolution ofsynchrotron emission in C-Mod. Initial results show good agreement between experimentand the new synthetic diagnostic SOFT (Synchrotron-detecting Orbit-Following Toolkit) [4].Finally, a first look at synchrotron polarization data is presented.

Introduction:• In a plasma, a particle’s collision frequency with the background plasma varies as n/v3

• Given a large enough [7,8] electric field, electrons can “run away” to relativistic speeds• In a magnetic field, relativistic REs emit synchrotron radiation

mcdp

dt= −𝐞𝐄 − 𝐅𝐜

𝐧

𝐯𝟐− 𝐅𝐫𝐚𝐝(𝐩, 𝐁, 𝐙𝐞𝐟𝐟, … )

Motivation:• Relativistic “runaway” electrons (REs) can have

energies of ~10s of MeV• Full conversion of plasma current to RE current has

been seen in TCV, JET [5,6] 15 MA for ITER!• Unconfined REs can seriously damage vacuum vessel

RE radiation can provide important information about the evolution of REs in momentum and real space

1160824028

Experimental setup:

Visible camera(spatial distribution)

Motional Stark Effect(polarization, spatial)

Absolutely-calibrated spectrometers (energy)

HXR Camera(energy, spatial)

Synchrotron Orbit-Following Toolkit (SOFT):• Synthetic diagnostic of camera AND spectrometer

Implements:• Camera viewing geometry• Magnetic geometry (toroidally symmetric)• Momentum space and radial distribution

MSE Polarization data:• MSE views tangency radii of r/a ~ -0.4 – 0.8• Polarization angle and fraction of polarized light

can give information about spatial distribution of runaways

• Soon to be implemented into SOFT

=

Summary:• RE synchrotron spectra were measured in plasmas with similar driving electric fields

(E/Ec) on axis, but three different magnetic fields: 2.7, 5.4, and 7.8 T.• From absolute brightness calibrations, the normalized brightness was computed to

compare the relative change in spectral shape over time.• A test particle model [1,8,11] and kinetic solver CODE [2,3] were used to simulate RE

evolution in momentum space, predicting RE energy, pitch, and density.• Using appropriate brightness formulas [8,11,12] and synthetic diagnostic SOFT [4], the

normalized brightness was calculated and compared to experiment. For B = 2.7 T, the TPM (also with SOFT) reproduces experimental data,

indicating that ~35-40 MeV REs are plausible near the plasma core. For B = 5.4 T, the experimental data are best reproduced by the TPM if

REs at the q = 3/2 and 2 surfaces contribute significantly to the spectra. For B = 7.8 T, the predicted normalized brightness from the TPM does

not match experiment, but the trend in data is captured.• First results using the SOFT synthetic camera are able to match experimental

synchrotron images well. More work is to be done implementing spatial distributions and spatial evolution of REs.

• MSE polarization data shows qualitative similarities with theory, but more work is needed for quantitative comparisons. This could provide important spatial information about REs.

SOFTCODE

Synchrotron emission from REs measured at three magnetic fields: 2.7 T, 5.4 T, 7.8 T• Similar RE parameters (E0/Ec) [8] on axis• For a fixed RE energy and pitch, increasing the

magnetic field will shift the synchrotron spectra toward shorter wavelengths

B = 2.7 T

REs strike limiter in C-Mod

+

The difficulties with synchrotron spectral measurements:• RE density affects amplitude and is hard to reproduce• No unique solution, i.e. many combinations of energy and pitch produce same

brightness spectraNormalize at one wavelength (675 nm) to observe relative change of spectral shape

magnetic axisq = 1, 3/2, 2

Methodology:For each flux surface (magnetic axis; q = 1, 3/2, 2)1. Obtain flux-surface-averaged plasma parameters (n, T, E, B)2. Run kinetic solver like CODE [2,3] OR3. Initiate a RE at many times within the discharge4. Evolve the test particle model (TPM) for momentum [1]

and density [8,11]5. Check for particle confinement (orbit loss) [12]6. Determine the maximum allowable RE momentum

(expected to dominate synchrotron spectra) 7. Compute the normalized synchrotron brightness from

[10,13,14] and using SOFT [4], and compare to experiment

B = 5.4 T

B = 7.8 T

The normalized brightness as calculated from the test particle model using [10,13,14] and normalized power from SOFT are in good agreement with experiment. As the TPM predicts, this suggests that most high energy REs are generated in the

plasma core, near the magnetic axis and q = 1 surface. There is little contribution from REs at the outer edge of the plasma. Data is consistent with RE energies of ~35-40 MeV.

The normalized brightness predicted from the CODE momentum-space distribution does not match as well with experiment. The reason for this is unclear.

The normalized brightness and power calculated from the TPM using [10,13,14] and SOFT match experiment best for the outer q surfaces of q = 2 and 3/2, respectively. REs off axis are contributing significantly to the

observed spectra. Data is consistent with RE energies of ~20-25 MeV. If only REs on axis were considered, the predicted

spectra would only match the measured spectra for Zeff = 1.

Neither the normalized brightness calculated from [10,13,14] nor that from SOFT can reproduce the relative change of spectral shape observed in experiment, even for the highest energy REs at the q = 3/2 and q = 2 surfaces. This could imply higher energy REs than predicted.

In fact, the q = 2 surface is just barely observed by the spectrometer. A slight shift in radial position (perhaps due to drifts) could improve the match with experiment.

However, the general trend of spectra shifting toward longer wavelengths in time is consistent with experiment.

[10]

[10]

[9]

References:[1] JR Martín-Solís, et al., PoP 5 (1998)[2] M Landreman, et al., CPC 185 (2014)[3] A Stahl, et al., NF 56 (2016)[4] M Hoppe, et al., Submitted to NF (2017)[5] G Papp, et al., EX/9-4, 26th IAEA (2016)[6] C Reux, et al., 30th MHD ITPA (2017)[7] H Dreicer, Phys Rev 115 (1959)

[8] J Connor and R Hastie, NF 15 (1975)[9] J Decker, et al., PPCF (2014)[10] IM Pankratov, Plasma Phys. Reports 25 (1999)[11] M Rosenbluth and S Putvinski, NF 37 (1997)[12] H Knoepfel and DA Spong, NF 19 (1979)[13] JH Yu, et al., PoP 20 (2013)[14] A Stahl, et al., PoP 20 (2013)[15] YM Sobolev, ISSN 1562-6016, BAHT (2013)

0 2.5 5.0 7.5 10

E (MeV) ~ p/2

0

1

2

t (s)

B = 2.7 T

2.7 T 5.4 T 7.8 T