Transcript 杨程 - 等离子体物理理论与数值模拟研究室

低杂波加热的射线追踪以及全
波解模拟
杨程
等离子体所五室
OUTLINE
• 理论以及数值模拟方法介绍
1）射线追踪（GENRAY）
2）全波解（TORLH）
• 两种方法对一个物理问题——低杂波的
spectral gap的研究和比较
理论以及数值模拟方法介绍——
射线追踪方法
Ray tracing
Ray tracing is based on WKB approximation: if the wavelength is short
compared to the scale length of tokamaks, the wave quantities such as
electric field is taken to be wary as


E 0 ( x ) exp  i  k x ( x ) dx   ik y y  ik z z  i t 
 x

Here the x axis is taken to be along the direction of the inhomogeneity. The
wave propagates with fixed ky, kz, and ω and with kx(x) determined through
the local dispersion relation, the parameters of which are functions of x.
The WKB method can be extended to three dimensions: E 0 exp( i  i t ) , where
is a function of position such that    k , with k satisfying the local dispersion
relation.
The direction of energy flow is given by the group velocity
vg   / k
The vector k is determined by a dispersion relation of the form:
and the fact that along the ray path:
D ( , k , r )  0
dD / dt  (  D /  k )( dk / dt )  (  D /  r )( dr / dt )  0
The ray tracing equations can be given by:
dr / dt   (  D /  k ) /(  D /   )
dk / dt  (  D /  r ) /(  D /   )
.
Along the ray path, based on the requirement of conservation of energy flux
and weak attenuation approximation, the evolution of the electrostatic wave
power is given by
𝑑𝑃
=−2𝛾𝑡 P, 𝛾𝑡 = 𝛾𝑒 + 𝛾𝑖 + 𝛾𝑐
𝑑𝑡
Where 𝛾𝑡 is calculated by writing:
Expending D for small 𝛾𝑡 and 𝐷(𝛼) with 𝐷0 (𝒙, 𝒌, 𝜔)=0, yields:
The expressions for 𝐷 (𝛼) are given by [*]
[*]
理论以及数值模拟方法介绍——
全波解方法
TORLH
The TORLH[J. C. Wright, 2009] is a full wave electromagnetic code, which solves a wave
equation of the form:
c2
𝛻
ω2
× 𝛻 × E = SE⊥ + iD(b ×E⊥ ) + PE∥ +
4πi 2
J
ω e
λ0 =
=
c2
− ω2 𝛻⊥
× λ0 𝛻⊥ × E⊥
1 ω2pe v2te
[−xe Z
2 Ω2ce c2
xe ],
+
4πi
2
(
J
e
ω
c2
i ω2 {𝛻⊥
ξ0 =
+ JANT ) ,
× [ξ0 b ∙ 𝛻 E∥ b + b b ⋅ 𝛻 ξ0 b ∙ 𝛻⊥ × E⊥ } ,
1 ω2pe v2te 2 ′
[xe Z
2 ωΩce c2
xe ]
where S, D, and P, are the Stix components of the dielectric tensor, expanded in the LH range of
frequencies (Ω𝑐𝑖 ≪ 𝜔 ≪ Ω𝑐𝑒 )
A semi‐spectral ansatz is used in TORLH to represent the electric field in terms of poloidal (m)
and toroidal (n) Fourier modes:
Em,n (Ψ)eimθ+inϕ
E x =
m,n
The radial dependence of the coefficients in the above equation is expressed in terms of cubic
Hermite interpolating polynomials. In the full wave simulations carried out for EAST we have
taken 2047 poloidal modes (−1023 < 𝑚 < +1023) and 480 radial elements in order to
resolve the short wavelength electrostatic LH wave (𝜆⊥ ∼ 1 − 3𝑚𝑚).
Spectral gap problem
Quasilinear calculations[*] have shown that damping of LH waves begins at a phase
velocity of 𝜔 𝑘∥𝑣𝑡𝑒 = 2.5, which can be restated as 𝑛∥ = 5.7 𝑇𝑒 (𝑘𝑒𝑉) . Current drivie
efficiency scales[#] as 𝑛𝑒 /𝑛∥ , thus the accessibility limit sets the maximum current drive
effciency. If the electron temperature is not sufficiently high, it requires high 𝑛∥ to have
sufficient Landau damping. The difference between the launched 𝑛∥ ,which tends to be
just above the accessibility limit, and the slower wave velocity (require higher 𝑛∥ ) at
which damping occurs is known as the spectral gap.
[*]P.T. Bonoli, IEEE Trans. Plasma Sci. 12,95(1984)
[#]N.J. Fisch, PRL, 41,873(1978)
Note that in toroidal geometry the poloidal mode number of LH wave is not a
conserved quantity and can vary considerably
𝑘∥ = 𝒌 ∙ 𝒃 = [
𝑚
𝑟
n
R
𝑩𝜽 + ( ) 𝑩𝝓 ]/ 𝑩
Toroidal effects move the value of 𝑘∥ up and down during propagation: A rise
helping absorption, and a fall retarding absorption.
The rise in 𝑘∥ is a stronger effect as density, poloidal field inside the plasma
increase and as aspect ratio (R/a) decreases [*]
Motivations of the following works:
The spectral gap for LHW in EAST is expected to be very large, and the LHW is in a
low single pass absorption regime, which poses a challenge for the simulation of
LHCD in EAST device with large aspect ratio.
[*] K. Kupfer and D. Moreau, \Wave chaos and the dependence of LHCD eciency on temperature," Nucl. Fusion 32 1845 (1992).
H. Takahashi, et al., \Accessibility for lower hybrid waves in PBX{M," in 1993 Conference on Controlled Fusion and Plasma
Physics (Proc. Conf. Lisbon, 1994) Vol 17C, Part III, European Physical Society, Geneva (1993) 901.
F. Paoletti, et al., \LHCD Accessibility Study with Reconstructed Magnetic Equilibria in PBX-M," Nucl. Fusion 34 771-776 (1994).
Study of the Spectral gap problem in EAST LHCD using ray-tracing (GENRAY)
LHW(2.45GHz, 𝑛∥ ~2.1 ) power deposition profile
modeled by the GENRAY with the parameters: Te
=1.5keV, 𝑛𝑒 = 1.0 × 1019 𝑚−3
LHW undergoes multiple passes before toroidal effects give rise to sufficient upshift in the parallel wavenumber
Simulation results shows that the spectral gap in EAST LHCD is large and the Landau damping
relies very much on the up-shift in the parallel wavenumber due to the toroidal effects, since
the power is not able to be absorbed on-axis where the temperature is higher but the poloidal
field ~0.
Comparisons between the GENRAY and the TORLH simulation for the LHW in EAST for
the Maxwellian case
Maxwellian electron distribution function assumed for both the GENRAY and the TORLH simulation:
S p e c tr u m o f R e ( E ) v e r s u s m
~ r/a
0 .1 0
0 .3 0
0 .5 0
0 .7 0
0 .9 0
1 .0 0
0
|F F T ( R e ( E ) ) |
10
10
2
10
4
10
6
8
10
1000
500
0
500
1 00 0
P o lo id a l m o d e n u m b e r - m
Poloidal spectrum modeled by TORLH
Parameters:
𝑛∥ = 2.1, 𝑇𝑒 = 1.5𝑘𝑒𝑉,
𝑛𝑒 = 1.0 × 1019
The power deposition related to the up-shift of 𝑛∥ due to the toroidal effect calculated by the GENRAY and the TORLH agree
well with each other. Besides, the full wave effects play another important role in the 𝑛∥ up-shifting.
Comparisons between the GENRAY and the TORLH with different temperature
1.5 keV
5.0 keV
电子温度越高，低杂波spectral gap越小，吸收加强，朗道阻尼对𝑘∥ 上移的依赖性越小，采用
不同的最大极向模数的全波结果就越接近，全波结果和射线追踪结果也越接近。
Comparisons between the GENRAY and the TORLH simulation for the LHW in EAST for
the Non-Maxwellian case
Using the same Non-Maxwellian electron distribution
function provided by the iteration between GENRAY and
CQL3D for both the GENRAY and the TORLH simulation:
Maxwellian results of wave
electric field modeled by TORLH:
Quasi-linear results of wave electric
field modeled by TORLH
For the non-Maxwellian case, the two deposition profiles look similar to
each other except for a slight difference between the positions of the
peaks, which indicates that the discrepancy of the two codes is diminished
as the results of the iteration between the GENRAY and the CQL3D being
used.
Summary
•
The spectral gap for LHCD in EAST has been studied and proved to be very large as
expected. The electron Landau damping of the LHW relies very much on the upshift in the parallel wavenumber, and the LHW is in a low single pass absorption
regime, which poses a challenge for the simulation of LHCD in EAST.
•
The GENRAY and TORLH agree well with each other on the calculation of the upshift of 𝑛∥ due to the toroidal effect. The full wave effects also play important role
in up-shifting in the parallel wavenumber with sufficient maximum poloidal mode
numbers. But the discrepancy of the two codes is diminished in the quasi-linear
case as the damping becomes stronger.
•
Simulations using Genray-CQL3D and TORLH models for LHCD in EAST would be
an important addition to the validation studies of those models being used in low
single pass absorption scenarios of tokamaks with large aspect ratio.