The numerical resolution of the Vlasov equation is usually performed by particle methods (PIC) which consist in approximating the plasma by a finite number of particles. The trajectories of these particles are computed from the characteristic curves given by the Vlasov equation, whereas selfconsistent fields are computed on a mesh of the physical space. This method allows to obtain satisfying results with a relatively small number of particles. However, it is well known that, in some cases, the numerical noise inherent to the particle method becomes too important to have an accurate description of the distribution function in phase space. To remedy to this problem, methods discretizing the Vlasov equation on a mesh of phase space have been proposed [2,5,6].The major drawback of Vlasov methods using a uniform and fixed mesh is that their numerical cost is high, which makes them rather inefficient when the dimension of phase-space grows. For this reason we are investigating a method using an adaptive mesh. The adaptive method is overlayed to a classical semi-Lagrangian method which is based on the conservation of the distribution function along particle trajectories. The phase-space grid is updated using a multiresolution technique.The model we consider throughout this paper is the nonrelativistic Vlasov equation coupled self-consistently with Poisson's equation. It readsthe self electric field E is computed from Poisson's equationsThe magnetic field is external and considered to be known. In the present work, we have chosen to introduce a phase-space mesh which can be refined or derefined adaptively in time. For this purpose, we use a technique based on multiresolution analysis which is in the same spirit as the methods developed in particular by S. Bertoluzza [1], A. Cohen et al. [3] and M. Griebel and F. Koster [4]. We represent the distribution function on a wavelet basis at different scales. We can then compress it by eliminating coefficients which are small and accordingly remove the associated mesh points. Another specific feature of our method is that we use an advection in physical and velocity space forward in time to predict the useful grid points for the next time step, rather than restrict ourselves to the neighboring points. This enables us to use a much larger time step, as in the semi-Lagrangian method the time step is not limited by a Courant condition. Once the new mesh is predicted, the semi-Lagrangian methodology is used to compute the new values of the distribution function at the predicted mesh points, using an interpolation based on the wavelet decomposition of the old distribution function. The mesh is then refined again by performing a wavelet transform, and eliminating the points associated to small coefficients. This paper is organized as follows : we shall first present the context of multiresolution analysis for adaptive simulations, then describe the adaptive semi-Lagrangian method and validate our method on typical simulation in beam physics.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.