full-wave simulations of lower hybrid wave …...57th annual meeting of the aps division of plasma...
TRANSCRIPT
![Page 1: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/1.jpg)
Full-wave Simulations of Lower
Hybrid Wave Propagation
in the EAST Tokamak
P. T. BONOLI, J. P. LEE, S. SHIRAIWA,
J. C. WRIGHT, MIT-PSFC,
B. DING, C. YANG, CAS-IPP, Hefei
57th Annual Meeting of the APS Division of Plasma Physics
Savannah, GA
November 18, 2015
Poster NP12.00050
![Page 2: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/2.jpg)
AbstractStudies of lower hybrid (LH) wave propagation have been
conducted in the EAST tokamak where electron Landau
damping (ELD) of the wave is typically weak, resulting in
multiple passes of the wave front prior to its being absorbed in
the plasma core. Under these conditions it is interesting to
investigate full-wave effects that can become important at the
plasma cut-off where the wave is reflected at the edge, as well
as full-wave effects such as caustic formation in the core. High
fidelity LH full-wave simulations were performed for EAST
using the TORLH field solver [1]. These simulations used
suffcient poloidal mode resolution to resolve the perpendicular
wavelengths associated with electron Landau damping of the
LH wave at the plasma periphery, thus achieving fully
converged electric field solutions at all radii of the plasma.
Comparison of these results with ray tracing simulations [2, 3]
will also be presented.
![Page 3: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/3.jpg)
This work is the subject of the 2017 Theory
and Performance Target
• Lower Hybrid current drive (LHCD) will be indispensable for driving
off-axis current during long-pulse operation of future burning plasma
experiments including ITER, since it offers important leverage for
controlling damaging transients caused by magnetohydrodynamic
instabilities.
• However, the experimentally demonstrated high efficiency of LHCD is
incompletely understood [in weak damping regimes].
• In FY 2017, massively parallel, high resolution simulations with 480
radial elements and 4095 poloidal modes will be performed using full-
wave radiofrequency field solvers and particle Fokker-Planck codes to
elucidate the roles of toroidicity and full-wave effects.
• The simulation predictions will be compared with experimental data
from the superconducting EAST tokamak.
• This work spans the research interests of the CSWPI SciDAC Center,
the MIT Theory Grant, the Alcator C-Mod Project and the EAST /
KSTAR Scenario Extension and Control Collaboration.
![Page 4: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/4.jpg)
LHCD in the EAST Tokamak
• LHCD experiments in EAST are in the weak damping
regime where full-wave effects and interference effects
can potentially be very important [4]:
Ray tracing / Fokker Planck
simulations show multi-pass
nature of absorption:
[2] C. Yang et al, PPCF (2014)
![Page 5: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/5.jpg)
Simulation set-up
• Use experimental profiles and EFIT equilibrium reconstruction for discharge 048888.05500:
– Created Fourier MHD equilibrium representation for TorLH (equigs.dat) and kinetic profiles data file (equidt.data):
– B0 = 2.31 T, Ip = 373 kA, a = 0.42 m, R0 = 1.85 m, Te(0) = 3.2 keV, ne(0) = 3.6 1019 m-3.
• RF parameters:
– PLH = 2 MW, f0 = 4.6 GHz, n// = -1.80
• Status of TorLH / CQL3D executables (10-2015):
– NERSC: TorLH and CQL3D available.
– MIT Loki Cluster: TorLH and CQL3D avaliable.
– IPP Shenma Cluster: TorLH and CQL3D are available (PGI build).
![Page 6: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/6.jpg)
• Semi-spectral ansatz is assumed for the electric field:
,
,
( ) ( ) im in
m n
m n
E x E e
– Spectral decomposition in the poloidal (m) and toroidal (n)
directions.
– Em,n() are represented by finite elements in the radial direction
(cubic Hermite interpolating polynomials).
• Using the ansatz above, the wave equation can be put
in a weak variational form (Galerkin method):
– Each toroidal mode (n) is solved separately assuming Nm
poloidal modes and Nr radial elements.
– This results in a block tri-diagonal matrix to invert.
Numerical implementation and mode
resolution requirements for TorLH
![Page 7: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/7.jpg)
Matrix inversion is performed using a new “3D”
parallel solver [5]
A
D U 0 0 0
L D U 0 0
0 L D U 0
0 0 L D U
0 0 0 L D
• Parallel decomposition of blocks is distributed over p2 p3 =
pcblock processors.
• Solver distributes groups of rows over p1 processors.
• Matrix inversion time scales like Nr (Nm)3 using p1 p2 p3
cores.
L, D, U are each dense
block matrices of size (2
3 Nm)2
Total of Nr block rows
![Page 8: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/8.jpg)
• Must resolve the shortest perpendicular wavelength in the system, which is given by the LH dispersion relation:
Poloidal mode resolution requirements for TorLH to
simulate LH wave propagation in EAST
/ /
pe mk k
r
• For an EAST discharge at r/a ~ 0.5 with B0 = 2.3 T, Te(0) ~ 1 keV, ne(0) ~ 2.0 1019 m-3, f0 = 4.6 GHz, we have:
1
/ / / /0
1
/ / / /
17 (for ~n =2) 372
and 2 1 745
46 (for ~n ~5.5) 1019
and 2 1 2038
m
ELD
m
k cm n m rk
N m
k cm n m rk
N m
• May need Nm ~ 4000 to resolve LH wave at edge ~ 45 cm.
![Page 9: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/9.jpg)
Convergence is only achieved on innermost flux
surface (r/a 0.1) at Nm = 255
![Page 10: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/10.jpg)
Convergence still limited to flux surfaces at r/a < 0.3
as Nm is increased from 255 511
![Page 11: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/11.jpg)
Convergence is finally achieved out to r/a 0.3 as Nm
is increased to 1023
![Page 12: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/12.jpg)
Convergence continues to improve on flux surfaces
out to r/a 0.5 as Nm is increased to 2047
• Simulation performed using 2.5 hours of wall clock time on Hopper at NERSC using 7168 cores.
![Page 13: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/13.jpg)
Find that convergence is achieved on all flux
surfaces as Nm is increased to 4095
• Simulation required ~7.5 hours of wall clock time on Hopper platform
at NERSC using 32,256 cores:– 1.5 hours for matrix inversion and ~ 6 hours for power reconstruction.
![Page 14: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/14.jpg)
• Broadening of the LH power deposition profile is consistent with adding higher k// components to the spectral solution in TorLH that correspond to higher m.
LH power deposition profile is peaked on-axis at Nm
= 255, but starts to broaden as Nm is increased to 511
![Page 15: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/15.jpg)
• Broadening of the LH power deposition profile is consistent with adding higher k// components to the spectral solution in TorLH that correspond to higher m.
Off-axis peak in LH power deposition profile
appears as Nm is increased from 511 1023 2047
![Page 16: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/16.jpg)
Parallel electric field of LH wave starts to exhibit
features of multiple reflections and weak single pass
absorption as Nm is increased
Nm = 255 Nm = 1023 Nm = 2047
![Page 17: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/17.jpg)
At Nm = 4095 the full-wave absorption profile is
clearly off-axis (for Maxwellian electron damping).
Nm = 4095
![Page 18: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/18.jpg)
• New 3D solver [5] implemented in power reconstruction
subroutine of TorLH to recover E from solution:
– Results in a decrease in wall clock time for the 4095 mode simulation on
Hopper from 7.5 hours to ~ 1 hour.
– Time for power reconstruction reduced from ~ 6.5 hours to ~ 500 sec.
• Necessary step for performing iterations between TorLH and
CQL3D.
• May use vector extrapolation technique to improve
convergence of iteration (see [1]), especially in weak damping
regimes.
• Final step: Compare hard x-ray profiles and LH power
deposition and current density profiles from TorLH / CQL3D
and GENRAY / CQL3D.
Dramatic reduction achieved in wall
clock time for TorLH
![Page 19: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/19.jpg)
Iteration Procedure for Simulation
a) Execute TorLH in “toric” mode using Maxwellian
electron Landau damping (ELD):
i. Perform a resolution scan to determine how many poloidal
modes are needed to resolve the LH wave in EAST.
b) Re-run TorLH in “qldce” mode to compute the RF
diffusion coefficients (D_ql) from the electric field
solutions computed in Step (a):
i. Remap D_ql from the TorLH (radial, velocity) space mesh to the
CQL3D (radial / velocity) space mesh.
c) Run CQL3D to obtain first iterate for the quasilinear
electron distribution fe(v, v//, r):
i. Create look-up table for Im{zz} due to ELD.
d) Repeat steps (a) – (c) until fe(v, v//, r) and D_ql (fe) are
self-consistent.
![Page 20: Full-wave Simulations of Lower Hybrid Wave …...57th Annual Meeting of the APS Division of Plasma Physics Savannah, GA November 18, 2015 Poster NP12.00050 Abstract Studies of lower](https://reader035.vdocument.in/reader035/viewer/2022071020/5fd3da6bab36ac577c3ba337/html5/thumbnails/20.jpg)
[1] J. C. Wright et al, Physics of Plasmas 16, 072502 (2009).
[2] C. Yang et al, Plasma Physics and Controlled Fusion 56
125003 (2014).
[3] S. Shiraiwa et al, 21st Topical Conference on Radio-frequency
Power in Plasmas, 27-29 April 2015, Lake ArrowHead, CA.
[4] J. C. Wright et al, Plasma Physics and Controlled Fusion 56,
045007 (2014).
[5] J. P. Lee and J. C. Wright, Computer Physics
Communications 185, 2598 (2014).
Work supported by the US DOE under Contract No. DE-
SC0010492 and DE-FC02-01ER54648.
References and Acknowledgements