ColDICE: a parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation

Abstract

Resolving numerically Vlasov-Poisson equations for initially cold systems can be reduced to following the evolution of a three-dimensional sheet evolving in six-dimensional phase-space. We describe a public parallel numerical algorithm consisting in representing the phase-space sheet with a conforming, self-adaptive simplicial tessellation of which the vertices follow the Lagrangian equations of motion. The algorithm is implemented both in six- and four-dimensional phase-space. Refinement of the tessellation mesh is performed using the bisection method and a local representation of the phase-space sheet at second order relying on additional tracers created when needed at runtime. In order to preserve in the best way the Hamiltonian nature of the system, refinement is anisotropic and constrained by measurements of local Poincaré invariants. Resolution of Poisson equation is performed using the fast Fourier method on a regular rectangular grid, similarly to particle in cells codes. To compute the density projected onto this grid, the intersection of the tessellation and the grid is calculated using the method of Franklin and Kankanhalli (1993) generalised to linear order. As preliminary tests of the code, we study in four dimensional phase-space the evolution of an initially small patch in a chaotic potential and the cosmological collapse of a fluctuation composed of two sinusoidal waves. We also perform a "warm" dark matter simulation in six-dimensional phase-space that we use to check the parallel scaling of the code.

Code and illustration movies available at: http://www.vlasix.org/index.php?n=Main.ColDICE - Article submitted to Journal of Computational Physics

References in corpus (6)

• N-body simulations of gravitational dynamics
• Simulation techniques for cosmological simulations
• An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics
• Vlasov versus N-body: the Hénon sphere
• Dark matter halo's and self similarity
• Comment on 'Discreteness Effects in Simulations of Hot/Warm Dark Matter' by J. Wang & S.D.M. White

Cited by in corpus (44)

• Large-scale dark matter simulations
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Reducing Noise in Cosmological N-body Simulations with Neutrinos
• Solving the Vlasov equation in two spatial dimensions with the Schrödinger method
• The Generation of Vorticity in Cosmological N-body Simulations
• Post-collapse perturbation theory in 1D cosmology -- beyond shell-crossing
• Phase-space structure of protohalos: Vlasov versus Particle-Mesh
• Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms
• Lagrangian cosmological perturbation theory at shell-crossing
• Spurious Small-Scale Structure & Discreteness-Driven Relaxation in Cosmological Simulations
• The median density of the Universe
• Phase-space structure of cold dark matter halos inside splashback: multi-stream flows and self-similar solution
• Cosmological Vlasov-Poisson equations for dark matter: Recent developments and connections to selected plasma problems
• Perturbation theory with dispersion and higher cumulants: framework and linear theory
• Long-lived transient structure in collisionless self-gravitating systems
• Finding closure: approximating Vlasov-Poisson using finitely generated cumulants
• Topology and geometry of the dark matter web: a multistream view
• SCALAR: an AMR code to simulate axion-like dark matter models
• Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations
• Cold dark matter protohalo structure around collapse: Lagrangian cosmological perturbation theory versus Vlasov simulations
• A "metric" semi-Lagrangian Vlasov-Poisson solver
• The features of the Cosmic Web unveiled by the flip-flop field
• Particle number dependence in the non-linear evolution of N-body self-gravitating systems
• Non-Halo Structures and their Effects on Gravitational Lensing
• Phase-space structure analysis of self-gravitating collisionless spherical systems
• Tetrahedral collapse: a rotational toy model of simultaneous dark-matter halo, filament and wall formation
• Fast and accurate collapse-time predictions for collisionless matter
• The gravitational force field of proto-pancakes
• Universal multi-stream radial structures of cold dark matter halos
• Collisionless Boltzmann equation approach for the study of stellar discs within barred galaxies
• Starting Cosmological Simulations from the Big Bang
• On the probabilistic approach to the N-body problem
• Statistical exploration of halo anisotropic clustering and intrinsic alignments with the mass-Peak Patch algorithm
• Discreteness effects in $N$-body simulations of warm dark matter
• Axion Structure Formation II: The Wrath of Collapse
• The Boosted Potential
• A Mean-Field Approach to Simulating the Merging of Collisionless Stellar Systems Using a Particle-Based Method
• Dark matter halo dynamics in 2D Vlasov Simulations: a self-similar approach
• Excursion Sets with a "Perfect" Collapse Model
• Simulating the complexity of the dark matter sheet II: halo and subhalo mass functions for non-cold dark matter models
• Comparison of simulations and semi-analytical model for WDM subhalo mass functions
• Post-collapse Lagrangian perturbation theory in three dimensions
• Non-halo structures and their effects on gravitationally lensed galaxies
• Tidal adaptive softening and artificial fragmentation in cosmological simulations