Overview of edge modeling efforts for advanced divertor configurations in NSTX-U with magnetic perturbation fields H. Frerichs, O. Schmitz, I. Waters, G. P. Canal, T. E. Evans, Y. Feng and V. Soukhanovskii University of Wisconsin - Madison ISTW - November 3-6, 2015 H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U
Motivation: control of divertor head loads Control of divertor heat loads (both steady state and transient) remains one of the key challenges for a fusion based reactor and for a compact fusion nuclear science development facility (FSNF) Resonant magnetic perturbations (RMPs) are a promising method for ELM control breaking of axisymmetry Advanced divertors (specialized divertor geometry): Snowflake, X-divertor) for steady state heat flux reduction Likely, both concepts will have to work together, and NSTX-U is a well suited device to study this: What is the impact on neutral fueling and exhaust (density control), and how does this affect high recycling and transition to detachment? H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 1
1 Introduction to magnetic perturbations 2 The snowflake divertor configuration 3 Edge transport modeling Introduction to EMC3-EIRENE First simulation results H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 2
Perturbation of the magnetic separatrix lobes The separatrix associated with X has 2 branches (set of field line trajectories): W + = W = { } p lim F p (l) X l { } p lim F p(l) X l W W + X point F p (l): field line through p Both branches overlap in the unperturbed configuration (arrows indicate field direction) H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 3
Perturbation of the magnetic separatrix lobes The separatrix associated with X has 2 branches (set of field line trajectories): W + = W = { } p lim F p (l) X l { } p lim F p(l) X l + W W (arrows indicate field direction) F p (l): field line through p Both branches overlap in the unperturbed configuration Magnetic perturbations result in a splitting, and both branches may intersect each other transversely This opens up a connection between the plasma interior and the divertor targets H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 3
Lobe density increases with lower aspect ratio -135-140 n = 3-80 -90 n = 3 Vertical position [cm] -145-150 -155-160 A = 1.7 NSTX U Vertical position [cm] -100-110 -120-130 A = 2.8 DIII D 40 45 50 55 60 65 70 75 80 85 Major radius [cm] -140 100 110 120 130 140 150 160 Major radius [cm] Lobes become smaller and more frequent towards the X-point lobe density: 1 L n B tor 2π R B pol n B 0 2π R 0 B pol ( 1 δ ) 2 X A Higher density may facilitate heat flux spreading between lobes Radial extension of lobes depends on perturbation strength H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 4
Radial size of lobes defines poloidal extend of the magnetic footprint, lobe density not affected by I c -138-140 NSTX U: ISP I c [ka] = 1 2 3 4-142 Vertical Position [cm] -144-146 -148-150 -152 0 50 100 150 200 250 300 350 Toroidal Angle [deg] H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 5
1 Introduction to magnetic perturbations 2 The snowflake divertor configuration 3 Edge transport modeling Introduction to EMC3-EIRENE First simulation results H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 6
The snowflake divertor configuration Second order null of the poloidal magnetic field the separatrix acquires a characteristic hexagonal form This results in a flux expansion near the null-point, and while B (snowflake) pol r 2 (standard divertor) B pol r: distance from the null-point r SF facilitates longer connection lengths and two additional strike points Generalization: two first order nulls in (close) proximity D. Ryutov et al, Phys. Plasmas 15, 092501 (2008) H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 7
The snowflake configuration allows for a variety of magnetic topologies at NSTX-U Nearly exact Snowflake (SF0) Normalized poloidal flux Nearly exact snowflake : approximation to classical snowflake H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 8
The snowflake configuration allows for a variety of magnetic topologies at NSTX-U Nearly exact Snowflake (SF0) Normalized poloidal flux Nearly exact snowflake : approximation to classical snowflake Here SF0 is actually snowflake minus, which is topologically equivalent to connected double null (CDN) H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 8
The snowflake configuration allows for a variety of magnetic topologies at NSTX-U Nearly exact Snowflake (SF0) Tilted Snowflake Minus (SF ) Snowflake Plus (SF+) 135111 Njw_08 116313 NfHz0+Qf 116313 NfHz0+Qg Normalized poloidal flux Nearly exact snowflake : approximation to classical snowflake Here SF0 is actually snowflake minus, which is topologically equivalent to connected double null (CDN) The secondary X-points in the (tilted) SF- and SF+ configurations here are outside the divertor targets topology is lower single null (LSN) H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 8
SF0 has increased connection length, but flux compression on target (with respect to SD) Connection Length [ m] 140 120 100 80 60 40 20 PFR Outer Strike Point SOL SD SF- SF+ SF0 Distance from OSP [cm] 25 20 15 10 5 Flux surface connection SD SF- SF+ SF0 0-2 0 2 4 6 8 10 12 Distance from Separatrix [cm] 0 0 0.5 1 1.5 2 2.5 3 3.5 4 Distance from Separatrix (LFS midplane) [cm] Both SF0 and SF- feature longer connection lengths then SD at the outer strike point: may result in lower divertor temperatures SF0 features flux compression while SF- features flux expansion with respect to standard divertor (SD) configuration: can impact heat flux spreading H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 9
SF0 has increased connection length, but flux compression on target (with respect to SD) Connection Length [ m] 140 120 100 80 60 40 20 PFR Outer Strike Point SOL SD SF- SF+ SF0 SF0 (div. adjust) Distance from OSP [cm] 25 20 15 10 5 Flux surface connection SD SF- SF+ SF0 SF0 (div. adjust) 0-2 0 2 4 6 8 10 12 Distance from Separatrix [cm] 0 0 0.5 1 1.5 2 2.5 3 3.5 4 Distance from Separatrix (LFS midplane) [cm] Both SF0 and SF- feature longer connection lengths then SD at the outer strike point: may result in lower divertor temperatures SF0 features flux compression while SF- features flux expansion with respect to standard divertor (SD) configuration: can impact heat flux spreading H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 9
Lobe size and extension depends on equilibrium Vertical Position [cm] -120-125 -130-135 -140-145 -150-155 BT = 1T -160 IP = 2 MA SD -120-125 -130-135 -140-145 -150-155 -160 Vertical Position [cm] 40 45 50 55 60 65 70 75 80 85 BT = 1T IP = 1 MA 40 45 50 55 60 65 70 75 80 85 Major Radius [cm] Major Radius [cm] SF0 Lobe density is much higher in SF0, i.e. lobe size is smaller (I p = 1 MA) Lobe size scales with r 2 near the X-point 4 3 in SF0, because B pol does 2 very small change in size between 1 neighboring lobes 0 Lobe spacing [cm] 9 8 7 6 5 SD SF- SF+ SF0 0 5 10 15 20 Distance from X-point along separatrix [cm] H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 10
1 Introduction to magnetic perturbations 2 The snowflake divertor configuration 3 Edge transport modeling Introduction to EMC3-EIRENE First simulation results H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 11
A 3D steady state transport model for the edge plasma Field lines are reconstructed from a 3D block-structured grid which can be adapted to the topology at hand H. Frerichs et al., Comp. Phys. Commun. 181 (2010) 61 Classical transport (Braginskii) along field lines Self-consistent solution by iterative application Magnetic field structure FLARE: 3D field line grid Plasma Transport EMC3: Fluid model Monte Carlo method Plasma Background n, u, Te, Ti Sources/Sinks Sp, Sm, See, Sei Transport coefficients Neutral Gas Transport EIRENE: Kinetic Model Monte Carlo method The following simulations are based on: Particle in-flux Γ in = 3.12 10 20 s 1 (50 A) Recycling coefficient c rec = 0.99 Edge input power P in = 2 MW Anomalous cross-field transport: (particles) D = 0.3 m 2 s 1 (energy) χ e = χ i = 2.0 m 2 s 1 Impurity transport and radiation is neglected at this point H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 12
No RMPs: SF0 operates at higher divertor density H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 13
RMPs: Lobes at the LFS, diffusion to strong at HFS H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 14
Heat load analysis favours SF- configuration Target Heat Load [MW m -2 ] 8 7 6 5 4 3 2 1 Heat load in axisymmetric configuration 0-5 0 5 10 15 20 25 Distance from Separatrix [cm] SD SF- SF+ SF0 Small reduction of peak heat load in SF+ Moderate reduction of peak heat load in SF-, peak is found at 10 20 cm distance from the separatrix Significant increase of peak heat load in SF0 due to compression of flux surfaces at the divertor target. H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 15
Heat load analysis favours SF- configuration Target Heat Load [MW m -2 ] 8 7 6 5 4 3 2 1 Heat load in axisymmetric configuration 0-5 0 5 10 15 20 25 Distance from Separatrix [cm] SD SF- SF+ SF0 Target Heat Load [MW m -2 ] 8 7 6 5 4 3 2 1 Heat load in RMP configuration 0-5 0 5 10 15 20 25 Distance from Separatrix [cm] SD SF- SF+ SF0 Small reduction of peak heat load in SF+ Moderate reduction of peak heat load in SF-, peak is found at 10 20 cm distance from the separatrix Significant increase of peak heat load in SF0 due to compression of flux surfaces at the divertor target. No significant impact of RMPs on toroidally averaged heat loads H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 15
Summary and conclusions Application of RMPs results in the formation of helical lobes which have a higher frequency at lower aspect ratio The snowflake family of divertor configurations allows a variety of magnetic topologies ranging from flux compression to flux expansion at the divertor targets First transport simulations favour the SF- configuration for heat load reduction Outlook: How do RMPs and advanced divertor configurations impact neutral fueling and exhaust, and how does this affect high recycling and transition to detachment? What is the impact of plasma response effects (use NIMROD, M3D-C1) H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 16
2w h R 0 w aspect ratio: A = R 0 /w elongation: κ = h/w triangularity: δ = /w 2h X H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 17
Radial lobe size depends on perturbation strength in relation to equilibrium field Vertical Position [cm] Vertical Position [cm] -120-120 -125 SD -125 SF -130-130 -135-135 -140-140 -145-145 -150 0000-150 Ψ = 1.024 B T = 1T -155-155 B T = 1T IP = 2 MA Ψ = 1.024-160 IP = 2 MA -160 000000000000000000000000000000000000000000 111111111111111111111111111111111111111111 000000000000000000000000000000000000000000 40 000000000000000000000000000000000000000000 111111111111111111111111111111111111111111-120 45 50 55 60 65 70 75 80 85 40-120 45 50 55 60 65 70 75 80 85 000000000000000000000000000000000000000000 000000000000000000000000000000000000000000 000000000000000000000000000000000000000000 111111111111111111111111111111111111111111 111111111111111111111111111111111111111111-125 Major Radius [cm] Major Radius [cm] SF+ -125 SF0-130 -130-135 -135-140 -140-145 -145 Ψ = 1.024 Ψ = 1.033-150 -150 Ψ = 1.024 B -155 T = 1T B T = 1T -155 I P = 2 MA IP = 1 MA -160-160 40 45 50 55 60 65 70 75 80 85 40 45 50 55 60 65 70 75 80 85 Major Radius [cm] Vertical Position [cm] Vertical Position [cm] Major Radius [cm] H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 18
Field lines can connect through the lobes from inside the separatrix to the divertor targets H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 19
Moiré pattern caused by small, similar sized lobes Lobe frequency is much higher in SF0, i.e. lobe size is smaller (I p = 1 MA) Lobe size scales with r 2 near the X-point 4 3 in SF0, because B pol does 2 very small change in size between 1 neighboring lobes 0 Lobe spacing [cm] 9 8 7 6 5 SD SF- SF+ SF0 0 5 10 15 20 Distance from X-point along separatrix [cm] H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 20
A 3D steady state fluid model for the edge plasma Particle balance (n: plasma density) [nu e D e e n ] = S p D : anomalous cross-field diffusion, S p: ionization of neutral particles Momentum balance (u : fluid velocity parallel to magnetic field lines) [ m i nu 2 e η e e u D e e ( m i nu ) ] = e n (T e + T i ) + S m η T 5/2 i : parallel viscosity, η = m i nd : cross-field viscosity, S m: interaction (CX) with neutral particles Energy balance (T e, T i : electron and ion temperature) [ 5 ( 2 Te nu e D e e n ) ( ] ) κ ee e + χ en e e Te [ 5 2 T ( i nu e D e e n ) ( ] ) κ i e e + χ i n e e Ti = +k (T i T e) + S ee = k (T i T e) + S ei κ e,i T 5/2 e,i : classical parallel heat conductivity, χ e, χ i : anomalous cross-field transport, k n 2 T 3/2 e : energy exchange between el. and ions, S ee, S ei : interaction with neutral particles and impurities (radiation) H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 21
Field lines are reconstructed from a finite flux-tube grid The 3D grid is generated by field lines tracing starting from 2D base grids Grid nodes Pre calculated field lines > 3D finite flux tube grid Reconstructed field line ϕ base EMC3: bilinear interpolation within flux-tubes (4 pre-calculated field lines) 100 50 0-50 -100 Separatrix LCFS 100 50 LSN 0 DDN 0 CDN -50-100 Inner sepx. Outer sepx. LCFS 100 120 140 160 180 200 220 240 100 120 140 160 180 200 220 240 150 100 50-50 -100-150 Separatrix LCFS 20 40 60 80 100 120 140 160 Discretization in the cross-field direction can be adapted to the magnetic configuration at hand. Unlike grids for 2D transport modeling, the actual information about the magnetic configuration is not stored in the 2D base grid(s) but in the 3D grid. H. Frerichs (hfrerichs@wisc.edu) 3D edge modeling for NTSX-U 22