Source author record

Z. X. Lu

Z. X. Lu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.plasm-ph nucl-th physics.comp-ph

Catalog footprint

What is connected

5works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

The development of an implicit full f method for electromagnetic particle simulations of Alfvén waves and energetic particle physics

In this work, an implicit scheme for particle-in-cell/Fourier electromagnetic simulations is developed and applied to studies of Alfvén waves in one dimension and three-dimensional tokamak plasmas. An analytical treatment is introduced to achieve efficient convergence of the iterative solution of the implicit field-particle system. First, its application to the one-dimensional uniform plasma demonstrates its applicability in a broad range of $β/m_e$ values. Second, toroidicity induced Alfvén eigenmodes (TAE) are simulated in a three-dimensional axisymmetric tokamak plasma, using the widely studied case defined by the International Tokamak Physics Activity (ITPA) Energetic Particle (EP) Topical Group. The real frequency and the growth (or damping) rate of the TAE with (or without) EPs agree with previous results reasonably well. The full f electromagnetic particle scheme established in this work provides a possible natural choice for EP transport studies where large profile variation and arbitrary particle distribution functions need to be treated in kinetic simulations.

preprint2020arXiv

Development and testing of an unstructured mesh method for whole plasma gyrokinetic simulations in realistic tokamak geometry

In this work, we have formulated and implemented a mixed unstructured mesh-based finite element (FE)-Fourier decomposition scheme for gyrokinetic simulations in realistic tokamak geometry. An efficient particle positioning (particle-triangle mapping) scheme for the charge deposition and field scattering using an intermediate grid as the search index for triangles has been implemented and a significant speed-up by a factor of $\sim30$ is observed as compared with the brute force scheme for a medium-size simulation. The TRIMEG (TRIangular MEsh based Gyrokinetic) code has been developed. As an application, the ion temperature gradient (ITG) mode is simulated using the simplified gyrokinetic Vlasov-Poisson model. Our simulation and that using the ORB5 code for the DIII-D Cyclone case show reasonable agreement. As an additional application, ITG simulations using an ASDEX Upgrade equilibrium have been performed with density and temperature gradient profiles similar to the Cyclone case. Capabilities of the TRIMEG code for simulations with realistic experimental equilibria in the plasma core and in the whole plasma volume with open field lines are demonstrated.

preprint2019arXiv

ORB5: a global electromagnetic gyrokinetic code using the PIC approach in toroidal geometry

This paper presents the current state of the global gyrokinetic code ORB5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The ORB5 code solves the electromagnetic Vlasov-Maxwell system of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization density and drift-kinetic electrons. Variants of the physical model can be selected for electrons such as assuming an adiabatic response or a ``hybrid'' model in which passing electrons are assumed adiabatic and trapped electrons are drift-kinetic. A Fourier filter as well as various control variates and noise reduction techniques enable simulations with good signal-to-noise ratios at a limited numerical cost. They are completed with different momentum and zonal flow-conserving heat sources allowing for temperature-gradient and flux-driven simulations. The code, which runs on both CPUs and GPUs, is well benchmarked against other similar codes and analytical predictions, and shows good scalability up to thousands of nodes.

preprint2014arXiv

Particle simulation of lower hybrid waves in tokamak plasmas

Global particle simulations of the lower hybrid waves have been carried out using fully kinetic ions and drift kinetic electrons with a realistic electron-to-ion mass ratio. The lower hybrid wave frequency, mode structure, and electron Landau damping from the electrostatic simulations agree very well with the analytic theory. Linear simulation of the propagation of a lower hybrid wave-packet in the toroidal geometry shows that the wave propagates faster in the high field side than the low field side, in agreement with a ray tracing calculation. Electromagnetic benchmarks of lower hybrid wave dispersion relation are also carried out. Electromagnetic mode conversion are observed in toroidal geometry, slow waves are launched at the plasma boundary and converts to fast waves at the mode conversion layer, which is consistent with linear theory.

preprint2011arXiv

Theoretical and numerical studies of wave-packet propagation in tokamak plasmas

Theoretical and numerical studies of wave-packet propagation are presented to analyze the time varying 2D mode structures of electrostatic fluctuations in tokamak plasmas, using general flux coordinates. Instead of solving the 2D wave equations directly, the solution of the initial value problem is used to obtain the 2D mode structure, following the propagation of wave-packets generated by a source and reconstructing the time varying field. As application, the 2D WKB method is applied to investigate the shaping effects (elongation and triangularity) of tokamak geometry on the lower hybrid wave propagation and absorbtion. Meanwhile, the Mode Structure Decomposition (MSD) method is used to handle the boundary conditions and simplify the 2D problem to two nested 1D problems. The MSD method is related to that discussed earlier by Zonca and Chen [Phys. Fluids B 5, 3668 (1993)], and reduces to the well-known "ballooning formalism" [J. W. Connor, R. J. Hastie, and J. B. Taylor, Phys. Rev. Lett. 40, 396 (1978)], when spatial scale separation applies. This method is used to investigate the time varying 2D electrostatic ITG mode structure with a mixed WKB-full-wave technique. The time varying field pattern is reconstructed and the time asymptotic structure of the wave-packet propagation gives the 2D eigenmode and the corresponding eigenvalue. As a general approach to investigate 2D mode structures in tokamak plasmas, our method also applies for electromagnetic waves with general source/sink terms, either by an internal/external antenna or nonlinear wave interaction with zonal structures.