Explicit finite-difference and direct-simulation-MonteCarlo method for the dynamics of mixed Bose-condensate and cold-atom clouds

We present a new numerical method for studying the dynamics of quantum fluids composed of a Bose-Einstein condensate and a cloud of bosonic or fermionic atoms in a mean-field approximation. It combines an explicit time-marching algorithm, previously developed for Bose-Einstein condensates in a harmonic or optical-lattice potential, with a particle-in-cell MonteCarlo approach to the equation of motion for the one-body Wigner distribution function in the cold-atom cloud. The method is tested against known analytical results on the free expansion of a fermion cloud from a cylindrical harmonic trap and is validated by examining how the expansion of the fermionic cloud is affected by the simultaneous expansion of a condensate. We then present wholly original calculations on a condensate and a thermal cloud inside a harmonic well and a superposed optical lattice, by addressing the free expansion of the two components and their oscillations under an applied harmonic force. These results are discussed in the light of relevant theories and experiments.Comment: 33 pages, 13 figures, 1 tabl

EUROfusion-theory and advanced simulation coordination (E-TASC) : programme and the role of high performance computing

This paper is a written summary of an overview oral presentation given at the 1st Spanish Fusion High Performance Computer (HPC) Workshop that took place on the 27 November 2020 as an online event. Given that over the next few years ITER24 will move to its operation phase and the European-DEMO design will be significantly advanced, the EUROfusion consortium has initiated a coordination effort in theory and advanced simulation to address some of the challenges of the fusion research in Horizon EUROPE (2021-2027), i.e. the next EU Framework Programme for Research and Technological Development. This initiative has been called E-TASC, which stands for EUROfusion-Theory and Advanced Simulation Coordination. The general and guiding principles of E-TASC are summarized in this paper. In addition, an overview of the scientific results obtained in the pilot phase (2019-2020) of E-TASC are provided while highlighting the importance of the required progress in computational methods and HPC techniques. In the initial phase, five pilot theory and simulation tasks were initiated: towards a validated predictive capability of the low to high transition and pedestal physics; runaway electrons in tokamak disruptions in the presence of massive material injection; fast code for the calculation of neoclassical toroidal viscosity in stellarators and tokamaks; development of a neutral gas kinetics modular code; European edge and boundary code for reactor-relevant devices. In this paper, we report on recent progress made by each of these projects.Peer reviewe

Hierarchy of second order gyrokinetic Hamiltonian models for particle-in-cell codes

International audienceThe reduced-particle model is the central element for the systematic derivation of the gyrokinetic Vlasov-Maxwell equations from first principles. Coupled to the fields inside the gyrokinetic field-particle Lagrangian, the reduced-particle model defines polarization and magnetization effects appearing in the gyrokinetic Maxwell equations. It is also used for the reconstruction of the gyrokinetic Vlasov equation from the particle characteristics. Various representations of reduced-particle models are available according to the choice of the gyrokinetic phase space coordinates. In this paper, the Hamiltonian representation of the reduced particle dynamics at an order suitable for the implementation in particle-in-cell simulations is explicitly derived from the general reduction procedure. The second-order (with respect to the fluctuating electromagnetic fields), full Finite Larmor Radius (FLR) Hamiltonian gyrokinetic particle model as well as the second-order model suitable specifically for the long-wavelength approximation (i.e., containing up to the second-order FLR corrections), are derived and compared to the model recently implemented in the particle-in-cell code ORB5. We show that the same long-wavelength approximate equations can also be derived by taking the proper limit of the full FLR model

Variational Framework for Structure-Preserving Electromagnetic Particle-in-Cell Methods

In this article we apply a discrete action principle for the Vlasov-Maxwell equations in a structure-preserving particle-field discretization framework. In this framework the finite-dimensional electromagnetic potentials and fields are represented in a discrete de Rham sequence involving general finite element spaces, and the particle-field coupling is represented by a set of projection operators that commute with the differential operators. With a minimal number of assumptions which allow for a variety of finite elements and shape functions for the particles, we show that the resulting variational scheme has a general discrete Poisson structure and thus leads to a semi-discrete Hamiltonian system. By introducing discrete interior products we derive a second type of space discretization which is momentum preserving, based on the same finite elements and shape functions. We illustrate our method by applying it to spline finite elements, and to a new spectral discretization where the particle-field coupling relies on discrete Fourier transforms

Global full-f gyrokinetic simulations of plasma turbulence

Critical physical issues can be specifically tackled with the global full-f gyrokinetic code GYSELA. Three main results are presented. First, the self-consistent treatment of equilibrium and fluctuations highlights the competition between two compensation mechanisms for the curvature driven vertical charge separation, namely, parallel flow and polarization. The impact of the latter on the turbulent transport is discussed. In the non-linear regime, the benchmark with the Particle-In-Cell code ORB5 looks satisfactory. Second, the transport scaling with p. is found to depend both on p. itself and on the distance to the linear threshold. Finally, a statistical steady-state turbulent regime is achieved in a reduced version of GYSELA by prescribing a constant heat source