Related papers: Vlasov versus N-body: the H\'enon sphere
Numerical simulations of self-gravitating systems are generally based on N-body codes, which solve the equations of motion of a large number of interacting particles. This approach suffers from poor statistical sampling in regions of low…
This work discusses the main analogies and differences between the deterministic approach underlying most cosmological N-body simulations and the probabilistic interpretation of the problem that is often considered in mathematics and…
We revisit in one dimension the waterbag method to solve numerically Vlasov-Poisson equations. In this approach, the phase-space distribution function $f(x,v)$ is initially sampled by an ensemble of patches, the waterbags, where $f$ is…
N-body simulations are essential for understanding the formation and evolution of structure in the Universe. However, the discrete nature of these simulations affects their accuracy when modelling collisionless systems. We introduce a new…
The effects of discreteness arising from the use of the N-body method on the accuracy of simulations of cosmological structure formation are not currently well understood. After a discussion of how the relevant discretisation parameters…
One of the fundamental aspects of statistical behaviour in many-body systems is exponential divergence of neighbouring orbits, which is often discussed in terms of Liapounov exponents. Here we study this topic for the classical…
In the mean field limit, isolated gravitational systems often evolve towards a steady state through a violent relaxation phase. One question is to understand the nature of this relaxation phase, in particular the role of radial…
N-body simulations were conducted of pairs of galaxies with a 3:1 mass ratio on parabolic orbits, in order to quantify the effect of dynamical friction. The effects of varying the ratio of the dark matter halo size to the distance of…
Using a generalized self-similar secondary infall model, which accounts for tidal torques acting on the halo, we analyze the velocity profiles of halos in order to gain intuition for N-body simulation results. We analytically calculate the…
The dynamics of sheared inelastic-hard-sphere systems are studied using non-equilibrium molecular dynamics simulations and direct simulation Monte Carlo. In the molecular dynamics simulations Lees-Edwards boundary conditions are used to…
We perform a quantitative comparison between N-body simulations and the Schr\"odinger-Poisson system in 1+1 dimensions. In particular, we study halo formation with different initial conditions. We observe the convergence of various…
The kinetic analyses of many-particle soft matter often employ many simulation studies of various physical phenomena which supplement the experimental limitations or compliment the theoretical findings of the study. Such simulations are…
Validity of fluid models breaks down for non-thermal or weakly collisional plasmas which often occur e.g. in the solar wind. In these regimes one has to resort to modelling through the first-principle Vlasov-Maxwell system, but its…
I propose a modification of the spherical infall model for the evolution of density fluctuations with initially Gaussian probability distribution and scale-free power spectra. I introduce a generalized form of the initial density…
Much of standard galaxy dynamics rests on the implicit assumption that the corresponding N-body problem is (near) integrable. This notion although leading to great simplification is by no means a fact. It is therefore important to develop…
We present the problematic of controlling the discreteness effects in cosmological N-body simulations. We describe a perturbative treatment which gives an approximation describing the evolution under self-gravity of a lattice perturbed from…
I present empirical measurements of the rate of relaxation in N-body simulations of stable spherical systems and distinguish two separate types of relaxation: energy diffusion that is largely independent of particle mass, and energy…
The determination of the resolution of cosmological N-body simulations, i.e., the range of scales in which quantities measured in them represent accurately the continuum limit, is an important open question. We address it here using…
Understanding the nature of dark matter is among the top priorities of modern physics. However, due to its inertness, detecting and studying it directly in terrestrial experiments is extremely challenging. Numerical N-body simulations…
A parallelized implementation of the Vlasov-Hybrid method [Nunn, 1993] is presented. This method is a hybrid between a gridded Eulerian description and Lagrangian meta-particles. Unlike the Particle-in-Cell method [Dawson, 1983] which…