Search:

New !

Codes and methods

Theory and analysis

Past events

A high-order adaptive scheme for six-dimensional convection simulations

Author of the code: Erwan Deriaz
Publication introducing and exploiting the code:

Introduction  Figure 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-3V Figure 2: Initial conditions for the interaction of two Plummer models. From Deriaz & Peirani (2018).   Figure 3: Three dimensional view in physical space of the collision of two Plummer spheres (IFRIT 3D Data Visualization). From Deriaz & Peirani (2018). Movie: anim40.mp4

Here, 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.

Conclusion

We 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

1. Yoshikawa K., Yoshida, N., Umemura, M., 2013, ApJ, 762, 116, 2013.
2. Springel V., 2005, MNRAS, 364, 1105
3. Springel V., Yoshida N., White S.D.M., 2001, NewA, 6, 51         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). Figure 5: Conservation analysis: plots of the mass, entropy and L2-norm. From Deriaz & Peirani (2018). Figure 6: Plot of the kinetic energy Ec obtained with the AMR code (solid line) compared to the one obtained in the `Gadget-2` simulation (dashed line). From Deriaz & Peirani (2018).                