Related papers: Some numerical aspects of the conservative PSM sch…
Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in physical space. The…
Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new,…
A model and code (TIFF) for isothermal gyrofluid computation of quasi-two-dimensional interchange and drift wave turbulence in magnetized plasmas with arbitrary fluctuation amplitudes (full-f) and arbitrary polarization wavelengths (full-k)…
Experience with core-collapse supernova simulations shows that accurate accounting of total particle number and 4-momentum can be a challenge for computational radiative transfer. This accurate accounting would be facilitated by the use of…
First-principles simulations of tokamak turbulence have proven to be of great value in recent decades. We develop a pseudo-spectral velocity formulation of the turbulence equations that smoothly interpolates between the highly efficient but…
Feedback stabilization of magnetohydrodynamic (MHD) modes is studied in a cylindrical model for a tokamak with resistivity, viscosity and toroidal rotation. The control is based on a linear combination of the normal and tangential…
We present a computational study for a family of discontinuous Galerkin methods for the one dimensional Vlasov-Poisson system that has been recently introduced. We introduce a slight modification of the methods to allow for feasible…
The connection between the poloidal and the toroidal rotation of plasma in tokamak is important for the high confinement regimes, in particular in reactor regime. The sudden onset of closed convection structures in the poloidal section, due…
Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct…
A computationally accurate and efficient numerical method under a unified framework is crucial to various multi-scale scientific and engineering problems. So far, many numerical methods have encountered various challenges in efficiently…
Recent development of structure-preserving geometric particle-in-cell (PIC) algorithms for Vlasov-Maxwell systems is summarized. With the arriving of 100 petaflop and exaflop computing power, it is now possible to carry out direct…
The Dimits shift, an upshift in the onset of turbulence from the linear instability threshold, caused by self-generated zonal flows, can greatly enhance the performance of magnetic confinement plasma devices. Except in simple cases, using…
The controlled ramp down of the toroidal plasma current in the ITER tokamak is simulated using a simple model that employs cylindrical geometry. The magnetohydrodynamical (MHD) stability of the plasma throughout the whole current ramp is…
Turbulence simulation codes can exploit the flute-like nature of plasma turbulence to reduce the effective number of degrees of freedom necessary to represent fluctuations. This can be achieved by employing magnetic coordinates of which one…
Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the SOL is highly collisional, an assumption that is valid in many present day experiments but is questionable in the high-power…
Most plasmas are only partially ionized. To better understand the dynamics of these plasmas, the behaviors of a mixture of neutral species and plasma in ideal magnetohydrodynamic states are investigated. The current approach is about the…
A paradigm based on the absolute equilibrium of Galerkin-truncated inviscid systems to aid in understanding turbulence [T.-D. Lee, "On some statistical properties of hydrodynamical and magnetohydrodynamical fields," Q. Appl. Math. 10, 69…
Low density plasmas are characterized by a large scale separation between the gyromotion of particles around local magnetic fields and the macroscopic scales of the system, often making global kinetic simulations computationally…
The pullback scheme is implemented in the global gyrokinetic particle-in-cell code ORB5 [S. Jolliet et al, Comp. Phys. Comm., 177, 409 (2007)] to mitigate the cancellation problem in electromagnetic simulations. The equations and the…
A high resolution fourth-order Pad\'e scheme is used to simulate locally isothermal 3D disk turbulence driven by the vertical shear instability (VSI) using 268.4 M points. In the early non-linear period of axisymmetric VSI, angular momentum…