Simulations and laboratory tests are used to design and optimize a quasi-optical system for cross-polarization scattering (CPS) measurements of magnetic turbulence on the DIII-D tokamak. The CPS technique uses a process where magnetic turbulence scatters electromagnetic radiation into the perpendicular polarization enabling a local measurement of the perturbing magnetic fluctuations. This is a challenging measurement that addresses the contribution of magnetic turbulence to anomalous thermal transport in fusion research relevant plasmas. The goal of the new quasi-optical design is to demonstrate the full spatial and wavenumber capabilities of the CPS diagnostic. The approach used consists of independently controlled and in vacuo aiming systems for the probe and scattered beams (55-75 GHz).

This paper reports on a new design quasi-optical system for cross-polarization scattering (CPS)1–12 measurements of magnetic turbulence on the DIII-D tokamak. Previous studies1–12 describe in detail the CPS theory and experimental method, experimental measurements, validity tests, etc., while this paper focusses on a new optimized system. The system presented here employs equipment and techniques not heretofore used in these types of measurements (e.g., internal and independent poloidal and toroidal steering). These improvements provide a new ability to: selectively detect CPS signals from all parts of the probe beam, access an expanded range of detected wavenumbers, and access higher wavenumbers. The design is based upon plasma and quasi-optical simulations along with laboratory tests. The goal of this new design is to demonstrate the full spatial and wavenumber capabilities of the CPS diagnostic.

The CPS technique is based upon the scattering of an incident microwave beam into an orthogonal polarization by magnetic fluctuations B̃. The scattered cross-polarized signal provides a measure of the internal magnetic fluctuations that induce the scattering. In the CPS process, the incident electric field Ei accelerates plasma electrons that interact with the magnetic perturbations B̃. The resulting induced current J̃ is perpendicular to Ei, J̃Ei×B̃.3 The induced current radiates an electric field Es with a polarization orthogonal to both Ei and B̃. It is this scattered field Es, proportional to B̃, that is detected. The CPS method was first used on Tore-Supra3,6 and was followed by experiments on GAMMA-10,7,8 FT-1,9 MAST,11 and DIII-D.10,12 The new design presented here is a significant improvement over that in Refs. 10 and 12 having larger radial and wavenumber ranges as well as significant aiming control. Multiple tests demonstrating the validity of the CPS data as well as the first data from the original CPS system can be found in Refs. 10 and 12.

In the initial design of the CPS technique on DIII-D,10 a Doppler backscattering (DBS) probe beam was chosen for its ability to improve the spatial localization and to provide simultaneous density turbulence and flow data. First introduced on TUMAN-3M13 and W7-AS,14 DBS is now a widely used technique that measures intermediate wavenumber (typically trapped electron modes to lower-k electron temperature gradient modes), density turbulence ñ, and flows.13–21 The DBS method consists of injecting a probe beam of a given frequency and polarization into a plasma which contains a cutoff for that frequency and polarization. The radiation is injected at an angle with respect to the cutoff layer (see Fig. 1) and scatters from internal density fluctuations. The beam propagates toward the cutoff and as the cutoff is approached refraction bends the beam upwards. Near the cutoff layer, a long-wavelength electric field pattern is formed roughly parallel to the flux surface (e.g., see Ref. 14). It is this field pattern that interacts strongly with the density fluctuations propagating in the poloidal direction, i.e., the ñ fluctuations just in front of the cutoff layer.

FIG. 1.

DBS probe launch (red beam) and scattered CPS receive (blue beam) geometry along with the wavevector matching conditions. The launch is assumed X-mode and the cross-polarized signal O-mode.

FIG. 1.

DBS probe launch (red beam) and scattered CPS receive (blue beam) geometry along with the wavevector matching conditions. The launch is assumed X-mode and the cross-polarized signal O-mode.

Close modal

CPS and DBS are scattering techniques and as such have well known physical properties and constraints.22 The most familiar constraint arises from conservation of momentum and energy. These conservation laws lead to a Bragg relation for the scattering process,2 viz.,

(1)
(2)

where ki and ks are the incident and scattered wavevectors, respectively, and kB̃ is the wavevector of the magnetic fluctuation (similarly for the frequencies ωi, ωs, and ωB̃). It is worth noting that, in the CPS process, ki and ks are of orthogonal polarizations (e.g., X-mode and O-mode) and, due to the different indices of refraction, can have very different wavelengths and propagation properties (e.g., see the vector diagram inset in Fig. 1). Since the Bragg relation is a vector relation, scattering only occurs if the kB̃ wavevector is both present (i.e., there needs to be turbulence with that wavevector) and aligned such that it scatters the incident wavevector into the desired scattered wavevector (Fig. 1 illustrates this desired matching condition). Note that the desired scattered wavevector ks should be in the physical direction and within the antenna pattern of the receiver otherwise little or no signal will be detected. The scattering process and the effect of this wavevector alignment have been examined by multiple authors11,22 and can be mathematically expressed as

(3)

where Ei and Es are the incident and scattered electric fields, aθ, ar are the incident beam electric field half-widths (the beams are assumed to be Gaussian shapes in the directions perpendicular to the axis of propagation), and the directions [θ, r] are perpendicular to the local magnetic field but within the flux surface and perpendicular to both the local magnetic field and the local flux surface, respectively. In this example, the wavenumbers parallel to the magnetic field kB,||, ki,||, and ks,|| are assumed to be negligibly small. From Eq. (3), it can be seen that if the differences in wavenumbers are not small (i.e., if ki,r + kB,rks,r and ki,θ + kB,θks,θ are not separately near zero), the resulting scattered field Es will be substantially reduced, possibly to effectively zero.

In Fig. 1, the DBS probe beam is X-mode and the cross-polarized scattered signal is O-mode. In this illustration, the CPS O-mode cutoff location is assumed to be far to the left in the figure, i.e., well away from the illustrated X-mode cutoff location. The X-mode probe beam is refracted upwards by the plasma, whereas the O-mode scattered beam has little refraction due to the assumed differences in indices of refraction. The wavenumber matching/scattering diagram illustrates the incident X-mode wavevector ki, the scattered O-mode wavevector ks, and the magnetic fluctuation scattering wavevector kB̃. In this illustration, the O-mode wavevector is much longer than the X-mode consistent with the larger O-mode index of refraction at this location, a common occurrence in lower density DIII-D plasmas.

In the following, Sec. II presents the scientific and engineering requirements, followed by 3-D raytracing calculations simulating cross-polarization scattering from an X-mode probe beam (DBS) to an O-mode scattered beam (CPS) in Sec. III. These simulations were run for a variety of plasmas to optimize the launch–receiver geometries. This is followed by Gaussian beam calculations using realistic optical elements in Sec. IV. These calculations were evaluated for a variety of focal lengths, distances, etc. Section V details how these elements come together into an integrated quasi-optical design that meets the scientific and engineering requirements.

The scientific requirements for the system were (1) poloidal and toroidal angular access of order ±20° and ±10°, respectively; this provides for DBS access to higher density turbulence wavevectors kθ (i.e., higher poloidal angle) as well as to match the toroidal pitch angle of these wavevectors which is known to be important at higher k;11,15 (2) independent polarization control of launch and receive signals, an important requirement as both CPS and DBS are polarization dependent; (3) independently aimed DBS and CPS antennas; this allows the CPS system to look at arbitrary locations along the probe beam, thus adjusting both the wavenumber probed as well as the spatial location; and (4) a controlled and optimized beam size within the plasma; the beam size at the reflection position is known to affect the wavenumber and spatial resolutions.23 Engineering constraints included (1) component setback from the plasma to avoid plasma heating and/or plasma contamination, and possible interaction with fast ions; (2) a reliable method of in vacuo mirror control; and (3) protection from ≥6 MW of 110 and 117 GHz electron cyclotron heating.

The design was arrived at by first examining the angular (poloidal and toroidal) requirements for both CPS and DBS. The radial, vertical, and antenna separations were varied to determine the optimal locations. The 3-D raytracing code GENRAY24 was utilized. The inputs to this code were measured density profiles (from Thomson scattering and profile reflectometry), and magnetic equilibria (from the EFIT code25) with toroidal symmetry are assumed. A variety of plasmas were examined for this optimization including L-mode, H-mode, and QH-modes, and the launch and receive locations were varied for each case. An L-mode discharge from this set of plasmas is displayed for illustration in Figs. 2 and 3. The magnetic equilibrium and cutoffs and resonances and the example of the 3-D raytracing results are shown in Figs. 2 and 3, respectively. The eight DBS cutoff locations are shown in Fig. 2 from which the 70 GHz frequency is selected to display in Fig. 3. In Fig. 3, the red curves are the launched X-mode DBS probe beam (70 GHz) and the blue curves are the scattered O-mode CPS rays. The solid curves are within the plasma, while the dashed lines are outside the last closed flux surface, where GENRAY uses an index of refraction N = 1. The antennas are indicated and the 70 GHz right-hand cutoff location within the plasma is shown. The DBS probe ray is launched upwards [Fig. 3(a)] and is refracted by the plasma, while the scattered O-mode ray is relatively unaffected by the plasma due to its larger index of refraction (i.e., the O-mode cutoff is much interior to that of the X-mode as can be seen in Fig. 2). Figure 3 illustrates two important aspects of this study. The first is the range of poloidal angles needed by the CPS antenna in order to sample the DBS probe beam. Table I shows how the detected magnetic wavenumber and scattering position vary with the receiver angle. Using these simulations, it was determined that the vertically displaced CPS/DBS antennas produced a superior wavenumber range as well as greater spatial control. Note that the toroidal propagation of the probe beam is shown in Fig. 3(b). This is partially due to launching the probe beam with a toroidal angle so as to maximize the scattering efficiency [i.e., maximize relation (3)] but also partly due to the natural tendency of a ray to travel toroidally (with the direction depending upon the ray polarization, O- or X-mode).

FIG. 2.

Example of cutoff and resonance profiles with the inset showing plasma shape. Filled disks indicate the cutoff locations of the eight DBS probe frequencies.

FIG. 2.

Example of cutoff and resonance profiles with the inset showing plasma shape. Filled disks indicate the cutoff locations of the eight DBS probe frequencies.

Close modal
FIG. 3.

Two views of DBS probe ray (red) launched from bottom and seven CPS scattering locations and rays (blue). Flux surfaces are indicated by solid black lines, with only midplane LCFS shown in (b).

FIG. 3.

Two views of DBS probe ray (red) launched from bottom and seven CPS scattering locations and rays (blue). Flux surfaces are indicated by solid black lines, with only midplane LCFS shown in (b).

Close modal
TABLE I.

Showing wavenumbers of the magnetic fluctuations detected at the seven locations shown in Fig. 3. ρ is the radial location in normalized flux surface coordinates, and wavenumber units are in cm−1.

ρksks,rks,⊥ks,θΔk = ks,⊥ks,θexp(Δk2a2/4)
0.85 22.4 22 22.3 −0.4 22.7 
0.75 22.5 22.3 22.5 21.5 
0.66 19.8 19.6 19.7 2.4 17.3 
0.61 12.7 12.3 12.6 4.1 8.4 
0.67 7.6 5.9 7.4 5.2 2.3 
0.75 7.1 4.1 5.9 1.1 0.06 
0.85 7.8 7.6 6.6 0.11 
ρksks,rks,⊥ks,θΔk = ks,⊥ks,θexp(Δk2a2/4)
0.85 22.4 22 22.3 −0.4 22.7 
0.75 22.5 22.3 22.5 21.5 
0.66 19.8 19.6 19.7 2.4 17.3 
0.61 12.7 12.3 12.6 4.1 8.4 
0.67 7.6 5.9 7.4 5.2 2.3 
0.75 7.1 4.1 5.9 1.1 0.06 
0.85 7.8 7.6 6.6 0.11 

The requirements to be met by the quasi-optical design were that it should (1) form a Gaussian profile beam of the desired width within the plasma over the range of frequencies used (i.e., 55–75 GHz), (2) be simple with as few moving parts as possible, (3) propagate arbitrary polarizations without significant polarization scrambling, and (4) fit into the available space. Based upon the previous experience at DIII-D as well as the published literature,17 a beam radius of ∼3 cm (1/e radius of magnitude electric field) was a reasonable starting point. While multiple focusing elements are possible, single element systems are less prone to alignment or other failure and generally have lower insertion losses. A metal lens26 was selected for its frequency dependent index of refraction which focusses higher frequencies further from the lens and thus deeper into the plasma [other options include dielectric lenses (e.g., quartz), focusing mirrors, either parabolic or ellipsoidal]. This better matches the focus to the reflection point being deeper in the plasma for higher frequencies. Other important features of the metal lens are its resistance to thermal damage, low differential expansion with respect to its metal mounting frame, and being unaffected by thin coatings of the various metals and compounds that internal components often suffer. Figure 4 shows a design drawing of a large metal lens. Although the shape is concave, the lens is a converging lens due to the index of refraction in the small circular metallic waveguides (i.e., the holes). The index of refraction for the metallic lens varies as n=[1(λ/λcutoff)2]1/2, where λ is the vacuum wavelength of the selected frequency, λcutoff = 1.706D, and D is the diameter of the holes in the metal lens. For this design, the hole diameter is chosen to be 4.41 mm, resulting in a cut-off frequency of 39.946 GHz. The radius of curvature of the lens is 17.89 cm. This results in focal lengths of 57.25 cm and 116.5 cm for 55 and 75 GHz frequencies, respectively.

FIG. 4.

Design drawing of a focusing metal plate lens. The hole diameter is 0.441 cm, the radius of curvature is 17.89 cm, and the hole coverage diameter is 9 cm.

FIG. 4.

Design drawing of a focusing metal plate lens. The hole diameter is 0.441 cm, the radius of curvature is 17.89 cm, and the hole coverage diameter is 9 cm.

Close modal

A Gaussian beam radius code (based upon a standard formulation27) was used to calculate the beam widths at various probe frequencies and distances from the lens. This was used to optimize the beam sizes for the frequency range 55-75 GHz and the distances of interest within the plasma. Note that this code does not take into account refraction due to the plasma. This will generally result in modest increases in the beam size with distance.28 This code was used to optimize the beam size in the plasma region of interest by varying the focal length and lens to antenna distance. The iteration of this process for selected frequencies covering the frequency range of interest allowed a selection of the best overall lens focal length and lens to antenna distance. Figure 5 shows the result of these optimizations for three representative frequencies in the band of interest (50, 60, and 70 GHz). The metallic lens is located at the circular waveguide aperture (see Sec. V for a presentation of the resulting design). The frequency dependence of the index of refraction results in the different focal lengths for different frequencies, as shown in Fig. 5. The beam radii minima are at ∼20, 27, and 35 cm from the lens surface (for 50, 60, and 70 GHz, respectively). The last closed flux surface (the approximate edge of the plasma) and normalized radial locations ρ = 0.5 and ρ = 0 are indicated. As can be seen, the radial beam widths are of order 3 cm for the radial range extending from the edge to ρ = 0.5. The higher frequency widths (≥70 GHz) remain near 3 cm interior to ρ = 0.5. Also shown in Fig. 5 are measurements of the beam profile widths for the frequencies shown. These measurements were made using a HDPE dielectric lens with a fixed focal length of 57 cm and serve to illustrate beam radius variation with frequency.

FIG. 5.

Calculated Gaussian beam radii and measured radii using optimized lens focal length and lens to antenna distance for three probe frequencies 50, 60, and 70 GHz. Locations of various plasma radii are indicated by vertical dashed black lines.

FIG. 5.

Calculated Gaussian beam radii and measured radii using optimized lens focal length and lens to antenna distance for three probe frequencies 50, 60, and 70 GHz. Locations of various plasma radii are indicated by vertical dashed black lines.

Close modal

The above elements were integrated into the design shown in Fig. 6, which illustrates the major elements. The DBS probe radiation enters the lower waveguide from the left via a quartz vacuum window (tilted 5° to avoid direct back reflections) and propagates down the waveguide to a 45° miter bend and then to a metal lens (forming the correct beam shape and size within the plasma, Sec. IV). The radiation is reflected from the adjustable flat mirror into the desired poloidal and toroidal angles. The corrugated circular waveguide (8.9 cm inner diameter) supports arbitrary polarization making it ideal for this application. The received CPS radiation enters the upper waveguide [Fig. 6(a)], passes through a metal lens and into the circular waveguide, and is then detected. The DBS retro-reflected signal retraces its outward propagation path back into the waveguide where it is also detected. The polarization, both launched and received, is set by large, external wire grid polarizers (not shown) similar to those in Refs. 10 and 12. The components closest to the plasma are recessed approximately 13 cm from the tile surface, thus avoiding interactions with either fast neutral beam particles (calculated to have approximately 3° deviation or maximum 1.8 cm penetration into the port box) or the plasma itself. The mirror movement is effected by the mirror control rod assembly shown in Fig. 6(b). The mirror is mounted on a “U” shaped yoke by a rotating hinge as shown by the arrow in Fig. 6(b). This movement is controlled by a sliding arm connected to the mirror providing one axis adjustment. The second movement axis is provided by a rotation about the long axis of the mirror assembly [arrow in Fig. 6(b)]. The sliding arm rotates with this long axis rotation allowing independent control in both directions. The mirror control design itself was strongly based upon the successful ECH aiming mechanism.29 

FIG. 6.

(a) Elevation and (b) plan views of the design showing major components with diagrammatic DBS launch and CPS receive beams. The X-mode probe beam is refracted by the plasma.

FIG. 6.

(a) Elevation and (b) plan views of the design showing major components with diagrammatic DBS launch and CPS receive beams. The X-mode probe beam is refracted by the plasma.

Close modal

This work was supported by the US Department of Energy under Nos. DE-FG02-08ER54984 and DE-FC02-04ER54698. Discussions with Dr. John Lohr and Dr. Mirela Cengher are gratefully acknowledged.

This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. DIII-D data shown in this paper can be obtained at https://fusion.gat.com/global/D3D_DMP.

1.
T.
Lehner
,
J. M.
Rax
, and
X. L.
Zou
,
Europhys. Lett.
8
,
759
(
1989
).
2.
L.
Vahala
 et al.,
Phys. Fluids B
4
,
619
(
1992
).
3.
X. L.
Zou
 et al.,
Phys. Rev. Lett.
75
,
1090
(
1995
).
4.
L. G.
Bruskin
 et al.,
Plasma Phys. Controlled Fusion
36
,
681
(
1994
).
5.
L. G.
Bruskin
 et al.,
Plasma Phys. Controlled Fusion
37
,
255
(
1995
).
6.
L.
Colas
 et al.,
Nucl. Fusion
38
,
903
(
1998
).
7.
A.
Mase
 et al.,
Rev. Sci. Instrum.
68
,
454
(
1997
).
8.
Y.
Kogi
 et al.,
Rev. Sci. Instrum.
70
,
99
(
1999
).
9.
A. D.
Gurchenko
 et al.,
Plasma Phys. Rep.
31
,
476
(
2005
).
10.
T. L.
Rhodes
 et al.,
Rev. Sci. Instrum.
85
,
11D838
(
2014
).
11.
J. C.
Hillesheim
 et al.,
Nucl. Fusion
55
,
073024
(
2015
).
12.
K.
Barada
 et al.,
Rev. Sci. Instrum.
87
,
11E601
(
2016
).
13.
V. V.
Bulanin
 et al.,
Plasma Phys. Rep.
26
,
813
(
2000
).
14.
M.
Hirsch
 et al.,
Plasma Phys. Controlled Fusion
43
,
1641
(
2001
).
15.
P.
Hennequin
 et al.,
Rev. Sci. Instrum.
75
,
3881
(
2004
).
16.
G. D.
Conway
 et al.,
Plasma Phys. Controlled Fusion
46
,
951
(
2004
).
17.
M.
Hirsch
 et al.,
Plasma Phys. Controlled Fusion
48
,
S155
(
2006
).
18.
L.
Schmitz
 et al.,
Rev. Sci. Instrum.
79
,
10F113
(
2008
).
19.
X.
Weiwen
 et al.,
Plasma Sci. Technol.
10
,
403
(
2008
).
20.
T.
Happel
 et al.,
Rev. Sci. Instrum.
80
,
073502
(
2009
).
21.
W. A.
Peebles
 et al.,
Rev. Sci. Instrum.
81
,
10D902
(
2010
).
22.
R. E.
Slusher
and
C. M.
Surko
,
Phys. Fluids
23
,
472
(
1980
).
23.
M.
Hirsch
and
E.
Holzhauer
,
Plasma Phys. Controlled Fusion
46
,
593
(
2004
).
24.
A. P.
Smirnov
 et al., in
Proceedings of 15th Workshop on ECE and ECRH
(
Singapore World Scientific
,
2009
), p.
301
, www.compxco.com/genray.
25.
L. L.
Lao
 et al.,
Nucl. Fusion
25
,
1611
(
1985
).
26.
P. F.
Goldsmith
,
IEEE Trans. Antennas Propag.
39
(
6
),
834
(
1991
).
27.
P. F.
Goldsmith
,
Proc. IEEE
80
(
11
),
1729
(
1992
).
28.
J. C.
Hillesheim
 et al.,
Rev. Sci. Instrum.
83
,
10E331
(
2012
).
29.
R.
Ellis
 et al., “
Design of the P2001 ECH launcher for DIII-D
,” in
Proceedings of the 19th IEEE/IPSS Symposium on Fusion Engineering
(
IEEE
,
2002
).