Turbulent transport reduction induced by transition on radial electric field shear and curvature through amplitude and cross-phase in torus plasma.

Spatiotemporal evolutions of radial electric field and turbulence are measured simultaneously in the H-mode transition, which is a prototypical example of turbulence structure formation in high-temperature plasmas. In the dynamical phase where transport barrier is established abruptly, the time-space-frequency-resolved turbulent particle flux is obtained. Here we report the validation of the mechanism of transport barrier formation quantitatively. It is found that the particle flux is suppressed predominantly by reducing density fluctuation amplitude and cross phase between density fluctuation and potential fluctuation. Both radial electric field shear and curvature are responsible for the amplitude suppression as was predicted by theory. Turbulence amplitude reduction immediately responds to the growth of the radial electric field non-uniformity and saturates, while cross phase continuously approaches zero.

Introduction

Structure of flows and turbulence in non-equilibrium plasmas has attracted much attention because of its great impact on the entire media dynamics. One of the prototypical examples can be seen in the solar tachocline[1], across which the transition from the differential rotation in the solar convective zone to the rigid rotation in the radiative interior occurs. The strong shear flow within the thin layer is believed to amplify the magnetic field as the solar dynamo. Another example is the edge transport barrier (ETB) in toroidal fusion plasmas[2], which is spontaneously formed in the Low-confinement mode to High-confinement mode transition (L-H transition)[3]. Because of its rich nature of nonlinear dynamics as well as its promising capability for achieving good plasma performance toward the fusion reactor, much attention has been paid for clarifying the underlying physics for decades. Numerous studies have shown essential roles of an edge-localized poloidal E × B flow structure on confinement improvement[4-7]. However, definitive conclusion regarding what aspect of the E × B flow suppresses the turbulent transport is still under debate[8]. In order to obtain an understanding based on first principles, interplay between the flow structure and the turbulence must be diagnosed with high spatiotemporal resolutions, which remains challenging.

In modeling works, roles of inhomogeneous E × B flow on the turbulence suppression are classified in two elements: shear and curvature. Effects of both non-uniformities on mode instability were discussed[9,10]. Nonlinear saturation level of turbulence was studied using a statistical approach, showing an essential role of both shear and curvature of the E × B flow[11]. In a newly developed model[12] described in “Method”, shear and curvature were simultaneously treated in a single formulae that provided turbulence reduction rate. In that model, responsible physics of shear and curvature were considered as the E × B flow shear decorrelation of turbulence eddy[13-15] and the modulational coupling causing energy transfer from turbulence to macroscopic E × B flow structure[2,12,16], respectively. Experimental examinations of those models have been performed[8,17-25]. Focusing on detailed physics of transport suppression, direct fluctuation measurement by use of electrostatic probes has been promoted. It was found that the turbulent particle flux is reduced not only by the density fluctuation amplitude suppression but also by the cross phase alternation between the density fluctuation and the potential fluctuation[19-22]. Although individual elements regarding the turbulent transport suppression by the inhomogeneous E × B flow have been raised, i.e., shear and curvature of E × B flow and amplitude and cross phase of fluctuations, the mutual relation of these elements remain unclear.

In this paper, we investigate response of turbulent particle flux on radial electric field non-uniformity by analyzing data from a heavy ion beam probe (HIBP). In particular, interrelations among shear and curvature of radial electric field and amplitude and cross phase of fluctuations are shown for the first time. Electron density and electrostatic potential measured with high spatiotemporal resolutions allow us to perform time-space-frequency-decomposition to the fluctuation induced particle flux. The particle flux is suppressed predominantly by reducing both the density fluctuation amplitude and the cross phase between density fluctuation and potential fluctuation. Both curvature and shear are responsible for transport reduction. Turbulence amplitude reduction immediately responds to the growth of the radial electric field non-uniformity and saturates, while cross phase continuously approaches zero. The time scales of these dynamics have an order of magnitude difference. As a result, turbulent transport reduction occurs with two different time scales.

Results

Time-space-frequency-resolved turbulence spectrum

Figure 1 (a) shows time evolutions of emission from divertor region, inverse density gradient length , where is the HIBP secondary beam current that reflects the local electron density, and radial electric field at cm, in which the bottom of the E -well structure emerges in H-mode. The global profile changes are represented by low-pass filtered signals with a cut-off frequency of 2 kHz. The cut-off frequency is chosen to eliminate dynamics of the limit-cycle oscillation at 4.5 kHz[26]. As clearly indicated, transition in and E occurs twice[27]. The first transition is triggered by the reach of a sawtooth crash-induced heat pulse at the edge region. The onset time is defined using the emission signal and a soft-x-ray signal. After the first transition within several hundred microseconds, the plasma experiences a short period quasi-stationary state of ∼2 ms, which we call the “MH-mode” (meta-stable H-mode). Then, the plasma reaches the final H-mode by completing the transport barrier with the second transition within ∼1 ms. A dedicated study revealed that the radial current induced by the neoclassical bulk viscosity[4,5], and the ion loss-cone loss[5] plays an important role for building E structure during the first transition[28]. Reynolds stress driven E was found to be much smaller.

Figure 1

Time traces of (a) emission, inverse density gradient length and radial electric field, wavelet power spectra of (b) relative density fluctuation and (c) potential fluctuation, (d) and (e) squared cross coherence and cross phase between density fluctuation and potential fluctuation, respectively, and (f) particle flux normalized by mean density at cm. Black curves on contour plots show the expected frequency Doppler shift in the laboratory frame.

Contour plots in Fig. 1 show time-frequency-resolved wavelet power spectra of (b) relative density fluctuation and of (c) potential fluctuation , (d) and (e) squared cross coherence and cross phase between and denoted as and , respectively, and (f) particle flux normalized by mean density . Here, wavelet transform is performed with a time interval of 50 μs, after which four sequential time samples are ensemble averaged. Squared cross coherence and cross phase are defined as and with the power spectrum of and and cross spectrum between them denoted as , , and , respectively. The particle flux is defined as [29]. The poloidal wavenumber of turbulence k is determined to be in kHz in a dedicated discharge, in which the sample volumes of HIBP are aligned in a magnetic surface at cm[30]. The poloidal wavenumber in kHz is much smaller compared to that in kHz so that the particle flux driven by those low frequency components play a minor role for confinement. In this paper, the poloidal direction is defined as the electron diamagnetic direction. At the central frequency kHz, phase velocity is approximately equal to the electron diamagnetic velocity in the plasma frame.

Properties of turbulence and turbulent particle flux change drastically through the transitions. In L-mode, a broadband turbulence spectrum in kHz arises both in and spectra. Cross phase is slightly positive in average and turbulence-driven particle flux is directed outward. During the L-MH transition, the central frequency of the turbulence spectrum rises sharply up to ∼100 kHz. Black curves superimposed on Fig. 1 (b)-(f) show the expected Doppler frequency shift for turbulence component at kHz in L-mode, where is assumed to remain unchanged through L-MH transition. Below the curve, amplitude of fluctuation is significantly reduced, in which fluctuation wavenumber is found to be substantially lower than that of the high frequency turbulence component[30]. In contrast, the higher frequency component suffers a moderate amplitude reduction. Squared cross coherence remains close to unity but cross phase decreases toward zero. A considerable reduction of the particle flux occurs, which is responsible for density pedestal formation. Further deepening of the E -well occurs in the final H-mode. Doppler shifted frequency exceeds measurable frequency band. Therefore we avoid investigating the final H-mode period in this paper.

Figure 2 shows radial profile of time-averaged turbulence spectrum. Radial profile of E is superimposed as black curves. In L-mode, turbulence has a broad spectrum both in frequency and in space. This broadband turbulence generates a widely distributed outward directed particle flux that is related to the confinement degradation. After the L-MH transition, the E -well structure with the full width at half maximum of ∼1.4 cm emerges at cm. In MH-mode, Doppler frequency shift occurs up to ∼100 kHz with a moderate amplitude reduction at the bottom of the E -well. The Doppler frequency shift is only visible within the E -well. While, the turbulence is strongly stabilized in the outer shear region ( cm), the inner shear region ( cm), and further inside (). As a result, the outward particle flux is entirely suppressed outside the E -well. Even inside the E -well, the particle flux is clearly reduced mainly because of change in the cross phase as discussed above.

Figure 2

Radial profiles of time averaged wavelet spectra of (a and c) relative density fluctuation and (b and d) particle flux normalized by mean density in L-mode and in MH-mode, respectively. Black curves in (a) and (c) are radial profile of radial electric field.

Response of turbulent particle flux on radial electric field non-uniformity

On frequency integrated turbulence properties shown in Fig. 3, we quantitatively show response of turbulent particle flux on radial electric field non-uniformity. Radial profiles of fluctuation amplitude for and [Fig. 3(a) and (b)] are given by integrating the spectra in . In L-mode, shapes of both amplitude profiles are similar and normalized turbulence amplitudes are approximately in the same level, , where an equivalent order of electron temperature and ion temperature is assumed[7], and eV is obtained by a charge exchange spectroscopy. A small but finite phase difference between density fluctuation and potential fluctuation exists. The dominant turbulence source is considered as the resistive drift wave[26]. Two different elements of radial electric field non-uniformity, shear and curvature, are parameterized as non-dimensional factors Z 1 and Z 2, respectively. See “Method” for the definition. Radial profiles of Z 1 and Z 2 are shown in Fig. 3(f) and (g). Spatial derivative is taken for fifth order polynomial fit of the E profile. The model[12] predicts that the non-uniformity observed here is large enough to suppress turbulence amplitude, i.e., , where the turbulence scale factor is given as . Note that the criterion for is equivalent to the well-known shearing rate criterion[13], i.e., , if the turbulence decorrelation rate can be approximated by the drift wave frequency . The criterion for Z 2 corresponds to [16]. At the location where is held, a substantial fraction of the microscopic turbulence energy is transferred into the macroscopic radial electric field. If the transferred energy is large enough to maintain the radial electric field structure, the H-mode transition can be described by an integrated system of physics, which is referred to as “single-step” transition in Ref.[31]. In the preset case, however, the turbulence driven radial electric field accounts for only a small fraction of the total one[28] (referred to as “two-step” transition[31]).

Figure 3

Radial profiles of (a) relative density fluctuation amplitude, (b) potential fluctuation amplitude, (c) and (d) squared cross coherence and cross phase between density fluctuation and potential fluctuation, respectively, (e) radial electric field, (f) shear parameter , (g) curvature parameter , and (h) particle flux normalized by mean density. Arrows labeled “A” and “B” in (h) indicate radii in which detailed time traces are shown in Fig. 4.

Figure 4

Time traces of (a) radial electric field , (b) relative density fluctuation amplitude, (c) cross phase between density fluctuation and potential fluctuation, and (d) particle flux at cm. (e–h) Those at cm. Blue-solid and red-dashed horizontal lines in (b–d) and (f–h) are mean values in L-mode and MH-mode, respectively.

Comparing the profiles in L-mode and MH-mode, Z 1 seems to be more effective for reducing turbulence amplitude than Z 2. Note that around the radius of the Z 2 peak, the density pedestal appears in which linear energy input to the turbulence is enhanced. The profile has a pivot point at cm, inside or outside which the turbulence amplitude decreases or increases, respectively. A sign dependence of E shear for suppression possibly exists as discussed in Ref.[19]. A mild reduction in is also seen. Cross phase approaches zero in , while outside the region becomes negative. Non-adiabatic response of electrons on potential perturbation that gives birth to the outward particle flux is weakened. Figure 3(h) shows profile of the particle flux normalized by mean density. The value in L-mode is in a similar order to that reported in Ref.[32] in which the turbulence particle flux is identified as the dominant loss channel of plasma density. As a result of changes mainly in the density fluctuation amplitude and the cross phase, the particle flux is drastically reduced.

Detailed time evolution during the transition is shown in Fig. 4. Left and right columns are time traces for cm (#90048) and for cm (#90055), respectively, in which either Z 1 or Z 2 predominantly varies. Time scale of change in E non-uniformity is considered to be equivalent as that in E . For both cases, time scale of the transition is the order of 100 μs. After the transition, a slow draft of E further deepens the transport barrier in the MH-mode state. The relative density fluctuation amplitude immediately responds to change in E with the time scale of the order of 100 μs. The change in the amplitude saturates prior to E . Meanwhile, the change in the cross phase is much slower. The time scale of the change is 500 μs to 1 ms. These tendencies are common for both cases. The difference of time scale may suggest different underlying physics for density fluctuation amplitude suppression and cross phase modification. Note that the potential fluctuation amplitude is less sensitive to the change in the radial electric field in particular at the E -well location as shown in Fig. 3. Models taking into account the different responses in density fluctuation and potential fluctuation are highly desirable in future.

For the case of #90055, the relative density fluctuation amplitude, the cross phase, and the particle flux are plotted as a function of E in Fig. 5. The relative density fluctuation amplitude decreases as E grows when kV/m. Correspondingly, the particle flux decreases, suggesting that the paradigm of turbulence transport suppression through amplitude effectively works when the E non-uniformity is relatively small. Throughout the change in E the cross phase gradually approaches zero. With larger value of the E non-uniformity the change in the cross phase dominates change in the particle flux.

Figure 5

(a) Relative density fluctuation amplitude and particle flux and (b) cross phase plotted as a function of radial electric field . Color of symbols corresponds to time, which is shown by the color bar.

Discussion and summary

Even inside radii where radial electric field non-uniformity are small, turbulence amplitude in both relative density fluctuation and potential fluctuation are reduced, as observed in CCT[18]. One of the possible mechanisms that can bring this disparate space turbulence suppression is “turbulence spreading”[33]. Dynamic inward transmission of a turbulence packet was observed in a limit cycle oscillation event[26]. This observation can provide a clue to address a long-standing mystery, that is, the fast improvement of core confinement by edge transport barrier[34].

In conclusion, we investigated response of turbulent particle flux on radial electric field shear and curvature measured by heavy ion beam probe. Particle flux was mainly suppressed by reducing density fluctuation amplitude and cross phase between density fluctuation and potential fluctuation. Both radial electric field shear and curvature were responsible for transport suppression. Turbulence amplitude reduction immediately responded to the growth of the radial electric field non-uniformity and saturates, while cross phase continuously approached zero.

Methods

JFT-2M

The experiments were conducted in the JFT-2M tokamak, which has a major radius of 1.3 m and a minor radius of 0.3 m. Plasma was auxiliary heated by a co-injected neutral beam (NB) with the power of 750 kW, which is just above the L-H transition threshold. Line averaged electron density was before the L-H transition. An upper single-null divertor configuration was employed with the drift of ion directed toward the X-point. Other operation parameters were as follows: toroidal magnetic field at the magnetic axis of 1.17 T, safety factor at the flux surface enclosing 95% of the total poloidal flux, , of 2.9, and plasma current I p of 190 kA.

Heavy Ion Beam Probe (HIBP)

In order to diagnose electrostatic potential ϕ, heavy ion beam is injected from the top-side of the torus, which is ionized doubly inside the confined plasma. By analyzing the secondary beam energy, ϕ is given at four sample volumes () with a sampling rate of 1 μs[30]. Radial distance between each sample volume projected in the outer mid-plane is ~2.5 mm. With precise tuning of the HIBP measurement conditions, such as the primary beam energy, the toroidal magnetic field, and the incident angle of the beam, measurement positions can be scanned in an edge region () on a shot-to-shot basis. Relative secondary beam current fluctuation is regarded to be equivalent to relative electron density fluctuation , since beam attenuation effect is estimated to be negligibly small at the edge[35]. Gradient of quantities is defined by the finite difference of the neighboring sampling volumes.

Turbulence suppression model

In the newly developed theoretical model[12], amplitude of turbulence having a perpendicular wavenumber k is reduced by the shear factor Z 1 and the curvature factor Z 2 as , where I and I 0 are the reduced turbulence amplitude and the intrinsic turbulence amplitude without E effects, respectively. Shear factor and curvature factor are defined as indicators of significance of shear decorrelation and of modulational coupling, respectively, where is the ion gyro-radius, is the diamagnetic velocity, and prime is the radial derivative. In the present case, is considered to hold[7] so that the toroidal velocity correction term for Z 2 is neglected. Recent investigation in JT-60U[8] showed an important role of not only Z 1 but also of Z 2 by referring to temperature gradient profile as an indicator of confinement intensity.