A system for studying the spatiotemporal dynamics of fluctuations in the boundary of the W7-X plasma using the “Gas-Puff Imaging” (GPI) technique has been designed, constructed, installed, and operated. This GPI system addresses a number of challenges specific to long-pulse superconducting devices, such as W7-X, including the long distance between the plasma and the vacuum vessel wall, the long distance between the plasma and diagnostic ports, the range of last closed flux surface (LCFS) locations for different magnetic configurations in W7-X, and management of heat loads on the system’s plasma-facing components. The system features a pair of “converging–diverging” nozzles for partially collimating the gas puffed locally ≈135 mm radially outboard of the plasma boundary, a pop-up turning mirror for viewing the gas puff emission from the side (which also acts as a shutter for the re-entrant vacuum window), and a high-throughput optical system that collects visible emission resulting from the interaction between the puffed gas and the plasma and directs it along a water-cooled re-entrant tube directly onto the 8 × 16 pixel detector array of the fast camera. The DEGAS 2 neutral code was used to simulate the Hα (656 nm) and HeI (587 nm) line emission expected from well-characterized gas-puffs of H2 and He and excited within typical edge plasma profiles in W7-X, thereby predicting line brightnesses used to reduce the risks associated with system sensitivity and placement of the field of view. Operation of GPI on W7-X shows excellent signal-to-noise ratios (>100 at 2 Mframes/s) over the field of view for minimally perturbing gas puffs. The GPI system provides detailed measurements of the two-dimensional (radial and poloidal) dynamics of plasma fluctuations in the W7-X edge and scrape-off layer and in and around the magnetic islands outside the LCFS that make up the island divertor configuration employed on W7-X.

The understanding of edge and scrape-off layer (SOL) transport in magnetically confined plasmas continues to be of crucial importance in the development of a viable fusion power plant. It affects critical aspects of energy and particle exhaust handling, e.g., divertor and first-wall lifetimes, and fuel recycling, as well as critical aspects of core plasma performance (via the pedestal) and fueling. It has long been recognized that turbulence is important to edge transport, and the past two decades have seen significant advances in our understanding of edge phenomena in tokamaks and stellarators. These advances have been made possible by the combination of advanced theoretical/computational models and detailed diagnostic information. In particular, a number of imaging diagnostic systems yielding two-dimensional information have been developed for the visualization of turbulent dynamics, including Beam Emission Spectroscopy (BES),1,2 Electron Cyclotron Emission Imaging (ECEI),3–5 Microwave Imaging Reflectometry (MIR),6,7 and (most importantly in the context of this report) Gas-Puff-Imaging (GPI).8–10 GPI is a mature diagnostic technique, which has been of particular importance for identifying the structure and dynamics of edge phenomena in tokamaks: Alcator C-Mod,11 NSTX,10 ASDEX-Upgrade,12 TEXTOR,13 EAST,14 and TCV,15 and in a reversed-field-pinch—RFX-Mod16—see the review of GPI by Zweben et al.17 This report describes the design, component characterization, and performance of a GPI system deployed on the Wendelstein 7-X (W7-X) stellarator.18–20 

GPI has proven to be valuable for understanding a number of key physics issues of the plasma boundary and the SOL in tokamaks, including but not limited to the following examples: the dynamics of highly intermittent fluctuations in the SOL and divertor, often referred to as filaments or blobs;12,21,22 the dominant role of these filaments in the far-SOL (see, e.g., Ref. 23); the development of models for filament intermittency and filament dynamics;24,25 dynamics during the L–H transition;26,27 SOL dynamics during ELMs;28 the role that filamentary radial transport plays in the tokamak density limit;29 the non-linear turbulent kinetic energy transfer from the background turbulence into sheared quasi-static flows;30 zonal flows and geodesic acoustic modes;31 and the characterization of coherent electromagnetic modes.11,32

While filaments are typically present in the W7-X SOL, their transport and intermittency statistics are very different from what is found in tokamaks.33,34 The presence of edge magnetic islands in W7-X appears to play a role in the filament dynamics, as does the field-line averaged curvature.33–35 2D imaging by GPI at the outboard island in W7-X will help to elucidate those dynamics and the effects of the islands on the filaments.36 Regarding edge coherent modes in W7-X, a number have been observed,33,37–40 and GPI will further characterize and investigate those modes.

The GPI technique is well-described in the GPI review publication.17 We reiterate the essentials and key components of the technique here. A “gas puff” component provides a source of hydrogenic or helium atoms with a limited toroidal extent at the edge of the plasma. (A vertical “sheet” of atoms is the ideal.) As the atoms from the puff penetrate the edge plasma, line emission is excited by the local plasma electrons. Under conditions41 typically met in the edge plasmas of stellarators and tokamaks, the line emission from Balmerα (when puffing H2 or D2) or the 33D − 23P transition of HeI at 587 nm (when puffing He) will represent excitation by unique ne and Te in the local plasma and therefore also the combined effects of fluctuations in those quantities. The line emission is collected from the region “illuminated” by the gas puff and conveyed, by the “light collection” component, to a fast-framing two-dimensional detector, the “fast camera” component. The viewing geometry is important. Ideally, the viewing chords should be aligned with the magnetic field lines at the gas puff location since the fluctuations are primarily field-aligned and the most interesting viewing field is therefore perpendicular to them. Any toroidal extent of the emission from the gas cloud will limit the spatial resolution in the images compared to that obtainable using the ideal “sheet” illumination. While the GPI concept is relatively simple, the realization of an actual system on a magnetic-confinement fusion-relevant device, such as W7-X, is quite challenging as will be discussed below. The CAD model of these GPI system components, as realized on W7-X, is shown in Fig. 1.

FIG. 1.

CAD drawing of the GPI system on W7-X. The plasma (not shown) occupies the bottom left of the drawing. The “gas puff” component is located in the horizontal vessel access port on the right. The “light collection” component is in the water-cooled re-entrant immersion tube in the port on the left. The “fast camera” component (in pink) is mounted at the end of the tube just beyond the large port flange (in yellow). The distance from this flange to the plasma vessel wall is ≈1.9 m, and from there to the LCFS of the “standard” W7-X magnetic field at the location of the nozzle is ≈0.14 m.

FIG. 1.

CAD drawing of the GPI system on W7-X. The plasma (not shown) occupies the bottom left of the drawing. The “gas puff” component is located in the horizontal vessel access port on the right. The “light collection” component is in the water-cooled re-entrant immersion tube in the port on the left. The “fast camera” component (in pink) is mounted at the end of the tube just beyond the large port flange (in yellow). The distance from this flange to the plasma vessel wall is ≈1.9 m, and from there to the LCFS of the “standard” W7-X magnetic field at the location of the nozzle is ≈0.14 m.

Close modal

W7-X is a large superconducting stellarator, which provides plasma parameters and engineering demonstrations relevant to a future stellarator power plant. Its design was based on optimizing seven criteria, one of which is reduced neoclassical transport of the thermal plasma.19 It is developing the capability to maintain plasma durations of thirty minutes, and plasmas of 8 minutes duration and over 1 GJ of input (and exhaust) energy have already been realized.20 W7-X can accommodate a variety of magnetic configurations and utilizes the magnetic island divertor concept for heat and particle exhaust handling.42 Elucidation of the islands’ role in SOL transport is a primary goal for GPI investigations. These realities present significant challenges for implementing GPI on W7-X. As noted above, the “gas puff” component must be in-vessel and close to the plasma boundary. Ports with tangential views of the plasma that fulfill the requirement of field-aligned viewing chords are too small for the light collection requirements. Therefore, the “light collection” component must enter from one of the large horizontal ports and have parts relatively close to the plasma. These port flanges are roughly 2 m radially from the plasma edge in order to accommodate the coils and cryostat. Furthermore, the variety of magnetic configurations created in W7-X results in a range of last closed flux surface (LCFS) locations and island geometries that affect the decisions about the placements of the nozzle and field of view (FOV). We summarize the design challenges for GPI on W7-X below. We believe that each represents a challenge not faced by any current or previously existing GPI system.

  • The variety of available SOL and island locations in the W7-X configuration space requires a large distance (∼0.1 m) between the “gas puff” component and the field of view in most of those configurations. This, in turn, makes it highly desirable that the flux of atoms/molecules from the nozzle be somewhat collimated in the toroidal dimension in order to minimize the loss in spatial resolution due to smearing of field-aligned fluctuation features brought about by chordal integration through the gas cloud.

  • The variety of LCFS locations and island locations makes it desirable to have the capability to change the FOV.

  • Longer pulse/steady-state higher-power operation is planned for W7-X, and plasma radiation from those discharges will result in large heat loads on in-vessel diagnostics with plasma-facing components. These are required to be designed to withstand sustained radiative heat fluxes from the plasma boundary of up to 100 kW/m2. For GPI, this means thoughtful heat load management for its in-vessel components, especially the vacuum window, plasma-facing components, and nozzle structure.

  • The heat load concerns precluded the use of in-vessel optical fibers that would transmit light gathered from imaging optics mounted (say) on the vessel wall to the ex-vessel “fast camera” component. This required that a re-entrant tube of ∼2 m radial extent be used to provide the view of the gas puff and transmit the image to the “fast camera.”

In addition to the challenges specific to W7-X, there are following challenges that must be overcome for any successful GPI realization:

  • Puff enough gas and gather enough light on turbulence timescales (1μs) for good signal-to-noise in the images.

  • Image with a spatial resolution good enough to discern the turbulence spatial scales of greatest interest, targeting in the W7-X case the filament dynamics and edge quasi-coherent modes, targeting kρs in the range 0.01 to 0.1. This includes minimizing the toroidal spread of the gas puff at the FOV of the light collection component as well as choosing the spot size of a detector pixel in the object plane of the collecting optics, i.e., in the focal plane in front of the “gas puff” component/nozzles. The spatial resolution (5 mm) is analyzed in detail in Sec. IV B.

  • Do not puff gas in amounts that will significantly perturb either the global plasma or ne and Te of the region whose turbulence characteristics are being measured. This is addressed in detail in Sec. IV A.

The challenges listed in Sec. I B formed the basis for the design criteria of the system that was ultimately realized. Those criteria are summarized in Table I, along with the final design decisions.

TABLE I.

Design criteria for the GPI system, the realized design decisions, and their motivations. The “standard” W7-X magnetic configuration is the most frequently run n/m = 5/5 island configuration, with n (m) as the toroidal (poloidal) mode number of the island chain. The “bean” shape of the plasma is the most vertically extended cross section of the W7-X plasma.

ParameterConstraintReason for constraintRealized value
Angle between sight lines and local B <15° Ideally parallel to reduce spatial smearing 8–11° 
Poloidal/radial location of FOV Centered between 0 and 60 mm out from LCFS and contain island O-pt. of “standard” mag. config Diagnostic should measure scrape-off layer fluctuations and dynamics at island Center of FOV is ≈25 mm outside the LCFS and contains O-pt 
Radial dist. between FOV and gas-puff nozzles ≤150 mm Gas cloud spreads toroidally and vertically reducing resolution and brightness with distance (see Sec. II E 2110 mm 
Dist. between puff plane and turning mirror 600 ≤ d ≤ 1000 mm Maximize practical etendue 778 mm 
Toroidal location of FOV Within 15° of “bean” shape Core ITG turbulence predicted to be greatest outboard of “bean” shape 13.8° 
ParameterConstraintReason for constraintRealized value
Angle between sight lines and local B <15° Ideally parallel to reduce spatial smearing 8–11° 
Poloidal/radial location of FOV Centered between 0 and 60 mm out from LCFS and contain island O-pt. of “standard” mag. config Diagnostic should measure scrape-off layer fluctuations and dynamics at island Center of FOV is ≈25 mm outside the LCFS and contains O-pt 
Radial dist. between FOV and gas-puff nozzles ≤150 mm Gas cloud spreads toroidally and vertically reducing resolution and brightness with distance (see Sec. II E 2110 mm 
Dist. between puff plane and turning mirror 600 ≤ d ≤ 1000 mm Maximize practical etendue 778 mm 
Toroidal location of FOV Within 15° of “bean” shape Core ITG turbulence predicted to be greatest outboard of “bean” shape 13.8° 

In this section, we describe the details of the various GPI components as realized on W7-X and provide, in some cases, reasons for the design choices. The realized system accommodates puffing H2 or He and imaging either Hα or HeI 587 nm line emission. D2 puffing with Dα imaging could also be accomplished with only a change in the gas bottle if W7-X were to run deuterium plasmas, which is foreseen in the future.

In order to maximize the chances of satisfying the crucial criterion of a good signal-to-noise ratio in the images with minimal plasma perturbation by the gas puff, we chose to image the collected light directly onto detectors with optimal sensitivity and an optimal signal-to-noise ratio at low light incidence. Avalanche photodiodes (APDs) fulfill these sensitivity criteria.43 Dynamic range flexibility is provided by changing the high-voltage bias applied to APDs. A camera that utilizes arrays of APDs and provides on-board high-voltage bias control, adequate magnetic field shielding (up to 100 mT), onboard Peltier cooling of the APD arrays, and onboard 10 Gbit optical Ethernet data-transfer capability is commercially available from Fusion Instruments, Kft.44 The Fusion Instruments APDCAM-10G camera was provided with a custom APD array of 8 × 16 pixels, utilizing four Hamamatsu S8550 4 × 8 pixel array chips45 in the 8 × 16 configuration shown in Fig. 2. All APD pixels are covered by plastic microlenses that increase the effective packing fraction from ∼50% to ∼90%, make the effective pixel size 2.3 × 2.3 mm2, and extend the full detector area to a 20.8 × 38.3 mm2 rectangle. The noise on the 14-bit (16 384 counts maximum) signals is ≃25–30 counts at the controlled detector temperature of 25 °C. The camera can be read out to the data archive at frame rates of up to 4 Msamples/s, although 2 Msamples/s has been the typical frame rate.

FIG. 2.

Photo of the APDCAM-10G detector’s 8 × 16 array of APDs covered by the microlenses. The long dashed rectangles indicate pixels connected to the four independent high-voltage bias generators. Also shown are the five high-intensity red LEDs (circled in red) and the two 1 mm diameter fibers in the detector plane (circled in white).

FIG. 2.

Photo of the APDCAM-10G detector’s 8 × 16 array of APDs covered by the microlenses. The long dashed rectangles indicate pixels connected to the four independent high-voltage bias generators. Also shown are the five high-intensity red LEDs (circled in red) and the two 1 mm diameter fibers in the detector plane (circled in white).

Close modal

A 50 × 50 mm-square bandpass interference filter is mounted 14 mm in front of the detector plane. Special high transmission (97%–99%) filters46 are used (for Hα: center-λ = 657 nm with 5.2 nm passband; for He I: center-λ = 587.4 nm with 4.4 nm passband).

Figure 2 also shows the five high-brightness LEDs and two 1-mm fibers that are located in the detector plane. These are used (during periods of in-vessel access) as light sources for back-illumination through the collection optics onto a target at the toroidal angle of the gas-puffing nozzles that defines the “FOV registration plane.” Subsection II B describes how this back-illumination procedure achieves an accurate registration of the FOV in W7-X vessel coordinates. The two fibers can also be used to couple collected light to spectrometers if visible spectra from the FOV region are desired.

The usefulness of GPI is limited to plasma regions where there is enough atomic line emission from the puffed gas to allow imaging at the appropriate timescales (1μs). In W7-X, this means GPI is limited to the SOL and regions just inside the LCFS, hence the second constraint listed in Table I that the FOV center be within 0 and 60 mm radially out from the LCFS. Following the choice of detector in the fast camera (pixel spacing ≃2.3 mm) and a desired spatial resolution (∼5 mm) compatible with the cross-field size scale expected of the SOL filamentary structures (∼10 mm),17 the demagnification of the collection optics was chosen to be ∼0.5 (actual value is 0.510—see Sec. II D). The size of the FOV in the object plane of the collection optics is thus 41 × 75 mm. The toroidal location of the GPI “FOV registration plane” is in front of the nozzle. The availability of horizontal ports and mounting choices along with the last three constraints listed in Table I determined the toroidal location of the nozzle at 274.2°, i.e., 13.8° clockwise from the “bean” cross section in the 5th of W7-X’s fivefold symmetric modules (M5). The cross section of the “standard” magnetic field configuration at that toroidal angle is shown in Fig. 3. The FOV-positioning decisions were made in this configuration because it is the one most frequently produced on W7-X. As can be seen in Fig. 3, two GPI fields-of-view can be selected for any given experimental run-day by rotating the camera by 90° about the optical axis, a feature built into the camera-mounting apparatus. In the “default” FOV orientation, the long dimension is tangential to the local flux surfaces and thus in the poloidal direction, yielding 8 radial columns of pixel views, each with 16 poloidal views. Alternatively, rotation of the fast camera about the optical axis into the “rotated” orientation provides 16 radial columns with 8 poloidal views each, covering a greater radial range and including views inside the LCFS. The manual camera rotations are reproducibly achieved in ∼30 min of work at the camera location.

FIG. 3.

(a) Poincaré plot showing the structure of the W7-X “standard” magnetic field configuration at the toroidal angle of the GPI nozzle and “FOV registration plane” of the GPI collection optics. Also shown is the GPI FOV in the “default” orientation. (b) Close-up of the GPI FOVs relative to the magnetic island there. View centers in the “default” (“rotated”) orientation are shown by the red circles (blue stars). The locations of the nozzles are indicated by the red dots at the right. Also shown color-coded are the connection lengths in this region.

FIG. 3.

(a) Poincaré plot showing the structure of the W7-X “standard” magnetic field configuration at the toroidal angle of the GPI nozzle and “FOV registration plane” of the GPI collection optics. Also shown is the GPI FOV in the “default” orientation. (b) Close-up of the GPI FOVs relative to the magnetic island there. View centers in the “default” (“rotated”) orientation are shown by the red circles (blue stars). The locations of the nozzles are indicated by the red dots at the right. Also shown color-coded are the connection lengths in this region.

Close modal

The “default” FOV, shown in red in Fig. 3, includes an island O-point in the “standard” configuration (one of the design criteria) and is entirely in the SOL. The choice of its radial location was also informed by (time-independent) modeling of emission resulting from the interaction of the SOL plasma with the gas puff (see Sec. III). The FOV was placed so that the radial location of the peak line emission predicted by the modeling was just inside the outer edge of the FOV under “typical” SOL conditions. Accurate registration of the FOV is achieved by mounting a target onto the nozzle structure whose location is known in W7-X vessel coordinates from in-vessel metrology. This target defines the “FOV registration plane,” which is orthogonal to the front face of the nozzle structure and is essentially in the ϕ = 274° (R, Z) plane. The geometry is illustrated in Fig. 4. The fast-camera LEDs and fibers (see Fig. 2) are back-illuminated through the collection optics and imaged onto the target by the optics. Since the pixel positions relative to the LEDs/fibers are known, the locations viewed by the camera pixels are accurately registered in W7-X coordinates, thereby defining the FOV in this plane. In the “standard” configuration, the radial/normal-to-flux-surface extent of the “default” FOV is from 5 to 44 mm outside of the LCFS. For this and the other common W7-X magnetic configurations, the FOV (relative to the LCFS) are listed in Table II, where it can be seen that there is good coverage in all but the “high iota” and “high mirror” configurations, for which the emission in the outer views is weak because they are so far out in the SOL.

FIG. 4.

CAD drawing of the nozzle structure viewed from above and in front with illustrations of the “FOV registration plane,” the object plane of the collection optics, the optical axis, and a field line, both piercing the center of the FOV. The optical axis is normal to the object plane. The gas cloud issues from the nozzles into the FOV.

FIG. 4.

CAD drawing of the nozzle structure viewed from above and in front with illustrations of the “FOV registration plane,” the object plane of the collection optics, the optical axis, and a field line, both piercing the center of the FOV. The optical axis is normal to the object plane. The gas cloud issues from the nozzles into the FOV.

Close modal
TABLE II.

Radial extents of the GPI FOV in the “default” and “rotated” camera orientations, calculated assuming plasma β = 0 and no toroidal plasma current.

ConfigurationRadial coverage relative to LCFS (mm)Comment
Standard 5 to 44 (“default”) Includes O-pt 
−12 to 65 (“rotated”) Includes O-pt and LCFS 
Low iota −18 to 21 (“default”) Spans LCFS and includes X-pt 
−35 to 42 (“rotated”) Spans LCFS and includes X-pt 
High iota 31 to 70 (“default”) Weak emiss. on outboard views 
13 to 90 (“rotated”)  
Low mirror −3 to 36 (“default”) Spans LCFS 
−20 to 56 (“rotated”) Spans LCFS 
High mirror 31 to 71 (“default”) Includes O-pt 
14 to 91 (“rotated”) Includes O-pt 
ConfigurationRadial coverage relative to LCFS (mm)Comment
Standard 5 to 44 (“default”) Includes O-pt 
−12 to 65 (“rotated”) Includes O-pt and LCFS 
Low iota −18 to 21 (“default”) Spans LCFS and includes X-pt 
−35 to 42 (“rotated”) Spans LCFS and includes X-pt 
High iota 31 to 70 (“default”) Weak emiss. on outboard views 
13 to 90 (“rotated”)  
Low mirror −3 to 36 (“default”) Spans LCFS 
−20 to 56 (“rotated”) Spans LCFS 
High mirror 31 to 71 (“default”) Includes O-pt 
14 to 91 (“rotated”) Includes O-pt 

The gas puff emission is viewed from the horizontal port adjacent to the one holding the nozzle structure (Fig. 1). The views are designed to be as close to B-field-aligned as is feasible. We set as a constraint (the first in Table I) that the sight lines be within 15° of the field lines local to the gas puff. A turning mirror is required to fulfill this requirement, which we decided to combine with the need to shutter the vacuum window of the re-entrant tube holding the collection optics. We mounted a polished stainless steel mirror (flatness of 1 λ at 630 nm and 60/40 scratch dig) on the back side of a hinged shutter plate and control the opening/closing of the shutter and mirror with a pneumatically actuated linear motion vacuum feed-through. The full-open angle of the shutter sets the location of the FOV, so it is critical that hard stops provide a very reproducible “open” condition. Repeated testing showed that the views in the “FOV registration plane” were reproducible to within ≤±1 mm.

The desire to field-align the sight lines needed to be balanced with the need to minimize/eliminate interaction between the shutter and the plasma, including fast ions from the neutral beams. This compromise resulted in the following: (1) we placed the center of the open turning mirror at W7-X coordinates [R, ϕ, Z] = [6363 mm, 4.900 rad, −210 mm]; (2) angles between the sight lines and the B-field local to the gas puff are therefore from 8 to 11°; (3) since we were free to choose the vertical location of the turning mirror, the sight lines and the local field lines piercing the “FOV registration plane” are at essentially the same angle (±2°) relative to horizontal planes; (4) the shutter is fully recessed into the port extension when closed; and (5) the part of the open shutter closest to the plasma is ∼115 mm from the LCFS of the largest of the most common magnetic configurations (i.e., the “low-iota” configuration). Reciprocating probe measurements have shown that the non-radiative power fluxes at locations >100 mm outside the LCFS surface are negligible.47,48 Nonetheless, radiative heat fluxes on the plasma-facing shutter of up to 100 kW/m2 had to be designed for, an issue discussed in Sec. II F.

The 98 mm clear-aperture vacuum window in the re-entrant tube is a critical component for machine safety. The two main threats to the window integrity are radiative heating by the plasma and stray ECRH power. These risks were mitigated by recessing the window by ∼660 mm up the water-cooled re-entrant tube, whose inner diameter is 100 mm up to that window location. The fused silica window has a broadband anti-reflection coating and is rated for temperatures up to 200 °C and for heating rates up to 25 °C/min.49 

In addition to the turning mirror, three lenses in the re-entrant tube on the atmosphere side of the vacuum window make up the “collection optics.” They collect light from sight lines passing through the region in front of the nozzles, transmit it up the re-entrant tube, and focus it onto the detectors. It was a primary design priority to collect as much light as was feasible, given the constraints of the size and availability of windows and lenses. This is accomplished using three 100 mm diameter plano-convex glass lenses with broadband anti-reflection coatings on all surfaces. The optical design is shown in Fig. 5 and is summarized here: the front lens, L1, is located just behind the vacuum window and renders light from the object plane parallel so that it passes down the tube. As illustrated in Fig. 4, the object plane intersects the vertical “FOV registration plane” (normal to the nozzle face) at their respective center points. The “FOV registration plane” is essentially the ϕ = 274° (R, Z) plane, while the object plane is rotated by 20° about a vertical axis and by 4.2° about a horizontal axis with respect to that plane. In other words, the optical axis strikes the “FOV registration plane” at a downward angle of 4.2 deg, while the total angle between the optical axis and the surface normal of this plane is 20°. The focal length of L1 is 1475 mm, equal to the distance from L1 to the “FOV registration plane” along the optical axis. Lenses L2 and L3 are placed together ∼665 mm beyond L1 and focus the parallel light onto the detector plane with an effective focal length of 753 mm, yielding a demagnification of 0.510. The 20° angle between the “object plane” and the “FoV registration plane” results in some differences in the optical properties across the FOV. The (geometrical) optics imaging properties of the system have been obtained using the ZEMAX-OpticStudio50 commercial software. The modeled images of points in the detector plane back-imaged into the “FoV registration plane” have diameters (90% enclosed energy) ranging from 0.08 to ≈0.3 mm, with the largest spot sizes occurring at largest and smallest Rmaj locations in the FOV. However, this imaging only increases the FWHM of back-imaged pixels from perfectly imaged 4.5–4.6 mm (minimum) and 4.7 mm (maximum). Thus, we will treat the optical resolution in the “FoV registration plane” as 4.65 mm. Thus, 2.3 × 2.3 mm2 pixels collect essentially all of their incident light from 4.65 × 4.65 mm2 areas in the “FoV registration plane.” The aperture stop in the system is the 96 mm diameter clear aperture of the L1 lens mount. The etendue for an on-axis pixel is 6.5 × 10−2 mm2 ster. Light from off-axis points is vignetted because it is at a small angle with respect to the central ray on being rendered parallel by L1, yielding an etendue for the most off-axis pixels (the 4 corner pixels) of 4.7 × 10−2 mm2 ster. The lenses are secured in the re-entrant tube using a removable stiff frame of angled bars that is secured at each end. This can be seen in the cut-away rendering in Fig. 1.

FIG. 5.

Optical layout of collection optics.

FIG. 5.

Optical layout of collection optics.

Close modal

The gas puffing component of the system consists of the nozzles, the nozzles’ housing, a capillary feed line, a gas vacuum feed-through, an injection valve, a 0.80 l plenum, and a gas control system. The first four of these items are shown in Fig. 1. The gas control system (not shown in Fig. 1) is mounted on a panel and connects to the injection valve through an ∼1 m long flexible stainless steel (SSTL) hose.

We wanted to locate the nozzle structure as close to the plasma as feasible while avoiding direct heating by plasma conduction and convection. As noted in Subsection II C, reciprocating probe measurements showed that conducted/convected power was quite small at locations ⪆100 mm into the outboard SOL. Considering the range of LCFS locations of the available W7-X magnetic configurations, we placed the front face of the nozzle housing ≈95 mm radially outward from the LCFS of the outermost of the more common configurations. This places the nozzles, which are recessed 2 mm into the housing, ≈135 mm from the LCFS in the “standard” configuration and ≈110 mm from the FOV center. The vertical center of the nozzles was set to align with the FOV center. Its toroidal location was set by mounting it onto the water-cooled port protection liner, thereby minimally affecting viewing access by other diagnostics utilizing that port. The in-vessel metrology activity mentioned in Sec. II B registered the location of the nozzle housing front face in W7-X plasma vessel coordinates at [R, ϕ, Z] = [6156.4 mm, 4.7856 rad, −267.7 mm].

We surrounded the SSTL nozzle body in a graphite housing made from the same material that is used for the W7-X divertor targets. A short section of 1 mm I.D. capillary that couples to the nozzle body via a press fit exits radially out the back of the housing and connects with a 1.9 m length of 1 mm capillary that runs along the side of the port liner to a vacuum feedthrough mounted on the port extension close to the port flange. The puff valve is mounted as close as possible to the feedthrough as part of the effort to minimize the volume between the nozzles and the puff valve since that volume of gas (Vmin = 3.02 ml) must necessarily flow into W7-X. Thus, the minimum gas load on W7-X is po × Vmin, where po is the plenum pressure backing the puff valve. The gas control system performs several functions in addition to opening the puff valve at programmed times and durations, which it can do up to four times during a single W7-X shot/program. The gas control hardware consists of a plenum volume, an absolute pressure gauge to measure the pressure in the plenum, a differential pressure gauge to measure the rate of change in the pressure occurring during each gas puff, valves to access or isolate various parts of the system, and a mechanical pump to evacuate the system into the W7-X exhaust line, when necessary. A schematic of the GPI gas control system is presented in Fig. 6. The timing for the opening and closing of the valves, as well as the monitoring and digitizing of the system pressures, is done by programmed Red Pitaya’s FPGAs/CPU/digitizers. It is mounted on the gas control panel and is interacted with using a Python graphical user interface program running on a network computer.

FIG. 6.

Schematic of the GPI gas control system.

FIG. 6.

Schematic of the GPI gas control system.

Close modal

The plenum pressure is manipulated by the valves V5 and V4 connecting it with either the regulated gas bottle or the mechanical pump while feeding back on the pressure of the absolute pressure gauge, AG.51 Measurements of the gas puff flow rate and total amount puffed are accomplished by digitizing the differential pressure from gauge DG52 after closing the normally open valve V3 for the duration of the gas puffs. With V3 closed, the time history of the pressure difference between the pre-puff pressure and the plenum pressure is measured, and the time derivative of that differential pressure times the plenum volume (Vplenum = 0.802 l) provides the gas flow rate out of the plenum. The pressure range of the differential gauge is 0–39 mbar, with an estimated error of <1% for typical pressures. Thus, the maximum measurable amount of gas input into a single W7-X plasma is limited to 31 mbar l, which is greater than the typical amount delivered in four puffs, i.e. a total of ∼12 mbar l, at ∼3 mbar l per puff.

1. Nozzle design and gas puff flow rate measurements

Having defined and registered the locations of the FOV and the nozzles, the radial distance between the nozzles and the center of the FOV is 110 mm, satisfying one of the design criteria listed in Table I. However, that large distance made improved collimation of the gas cloud a high priority. Another factor in choosing the nozzle design was the requirement of minimizing the perturbation of the plasma by the gas puff, which is essentially a constraint on the gas flow rate. In this subsection, we will present (1) the details of the “converging–diverging” (C–D) nozzle design, which was driven by the collimation and perturbation constraints, (2) its fabrication, and (3) the measured flow rates from the nozzles and comparison to predictions for C–D nozzles.

We characterize the degree of collimation using the half angle of the gas cloud cone as it expands from a single nozzle into the relative vacuum, defined here as α1/2 = tan−1(HWHM/Δ), where HWHM is the half-width at half-maximum of the gas pressure distribution in a normal plane at a distance Δ from the nozzle. The types of nozzles we considered were as follows: unshaped 1 mm diameter capillary tubes, de Laval nozzles, providing, in theory, optimal collimation,53 and C–D nozzles54 with good collimation. The de Laval shape was too difficult to fabricate at the size estimated for the desired flow rate. C–D nozzles were difficult, but possible, to fabricate. Unshaped capillaries were a last resort since our only guidance was a measurement from Ref. 55 of α1/2 ≈ 25° for D2 at 600 mbar backing pressure, and we desired better collimation than that.

We chose to use two nozzles, each with a C–D shape, separated vertically by 35 mm. The two nozzles with that separation were necessary in order to properly “illuminate” the full poloidal FOV in the “default” camera orientation. For the C–D nozzle, one specifies the converging angle of the nozzle shape on the high gas pressure side, the diverging angle on the outflow side, and crucially the “throat” diameter where the cones meet. Under isentropic conditions and with a high enough pressure difference between the input and output sides, Mach 1 flow is achieved at the throat with a flow rate there that is given by56 
(1)
where Nṗ is the gas flow rate [No./s], nthroat is the gas particle density (at the throat), vthroat is the Mach 1 speed, Athroat is the throat area, C(γ) = 0.58 for H2 (and = 0.57 for He) is a constant that relates the gas density and temperature at the throat to the gas density and temperature of the backing reservoir and only depends on γ, the ratio of specific heats for the specific gas (1.41 for H2 and 1.67 for He), po is the backing reservoir pressure in Pa, To is the reservoir gas temperature in K, k is Boltzmann’s constant, and R is the ideal gas constant for that gas, e.g., 4124 J/kg/K for H2.

Measurements of α1/2 on a C–D nozzle with a 350 µm throat diameter54 showed that α1/2 decreased significantly as the backing pressure increased over the range of interest (∼0.4 to ∼1.3 bar). Since the flow rate is proportional to the backing pressure, we had to balance the improved collimation at the higher backing pressures with the risk of perturbing the W7-X plasma. According to Eq. (1), the throat diameter of the C–D nozzle design and the backing pressure determine the gas flow rate. The maximum acceptable gas flow rates were determined by modeling a hypothetical gas puff that perturbs the W7-X global plasma density by 10%. This modeling therefore informs the choice of throat diameter for the C–D nozzle shape.

We modeled the total electron inventory in the W7-X plasma, Netot, as being sustained by a recycling source, Φrec, that is proportional to the total inventory Netot and balanced by a sink defined by an effective particle confinement time, τp. We defined another source, Φext(t), with time-dependence similar to that predicted from finite-element analysis of gas transport occurring during a 50 ms opening of the puff valve to our plenum and accounting for flow through the 1.9 m long feed capillary to two C–D nozzles with throat diameters of roughly 100 μm. The time-profile of the external source was an ∼40 ms rise time, a 50 ms flattop, and a 0.5 s exponential decay. This external puff was assumed to fuel the plasma with an efficiency, ɛfuel, of 0.8, as will be discussed below. The model equation was thus given by
(2)
τp was taken to be 300 ms (using Fig. 8 from Ref. 57). Taking Netot=VplasmaCnel, where Vplasma is the plasma volume (30 l), nel is the line-integrated density, and C is the ratio between the volume-averaged and line-integrated density (≈0.8 m−1), Φrec is evaluated using the measured nel with Φext = 0. We note that using essentially the same model and actual W7-X discharges (produced in 2018) into which He gas was puffed using a divertor gas puff system, the authors of Ref. 58 evaluated τp as (0.258 ± 0.124 s) and the fueling efficiency of those divertor puffs to be in the range of 0.3 to 0.44. Using other 2018 discharges with nel values of 2.2 × 1019, 5.5 × 1019, and 9 × 1019 m−2, we modeled the response to an external gas puff with the time-profile described above and the conservative guess for ɛfuel of 0.8. The peak flow rates that yielded a 10% increase in nel’s of those three modeled discharges are listed in Table III. Also shown there are the backing pressures that according to Eq. (1) would provide those rates through two C–D nozzles, each with a 70 µm throat diameter. A 70 µm throat diameter was the smallest that could reasonably be fabricated.
TABLE III.

Peak H2 flow rates yielding a modeled 10% density perturbation and corresponding required backing pressures for two 70 µm throat diameter C–D nozzles.

W7-X nelPeak flow rate [No./s] (10% perturbation)Backing pressure (bar)
2.2 × 1019 m−2 4 × 1019 0.34 
5.5 × 1019 m−2 14 × 1019 1.17 
9 × 1019 m−2 30 × 1019 2.5 
W7-X nelPeak flow rate [No./s] (10% perturbation)Backing pressure (bar)
2.2 × 1019 m−2 4 × 1019 0.34 
5.5 × 1019 m−2 14 × 1019 1.17 
9 × 1019 m−2 30 × 1019 2.5 

Thus, we specified a 70 μm throat diameter. The half-angle for the diverging cone was specified as 25° since that was seen to be in the optimal range in Ref. 59. A 45° half-angle was chosen for the converging cone. The two nozzles were cut into a trough milled down to 1 mm thickness in a SSTL blank with an overall thickness of 3 mm. This trough provides a 0.62 ml manifold behind the nozzles.

Two vendors willing to attempt fabrication to those specifications were identified. One vendor used fine drills with tips cut at the desired angles; the other used a laser to sculpt the shape by varying the focal spot. Both fabrications were to be electro-polished to a surface finish with a roughness average Ra better than 0.4 μm. Examination of results with a microscope showed that the mechanical drilling produced far superior results. The cone walls were smooth, and the throat was well-defined, albeit with a diameter of 83μm, larger than the specification. The walls of laser-drilled nozzles were rippled; the throat hole was ragged, and this fabrication was not used. The nozzle plate was then vacuum welded onto a back cover plate that holds a press-fit length of 1 mm I.D. capillary protruding out of the backside of the back plate. The assembly was vacuum leak-tested by pumping on the capillary after temporarily blocking the nozzle holes. An exploded view of the nozzle structure is shown in Fig. 7, and an assembled rendering is shown in Fig. 4.

FIG. 7.

CAD rendering of the exploded nozzle structure. From left to right: left clamp, bracket for mounting to the port extension fixture, Sigraflex pad60 for better heat conduction, left half of the graphite housing, the welded SSTL piece with two C–D nozzles (separated by 35 mm vertically) and feed capillary, right half of the graphite housing, and right clamp. The short feed capillary exiting the back of the nozzle structure connects to the longer 1.9 m feed capillary. The clamping and mounting bolts are not shown.

FIG. 7.

CAD rendering of the exploded nozzle structure. From left to right: left clamp, bracket for mounting to the port extension fixture, Sigraflex pad60 for better heat conduction, left half of the graphite housing, the welded SSTL piece with two C–D nozzles (separated by 35 mm vertically) and feed capillary, right half of the graphite housing, and right clamp. The short feed capillary exiting the back of the nozzle structure connects to the longer 1.9 m feed capillary. The clamping and mounting bolts are not shown.

Close modal

The flow rates of H2 and He, the gases of interest for GPI, as well as N2 and Ar were measured multiple times over 4 years using the gas control system. The 1.9 m length of feed capillary was present in the system for these measurements. The results are shown in Fig. 8. For one set of H2, He, and N2 measurements, the puffs were into the W7-X vacuum vessel with the gate valves to the pumps closed and with the vessel pressure measured by a calibrated ASDEX-Upgrade Baratron gauge. Comparisons of the total amount of gas puffed during each sequence at each plenum pressure for each gas showed that (ΔpGPI × Vplenum) and (ΔpW7-X × VW7-X) were the same to within 5%. This provided independent confirmation of the accuracy of the GPI differential pressure measurements. Regarding the flow rate measurements, we note the following:

  • The flow rates are linear with backing pressure as predicted by Eq. (1).

  • The flow rates are ≈25% larger than what is predicted by Eq. (1) using the measured throat diameter. We do not know the reason(s) for this. Note that, in Fig. 8, the 93 μm throat diameter was chosen in order to provide a good match to the experimental flow rates.

  • The scaling among the four different gases is approximately what is predicted by Eq. (1). It is good for H2 and He, but somewhat worse for N2 and Ar relative to H2 and He.

FIG. 8.

Flow rate measurements of four gases as a function of the measured plenum pressure. The dashed lines are the calculated flow rates using Eq. (1) and assuming a 93 μm diameter throat for each of the two nozzles. The solid circles are from puffs into W7-X, while +’s and △’s are from lab measurements.

FIG. 8.

Flow rate measurements of four gases as a function of the measured plenum pressure. The dashed lines are the calculated flow rates using Eq. (1) and assuming a 93 μm diameter throat for each of the two nozzles. The solid circles are from puffs into W7-X, while +’s and △’s are from lab measurements.

Close modal

2. Measurements of gas cloud extent

As noted in Subsection II E 1, we had reason to believe that C–D shaped nozzles would provide some degree of collimation of the outflowing gas and that it should improve with increasing plenum pressure.54 We chose to measure this for the nozzles as fabricated. In this subsection, we present the measurements of the spatial distributions of the gas emerging from the nozzles. These distributions were measured using the vacuum chamber and the computer-controlled positioning system capable of scanning a probe in a 2D plane of the VINETA II linear plasma device61 at IPP. By mounting a custom miniaturized hot filament ionization gauge onto VINETA’s positioning system, we measured the spatial distributions of the local pressure within the puffed gas cloud as functions of gas backing pressure.

The vacuum chamber of VINETA II is roughly cylindrical, ∼1.6 m long and ∼0.97 m in diameter, with a typical base pressure of 2×104 Pa. The custom ionization gauge is based on a design from Ref. 62. It consists of three tips protruding ∼3 mm from channels in a long 5 mm diameter ceramic cylinder. One tip is a U-shaped loop of tungsten wire that is heated to thermionic emission and acts as an electron source (the cathode at V = 0, i.e., ground) for a second probe (the anode) that is biased (at +90 V) to draw ∼1 mA of electron current. The third probe (the collector) is biased (at −60 V) to collect ions of the gas local to the probe tips that have been ionized by the cathode-to-anode electrons. The anode and collector are separated by ≈1 mm. The collector current is amplified and digitized. The electrical circuit is shown in Fig. 9.

FIG. 9.

Circuit diagram for the miniaturized ion gauge probe. Inset: an example of the calibration of the ion gauge for H2 gas.

FIG. 9.

Circuit diagram for the miniaturized ion gauge probe. Inset: an example of the calibration of the ion gauge for H2 gas.

Close modal

The measured ion currents are small (0–0.25 μA), non-linear (but monotonic) functions of the gas density at the probe and dependent upon the gas species. Calibrations of the ion current response to the pressure of each puffed-gas species were performed by introducing increasing amounts of gas into VINETA II via a needle valve and logging the absolute pressure measured by a Baratron capacitive pressure gauge along with the measured ion current. The calibration range covered the pressure range encountered in the measurements on the GPI puffs into VINETA II (from about 0.005 to 0.5 Pa), and calibrations were performed for each gas species and after each vacuum break. An example of a calibration curve for H2 is shown in the inset of Fig. 9.

The miniaturized ion gauge probe was scanned spatially within planes that were 54 and 102 mm in front of and normal to the 2-nozzle plate. Additionally, time histories at a single location 54 mm directly in front of the top nozzle for puffs with different valve-open durations were measured, shown in Fig. 10. The rise time is ≈40 ms, and the decay time is ≈350 ms. The long decay time is consistent with our finite-element modeling of the gas transport through the system, and indicative of the gas in the volume between the valve and the nozzles (Vmin) draining out of the nozzles after the valve is closed. Also evident in Fig. 10 are the “flattop” durations that roughly match the valve-open durations 0.1 s.

FIG. 10.

Time histories of pressure in front of the top nozzle for different valve-open durations. The backing pressure was 1.333 bar of He gas.

FIG. 10.

Time histories of pressure in front of the top nozzle for different valve-open durations. The backing pressure was 1.333 bar of He gas.

Close modal

2D time-resolved pressure measurements were made for He puffs at eight backing pressures in the 54 mm plane and at 3 backing pressures in the 102 mm plane and for H2 puffs at two backing pressures in the 54 mm plane. For each gas puff, the probe responded to the rapid direct pressure rise within the gas cloud as well as the less swift “base pressure” rise in VINETA II. The latter was therefore measured with the probe well away from the nozzle during the puff and subtracted from the overall signal at each measurement location within the gas cloud. Full 2D scans showed left–right (toroidal) and up–down (vertical) symmetry of the gas distribution around the center of the nozzle structure. An example of a pressure distribution along a cut in the toroidal dimension at the height of the top nozzle is shown in Fig. 11. The distribution is well described by a Lorentzian whose HWHM is equal to α1/2ϕ. The pressure distribution along a cut in the vertical dimension centered at the peak of the pressure in the toroidal dimension is shown in Fig. 12. It is well described by two Lorentzian distributions separated by 35 mm with equal HWHMs (=α1/2Z). The results for the α1/2ϕ evaluations vs backing pressure are summarized in Fig. 13. The measurements were evaluated during the puff “flattop,” during which time α1/2s are smallest. We plot only α1/2ϕ since they are the ones that matter for constraining the spatial resolution, but note that the α1/2Z values were very close to the α1/2ϕ values. We note the following:

  • The degree of collimation improves with the increasing backing pressure, as was also observed in Ref. 54. The minimum α1/2 appears to be 12°.

  • Also plotted are the half angles measured for the gas distribution from a single unshaped capillary with a 1 mm I.D. It is apparent that our C–D nozzles improve the collimation by roughly a factor of 2.

  • The half angles in the plane 54 mm from the nozzles are essentially the same as those in the plane at 102 mm. Recall that the center of the FOV is in the plane 110 mm away from the front surface of the nozzle plate.

FIG. 11.

Measured gas puff pressure profile in the toroidal direction 102 mm in front of the nozzles for a plenum pressure of 1.07 bar.

FIG. 11.

Measured gas puff pressure profile in the toroidal direction 102 mm in front of the nozzles for a plenum pressure of 1.07 bar.

Close modal
FIG. 12.

Measured gas puff pressure profile measured in the Z direction 102 mm in front of the nozzles for a plenum pressure of 1.07 bar.

FIG. 12.

Measured gas puff pressure profile measured in the Z direction 102 mm in front of the nozzles for a plenum pressure of 1.07 bar.

Close modal
FIG. 13.

Summary of the gas puff’s toroidal half-angle vs plenum pressure measurements. The black and red circles are for He puffs measured in the 54 and 102 mm planes, respectively. The blue diamonds are for H2 puffs. The blue stars are α1/2 for H2 distributions from a single unshaped 1 mm I.D. capillary and are included for comparison with those from the C–D nozzles.

FIG. 13.

Summary of the gas puff’s toroidal half-angle vs plenum pressure measurements. The black and red circles are for He puffs measured in the 54 and 102 mm planes, respectively. The blue diamonds are for H2 puffs. The blue stars are α1/2 for H2 distributions from a single unshaped 1 mm I.D. capillary and are included for comparison with those from the C–D nozzles.

Close modal

We placed the diagnostic’s plasma-facing components in locations where heat loading due to plasma conduction and convection was likely to be small. However, radiative loads on those components cannot be avoided, and, as noted in the Introduction, all in-vessel components on W7-X must be designed to withstand long radiative emission from the plasma boundary of up to 100 kW/m2. W7-X has infrastructure for extensive water cooling that we used for active cooling of GPI’s re-entrant tube. In addition, the nozzle structure was mounted to the water-cooled port liner assembly. Since the nozzle housing is made from the same graphite that is used for the W7-X divertor tiles and the nozzles are SSTL, we depended upon heat conduction to the port liner and radiative loss to keep the peak finite-element modeled temperature of the nozzle structure 700°C when subjected to a 100 kW/m2 steady-state power flux emanating from the plasma.

Active cooling of the re-entrant tube was provided by 12 mm I.D. tubes that are welded to the outside of the re-entrant tube and carry flowing water. There are six channels distributed around the larger diameter tube that holds the optics and four channels distributed around the 100 mm I.D. tube on which the turning mirror/shutter is mounted. These two cooling circuits are independently fed by the W7-X cooling water at pressures of ∼5 bar and a temperature of ∼20 °C. Tests and modeling of the six-channel circuit showed that the realized flow rates (∼1.5 m/s) in the different channels are similar, with a full exchange of water every 2–3 s. The cooling tubes run lengthwise and are combined in a collector ring around the vacuum window position on the larger tube and in a “cold foot” on the smaller tube to which the mirror/shutter attaches. The collector ring provides additional local cooling to the window, which is a critical machine safety component. Modeling shows that a 100 kW/m2 steady-state power flux emanating from the plasma and intercepted by the re-entrant tube in the actual geometry increases the water temperature in the circuit by ∼5 °C for an assumed flow rate of 2 m/s. The main design difficulty arises from the fact that the component intercepting the most heat flux, i.e., the movable shutter/turning mirror, could not be actively cooled.

The design solution for keeping the shutter/mirror within an acceptable temperature bound was to use flexible “straps” that improve the thermal conductivity between the movable shutter and the water-cooled “cold foot” described above. The flexible “straps” are made of stacks of 425 very thin (25 μm) Cu foils.63 A labeled photo of the shutter/mirror component is shown in Fig. 14. The “front” ends of each of the two straps are coupled to a Cu plate sandwiched between the SSTL mirror on the inside of the shutter and the SSTL plasma-facing side of the shutter. Sigraflex pads60 (thin, compressible, high-vacuum-compatible graphite pads with good thermal conduction) are used at each coupling interface. The measured thermal conductivity from the sandwiched Cu plate surface to the cold foot was ≈1.0 W/K. With the straps, the finite-element modeled response of the shutter temperature to plasma boundary emission of 100 kW/m2 yields the temperature distribution on the shutter assembly shown in Fig. 15, where the shield covering the straps (not shown in Fig. 14) reaches a peak steady-state value of 700°C, and the plasma-facing surface of the shutter reaches 500°C in the modeling. The straps reduce the steady-state temperatures by 20%–30% at high power loading but result in a much greater reduction through heat conduction at lower loads where the radiative loss from the components is small.

FIG. 14.

Photo of the front end of the re-entrant tube showing the shutter and the two flexible Cu straps. The cooling water tubes and cold foot are also apparent. The shield covering for the Cu straps (whose inside surface was Cu-plated) has been removed for the photo.

FIG. 14.

Photo of the front end of the re-entrant tube showing the shutter and the two flexible Cu straps. The cooling water tubes and cold foot are also apparent. The shield covering for the Cu straps (whose inside surface was Cu-plated) has been removed for the photo.

Close modal
FIG. 15.

The temperature distribution on the shutter structure when subjected to a 100 kW/m2 radiative heat flux from the plasma boundary as modeled using a finite-element heat transfer code.

FIG. 15.

The temperature distribution on the shutter structure when subjected to a 100 kW/m2 radiative heat flux from the plasma boundary as modeled using a finite-element heat transfer code.

Close modal

In order to assess the actual temperatures on the shutter and on the vacuum window, thermocouples were placed on the front face of the shutter and on the quartz window surface just inside the window edge. These temperatures are monitored at all times during W7-X operation.

It was important to model the expected interaction between the W7-X edge plasma and the gas puff to reduce the risk that we would not collect enough light to perform the imaging and to guide the radial placement of the FOV. The modeling tools that were used were the 3D Monte Carlo neutral code DEGAS 264 and the kinetic neutral 1D transport code KN1D.65 Both codes compute steady-state solutions. DEGAS 2 uses user-input profiles of Te and ne with Monte Carlo neutrals (H2 or He) launched in a user-defined geometry and calculates, among other things, the 3D emissivities of Hα or He I (587 nm) line radiation. A viewing geometry is also input, and chord brightnesses in that geometry are computed from the emissivities. Even though DEGAS 2 is 3D for the neutrals and the emissivities, it uses an axisymmetric target plasma. To adapt this constraint to the non-axisymmetric W7-X geometry, we assumed that the relevant GPI region was small enough that it could be approximated by a greatly expanded axisymmetric Alcator C-Mod equilibrium, whose equilibrium files were readily available. The LCFS of the expanded equilibrium was placed at a major radius R of 6.015 m, and the gas puff was centered vertically at the midplane. Thus, the plasma there was up–down symmetric and not tilted as in the W7-X case (see Figs. 3 and 21). We specified the nozzle spatial locations relative to the nominal LCFS of the equilibrium and the Te and ne profiles relative to the same LCFS. We specified the viewing geometry by maintaining the W7-X viewing angles relative to the field lines local to the design FOV at the gas cloud as well as the path lengths for integration of the emission through the gas cloud as they would be on W7-X. The Te and ne profiles were approximated using measurements from the W7-X reciprocating probe system in the “standard” magnetic configuration48,66 made on discharges produced during the 2018 experimental campaign, OP 1.2. Simulations of H2 and He gas puffs with cloud half-angles of 18° and flow rates of 4.3 × 1019 No./s were performed. The resulting chord brightnesses in the simulated viewing geometry are shown in Fig. 16 for Hα and Fig. 17 for the He I (587 nm) line.

FIG. 16.

Results of DEGAS 2 and KN1D simulations of H2 flowrates of 4.3 × 1019 molecules/s into a SOL with “standard” configuration ne and Te profiles. The black circles show the DEGAS 2-predicted Hα brightness profile in units of mW/cm2/ster. The red curve shows the KN1D Hα emissivity profile scaled to match the DEGAS 2 brightness profile. The blue crosses show the measured time-averaged brightness profile during the first H2 puff into W7-X plasma 20230323.049 scaled by the ratio of the simulation gas flow rate over the actual flow rate for that gas puff (see Sec. IV B). The shaded regions show the extents of the viewed profiles with the FOV with the camera in the “default” and “rotated” orientations.

FIG. 16.

Results of DEGAS 2 and KN1D simulations of H2 flowrates of 4.3 × 1019 molecules/s into a SOL with “standard” configuration ne and Te profiles. The black circles show the DEGAS 2-predicted Hα brightness profile in units of mW/cm2/ster. The red curve shows the KN1D Hα emissivity profile scaled to match the DEGAS 2 brightness profile. The blue crosses show the measured time-averaged brightness profile during the first H2 puff into W7-X plasma 20230323.049 scaled by the ratio of the simulation gas flow rate over the actual flow rate for that gas puff (see Sec. IV B). The shaded regions show the extents of the viewed profiles with the FOV with the camera in the “default” and “rotated” orientations.

Close modal
FIG. 17.

Same as Fig. 16 but with the HeI (587 nm) line profiles predicted for He flowrates of 4.3 × 1019 atoms/s. The blue crosses show 3× the measured time-averaged brightness profile of HeI (587 nm) line at the peak of a He gas puff into W7-X He-majority plasma 20230330.031 scaled by the ratio of the simulation gas flow rate over the actual flow rate for that gas puff (see Sec. IV B).

FIG. 17.

Same as Fig. 16 but with the HeI (587 nm) line profiles predicted for He flowrates of 4.3 × 1019 atoms/s. The blue crosses show 3× the measured time-averaged brightness profile of HeI (587 nm) line at the peak of a He gas puff into W7-X He-majority plasma 20230330.031 scaled by the ratio of the simulation gas flow rate over the actual flow rate for that gas puff (see Sec. IV B).

Close modal

The 1D kinetic neutrals code, KN1D, was run for hydrogen, where it computes the neutral transport, as well as neutral–plasma interactions (dissociation, ionization, etc.) with stationary 1D radial plasma profiles, and yields, among other things, radial profiles of H, H2, Lymanα emission, and Hα emission. It considers both the atoms and the molecules. It is run with a boundary condition specifying the radially inward neutral particle flux at a boundary far from the plasma. It was also run for He, yielding radial profiles of He I atoms and HeI 587 nm line emission. We used the same Te and ne profiles that were used for the DEGAS 2 simulations. Since KN1D is a 1D simulation and does not consider the cloud shape and viewing geometry, we were primarily interested only in the radial profile of the line emissivities, and those profiles are compared with the brightness profiles from the DEGAS 2 simulations in Figs. 16 and 17. The KN1D emissivity profiles have been scaled to roughly match the DEGAS 2 brightness profiles.

We primarily learned two important things from these simulations:

  • The predicted brightnesses over the FOV are in the range ∼0.2 to ∼2.5 mW/cm2/ster for puff flow rates of 4×1019 No./s into “standard” configuration plasmas. Recall (Table III) that modeling indicated that this flow rate would produce a 10% density perturbation of a low-density W7-X plasma.

  • Locating the FOV in its design location places the maximum in the predicted brightness profile near the outside edge of the “default” orientation FOV (see Fig. 3).

Using this brightness range, we estimated the expected signal levels on the detectors in order to ensure that the system is sensitive enough to perform its imaging function. The estimated photon flux onto an on-axis pixel is
(3)
where Bph is the brightness in units of [photons/s/unit area/ster], BmW is the brightness in units of [mW/cm2/ster]; Rm is the reflectivity of the turning mirror; the Ts are the transmissions of the window, lenses, and filter, respectively, and AdetΩimage is the on-axis etendue. The estimate (supplied by the manufacturer) of the APDCAM detector sensitivity at a mid-range bias-voltage of 360 V was as follows: 7 × 106 ph/s corresponds to 1 digitization count. Thus, an Hα brightness of 1 mW/cm2/ster with a bias-voltage of 360 V would result in a 170 000 count signal level. Since this is 10.4× greater than the 14-bit full-scale level, we were assured that the system sensitivity was more than adequate and that a neutral density filter may be required in the system to reduce the light levels. This is addressed in Subsection IV B.

The GPI system was installed on W7-X during the 2021–2022 time period between experimental campaigns OP 1.2 and OP 2.1 and commissioned just prior to the start of OP 2.1. Before describing the performance details of the key GPI components and the images obtained, we note that GPI measurements were made on over 80% of the viable W7-X shots/programs during the OP 2.1 campaign (November 2022 through March 2023), i.e., ∼1800 gas puffs distributed over ∼800 discharges.

We address the machine safety issues first. Operationally, we chose to open the shutter at the beginning of a run day, i.e., not cycling it (closed–open–closed) for each program. The maximum temperature increase registered by the shutter thermocouple due to a single plasma program was 55 °C in response to a 113 s long Pin = 5 MW plasma with a 100% radiation fraction. An initial rapid decay from that temperature increase, followed by a much longer decay time, indicates that the longer decay time is more representative of the bulk temperature of the shutter. Thus, the maximum bulk temperature rise for the shutter is probably significantly lower than 55 °C. In any case, heat loads from OP 2.1 plasmas were not a problem for the shutter/mirror structure, implying as well that the radiative fluxes from the OP 2.1 plasmas are much less than 100 kW/m2. This includes the plasma program with 1.3 GJ of input power, a W7-X record to date.20 Furthermore, the maximum temperature rise of the vacuum window from any W7-X plasma during the campaign, as measured by the window thermocouple, was 1.4 °C and, therefore, also appears to be robust to operational thermal loading.

Since a key heating system for W7-X is the Electron Cyclotron Resonance Heating (ECH) system (providing up to 7.5 MW) and since the GPI re-entrant tube is in a port adjacent to an ECH launcher port, we were also concerned about a heat load on the fast camera from stray ECH radiation that propagates up the re-entrant tube to the camera. We were not confident in the results from attempts to model this possible heat load on the camera. Thus, during OP 2.1, we measured the ECH heat flux through the aperture immediately in front of the camera using an ECH bolometer provided by the ECH group.67 With the GPI shutter open and using the same 113-s PinECH = 5 MW plasma discharge that resulted in the 55 °C temperature increase on the shutter thermocouple, the stray ECH power flux passing through the aperture just in front of the camera was found to be 70 W/m2, a level much smaller than one that would damage the camera.

In Subsections IV A and IV B, we review the performance of the GPI system in three crucial areas: the perturbation of the gas puff on the global W7-X plasma, the signal-to-noise ratio in the images, and the spatial resolution of the images.

We quantify the perturbation to the W7-X plasma from the GPI gas puff by examining the line-integrated density obtained from the W7-X single channel dispersion interferometer system.68 The radially directed beam of the interferometer passes through the plasma center and integrates through all (non-island) flux surfaces. It is over 10 m from the gas puff and as such measures only the global plasma response to the puff. The perturbation analysis is complicated by W7-X’s excellent density feedback control that uses the line integrated density, nel(t), as the sensor. The results of the analysis are summarized in Fig. 18, where the maximum change in nel(t) due to the puff is plotted vs the total amount of (H2) gas puffed for both 0.05 and 0.1 s valve-open durations. The red circles are the results of a GPI-plenum pressure scan accomplished during the commissioning phase of the OP 2.1 campaign before the density control feedback system was operational. These feedback-off perturbations are <1.2×1018m2 when the puff amounts are in the range typically used during the campaign, indicated by the blue shaded region. Recalling the perturbation modeling discussed in Subsection II E 1, we note that this perturbation is <6% of the lowest density plasmas (nel ≈ 2 × 1019 m−2) typically run in W7-X and implies a fueling efficiency, ɛfuel, of roughly 0.4.

FIG. 18.

The maximum change in the line-integrated density, Δnel, measured as a function of the amount of H2 puffed in a GPI single puff. Red circles: perturbation with the density control feedback off; black circles: feedback on. Blue band: typical range of puffed amounts. Light green band: minimum detectable change in nel. These evaluations were made on plasmas with nel’s that ranged from 2.2 to 10 × 1019 m−2.

FIG. 18.

The maximum change in the line-integrated density, Δnel, measured as a function of the amount of H2 puffed in a GPI single puff. Red circles: perturbation with the density control feedback off; black circles: feedback on. Blue band: typical range of puffed amounts. Light green band: minimum detectable change in nel. These evaluations were made on plasmas with nel’s that ranged from 2.2 to 10 × 1019 m−2.

Close modal

The vast majority of OP 2.1 GPI measurements were done with the density control feedback on. A sub-sample of those for which a perturbation assessment could be made is shown by the black circles in Fig. 18. Most of those feedback-on points show a maximum nel(t) perturbation <0.4×1018m2, indicating that the system successfully corrects for the GPI injection by reducing the fueling flow rate, with the correction occurring within roughly 0.3 s. The feedback-on points above that value are from plasmas very early in the campaign when the control system was being tuned up. We do not know the reason for the two feedback-on outliers at 5 mbar l. They are from the same run day, and the feedback response looks uncharacteristically slow. In any case, the global perturbations are small, typically ≪3%, and this is a primary reason why the GPI measurements were made on such a large percentage of W7-X shots/programs.

It was made clear in Subsections I B and II A that getting images with excellent signal-to-noise ratio (SNR) was a primary priority. The DEGAS 2 modeling indicated that this would be the case for expected gas puff flow rates, and this proved to be true.

Almost all of the OP 2.1 results were for H2 puffs into H-majority plasmas, although for two run days, helium was puffed into He-majority plasmas. Since the puff perturbations were small, we chose to maintain good spatial resolution by keeping the backing pressures/flow rates in mid-range in order to keep the toroidal spread of the gas cloud small (i.e., with α1/2 ≈ 14°—see Fig. 13). Under these conditions, we had to decrease the light flux onto the detectors to keep them from saturating. This was typically accomplished using a 5% neutral density filter in front of the Hα [or He (587 nm)] interference filter. This also allowed us to maintain high signal levels by changing the APD gains in response to changes in light levels due to different plasma conditions. The gain curve from one of the detector pixels is shown in Fig. 19; a dynamic range of ≈6 is evident.

FIG. 19.

The relative gain measured as a function of applied bias voltage for one of the 8 × 16 APD detector elements. The gain curves for the other APD elements are very similar. A normalized gain of 1.0 has been chosen for the 360 V bias voltage. The solid line is a fit to the measured values (red dots).

FIG. 19.

The relative gain measured as a function of applied bias voltage for one of the 8 × 16 APD detector elements. The gain curves for the other APD elements are very similar. A normalized gain of 1.0 has been chosen for the 360 V bias voltage. The solid line is a fit to the measured values (red dots).

Close modal

Typically, the high voltage for each of the four regions on the APD array was set independently throughout an experiment day so that signals from each region would be as high as possible without incurring any saturation of the 14-bit (16 384 count) digitization range. A typical single pixel signal time history is shown in Fig. 20. The electronic noise on the digitized signals is 25–30 counts. During the “physics-measurement” portion of the time history (the cyan-shaded region), the average signal level in Fig. 20 example is ∼10 000 counts (about 60% of full-scale), yielding an SNR of 370. For plasmas in the “standard” configuration, typical SNRs vary from ∼50 to ∼500 and can vary significantly over a single image. Nonetheless, these typical values are certainly enough to provide very low-noise images. The “hash” on the signal during the puff duration is from emission fluctuations due to plasma fluctuations and is exactly what GPI is supposed to measure. For actual videos of the images, see links noted in Ref. 69.

FIG. 20.

The signal time history from the pixel viewing the labeled [Rmaj, Z] location in the FOV. The absolute Hα brightness is given by the y axis values on the right. The gas from a puff with a 0.05 s valve-open duration and a 1.002 bar backing pressure entered the view at about 2.08 s. The total amount of gas puffed was 3.47 mbar l. The 0.07 s interval (cyan) around the maximum brightness is when the “physics measurements” are typically made. The interval before the puff (pink) gives the level of intrinsic Hα emission.

FIG. 20.

The signal time history from the pixel viewing the labeled [Rmaj, Z] location in the FOV. The absolute Hα brightness is given by the y axis values on the right. The gas from a puff with a 0.05 s valve-open duration and a 1.002 bar backing pressure entered the view at about 2.08 s. The total amount of gas puffed was 3.47 mbar l. The 0.07 s interval (cyan) around the maximum brightness is when the “physics measurements” are typically made. The interval before the puff (pink) gives the level of intrinsic Hα emission.

Close modal

Brightness levels averaged over the FOV (in the “standard” configuration) are typically in the range of 1–2.5 mW/cm2/ster, with an example of the spatial variation of the time-averaged brightness shown in Fig. 21. The absolute brightness calibration comes from a full system calibration on the installed diagnostic using an absolute continuum source placed in the FOV during in-vessel access in 2023 and is not based on the estimates used when evaluating Eq. (3). The long dimension of the FOV in this camera orientation is essentially poloidal, while the short dimension is radial, i.e., normal to the LCFS. The poloidal distribution of the average brightness is roughly centered in the FOV and varies by about a factor of 2 for a constant radial coordinate. The radial profiles have brightness maxima somewhat inside the outermost edge of the FOV and vary by about a factor of three across the FOV at a given poloidal coordinate. Thus, the emission “coverage” throughout the FOV is quite good for the “standard” configuration, as was desired in the design. We have investigated the quality of our DEGAS 2 and KN1D modeling predictions by comparing them with actual measured profiles even though the predictions were based on ne and Te profiles from 2018 W7-X plasmas. The comparisons (made after scaling the measured brightnesses by the ratio of the simulated to measured flow rates) are shown in Figs. 16 and 17, with the Hα profile measurement coming from the time-average shown in Fig. 21. The predictions for H2 puffs are quite good, given that the predictions were made using profiles from a different plasma. The prediction for the He puff is 3× to 4× larger than what was measured although the profile shape is reasonable. We do not know the reason for this overestimation; it cannot be blamed on the different ne and Te profiles.

FIG. 21.

The spatial distribution of the brightness time-averaged over the “physics-measurement” interval described in Fig. 20. The x axis/y axis is the Z (vertical)/R (major radial) coordinate [in m] in the W7-X coordinate system. The GPI system records the images at 2 Msamples/s. A minor radial lineout through the poloidal center of the FOV is the profile used for the comparison made in Fig. 16.

FIG. 21.

The spatial distribution of the brightness time-averaged over the “physics-measurement” interval described in Fig. 20. The x axis/y axis is the Z (vertical)/R (major radial) coordinate [in m] in the W7-X coordinate system. The GPI system records the images at 2 Msamples/s. A minor radial lineout through the poloidal center of the FOV is the profile used for the comparison made in Fig. 16.

Close modal

Of course, the radial distributions of the emission change depending on the magnetic configuration, the position of its LCFS relative to the FOV (see Table II), and the local plasma parameters. There is sufficient emission in all of the common configurations except “high iota,” where emission from the outer part of the FOV is much reduced and necessitates the removal of the 5% neutral density filter.

It is evident in Fig. 20 that there is a non-zero light signal before the gas puff, which is due to intrinsic Hα emission along the line of sight. It is best for imaging local fluctuations that this emission be small compared to the emission due to the gas puff, as it is in the example shown in Fig. 20. We have examined the conditions under which the intrinsic emission may not be sufficiently small. The intrinsic emission is approximately linear with the plasma density but is still small compared to the puff emission as long as the plasma radiative fraction is less than about 60% and the line-integrated density is less than 8×1019 m−2. Above that radiative fraction, the plasmas tend to be detached from the divertor targets, and the puff emission can be small since the radiation front moves radially inwards, out of the FOV. Under those circumstances, the intrinsic emission can be 10%–50% of the puff emission and the images will be somewhat compromised. To investigate those plasmas, it is best to puff He gas into H2 majority plasmas, thereby avoiding the intrinsic emission issue completely. However, in the large majority of OP 2.1 W7-X plasmas, this intrinsic emission was very small and did not compromise the images.

We now examine the issue of spatial resolution in the images, estimating the wavenumber range, the minimum distance between two features that can be resolved, and the minimum size of the image of a very thin field-aligned emission filament when viewed in the actual geometry, which we define as the “instrumental resolution.” Ideally, this latter quantity would be the size of the region viewed by a single pixel, but as will be shown, there is smearing of filament images due to the toroidal extent of the gas cloud.

With the camera in the “default” orientation, 16 pixel views are distributed poloidally. The view centers are separated by Δpol ≈ 4.5 mm, resulting in a kpol range of 2π/(2L) = 0.45 cm−1kpol ⪅ 6.9 cm−1 = 2π/(2Δpol), where the lower limit is based on the assumption that we can discern one-half of a wavelength using the full poloidal length, L, of the FOV. The pixel-view spacing in the radial dimension is non-uniform, but has an average spacing of 5.3 mm and an estimated krad range of 0.84 cm−1krad ⪅ 5.9 cm−1. To resolve two features of equal brightness as distinct, at least three contiguous pixel views must be examined, with the brightness of the center view being less than that of its neighbors. In our case, this spatial resolution criterion gives a poloidal (radial) resolution of 2Δpol ≈ 9 mm (2Δrad ≈ 11 mm), respectively, as long as the “instrumental resolution” allows it, as we will examine next.

If the gas cloud were a thin sheet in the “FOV registration plane,” then the sizes, as well as [Rmaj, Z] locations, assigned to features in the images would be “exact,” limited only by the size of the area viewed by each camera pixel and the accuracy of the spatial calibration. However, the toroidal spread of the gas cloud and the fact that our sight lines are not quite field-aligned within the gas cloud lead to the smearing of field-aligned emission structures in the images. Our analysis of this “instrumental” effect relies on our knowledge of the gas puff’s double-Lorentzian distribution in planes normal to the nozzle face and knowledge of the sight lines over which the light is collected. Since we are most interested in imaging turbulent structures that are largely field-aligned with kk, we model field-aligned filaments. We assume the filament emissivities are proportional to the local gas pressure in the plane parallel to the nozzle face, i.e., the planes in which the distributions were measured. This ignores the reality that those planes are not tangent to the local fluxes surfaces (there is an 9° angle between them), and we presume that this introduces only second-order errors in the analysis. Each modeled filament is assigned a puncture point in the “FOV registration plane.” We examine brightnesses of a large number of sight lines that image to the detector plane, constructing a high-resolution image of the “FOV registration plane.” An example using a set of ten simulated filaments, each with a 2.4 mm diameter circular cross section, is shown in Fig. 22. It shows the “instrumental broadening,” with the images of the thin filaments strongly elongated dominantly in the major radial (Rmaj) dimension, leading to significant radial smearing, especially in the top half of the FOV. Also shown in Fig. 22 are the centers of the views of the actual detector pixels (red dots) along with the viewed area of a single pixel (the white rectangle). It is clear that the modeled smearing can be larger than the radial size of the area viewed by a single pixel. Combining the extent of this smearing with that of the areas viewed by actual pixels enables us to estimate (below) the “instrumental resolution,” which must always be greater than or equal to the area viewed by a single pixel.

FIG. 22.

Modeling results of ten 2.4 mm diameter field-aligned filaments illuminated by the gas cloud mapped onto the “FOV registration plane” viewed by the camera in its “default” orientation. The x axis (y axis) is the distance into the radial (poloidal) dimension of the FOV. The red dots are the centers of the camera pixel views in this plane. The white rectangle is the area viewed by a single pixel.

FIG. 22.

Modeling results of ten 2.4 mm diameter field-aligned filaments illuminated by the gas cloud mapped onto the “FOV registration plane” viewed by the camera in its “default” orientation. The x axis (y axis) is the distance into the radial (poloidal) dimension of the FOV. The red dots are the centers of the camera pixel views in this plane. The white rectangle is the area viewed by a single pixel.

Close modal

We summarize the key findings from our spatial resolution analyses and modeling:

  • The detector-limited wavenumber resolutions are 0.45 cm−1 < kpol < 6.9 cm−1 and 0.84 cm−1 < krad < 5.9 cm−1 (with the camera in the “default” orientation).

  • To resolve two features of equal brightness as distinct in the poloidal (radial) dimension, they must be separated by at least 9 (11) mm (in the “default” orientation). However, the “instrumental resolution” in the top portion of the FOV is such that the 11 mm spatial resolution in the radial dimension is marginal there—see the following.

  • Modeled images of field-aligned structures show elongation dominantly in the Rmaj dimension. The FWHM along the long axes of these elongated projections range from ≈6 to 10 mm. The presence of such elongated structures in the images may complicate any cross-correlation analysis if they are moving in a direction other than one parallel to either the major or minor axis.

  • In the bottom half of the FOV, the radial and poloidal FWHM of the modeled images range from 2.5 to 4.7 mm and 2 to 3 mm, respectively, and thus are essentially the same or smaller than the 4.65 mm dimension of the area viewed by a single pixel in that part of the FOV.

  • Modeled images of field-aligned structures in the top half of the FOV show mostly radial elongation with radial FWHM of 4 to 9 mm (in the “default” orientation).

  • Combining the optical resolution with the modeled smearing yields values of “instrumental resolution” in the range 4.6 to 5.5 mm in the poloidal dimension and 6 to 10 mm in the radial dimension (in the “default” orientation). The resolutions are better closer to the nozzle, i.e., further out into the SOL.

  • From the modeling, a field-aligned structure piercing the “FOV registration plane” at [Rmaj, Z] is imaged with a brightness centroid within the pixel viewing that [Rmaj, Z] point in the “FOV registration plane.” In other words, the [Rmaj, Z] location assigned to a feature in the images will have pierced that viewed region at the toroidal angle of the “FOV registration plane.”

We now provide an example illustrating the measurement capabilities of the GPI system for characterizing 2D dynamics of plasma fluctuations in the SOL island region. We investigate a quasi-coherent mode (QCM) present in some W7-X discharges.33,38 In this instance, the OP 2.1 shot/program is 20221214.013. We analyze images from a GPI puff at 7.0 s when PECH = 2.0 MW, nel = 4.5 × 1019 m−2, and Ip = 2.5 kA. A QCM with a center frequency of ≈18 kHz is present in some pixel views but definitely absent in others. The spatial arrangement of the magnitudes of the signals’ power spectral densities (PSDs) in the mode is shown in Fig. 23(a). The mode is present in two radially separated regions. The W7-X magnetic configuration is “standard,” so the view of the SOL and island is similar to the one shown in Fig. 3, but the radial width of the island is significantly larger in this specific case. The island location and size are illustrated in Fig. 23(a), where the white solid circle indicates the location of the island O-point, and the white dashed lines show the contour for a 2 km connection length. It is a plausible explanation that the island plays a key role in the mode location. Examining the dynamics of fluctuations both visually in the videos and analytically, we find the motion is very dominantly in the poloidal direction, which allows for accurate evaluations of poloidal wavenumbers and phase velocities using Fourier analysis in the poloidal dimension and time. Frequency-normalized (kpol|f) Fourier spectrograms, two of which are shown in Figs. 23(b) and 23(c), provide a wealth of information. There is a layer of strong radial shear in the poloidal phase velocities, with the propagation direction changing sign from downward (in the two innermost columns of views) to upward (in the radially outboard columns of views where the QCM reappears). Multiple shear layers within the radial extent of the GPI FOV are not uncommon. Focusing on the QCM, we find it propagating poloidally in the innermost column of views (ρ = rviewrLCFS = 16.4 mm) at 2.6 km/s (in the lab frame), while in the ρ = 43.8 mm column of views, the phase velocity is 3.5 km/s. The poloidal wavenumbers of the mode are kpol ≈ ± 0.4 cm−1, i.e., approximately GPI’s minimum resolvable kpol. The mode is moving poloidally along with the broadband fluctuations. The frequency-FWHM of the mode, ΔFWHM, is 7 kHz; thus, ΔFWHM/f0 ≈ 0.4. Finally, we note that there is another QCM present in these data—the ∼2 kHz mode present in most W7-X plasmas.40 It is apparent as the strong very low-frequency feature at kpol ≈ 0 in spectrograms of Fig. 23. This brief analysis illustrates some of the phenomena being investigated by GPI on W7-X.

FIG. 23.

(a) The spatial distribution of the 18 kHz QCM. It is present on the innermost column of views (ρ = 16.4 mm) and on the columns of views with ρ between 36 and 46 mm. The color-coded “QCM magnitude” is proportional to the integral over frequency of the mode’s PSD. (b) A frequency-normalized (kpol|f) Fourier spectrogram of the fluctuations in the ρ = 16.4 mm column of views. (c) Same as in (b) but in the ρ = 43.8 mm column of views.

FIG. 23.

(a) The spatial distribution of the 18 kHz QCM. It is present on the innermost column of views (ρ = 16.4 mm) and on the columns of views with ρ between 36 and 46 mm. The color-coded “QCM magnitude” is proportional to the integral over frequency of the mode’s PSD. (b) A frequency-normalized (kpol|f) Fourier spectrogram of the fluctuations in the ρ = 16.4 mm column of views. (c) Same as in (b) but in the ρ = 43.8 mm column of views.

Close modal

In this report, we describe the considerations that went into the design of a new GPI system for the W7-X stellarator. We also describe in detail the different components that make up the W7-X GPI system and, finally, the performance of the realized system. Most important is the result that the system operated reliably during the OP 2.1 experimental campaign, acquiring excellent images for a large fraction of W7-X OP 2.1 plasmas. Detailed analyses of the images will provide an assortment of new and valuable scientific information. We reiterate the key performance results:

  • Excellent SNR in the 2D (radial and poloidal) images acquired at 2 Msamples/s.

  • Minimal perturbation to the W7-X plasma by the gas puff.

  • Spatial coverage over an 40×75 mm FOV with a spatial resolution of ≈5 mm poloidally and ≈8 mm radially.

  • A FOV that is typically in the SOL and includes the O-point of an outboard magnetic island in some W7-X configurations and an X-point in another configuration.

It is also important to emphasize some of the unique features of this GPI system:

  • A “converging–diverging” nozzle design and fabrication that provides significant collimation of the gas cloud and Mach 1 flow at the nozzle “throat.”

  • Confirmed collimation of the gas cloud that facilitated placing the nozzle at a safe 110 mm from the center of the FOV. To the best of our knowledge, this separation is the largest of any GPI system to date.

  • A re-entrant tube of ∼2 m length holding a shutter, a turning mirror, a vacuum window, and collection optics that traverses the radial extent of W7-X’s superconducting field coils and their cryostat.

  • A realized design that features robustness to heat loading on plasma-facing components in a long-pulse/steady-state fusion-relevant device. This includes water cooling of the re-entrant tube. Also implemented are a pair of flexible Cu straps that conduct heat from the movable shutter to a water-cooled “cold foot”. We cite the flexible straps as a possible solution for similar future situations on fusion devices.

  • Direct imaging of the light collected via the turning mirror and lenses onto the fast camera sensor. By building in the capability to rotate the camera about the optical axis of the collection optics, we have registered the locations of two fields-of-view and can switch between them as desired. The “default” FOV has the long dimension of the view oriented poloidally, while the long dimension of the “rotated” FOV is radial, i.e., normal to the LCFS.

The authors thank Dr. Timo Schröder for crucial contributions to the data acquisition process, the entire W7-X staff at IPP Greifswald, and engineer Rui Vieira at the MIT-PSFC. This work was supported, in part, by the US Department of Energy, Fusion Energy Sciences, Award No. DE-SC0014251. This work has been carried out within the framework of the EUROfusion Consortium, funded by the European Union via the Euratom Research and Training Program (Grant Agreement No. 101052200—EUROfusion). Views and opinions expressed are, however, those of the author(s) only and do not necessarily reflect those of the European Union or the European Commission. Neither the European Union nor the European Commission can be held responsible for them.

The authors have no conflicts to disclose.

J. L. Terry: Conceptualization (lead); Data curation (equal); Formal analysis (equal); Funding acquisition (lead); Investigation (equal); Methodology (lead); Project administration (lead); Resources (lead); Software (equal); Supervision (lead); Validation (equal); Visualization (equal); Writing – original draft (lead); Writing – review & editing (lead). A. von Stechow: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – review & editing (equal). S. G. Baek: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Software (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – review & editing (equal). S. B. Ballinger: Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Software (equal); Validation (equal); Visualization (equal); Writing – review & editing (equal). O. Grulke: Conceptualization (equal); Funding acquisition (equal); Project administration (equal); Resources (equal); Supervision (equal); Writing – review & editing (equal). C. von Sehren: Formal analysis (equal); Investigation (equal); Methodology (equal). R. Laube: Formal analysis (equal); Investigation (equal); Methodology (equal). C. Killer: Data curation (supporting); Investigation (supporting); Methodology (supporting); Validation (supporting); Writing – review & editing (equal). F. Scharmer: Formal analysis (supporting); Investigation (supporting); Methodology (supporting); Writing – review & editing (equal). K.J. Brunner: Data curation (supporting); Investigation (supporting). J. Knauer: Data curation (supporting); Investigation (supporting). S. Bois: Data curation (supporting); Formal analysis (supporting); Investigation (supporting); Validation (supporting); Writing – review & editing (equal).

The data that support the findings in this study are openly available in Harvard Dataverse at https://doi.org/10.7910/DVN/QAUPRJ, Ref. 70.

1.
S. F.
Paul
and
R. J.
Fonck
, “
Neutral beam emission spectroscopy diagnostic for measurement of density fluctuations on the TFTR tokamak
,”
Rev. Sci. Instrum.
61
,
3496
3500
(
1990
).
2.
G.
McKee
,
R.
Ashley
,
R.
Durst
,
R.
Fonck
,
M.
Jakubowski
,
K.
Tritz
,
K.
Burrell
,
C.
Greenfield
, and
J.
Robinson
, “
The beam emission spectroscopy diagnostic on the DIII-D tokamak
,”
Rev. Sci. Instrum.
70
,
913
916
(
1999
).
3.
G.
Cima
,
B.
Deng
,
C. W.
Dormer
,
W. R.
Geck
,
R. P.
Hsia
,
C.
Liang
,
F.
Jiang
,
N. C.
Luhmann
, Jr.
, and
D.
Brower
, “
ECE imaging on TEXT-U
,”
Fusion Eng. Des.
34–35
,
515
518
(
1997
).
4.
B. H.
Deng
,
C. W.
Domier
,
N. C.
Luhmann
, Jr.
,
D. L.
Brower
,
A. J.
Donné
,
T.
Oyevaar
, and
M. J.
van de Pol
, “
Electron cyclotron emission imaging diagnostic of Te profiles and fluctuations
,”
Phys. Plasmas
8
,
2163
2169
(
2001
).
5.
B.
Tobias
,
C. W.
Domier
,
T.
Liang
,
X.
Kong
,
L.
Yu
,
G. S.
Yun
,
H. K.
Park
,
I. G.
Classen
,
J. E.
Boom
,
A. J.
Donné
et al, “
Commissioning of electron cyclotron emission imaging instrument on the DIII-D tokamak and first data
,”
Rev. Sci. Instrum.
81
,
10D928
(
2010
).
6.
E.
Mazzucato
, “
Microwave imaging reflectometry for the visualization of turbulence in tokamaks
,”
Nucl. Fusion
41
,
203
(
2001
).
7.
T.
Munsat
,
E.
Mazzucato
,
H.
Park
,
B. H.
Deng
,
C. W.
Domier
,
N. C.
Luhmann
, Jr.
,
J.
Wang
,
Z. G.
Xia
,
A. J.
Donné
, and
M.
van de Pol
, “
Microwave imaging reflectometer for TEXTOR (invited)
,”
Rev. Sci. Instrum.
74
,
1426
1432
(
2003
).
8.
J. L.
Terry
,
R.
Maqueda
,
C. S.
Pitcher
,
S. J.
Zweben
,
B.
LaBombard
,
E. S.
Marmar
,
A. Y.
Pigarov
, and
G.
Wurden
, “
Visible imaging of turbulence in the sol of the Alcator C-Mod tokamak
,”
J. Nucl. Mater.
290–293
,
757
762
(
2001
).
9.
R. J.
Maqueda
,
G. A.
Wurden
,
S.
Zweben
,
L.
Roquemore
,
H.
Kugel
,
D.
Johnson
,
S.
Kaye
,
S.
Sabbagh
, and
R.
Maingi
, “
Edge turbulence measurements in NSTX by gas puff imaging
,”
Rev. Sci. Instrum.
72
,
931
934
(
2001
).
10.
S. J.
Zweben
,
R. J.
Maqueda
,
D. P.
Stotler
,
A.
Keesee
,
J.
Boedo
,
C. E.
Bush
,
S. M.
Kaye
,
B.
LeBlanc
,
J. L.
Lowrance
,
V. J.
Mastrocola
,
R.
Maingi
,
N.
Nishino
,
G.
Renda
,
D. W.
Swain
,
J. B.
Wilgen
, and
NSTX Team
, “
High-speed imaging of edge turbulence in NSTX
,”
Nucl. Fusion
44
,
134
153
(
2004
).
11.
I.
Cziegler
,
J. L.
Terry
,
J. W.
Hughes
, and
B.
LaBombard
, “
Experimental studies of edge turbulence and confinement in Alcator C-Mod
,”
Phys. Plasmas
17
(
5
),
056120
(
2010
).
12.
G.
Fuchert
,
G.
Birkenmeier
,
D.
Carralero
,
T.
Lunt
,
P.
Manz
,
H. W.
Müller
,
B.
Nold
,
M.
Ramisch
,
V.
Rohde
,
U.
Stroth
, and
ASDEX Upgrade Team
, “
Blob properties in L- and H-mode from gas-puff imaging in ASDEX upgrade
,”
Plasma Phys. Controlled Fusion
56
,
125001
(
2014
).
13.
I.
Shesterikov
,
Y.
Xu
,
M.
Berte
,
P.
Dumortier
,
M.
Van Schoor
,
M.
Vergote
,
B.
Schweer
, and
G.
Van Oost
, “
Development of the gas-puff imaging diagnostic in the TEXTOR tokamak
,”
Rev. Sci. Instrum.
84
,
053501
(
2013
).
14.
S. C.
Liu
,
L. M.
Shao
,
S. J.
Zweben
,
G. S.
Xu
,
H. Y.
Guo
,
B.
Cao
,
H. Q.
Wang
,
L.
Wang
,
N.
Yan
,
S. B.
Xia
et al, “
New dual gas puff imaging system with up-down symmetry on experimental advanced superconducting tokamak
,”
Rev. Sci. Instrum.
83
,
123506
(
2012
).
15.
N.
Offeddu
,
C.
Wüthrich
,
W.
Han
,
C.
Theiler
,
T.
Golfinopoulos
,
J. L.
Terry
,
E.
Marmar
,
C.
Galperti
,
Y.
Andrebe
,
B. P.
Duval
et al, “
Gas puff imaging on the TCV tokamak
,”
Rev. Sci. Instrum.
93
,
123504
(
2022
).
16.
M.
Agostini
,
R.
Cavazzana
,
P.
Scarin
, and
G.
Serianni
, “
Operation of the gas-puff imaging diagnostic in the RFX-mod device
,”
Rev. Sci. Instrum.
77
,
10E513
(
2006
).
17.
S. J.
Zweben
,
J. L.
Terry
,
D. P.
Stotler
, and
R. J.
Maqueda
, “
Invited review article: Gas puff imaging diagnostics of edge plasma turbulence in magnetic fusion devices
,”
Rev. Sci. Instrum.
88
,
041101
(
2017
).
18.
T.
Klinger
,
T.
Andreeva
,
S.
Bozhenkov
,
C.
Brandt
,
R.
Burhenn
,
B.
Buttenschön
,
G.
Fuchert
,
B.
Geiger
,
O.
Grulke
,
H. P.
Laqua
et al, “
Overview of first Wendelstein 7-X high-performance operation
,”
Nucl. Fusion
59
,
112004
(
2019
).
19.
R. C.
Wolf
,
A.
Alonso
,
S.
Äkäslompolo
,
J.
Baldzuhn
,
M.
Beurskens
,
C. D.
Beidler
,
C.
Biedermann
,
H.-S.
Bosch
,
S.
Bozhenkov
,
R.
Brakel
et al, “
Performance of Wendelstein 7-X stellarator plasmas during the first divertor operation phase
,”
Phys. Plasmas
26
,
082504
(
2019
).
20.
O.
Grulke
,
C.
Albert
,
J. A.
Alcuson Belloso
,
P.
Aleynikov
,
K.
Aleynikova
,
A.
Alonso
,
G.
Anda
,
T.
Andreeva
,
M.
Arvanitou
,
E.
Ascasibar
et al, “
Overview of the first Wendelstein 7-X long pulse campaign with fully water-cooled plasma facing components
,”
Nucl. Fusion
64
,
112002
(
2024
).
21.
S. J.
Zweben
,
J. R.
Myra
,
W. M.
Davis
,
D. A.
D’Ippolito
,
T. K.
Gray
,
S. M.
Kaye
,
B. P.
LeBlanc
,
R. J.
Maqueda
,
D. A.
Russell
,
D. P.
Stotler
et al, “
Blob structure and motion in the edge and SOL of NSTX
,”
Plasma Phys. Controlled Fusion
58
,
044007
(
2016
).
22.
C.
Wüthrich
,
C.
Theiler
,
N.
Offeddu
,
D.
Galassi
,
D.
Oliveira
,
B.
Duval
,
O.
Février
,
T.
Golfinopoulos
,
W.
Han
,
E.
Marmar
,
J.
Terry
,
C.
Tsui
, and
TCV Team
, “
X-point and divertor filament dynamics from gas puff imaging on TCV
,”
Nucl. Fusion
62
,
106022
(
2022
).
23.
F.
Militello
and
J. T.
Omotani
, “
Scrape off layer profiles interpreted with filament dynamics
,”
Nucl. Fusion
56
,
104004
(
2016
).
24.
O. E.
Garcia
,
S. M.
Fritzner
,
R.
Kube
,
I.
Cziegler
,
B.
LaBombard
, and
J. L.
Terry
, “
Intermittent fluctuations in the Alcator C-Mod scrape-off layer
,”
Phys. Plasmas
20
,
055901
(
2013
).
25.
D. A.
D’Ippolito
,
J. R.
Myra
, and
S. J.
Zweben
, “
Convective transport by intermittent blob-filaments: Comparison of theory and experiment
,”
Phys. Plasmas
18
,
060501
(
2011
).
26.
A.
Diallo
,
S.
Banerjee
,
S. J.
Zweben
, and
T.
Stoltzfus-Dueck
, “
Energy exchange dynamics across L–H transitions in NSTX
,”
Nucl. Fusion
57
,
066050
(
2017
).
27.
G. S.
Xu
,
L. M.
Shao
,
S. C.
Liu
,
H. Q.
Wang
,
B. N.
Wan
,
H. Y.
Guo
,
P. H.
Diamond
,
G. R.
Tynan
,
M.
Xu
,
S. J.
Zweben
et al, “
Study of the L–I–H transition with a new dual gas puff imaging system in the east superconducting tokamak
,”
Nucl. Fusion
54
,
013007
(
2013
).
28.
M.
Lampert
,
A.
Diallo
,
J. R.
Myra
, and
S. J.
Zweben
, “
Dynamics of filaments during the edge-localized mode crash on NSTX
,”
Phys. Plasmas
28
,
022304
(
2021
).
29.
J. L.
Terry
,
N. P.
Basse
,
I.
Cziegler
,
M.
Greenwald
,
O.
Grulke
,
B.
LaBombard
,
S. J.
Zweben
,
E. M.
Edlund
,
J. W.
Hughes
,
L.
Lin
et al, “
Transport phenomena in the edge of Alcator C-Mod plasmas
,”
Nucl. Fusion
45
,
1321
(
2005
).
30.
I.
Cziegler
,
G. R.
Tynan
,
P. H.
Diamond
,
A. E.
Hubbard
,
J. W.
Hughes
,
J.
Irby
, and
J. L.
Terry
, “
Zonal flow production in the L–H transition in Alcator C-Mod
,”
Plasma Phys. Controlled Fusion
56
,
075013
(
2014
).
31.
I.
Cziegler
,
P. H.
Diamond
,
N.
Fedorczak
,
P.
Manz
,
G. R.
Tynan
,
M.
Xu
,
R. M.
Churchill
,
A. E.
Hubbard
,
B.
Lipschultz
,
J. M.
Sierchio
,
J. L.
Terry
, and
C.
Theiler
, “
Fluctuating zonal flows in the I-mode regime in Alcator C-Mod
,”
Phys. Plasmas
20
,
055904
(
2013
).
32.
S.
Banerjee
,
A.
Diallo
, and
S.
Zweben
, “
Observation of quasi-coherent edge fluctuations in ohmic plasmas on National Spherical Torus Experiment
,”
Phys. Plasmas
23
,
044502
(
2016
).
33.
S.
Zoletnik
,
G.
Anda
,
C.
Biedermann
,
A. D.
Carralero
,
G.
Cseh
,
D.
Dunai
,
C.
Killer
,
G.
Kocsis
,
A.
Krämer-Flecken
,
M.
Otte
,
B.
Shanahan
,
T.
Szepesi
,
M.
Vecsei
, and
L.
Zsuga
, “
Multi-diagnostic analysis of plasma filaments in the island divertor
,”
Plasma Phys. Controlled Fusion
62
,
014017
(
2019
).
34.
C.
Killer
,
B.
Shanahan
,
O.
Grulke
,
M.
Endler
,
K.
Hammond
, and
L.
Rudischhauser
, “
Plasma filaments in the scrape-off layer of Wendelstein 7-X
,”
Plasma Phys. Controlled Fusion
62
,
085003
(
2020
).
35.
P.
Huslage
,
G.
Birkenmeier
,
B.
Shanahan
, and
C.
Killer
, “
Coherence of turbulent structures with varying drive in stellarator edge plasmas
,”
Plasma Phys. Controlled Fusion
66
,
015006
(
2023
).
36.
C.
Killer
,
Y.
Narbutt
,
O.
Grulke
, and
W7-X Team
, “
Turbulent transport in the scrape-off layer of Wendelstein 7-X
,”
Nucl. Fusion
61
,
096038
(
2021
).
37.
S. C.
Liu
,
Y.
Liang
,
P.
Drews
,
A.
Krämer-Flecken
,
X.
Han
,
D.
Nicolai
,
G.
Satheeswaran
,
N. C.
Wang
,
J. Q.
Cai
,
A.
Charl
et al, “
Observations of the effects of magnetic topology on the sol characteristics of an electromagnetic coherent mode in the first experimental campaign of W7-X
,”
Nucl. Fusion
58
,
046002
(
2018
).
38.
X.
Han
,
A.
Krämer-Flecken
,
T.
Windisch
,
M.
Hirsch
,
G.
Fuchert
,
J.
Geiger
,
O.
Grulke
,
S.
Liu
, and
K.
Rahbarnia
, “
Experimental characterization of a quasi-coherent turbulent structure in the edge plasma in Wendelstein 7-X
,”
Nucl. Fusion
60
,
016011
(
2019
).
39.
A.
Krämer-Flecken
,
X.
Han
,
M.
Otte
,
G.
Anda
,
S. A.
Bozhenkov
,
D.
Dunai
,
G.
Fuchert
,
J.
Geiger
,
O.
Grulke
,
E.
Pasch
et al, “
Investigation of turbulence rotation in the SOL and plasma edge of W7-X for different magnetic configurations
,”
Plasma Sci. Technol.
22
,
064004
(
2020
).
40.
S. B.
Ballinger
,
J. L.
Terry
,
S. G.
Baek
,
M.
Beurskens
,
K. J.
Brunner
,
G.
Fuchert
,
J.
Knauer
,
C.
Killer
,
E.
Pasch
,
K.
Rahbarnia
et al, “
Dynamics and dependencies of the configuration-dependent 1–2 kHz fluctuation in W7-X
,”
Nucl. Mater. Energy
27
,
100967
(
2021
).
41.
A.
Mathews
,
J. L.
Terry
,
S. G.
Baek
,
J. W.
Hughes
,
A. Q.
Kuang
,
B.
LaBombard
,
M. A.
Miller
,
D.
Stotler
,
D.
Reiter
,
W.
Zholobenko
, and
M.
Goto
, “
Deep modeling of plasma and neutral fluctuations from gas puff turbulence imaging
,”
Rev. Sci. Instrum.
93
,
063504
(
2022
).
42.
Y.
Gao
,
J.
Geiger
,
M. W.
Jakubowski
,
P.
Drewelow
,
M.
Endler
,
K.
Rahbarnia
,
S.
Bozhenkov
,
M.
Otte
,
Y.
Suzuki
,
Y.
Feng
,
H.
Niemann
,
F.
Pisano
,
A.
Ali
,
A.
Puig Sitjes
,
M.
Zanini
,
H.
Laqua
,
T.
Stange
,
S.
Marsen
,
T.
Szepesi
,
D.
Zhang
,
C.
Killer
,
K.
Hammond
,
S.
Lazerson
,
B.
Cannas
,
H.
Thomsen
,
T.
Andreeva
,
U.
Neuner
,
J.
Schilling
,
A.
Knieps
,
M.
Rack
,
Y.
Liang
, and
W7-X Team
, “
Effects of toroidal plasma current on divertor power depositions on Wendelstein 7-X
,”
Nucl. Fusion
59
,
106015
(
2019
).
43.
D.
Dunai
,
S.
Zoletnik
,
J.
Sárközi
, and
A. R.
Field
, “
Avalanche photodiode based detector for beam emission spectroscopy
,”
Rev. Sci. Instrum.
81
,
103503
(
2010
).
44.
See https://fusioninstruments.com (accessed 26 February 2024).
46.
See https://www.chroma.com/products (accessed 26 February 2024).
47.
P.
Drews
,
Y.
Liang
,
S.
Liu
,
A.
Krämer-Flecken
,
O.
Neubauer
,
J.
Geiger
,
M.
Rack
,
D.
Nicolai
,
O.
Grulke
,
C.
Killer
et al, “
Measurement of the plasma edge profiles using the combined probe on W7-X
,”
Nucl. Fusion
57
,
126020
(
2017
).
48.
C.
Killer
,
O.
Grulke
,
P.
Drews
,
Y.
Gao
,
M.
Jakubowski
,
A.
Knieps
,
D.
Nicolai
,
H.
Niemann
,
A. P.
Sitjes
,
G.
Satheeswaran
, and
W7-X Team
, “
Characterization of the W7-X scrape-off layer using reciprocating probes
,”
Nucl. Fusion
59
,
086013
(
2019
).
49.
See https://www.vacom.net/de/vpcf100uvq-l-nm.html (accessed 26 February 2024).
53.
V. A.
Soukhanovskii
,
H. W.
Kugel
,
R.
Kaita
,
R.
Majeski
, and
A. L.
Roquemore
, “
Supersonic gas injector for fueling and diagnostic applications on the National Spherical Torus Experiment
,”
Rev. Sci. Instrum.
75
,
4320
4323
(
2004
).
54.
S.
Collis
,
J.
Howard
,
B. D.
Blackwell
,
P.
Carlsson
,
M.
Abelsson
, and
B.
Powell
, “
A supersonic gas injection system for fuelling and probing fusion plasmas
,”
Plasma Sources Sci. Technol.
15
,
797
(
2006
).
55.
M.
Griener
,
O.
Schmitz
,
K.
Bald
,
D.
Bösser
,
M.
Cavedon
,
P.
De Marné
,
T.
Eich
,
G.
Fuchert
,
A.
Herrmann
,
A.
Kappatou
,
T.
Lunt
,
V.
Rohde
,
B.
Schweer
,
M.
Sochor
,
U.
Stroth
,
A.
Terra
,
E.
Wolfrum
, and
A. U.
Team
, “
Fast piezoelectric valve offering controlled gas injection in magnetically confined fusion plasmas for diagnostic and fuelling purposes
,”
Rev. Sci. Instrum.
88
,
033509
(
2017
).
56.
J. D.
Anderson
,
Modern Compressible Flow: With Historical Perspective
(
McGraw-Hill
,
New York
,
1990
), Vol.
12
, Chap. V.
57.
O.
Schmitz
,
Y.
Feng
,
M.
Jakubowski
,
R.
König
,
M.
Krychowiak
,
M.
Otte
,
F.
Reimold
,
T.
Barbui
,
C.
Biedermann
,
S.
Bozhenkov
,
S.
Brezinsek
,
B.
Buttenschön
,
K.
Brunner
,
P.
Drewelow
,
F.
Effenberg
,
E.
Flom
,
H.
Frerichs
,
O.
Ford
,
G.
Fuchert
,
Y.
Gao
,
D.
Gradic
,
O.
Grulke
,
K. C.
Hammond
,
U.
Hergenhahn
,
U.
Höfel
,
J.
Knauer
,
P.
Kornejew
,
T.
Kremeyer
,
H.
Niemann
,
E.
Pasch
,
A.
Pavone
,
V.
Perseo
,
L.
Rudischhauser
,
G.
Schlisio
,
T.
Sunn Pedersen
,
U.
Wenzel
,
V.
Winters
,
G.
Wurden
,
D.
Zhang
, and
W7-X Team
, “
Stable heat and particle flux detachment with efficient particle exhaust in the island divertor of Wendelstein 7-X
,”
Nucl. Fusion
61
,
016026
(
2020
).
58.
T.
Kremeyer
,
R.
König
,
S.
Brezinsek
,
O.
Schmitz
,
Y.
Feng
,
V.
Winters
,
L.
Rudischhauser
,
B.
Buttenschön
,
K.
Brunner
,
P.
Drewelow
,
E.
Flom
,
G.
Fuchert
,
Y.
Gao
,
J.
Geiger
,
M.
Jakubowski
,
C.
Killer
,
J.
Knauer
,
M.
Krychowiak
,
S.
Lazerson
,
F.
Reimold
,
G.
Schlisio
,
H.
Viebke
, and
W7-X Team
, “
Analysis of hydrogen fueling, recycling, and confinement at Wendelstein 7-X via a single-reservoir particle balance
,”
Nucl. Fusion
62
,
036023
(
2022
).
59.
W. F.
Louisos
,
A. A.
Alexeenko
,
D. L.
Hitt
, and
A.
Zilic
, “
Design considerations for supersonic micronozzles
,”
Int. J. Manuf. Res.
3
,
80
113
(
2008
).
61.
H.
Bohlin
,
A.
Von Stechow
,
K.
Rahbarnia
,
O.
Grulke
, and
T.
Klinger
, “
Vineta II: A linear magnetic reconnection experiment
,”
Rev. Sci. Instrum.
85
,
023501
(
2014
).
62.
T.
Schroeder
, “
Dynamic structures in expanding plasmas
,” Ph.D. thesis,
Ernst-Moritz-Arndt Universität
,
2014
.
63.
See https://thermotive.com/standard_series.html (accessed 26 February 2024).
64.
D.
Stotler
and
C.
Karney
, “
Neutral gas transport modeling with DEGAS 2
,”
Contrib. Plasma Phys.
34
,
392
397
(
1994
).
65.
B.
LaBombard
, “
KN1D: A 1-D space, 2-D velocity, kinetic transport algorithm for atomic and molecular hydrogen in an ionizing plasma
,”
PSFC Report RR-01-3
,
2001
, https://www1.psfc.mit.edu/∼Labombard/KN1D_Source_Info.html.
66.
D.
Nicolai
,
V.
Borsuk
,
P.
Drews
,
O.
Grulke
,
K.
Hollfeld
,
T.
Krings
,
Y.
Liang
,
C.
Linsmeier
,
O.
Neubauer
,
G.
Satheeswaran
,
B.
Schweer
, and
G.
Offermanns
, “
A multi-purpose manipulator system for W7-X as user facility for plasma edge investigation
,”
Fusion Eng. Des.
123
,
960
964
(
2017
), part of Special Issue: Proceedings of the 29th Symposium on Fusion Technology (SOFT-29) Prague, Czech Republic, September 5–9, 2016.
67.
J. W.
Oosterbeek
et al, “
Microwave stray radiation measurement techniques
,” in
The 33rd Symposium on Fusion Technology (SOFT 2025)
,
Dublin
,
Sept, 2024
.Submitted.
68.
K. J.
Brunner
,
T.
Akiyama
,
M.
Hirsch
,
J.
Knauer
,
P.
Kornejew
,
B.
Kursinski
,
H.
Laqua
,
J.
Meineke
,
H.
Trimiño Mora
, and
R. C.
Wolf
, “
Real-time dispersion interferometry for density feedback in fusion devices
,”
J. Instrum.
13
,
P09002
(
2018
).
69.

A 400-µs raw-data image-sequence (24ja005_VIDEO_230323049_2.100300_2.100700s_raw_pixels.mp4) from the puff shown in Fig. 20 can be downloaded from https://doi.org/10.7910/DVN/QAUPRJ. A video of the same sequence, but with a 50-µs running minimum subtracted from each pixel signal (24ja005_VIDEO_230323049_2.100300_2.100700s_50us_minsub.mp4), can be downloaded from the same URL.

70.
J. L.
Terry
,
A.
von Stechow
,
S. G.
Baek
,
S. B.
Ballinger
,
O.
Grulke
,
C.
von Sehren
,
R.
Laube
,
C.
Killer
,
F.
Scharmer
,
K. J.
Brunner
,
J.
Knauer
,
S.
Bois
, and The W7-X Team (2024). “
The Realization of a Gas Puff Imaging System on the Wendelstein 7-X Stellarator
,”
Harvard Dataverse
, V1.
Published open access through an agreement with Massachusetts Institute of Technology