L-H Transition Dynamics and Power Threshold Minimum M.A. Malkov 1, P.H. Diamond 1,2, K. Miki 2,3 and G.R. Tynan 1 1 University of California, San Diego, USA; e-mail: mmalkov@ucsd.edu 2 WCI, NFRI Daejeon, Korea, 3 JAEA, Kashiwa, Japan APS 2014 CP8.00048 1 / 23
Abstract A physical model to explain the minimum in the power threshold of the L H transition is suggested. Its crucial new elements, compared to the prototype due to (Miki & Diamond 2012), are the separately evolved electron and ion heat transfers and the independent power supply in each chanel. The power threshold minimum may be explained as follows. The low-density branch, where the heat initially goes to the electrons, is associated with growing eciency of e i heat transfer and the ow generation. This suppresses the transport and lowers the threshold P th.. The high-density branch is associated with increased collisional damping of the ow. Model studies also reveal: (a) an increase in threshold power for o-axis electron heat deposition and (b) the absence of a clear P thr (n) minimum for pure ion heating. 2 / 23
Outline 1 Basic physics of LH transition 2 Recent Experiments and Shortcomings of Available Models 3 Predecessor Model and its Extension to Studies of P th Minimum 4 Model Equations 5 Results and analysis 6 Conclusions 3 / 23
Mechanism and occurrence of LH transition originates via coupling of turbulence to low frequency shear ows by Reynolds work Reynolds work causes collapse of turbulence and then turbulent transport diamagnetic electric eld grows, associated with P LH transition can occur via a protracted I-phase or in a single burst of shear ow 4 / 23
Objectives of this work (ongoing) establish link between microscopics and macroscopics in power threshold scaling reproduce and understand observed threshold P th (n) minimum explore P th in terms of other parameters, such as e-i thermal coupling eciency, noise... investigate the role of heating prole in LH transition e-i heating branching ratio role of mean shear in locking-in of the transition 5 / 23
Observations of power threshold minimum Ryter et al 2013 [1] ion heat ux plays a dominant role in LH transition P electron channel thought to be ignorable But: in low-density regimes with dominant EC heating electrons must transfer energy to ions to steepen P electron description must be separated from that of ions temperature and density dependence of collision rate need to be included (average values do not suce) 6 / 23
1-D Numerical Model Based on 1-D numerical 5-eld model (Miki & Diamond 2012,13+) [2] signicantly extends Kim & Diamond 2003 [3] 0D model MD 2012 captures transition layer evolution but does not separate species modify MD 2012 by adding separate electron heat transport equation include e-i thermal coupling depending on locally evolved temperatures and density include these parameters in ZF damping description include trapped electron growth 7 / 23
Predator-Prey Model Equations Heat transport i,e: [ ] P i,e + 1 t r r rγ(p) = ± 2m (r a) 2 i,e Mτ (P e P i ) + Q exp 2 r 2 Density χ, D t = Γ = (χ neo + χ t ) P r [ ] n t + 1 a r (a r) 2 r r rγ(n) = Γ a L 2 exp 2L 2 dep dep Γ (n) = (D neo + D t ) n r τ c C 2 s I 1 + α t V E, 2 V E ( = ρ i C s L 1 p L 1 p ) L 1 n Vϑ 8 / 23
Equations cont'd DW turbulence with ITG and TEM instabilities ( ) I t = γ L ωi α 0 E 0 α V V E 2 I +χ N r I I r, χ N ω C 2 s ( ) C s R γ L = γ 0i R L ( ) + γ 0e C s L 1 T R L p L n R e + L 1 n ZF energy E 0 t crit α 0 E 0 I = 1 + ζ 0 V ϑ γ 2 dampe 0 mean ow shear ( ) V ϑ = α 5 γ 0i C 2 a R s R L t R L p L n R µ neo ν ( ) ii q 2 R 2 V ϑ 1.17C s ρ i L 1 T crit I r 9 / 23
Control parameters for transition plasma density - n ref -reference (xed) density - n -center line averaged (current) density Heat mix parameter Q i H mix = Q i + Q e Q i Q i.e., Q i = H mix Q, Q e = (1 H mix ) Q, where Q denotes the total power deposited into the plasma. widths of the heat sources r e = r i = 0.15a, heat deposition radii a e,i = 0.3a 10 / 23
Transition morphology L H transition event shown in four characteristic variables center line averaged: density, ZF energy, MF E B velocity, DW energy. shown as functions of heating rate Q (t). Data points are taken at equal time intervals so their density indicates both the rate at which Q is changing and how quickly the changes in the variables occur. 11 / 23
Transition morphology <n>, density 2.4 2.2 2.0 1.8 <V θ >, Mean Flow 0.037 0.035 0.033 0.031 ln<(e0)>, ZF Energy 1.6 0.0065 0.0075 0.0085-10 -11-12 -13-14 -15-16 Q(t) -17 0.0065 0.0075 0.0085 Q(t) ln<i>, DW Energy 0.029 0.0065 0.0075 0.0085-10.3-10.5-10.7 Q(t) -10.9 0.0065 0.0075 0.0085 Q(t) 12 / 23
Transition morphology Top row: - an example of LH transition with an extended pre-transition I-phase shown for the electron pressure P e, ZF energy and the MF velocity - strong, edge-localized MF is a marker of the H-mode. Bottom row: - an example of a failed transition with inward propagation - the edge MF jet starts to form but then merge with the large scale MF 13 / 23
Transition morphology 14 / 23
Identifying transitions take half-way to clean pedestal cross-check with DW,ZF,MF channels need transition criterion to scan P th ( n, He,i, L dep,... ) - very weak transitions occur and are hard to detect 15 / 23
Spatio-temporal dynamics of transition I-phase persists before transition clearly spatio-temporal behavior beyond 0-D model (similar to MD 2012) ZF signicantly advances into the core before transition 16 / 23
Spatio-temporal dynamics of transition slight temperature attening in the core due to enhanced turbulent transport I-phase in density 17 / 23
P th (n) scans Pth 0.008 0.007 0.006 0.005 0.004 Squares indicate strong transitions with the density jumps 0.1 Circles indicate weaker transitions. as a function of reference density 0.003 0.002 0.001 0.000 0.8 1.2 1.6 2.0 n, density Pth 0.008 0.006 0.004 P thr, shown against the center-line averaged n 0.002 0.000 1.0 1.5 2.0 2.5 3.0 Nref, reference density 18 / 23
P th (n, H mix ) scans 0.0 0.2 0.4 0.6 Hmix 0.8 1.0 0.0 0.5 1.0 1.5 n 2.0 0.006 0.004Pth 0.002 0.000 2.5 0.008 P thr in heating mix -density variables, H mix and n. n, H mix choice: electron/ion biased heating at lower/higher densities extended sub-sample at H mix = 1 is also included Monotonic dependence of H mix (n), arbitrarily chosen from the sample shown in the previous Figure and the resulting P thr (n) constrained by the above relation. Pth 0.010 0.008 0.006 0.004 0.002 0.000 Hmix (n) Pth (n) 0.0 0.5 1.0 1.5 2.0 n/10 20 m -3 1.0 0.8 0.6 0.4 0.2 Hmix 19 / 23
Model Findings no clear threshold minimum for pure ion heat deposition minimum consistend with observations [1, 4] found for mixed e i heating on low-density branch of P thr (n) electrons absorb most of the heat initially and as they transfer it to ions more eciently with growing density, P thr (n) decreases on high-density branch ions are heated and as ZF damping is growing with n, P thr must also grow 20 / 23
Model Findings the overall picture is consistent with the following two premises: - L-H transition is locked in by V E ( P i/n) - DW turbulence coupling to ow is a key trigger P th increases for o-axis electron heat deposition (reduced electron-ion coupling) shallow minimum power is predicted for heating mix scan as well as for density scan from above ndings a global minimum in multi-parameter space is predicted Ongoing work: quantifying strength of hysteresis in terms of multiple macroscopic parameters; relating this to observed back-transition shear ow and turbulence dynamics 21 / 23
Conclusions an extended 6-eld 1-D PDE model is developed (P e, P i, n, DW, ZF, Mean Flow) link between microscopics (e-i collisional heat exchange, turbulence) and macroscopics (transport barrier, P-n proles) in power threshold scaling is established threshold P th (n) minimum is reproduced and understood using a simple model of e-i heat transfer P th ( n, Ldep,.... ) is explored in terms of its dependence on other parameters, such as e-i thermal coupling eciency the role of heating prole in LH transition is investigated the role of e-i heating split ratio is studied, minimum of P th predicted role of mean shear in locking-in of transition is signicant 22 / 23
References [1] F. Ryter et al., Nucl. Fusion 53, 113003 (2013). [2] K. Miki et al., Phys. Plasmas 19, 092306 (2012). [3] E.-J. Kim and P. H. Diamond, Physical Review Letters 90, 185006 (2003). [4] C. F. Maggi et al., Nuclear Fusion 54, 023007 (2014). 23 / 23