Related papers: ORB5: a global electromagnetic gyrokinetic code us…
Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently has it become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison…
The first three-dimensional, nonlinear gyrokinetic simulation of plasma turbulence resolving scales from the ion to electron gyroradius with a realistic mass ratio is presented, where all damping is provided by resolved physical mechanisms.…
We present the TRIDK code, which is written in Python, that solves Triad Reduced Interactions for a Drift Kinetic system of equations. The four-dimensional drift kinetic system of equations captures nonlinear electromagnetic (or…
We present a new, high-performance coupled electrodynamics-micromagnetics solver for full physical modeling of signals in microelectronic circuitry. The overall strategy couples a finite-difference time-domain (FDTD) approach for Maxwell's…
We describe an implementation of a particle physics module available for the PLUTO code, appropriate for the dynamical evolution of a plasma consisting of a thermal fluid and a non-thermal component represented by relativistic charged…
The paper presents a 1D2V electrostatic PIC model with a drift-kinetic description of all particle types aiming at simulating classical longitudinal plasma transport in axially symmetric open traps. The model generalizes the semi-implicit…
We construct a relativistic resistive magneto-hydrodynamic (RRMHD) numerical simulation code for high-energy heavy-ion collisions. We split the system of differential equations into two parts, a non-stiff and a stiff part. For the non-stiff…
Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code…
We study the electrokinetic flow of viscoelastic fluids subjected to an oscillatory pressure gradient, and particularly focus on the resonance behaviors in the flow. The governing equations are restricted to linear regime so that the…
We consider a three-dimensional kinetic model for a two species plasma consisting of electrons and ions confined by an external nonconstant magnetic field. Then we derive a kinetic-fluid model when the mass ratio $m_e/m_i$ tends to zero.…
Our contributions are motivated by fusion reactors that rely on maintaining magnetohydrodynamic (MHD) equilibrium, where the balance between plasma pressure and confining magnetic fields is required for stable operation. In axisymmetric…
We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a…
This paper proposes a charge-conserving, variational, spatio-temporal discretization for the drift-kinetic Vlasov-Maxwell system, utilizing finite-elements for the electromagnetic fields and the particle-in-cell approach for the Vlasov…
The hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell…
We study a finite Larmor radius model used to describe the ions distributions in the core of a toka-mak plasma, that consist in a gyro-kinetic transport equation, coupled with an electro-neutrality equation. Since the last equation do not…
Parallel-in-time methods, such as multigrid reduction-in-time (MGRIT) and Parareal, provide an attractive option for increasing concurrency when simulating time-dependent PDEs in modern high-performance computing environments. While these…
We investigate tearing modes (TM) driven by current density gradient in collisionless tokamak plasmas by using the electromagnetic gyrokinetic simulation code ORB5. We elucidate the TM width by simulations for flat profiles, as the absence…
A recent proof-of-principle study proposes a nonlinear electrostatic implicit particle-in-cell (PIC) algorithm in one dimension (Chen, Chacon, Barnes, J. Comput. Phys. 230 (2011) 7018). The algorithm employs a kinetically enslaved…
The simulation of problems in kinetic plasma physics are often challenging due to strongly coupled phenomena across multiple scales. In this work, we propose a wavelet-based coarse-grained numerical scheme, based on the framework of…
Controllable arrays of ions and ultra-cold atoms can simulate complex many-body phenomena and may provide insights into unsolved problems in modern science. To this end, experimentally feasible protocols for quantifying the buildup of…