Codes and methods
Theory and analysis Past events |
A high-order adaptive scheme for six-dimensional convection simulationsAuthor of the code: Erwan Deriaz Six-dimensional Adaptive Simulation of Vlasov Equations using a Hierarchical Basis, Deriaz E., Peirani S., 2018, Multiscale Modeling and Simulation 16, 583 (hal-01419750 preprint)
Table of ContentsIntroductionFigure 1: Refinement of a two dimensional grid: on the left the refinement process, on the right the resulting adaptive grid. Even though computer power increases exponentially with time, resolving Vlasov-Poisson equations in full six-dimensional phase-space remains challenging nowadays and requires special computing techniques. Here we explore the use of Adaptive Mesh Refinement (AMR), which allows for a substantial gain in terms of memory and computational time. In six dimensions, the operations should have a complexity independent of the number of dimensions. In order to avoid excessive dissipation of the solution, we use a third-order scheme. Regarding refinement structure, we propose to use an adaptive point-centered dyadic grid, as illustrated by Figure 1. Let us check the result of a straightforward recipe with the simplest numerical tools: Eulerian discretization, polynomial interpolation and finite differences for hyperbolic partial differential equations. Interaction of two Plummer spheres in 3D-3VHere, we perform a simple test, similarly as in Yoshikawa et al. (2013), consisting in simulating the collision between two Plummer spheres (Figure 2). When isolated, these systems should remain stable. Thanks to the AMR technique, we are able to simulate their interaction in full six-dimensional phase-space (Figures 3 and 4), with a reasonable level of accuracy (Figures 5 and 6) given the high dimensionality of the problem. ConclusionWe demonstrated the feasibility of simulations of Vlasov-Poisson systems in six-dimensional phase-space thanks to an innovative high-order adaptive scheme. It is characterized by a fourth-order interpolation in time, third-order in (phase-)space for the convection and second-order in space for the calculation of the force field. Thanks to its weaker CFL condition compared to standard Eulerian schemes with uniform grids, it can compete with semi-Lagrangian solvers. While the current version of the code is operational on shared memory architectures with a 80% efficient parallelization with OpenMP, it remains to be optimized to run on large configurations. The next steps are thus MPI parallelization in addition to the implementation of compact and conservative schemes in order to obtain substantial gains regarding the conservations shown on Figure 5. Bibliography
Figure 4: View in phase space of the collision of two Plummer spheres: a slice in the (z,w) plane of six-dimensional phase space at times t=15.2, t=21.6 and t=33.6. First row: the adaptive grid. Second row: the density distribution function obtained with the hierarchical basis adaptive scheme. Third row: the result obtained with a N-body simulation performed with the public code Gadget-2 (Springel, Yoshida & White 2001; Springel 2005). From Deriaz & Peirani (2018). |