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.

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][5][6][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][14][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][18][19][20][21][22][23][24][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][20][21][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  , where I HIBP is the HIBP secondary beam current that reflects the local electron density, and radial electric field E r at − ∼ − . r a 0 8 cm, in which the bottom of the E r -well structure emerges in H-mode. The global profile changes are represented by lowpass filtered signals with a cut-off frequency of 2 kHz. The cut-off frequency is chosen to eliminate dynamics of and E r 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 t LH is defined using the α D 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 r structure during the first transition 28 . Reynolds stress driven E r was found to be much smaller.
Contour plots in Fig. 1 show time-frequency-resolved wavelet power spectra of (b) relative density fluctuation n n /  and of (c) potential fluctuation φ  , (d) and (e) squared cross coherence and cross phase between n  and φ  denoted as n 2 γ φ and n α φ , respectively, and (f) particle flux normalized by mean density n / r Γ . 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 γ = | | , with the power spectrum of  n and  φ and cross spectrum between them denoted as P n , φ P , and P n,φ , respectively. The particle flux is defined as The poloidal wavenumber of turbulence k θ is determined to be ∼ . θ − k 0 75 cm 1 in > f 20 kHz in a dedicated discharge, in which the sample volumes of HIBP are aligned in a magnetic surface at r a 1 − ∼ − cm 30 . The poloidal wavenumber in f 20 < kHz is much smaller compared to that in f 20 > 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 ∼ f 50 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 < f 80 kHz arises both in  n n / and  φ spectra. Cross phase α φ n 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 show the expected Doppler frequency shift for turbulence component at f 20 = kHz in L-mode, where θ k 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 γ φ n 2 remains close to unity but cross phase α φ n decreases toward zero. A considerable reduction of the particle flux occurs, which is responsible for density pedestal formation. Further deepening of the E r -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 r 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 r -well structure with the full width at half maximum of ∼1.4 cm emerges at r a 0 8 − = − . cm. In MH-mode, Doppler frequency shift occurs up to ∼100 kHz with a moderate amplitude reduction at the bottom of the E r -well. The Doppler frequency shift is only visible within the E r -well. While, the turbulence is strongly stabilized in the outer shear region ( − > − . r a 0 3 cm), the inner shear region ). As a result, the outward particle flux is entirely suppressed outside the E r -well. Even inside the E r -well, the particle flux is clearly reduced mainly because of change in the cross phase n α ϕ as discussed above.
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 n n /  and φ  [ Fig. 3(a) and (b)] are given by integrating the spectra in ≤ ≤ f 20 kHz 110 kHz. In L-mode, shapes of both amplitude profiles are similar and normalized turbulence amplitudes are approximately in the same level, ~~ẽ T n n / / eφ , where an equivalent order of electron temperature and ion temperature ∼ T T e i is assumed 7 , and ∼ T 130 i 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 r profile. The model 12  ρ > − 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 . 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 Z 1 , (g) curvature parameter Z 2 , 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. 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 ).
Comparing the  n n / 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 r a 0 4 − ∼ − . cm, inside or outside which the turbulence amplitude decreases or increases, respectively. A sign dependence of E r shear for  φ suppression possibly exists as discussed in Ref. 19 . A mild reduction in γ φ n 2 is also seen. Cross phase α φ n approaches zero in − ≤ − ≤ − . r a 3cm 0 3cm, while outside the region n α φ 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 − ∼ − . r a 1 48 cm (#90048) and for − ∼ − . r a 0 90 cm (#90055), respectively, in which either Z 1 or Z 2 predominantly varies. Time scale of change in E r non-uniformity is considered to be equivalent as that in E r . For both cases, time scale of the E r transition is the order of 100 μs. After the transition, a slow draft of E r further deepens the transport barrier in the MH-mode state. The relative density fluctuation amplitude immediately responds to change in E r with the time scale of the order of 100 μs. The change in the amplitude saturates prior to E r . 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 r -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 r in Fig. 5. The relative density fluctuation amplitude decreases as E r grows when E 7 r > − kV/m. Correspondingly, the particle flux decreases, suggesting that the paradigm of turbulence transport suppression through amplitude effectively works when the E r non-uniformity is relatively small. Throughout the change in E r the cross phase gradually approaches zero. With larger value of the E r non-uniformity the change in the cross phase dominates change in the particle flux.

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.

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 . × − m 1 1 10 19 3 before the L-H transition. An upper single-null divertor configuration was employed with the ∇B drift of ion directed toward the X-point. Other operation parameters were as follows: toroidal magnetic field at the magnetic axis B t of 1.17 T, safety factor at the flux surface enclosing 95% of the total poloidal flux, q 95 , 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 ( × 6 mm 2 mm) 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 (− < − < r a 5 cm 0 cm) on a shot-to-shot basis. Relative secondary beam current fluctuation I I / HIBP HIBP  is regarded to be equivalent to relative electron density fluctuation  n n / , 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 ρ = + + are defined as indicators of significance of shear decorrelation and of modulational coupling, respectively, where ρ i is the ion gyro-radius, ≡ | ′| V T n enB / d is the diamagnetic velocity, and prime is the radial derivative. In the present case,  θ E V B r tor 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.