On how structures convey non-diffusive turbulence spreading

We report on comprehensive experimental studies of turbulence spreading in edge plasmas. These studies demonstrate the relation of turbulence spreading and entrainment to intermittent convective density fluctuation events or bursts (i.e. blobs and holes). The non-diffusive character of turbulence spreading is thus elucidated. The turbulence spreading velocity (or mean jet velocity) manifests a linear correlation with the skewness of density fluctuations, and increases with the auto-correlation time of density fluctuations. Turbulence spreading by positive density fluctuations is outward, while spreading by negative density fluctuations is inward. The degree of symmetry breaking between outward propagating blobs and inward propagating holes increases with the amplitude of density fluctuations. Thus, blob-hole asymmetry emerges as crucial to turbulence spreading. These results highlight the important role of intermittent convective events in conveying the spreading of turbulence, and constitute a fundamental challenge to existing diffusive models of spreading.

Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence.Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Edge plasmas are usually strongly turbulent, and manifest intermittency in linear plasma devices and toroidal fusion devices [1][2][3][4].Spreading of turbulence occurs because inhomogeneous turbulence tends to relax and entrain laminar or more weakly turbulent regions [5].For example, a wake formed downstream of a moving object expands due to the turbulence spreading from the fully turbulent core [6].Another example is the detection of strong temperature fluctuations inside a magnetic island, due to transient spreading events across it [7].Generally, turbulence spreading refers to the spatial propagation of turbulence intensity or energy due to nonlinear interactions [5].Spreading can decouple the intensity field from the local instability growth rate.The observed intensity flux caused by spreading can break Fick's law for the local flux-gradient relation (Γ = −D∇P) [8].
In this letter, we present comprehensive experimental studies of turbulence spreading dynamics at the tokamak plasma edge.We focus on the non-diffusive characteristics of turbulence spreading and their relation to intermittent convective density fluctuation events (blobs/holes).Turbulence spreading is significantly enhanced at high collisionality or low adiabaticity.Characterizing turbulence spreading as a Fickian diffusion yield an unphysically large or singular diffusivity χ I , which departs radically from the turbulent particle diffusivity.Turbulence spreading manifests non-diffusive characteristics and thus the usual models are dubious.The turbulence spreading speed (named 'mean jet velocity') correlates linearly with the skewness of density fluctuations.Spreading induced by positive density fluctuation events is outward toward larger radii, and spreading induced by negative density fluctuation events is inward toward smaller radii, as shown by the diagram in figure 1. Increasing symmetry breaking between outgoing blobs and incoming holes causes enhanced turbulence spreading.These results are essential to understand the basic physics of turbulence spreading and to build a foundation for the prediction and interpretation of turbulence transport behaviors.
The experiments were carried out in Ohmic hydrogen discharges of the J-TEXT tokamak with a limiter configuration [40].Its major radius and minor radius are R = 1.05 m and a = 0.255 m, respectively.Edge plasmas with different th,e /ν ei ω * .Here, ν ii is the ion collision rate, ε is the inverse aspect ratio, v th,e = √ K B T e /m e is the electron thermal speed, k ∥ ∼ 1/qR is the parallel wavenumber, ν ei is the electron-ion collision rate, and ω * is the electron diamagnetic frequency [41].Since turbulence spreading refers to the spatial propagation of turbulence intensity or energy, it is usually defined as the spatial cross-field flux of turbulence intensity or energy.For turbulence intensity as measured by density fluctuations, ⟨ ṽr ñ2 ⟩ /2 is the relevant flux.It appears in the evolution equation of turbulence intensity, i.e.
) is the relevant flux [16].For turbulence kinetic energy, ⟨ ṽr ṽ2 ⊥ ⟩ is the relevant flux.Kinetic energy flux is observed to be about two orders of magnitude smaller than internal energy flux in the experiments [15,16], and so is negligible here.For simplicity, the turbulence spreading flux is calculated as ⟨ ṽr ñ2 ⟩ /2.This representation is also widely used in theory and simulations [8,23].
⟨ ṽr ñ2 ⟩ /2 at the plasma edge is measured by a Langmuir probe array [16,37].Here, ṽr represents the radial velocity fluctuation and ñ stands for the density fluctuation.Electron density n e , temperature T e , plasma potential ϕ p and ion saturation current I sat are measured by the triple probe.Ion temperature T i is assumed to be approximately equal to T e .The mean E × B poloidal flow v θ is calculated as ∇ r ⟨ϕ p ⟩/B t .Here, ⟨•⟩ indicates a time average.The fluctuating density ñ is estimated from ion saturation current fluctuations, i.e.
Ĩsat ⟨n e ⟩/⟨I sat ⟩.The fluctuating radial velocity ṽr is inferred from the difference between two poloidally separated floating potential probes, i.e. −∇ θ Ṽf /B t .The sampling rate of probes is 2 MHz.A digital finite impulse response filter is used to obtain the fluctuations with 2 < f < 100 kHz.From the auto-spectra and the wavenumber-frequency spectra of potential fluctuations or density fluctuations, the fluctuations above 100 kHz are at least 2-3 orders of magnitude weaker than those below 100 kHz, and thus are negligible and filtered out.
Figure 2 show the radial profiles of plasma density, electron temperature, electron collision rate ν ei , turbulence spreading flux ⟨ ṽr ñ2 ⟩ /2, skewness of density fluctuations and E × B poloidal flow velocity, for discharges with different collisionality ν * and electron adiabaticity parameter α near the last closed flux surface (LCFS).Curves with red diamond, green cross, and blue inverted triangle represent edge plasmas with lower ν * and higher α, medium ν * and medium α, higher ν * and lower α, respectively.The black dotted line marks the position of the LCFS at r = 25.5 cm.We observe that the edge turbulence spreading flux increases as collisionality ν * increases or adiabaticity α decreases, as shown in figure 2(d).For the hydrodynamic regime with α < 1, density fluctuations ñ and potential fluctuations φ are not proportional, as in the adiabatic regime with α > 1.The electron response thus resembles that of a convective cell [42].The result shows that as collisions become more frequent and the hydrodynamic response is stronger, edge turbulence spreading grows significantly and becomes more prominent.
Turbulent diffusion of intensity has long been used as a working model of turbulence spreading.This model has not been tested critically.As the convective properties of edge turbulence increase for higher collisionality in the hydrodynamic regime, the diffusive description of turbulence spreading should be reconsidered.Actually, it has been hypothesized that turbulence kinetic energy spreads via a process independent of turbulent viscosity in the formation of a turbulent wake [6].The hope that the turbulence energy diffusivity was the same as the turbulent viscosity or particle diffusivity seems unlikely to be realized.Heuristically, we can write: where Γ residual is the residual, non-diffusive flux of turbulence intensity.Extracting both χ * I and Γ residual from data likely requires dynamics perturbation studies, rather like gas puffing modulation experiments to determine diffusivity D and convection V contributions to the particle flux.These would be challenging new experiments, well beyond the scope of this paper.Here, then, we can argue by reductio ad absurdum.So we then the calculated χ I is likely unphysically large (This is a presumption, of course).So that (4) Γ residual ̸ = 0, and the residual spreading flux must be retained, i.e. the spreading has a significant non-diffusive component.
The diffusion coefficients of a Fickian model for turbulence spreading flux are calculated and compared to the turbulent particle diffusivity, to test whether the use of diffusive model to describe turbulence spreading is justified or sensible.The turbulence spreading diffusivity is defined as , while the turbulent particle diffusivity is defined as D n = −⟨ñṽ⟩ r /∂ r ⟨n⟩.As illustrated in figure 3, it is obvious that the turbulence spreading diffusivity χ I is significantly different from the turbulent particle diffusivity D n .For either the case of lower ν * and higher α (red diamond) or the case of higher ν * and lower α (blue inverted triangle), the profiles of D n vary gently, and are bounded between 0-5 m 2 s −1 , as shown in figures 3(c) and (f ).However, extremely large turbulence spreading diffusivity χ I (10-30 m 2 s −1 ) is observed in the region at 24.0 cm < r < 24.5 cm, as shown in figure 3(e).Besides, sharp spikes in the profiles of χ I exist in the region at 24.5 cm < r < 25.0 cm, as shown in figures 3(b) and (e).These singular χ I values are well over 10 m 2 s −1 , up to ∼2000 m 2 s −1 .This is because the turbulence spreading flux ⟨ ṽr ñ2 ⟩ /2 is considerable in these regions, while the local intensity gradient ∂ r ⟨ ñ2 ⟩ is close to zero, as shown in Γ residual ̸ = 0), and that Fick's law is not sensible for describing turbulence spreading in these regions.It is clear that turbulence spreading is strongly non-diffusive.Besides, when the collisional and hydrodynamic response intensify (i.e.ν * goes up and α goes down), D n increases by ∼2 times, as shown in figures 3(c) and (f ), while χ I increases by more than 10, as shown in figures 3(b) and (e).χ I departs radically from D n .These suggest that the usual diffusive model frequently used to describe the turbulence spreading (with The apparently non-diffusive character of turbulence spreading brings to mind the intermittency (bursty events) of density fluctuations due to the presence of coherent structures, such as blobs and holes, within ambient turbulence [43][44][45][46][47][48][49][50][51][52][53][54][55][56][57].Thus, we study the relation between turbulence spreading and blobs/holes.Regarding the non-diffusive character of spreading, it is useful to define the mean jet velocity of turbulence spreading V I = ⟨ṽ r ñ2 ⟩/⟨ñ 2 ⟩, which can be regarded as the characteristic velocity of spatial propagation of turbulence internal energy [6].V I is effectively an average of radial velocity over the ensemble of density fluctuation intensity.The positive/negative skewness of density fluctuations ñ was widely used to identify the existence of blobs/holes, respectively [58].
Figure 4(a) shows the variation of mean jet velocity of turbulence spreading V I versus the skewness of ñ in the edge region at 24.0 cm < r < 24.5 cm.There are several measurement data points (typically [3][4][5]  Our experimental results reveal that the mean jet velocity of turbulence spreading V I manifests a linear correlation with the skewness of density fluctuations.As shown in figure 4(a), V I increases linearly with Skew (ñ).V I ranges from −100 m s −1 to 300 m s −1 and the skewness ranges from −0.25 to 0.60.It is worth noting that the points of zero skewness and the zero mean jet velocity are approximately coincident.The more positive mean jet velocity of turbulence spreading coincides with the more positive skewness of density fluctuations, while the more negative mean jet velocity coincides with the more negative skewness.This indicates that the positive-dominant density fluctuation events coincide with a radially outward velocity of turbulence spreading, while the negative-dominant density perturbations events coincide with a radially inward velocity of turbulence spreading.The mean jet velocity V I also increases as the auto-correlation time τ ac of ñ increases, as shown by figure 4(b).Here, τ ac is determined from the efolding width of the envelope of the auto correlation function of density fluctuations.The results above demonstrate a strong relation between turbulence spreading dynamics and coherent structures, like blobs and holes.
The particle transport velocity V T = ⟨ṽ r ñ⟩/⟨ñ 2 ⟩ 1/2 , which can be seen as the characteristic velocity of turbulent particle transport, has no obvious correlation with both skewness and auto-correlation time, as shown in figures 4(c) and (d).This suggests that the transport of turbulence energy does not simply follow along with turbulent particle transport.Note that the turbulent particle flux is positive (i.e.outward) although the skewness of density fluctuations is negative at 24.0 cm < r < 24.5 cm for the discharge marked with red diamonds in figure 4. Actually, figures 3(c) and 2(e) have shown this result.For all the discharges in our experiments, turbulent particle flux ⟨ñṽ⟩ r > 0 and −∂ r ⟨n⟩ > 0 lead to positive (i.e.outward) particle diffusivity D n = −⟨ñṽ⟩ r /∂ r ⟨n⟩ in the edge region across [24.0, 26.0] cm.There is no necessary relationship between the sign of skewness of ñ and the direction of the particle flux.In contrast to the linear correlation between mean jet velocity V I and Skew (ñ), a similar correlation between V I and the kurtosis of ñ is not observed.Figure 5(a) shows the variation of mean jet velocity of turbulence spreading V I versus the skewness of ñ in the edge region at 24.5 cm < r < 25.0 cm.It shows results consistent with those shown in figure 4. The mean jet velocity V I exhibits   Note that skewness measures the asymmetry of probability distribution function.The strong correlation between mean jet velocity and skewness suggests the importance of blobhole asymmetry to turbulence spreading.Figure 6(a Figure 6(a) shows that, the turbulence spreading flux induced by positive density fluctuations is positive, i.e. the turbulence intensity contributed by positive ñ propagates outward.The turbulence spreading flux induced by negative density fluctuations is negative, i.e. the turbulence intensity contributed by negative ñ propagates inward.The net total turbulence spreading flux ⟨ ṽr ñ2 ⟩ /2 is obtained from the sum of positive and negative components, as given by the values in the legend of figure 6(a).Note that for the case of weak spreading, marked by red diamonds, the distribution function is nearly symmetric.For increasing outward spreading, marked by green crosses and blue inverted triangles, the distribution function is strongly asymmetric.The positive contribution of high-amplitude positive density fluctuations increases while the negative contribution of high-amplitude negative density fluctuations decreases.This is especially true for the contribution of blobs and holes with absolute amplitudes higher than 2σ (ñ).This can be directly seen from figure 6(b).The blobs-dominant events with ñ ⩾ 2σ (ñ) cause radially outward turbulence spreading, while the holes-dominant events with ñ ⩽ −2σ (ñ) cause a radially inward turbulence spreading.For increasing total turbulence spreading ⟨ ṽr ñ2 ⟩ /2, the positive components contributed by blobs, i.e.
Figure 7(a) shows the distribution of turbulence spreading flux P spreading (A) as a function of density fluctuation amplitude A = |ñ| /σ ñ in the region of extremely large turbulence spreading diffusivity χ I at 24.0 cm < r < 24.5 cm, as shown in figure 3(e).It also shows that the turbulence spreading flux induced by positive density fluctuations is positive, and the turbulence spreading flux induced by negative density fluctuations is negative.For enhanced outward spreading, marked by green crosses and blue inverted triangles, the distribution function is more asymmetric.This is especially true for the contribution of blobs and holes with absolute amplitudes higher than 2σ (ñ).It is noted that for the increased outward spreading, the contribution of blobs (> 2σ (ñ)) to spreading even decreases slightly, as shown by the green crosses and blue inverted triangles in figure 7(b).The enhanced net total outward spreading is not related to the increased contribution of blobs (>2σ (ñ)) but the decreased contribution of holes (<2σ (ñ)) in this region.The increased symmetry breaking between positive and negative density fluctuation events-with ñ > 0 dominant-causes enhanced outward turbulence spreading.We also note that the symmetry breaking between positive and negative density fluctuation events-with ñ < 0 dominant-causes the total net inward turbulence spreading, as shown by the red diamond symbol.Recall that the negative (or inward) turbulence spreading flux in the region at 24.0 cm < r < 24.5 cm for the case with lower ν * and higher α shown by figure 2(d), and that the inward mean jet velocity of spreading coincides with the negative skewness of density fluctuations shown by figure 4(a).
To further explore the role of blob-hole dynamics in turbulence spreading, conditional average [47,62] methods, which are routinely used to isolate blobs and holes from ambient turbulence, are adopted in our work.The conditional average results of ñ in the region of spikey, unphysical spreading diffusivity χ I at 24.5 cm < r < 25.0 cm for discharges with weak and strong turbulence spreading are shown in figures 8(a1) and (a2).The time window of counting blobs and holes is 3 ms corresponding to 6000 sampling points.The conditional average curve is obtained by averaging the accumulation of burst events with a peak value exceeding the threshold of ±2σ in a time interval of ±15µs around the peak.The time coordinates of burst events in ñ are used to select the time series of ṽr .The conditional average results of ṽr are shown in figures 8(b1) and (b2).The temporal and spatial features of blobs and holes can be obtained by analyzing the conditional average results.The burst rate γ b of blobs and holes, which is the number of burst events in ñ divided by the total time they are counted, characterizes the intermittency of density fluctuations.The life time τ life of blobs and holes is determined from the full width at 1/e of the peak value in the conditional average of ñ.The radial propagation speed v prop of blobs and holes is estimated from The burst rate γ b , life time τ life , radial propagation speed v prop and characteristic radial length l r of blobs and holes for weak/strong turbulence spreading are given in table 1.We also give the values of some familiar relevant quantities to facilitate the appreciation of the temporal and spatial features of blobs/holes.As turbulence spreading becomes stronger: (i) the burst rate of blobs increases from ∼7.9 kHz to ∼12 kHz, while the burst rate of holes decreases prominently from ∼9.1 kHz to ∼0.3 kHz.Here, the electron diamagnetic drift frequency f * is ∼45 kHz for weak spreading and ∼12 kHz for strong spreading.(ii) the life time τ life of blobs and holes increases simultaneously from ∼8 µs to ∼10 µs, which is reasonable since τ ac of ñ is ∼9 µs for weak spreading and ∼14 µs for strong spreading.(iii) the radial propagation speed v prop of blobs increases from 70 m s −1 to 630 m s −1 and v prop of holes increases from −18 m s −1 to −360 m s −1 .Here, the mean jet velocity of turbulence spreading V I increases from ∼55 m s −1 to ∼350 m s −1 , while the electron diamagnetic drift velocity v * decreases from ∼1550 m s −1 to ∼570 m s −1 .Note that the mean jet velocity is somewhere between the radial propagation speed of blobs and holes, for either weak spreading or strong spreading.This outcome is reasonable, as the mean jet velocity is an ensemble averaged quantity.(iv) the characteristic radial length l r of blobs increases (from 0.6 mm to 6.6 mm) and l r of holes also increases (from 0.1 mm to 3.7 mm).For weak spreading, the characteristic radial lengths l r of blobs and holes are close to the ion gyro radius ρ i ∼ 0.3 mm.For strong spreading, l r of blobs and holes are close to the radial correlation length of density fluctuations l corr ∼ 4.4 mm.This is reasonable as figure 8 implies that regimes of strong spreading are necessarily blob dominated.Note that v prop of blobs are positive and v prop of holes are negative, indicating blobs move outward while holes move inward.This is consistent with previous studies [58].This again shows that positive density fluctuation events cause radially outward spreading of turbulence intensity (ñ 2 ), while negative density fluctuation events cause radially inward spreading of turbulence intensity (ñ 2 ), as presented in figure 8.For the case with weak turbulence spreading, the burst rate of blobs and holes are very close to one another, which corresponds to the nearly symmetrical distribution function of spreading flux in figure 6(a).For the case with strong turbulence spreading (where the collisionality ν * is high and adiabaticity α is low), the burst rate of blobs is about 40 times larger than that of holes, which corresponds to a remarkable symmetry breaking in the distribution function of the spreading flux.Therefore, the difference in burst rates of blobs and holes plays a significant role in the turbulence spreading dynamics.For strong spreading, the width of the singular turbulence spreading diffusivity (χ I ) region with ∂ r ñ2 ∼ 0 is ∼5 mm, shown by figures 3(d) and (e).This width falls between the characteristic radial lengths l r of blobs and holes, i.e. 3.7-6.6mm.This clarifies the importance of intermittent convective coherent events to turbulence spreading.
To summarize, an in-depth study of turbulence spreading and its relation to convective coherent structures was carried out by using direct experimental measurements of edge plasmas.Results indicate that turbulence spreading is significantly enhanced as collisions become more frequent and the electron response becomes hydrodynamic.Turbulence spreading is shown to be strongly non-diffusive.Attempts to characterize turbulence spreading as a Fickian diffusion process yield a singular (spikey, unphysical) profile of turbulence spreading diffusivity χ I , which departs significantly from the turbulent particle diffusivity D n .Thus, turbulence spreading is not amenable to characterization as a diffusion process with a diffusivity closely related to the turbulent diffusivity, as often assumed.To quantify the strength of spreading, it is useful to define the mean jet velocity of turbulence spreading V I .The mean jet velocity V I manifests a linear correlation with the skewness of density fluctuations.The zero point of skewness and the zero point of the mean jet velocity approximately intersect.The more positive mean jet velocity of turbulence spreading coincides with the more positive skewness of density fluctuations, while the more negative mean jet velocity coincides with the more negative skewness.V I increases with the auto-correlation time of ñ while the particle transport velocity V T does not.This suggests that the transport of turbulent internal energy does not simply follow along with turbulent particle transport.The distribution of turbulence spreading flux exhibits interesting trends in its dependence on amplitudes of positive density fluctuations and negative density fluctuations.Turbulence spreading flux induced by positive density fluctuations is positive, i.e. propagating outward.Turbulence spreading flux induced by negative density fluctuations is negative, i.e. propagating inward.For the case with weak turbulence spreading, the burst rate of blobs and holes are very close to each other and the distribution function is symmetrical.For the case with strong turbulence spreading at higher collisionality or lower adiabaticity, the burst rate of blobs is much larger than that of holes and the distribution function of spreading flux is strongly asymmetric.The increased symmetry breaking between convective coherent structures, i.e. asymmetry between the populations of outgoing blobs and incoming holes, causes enhanced net spreading.The width of the region with singular turbulence spreading diffusivity (χ I ) and ∂ r ⟨ ñ2 ⟩ ∼ 0 falls between the characteristic radial lengths of blobs and holes.These results clarify the importance of intermittent convective coherent events in the basic structure of turbulence spreading.Spreading is seen to be fundamentally intermittent and non-diffusive.These results present a significant challenge to conventional models of turbulence spreading.
Regarding future work, we note that recently, there has been a resurgence of interest in turbulence spreading and related phenomena in the context of SOL heat load width or power scrape off width (λ q ) broadening.In many cases, especially in H-mode conditions, local SOL turbulence is quenched, so λ q collapses to the neoclassical prediction, which is unacceptably small.The hope is that turbulence spreading from the pedestal may sufficiently energize SOL turbulence, so as to broaden λ q with enhanced cross-field transport [13].More generally, some recent studies indicates that even in Lmode, spreading from the edge has a significant or even dominant impact on SOL turbulence levels, and thus on the SOL width.The findings of this letter indicate that such edge → SOL spreading is strongly intermittent and carried by structures (blobs).Thus, mean field models are destined to fail.The natural characteristic velocity for this spreading process is the jet or spreading velocity V I = ⟨ṽ r ñ2 ⟩/⟨ñ 2 ⟩. Figure 4(a) shows that V I correlates well with the skewness (indicative of structures) and with spreading.Thus, we conjecture that the SOL width resulting from spreading-induced broadening is λ SP ∼ V I τ d .Here, τ d is the SOL dwell time constrained by either the heat residence time or the auto-correlation time of turbulence, whichever is shorter.The detailed dependence of the heat residence time varies with collisionality.The relevant comparison is then λ SP vs λ HD , where λ HD ∼ V d τ d , and V d is the magnetic drift velocity.In this case, τ d ∼ Rq/c s .Thus, the natural figure of merit for comparing spreading and neoclassical processes is based on comparison of jet and drift velocities, i.e.V I and V d .This should be contrasted to the conventional approach, which is to compare ⟨ṽ 2 r 1/2 ⟩ to V d .
In particular, the conventional wisdom does not account for the coherence and cross phase of ṽr with ñ2 .This cross phase encapsulates the key physics of spreading.Similarly, criteria tied to turbulent transport are necessarily based upon the transport cross phase in ⟨ṽ r ñ⟩.This differs fundamentally from that which appears in ⟨ ṽr ñ2 ⟩ .In another word, turbulence spreading flux depends upon the coherence and cross phase of ṽr with ñ2 , while the turbulent particle flux depends upon the coherence and cross phase of ṽr with ñ.Increasing ñ will likely be beneficial to increase spreading flux and particle flux, but the actual spreading flux and particle flux need not increase in direct proportion, or even in the same direction simultaneously.Inward turbulence spreading but outward particle transport in the region of 24.0 cm < r < 24.5 cm has been observed for the discharge marked by red diamonds in figure 4 (also figures 2 and 3).Future work will focus on testing the conjecture set forth above and on quantifying the parametric dependences of the jet velocity.
We should also point out that, in our experiments, the collisionality near the LCFS is higher than 1 and the adiabaticity near the LCFS is lower than 1 (corresponding to the hydrodynamic regime).However, collisional plasmas are not always hydrodynamic.Theoretically, while the plasma edge is usually collisional, it can be either adiabatic or hydrodynamic, depending upon α = k 2 ∥ v 2 th,e /ν ei ω * .Turbulence in the adiabatic regime is observed to support avalanches.Turbulence avalanching involves turbulence spreading.Avalanches manifest 'bumps' (i.e.local excesses) or 'holes' (i.e.local deficits) in simulation works [8].A further comment is given in the appendix.

Appendix
One can distinguish two types of classification regarding the regimes for electron response: (i) collisional vs collisionlessdetermined by collision mean free path vs connection length, i.e. λ mfp < Rq vs λ mfp > Rq, where λ mfp ∼ v th,e /ν ei and Rq ∼ 1/k ∥ ; (ii) hydrodynamic vs adiabatic-depending upon α = k 2 ∥ v 2 th,e /ν ei ω * < 1 vs > 1.In practice, here k ∥ ∼ 1/Rq is used and all other quantities are measured.The adiabatic regime has electrons responding along the magnetic field at a rate exceeding the wave frequency).Then, ñ/n ∼ e φ /T and the fluctuations are drift wave like.The hydrodynamic regime is one of convective cells, where electrons oscillate faster than their response along the magnetic field.The collisionless electron response (λ mfp > Rq) is almost always adiabatic (k ∥ v th,e /ω * > 1) unless k ∥ → 0. While the plasma edge is usually collisional, it can be either adiabatic or hydrodynamic, depending upon k 2 ∥ v 2 th,e /ν ei ω * .In a burning plasma, the edge is more likely (but not definitely) to be adiabatic.Turbulence spreading can surely persist in adiabatic electron regimes.For wavelike fluctuations, a piece of turbulence spreading is simply the wave energy density flux [63], i.e.Γ * spread ∼ ∑ k V gr (k) ε k .Here, ε k is the wave energy density and V gr is the radial group velocity.This coexists with the eddy mixing process with Γ spread ∼ ⟨ṽ r ε⟩.The ratio of these two is closely related to the Kubo number, i.e.Ku ∼ ⟨ṽ r ε⟩ / ∑ k V gr (k) ε k .Turbulence spreading of wave energy can be alive and well in the adiabatic regime, when Ku is small.There are numerous gyrokinetic simulations which study turbulence spreading in collisionless or nearly collisionless regimes with adiabatic electrons.These are summarized in the paper by T.S. Hahm and P.H. Diamond [8].Simulation results show that turbulence spreading in the collisionless regime is closely related to avalanches.These avalanches manifest 'bumps' (i.e.local excesses) or 'holes' (i.e.local deficits), somewhat akin to 'blobs' and 'holes'.A precise comparison and contrast between these two sets of classification has not yet been undertaken.Whether the turbulence spreading in the adiabatic regime is related to the dynamics of 'blobs and holes' is still unresolved.

Figure 1 .
Figure 1.A simple diagram of turbulence spreading induced by positive (P) and negative (N) density fluctuation events at edge plasma.

Figure 3 .
Figure 3. Turbulence intensity ñ2 , turbulence spreading diffusivity χ I and turbulence particle diffusivity Dn for discharge with lower ν * and higher α (left column, red diamond) and for discharge with higher ν * and lower α (right column, blue inverted triangle).(Curves with red diamond correspond to the discharge with Bt ∼ 2.2 T, Ip ∼ 187 kA and n ∼ 4.1 × 10 19 m −3 ; curves with blue inverted triangle correspond to the discharge with Bt ∼ 1.6 T, Ip ∼ 126 kA and n ∼ 3.9 × 10 19 m −3 .).

Figure 4 .
Figure 4. (a) Mean jet velocity of turbulence spreading V I VS skewness of density fluctuations ñ; (b) mean jet velocity of turbulence spreading V I VS auto-correlation time of ñ; (c) particle transport velocity V T VS skewness of ñ; and (d) particle transport velocity V T VS auto-correlation time of ñ, in the edge region of 24.0 cm < r < 24.5 cm.(Data points with different symbols correspond to various discharge shot numbers.).

Figure 5 .
Figure 5. (a) Mean jet velocity of turbulence spreading V I VS skewness of density fluctuations ñ; (b) mean jet velocity of turbulence spreading V I VS auto-correlation time of ñ; (c) particle transport velocity V T VS skewness of ñ; and (d) particle transport velocity V T VS auto-correlation time of ñ, in the edge region of 24.5 cm < r < 25.0 cm.(Data points with different symbols correspond to various discharge shot numbers.).
) shows the distribution of turbulence spreading flux as a function of density fluctuation amplitude in the region of spikey, unphysical spreading diffusivity χ I at 24.5 cm < r < 25.0 cm.The horizontal coordinate is the density fluctuation amplitude normalized to a multiple of the standard deviation, i.e.A = |ñ| /σ ñ.The vertical coordinate is the accumulated and averaged turbulence spreading contributed by the density fluctuations within a certain range of amplitude, i.e.P spreading (A) = [ ∑ (A−0.05)σñ<ñ<(A+0.05)σñ(ṽ r ñ2 /2)]/M for positive ñ and P spreading (A) = [ ∑ −(A+0.05)σñ<ñ<−(A−0.05)σñ(ṽ r ñ2 /2)]/M for negative ñ.The radial velocity fluctuations with the same temporal sampling interval corresponding to the ñ amplitude range are taken for the calculations.M = 6000, which is the total number of sampling points of ñ or ṽr .The total turbulence spreading flux is the sum of distribution function, i.e.

Figure 8 .
Figure 8. Conditional average results of ñ and ṽr for discharges with weak turbulence spreading and strong turbulence spreading.

Table 1 .
The burst rate, life time, radial propagation speed and characteristic radial length of blobs/holes for weak/strong turbulence spreading.