Related papers: An adaptively refined phase-space element method f…
We calculate the cosmological evolution of the 1-point probability distribution function (PDF), using an analytic approximation that combines gravitational perturbation theory with the Edgeworth expansion of the PDF. Our method applies…
We present an energy-preserving mechanic formulation for dynamic quasi-brittle fracture in an Eulerian-Lagrangian formulation, where a second-order phase-field equation controls the damage evolution. The numerical formulation adapts in…
We introduce the N-body simulation technique to follow structure formation in linear and nonlinear regimes for the extended quintessence models (scalar-tensor theories in which the scalar field has a self-interaction potential and behaves…
We present a new formalism which allows to derive the general Lagrangian dynamical equations for the motion of gravitating particles in a non--flat Friedmann universe with arbitrary density parameter $\Omega$ and no cosmological constant.…
In these lecture notes, we describe the current state-of-the-art for numerical simulations of large-scale structure and galaxy formation. Numerical simulations play a central role in the preparation and the exploitation of large-scale…
The one-point probability distribution function (pdf) of the large-scale density field is an important tool to follow the evolution of cosmological structures. In this paper we present a new model for this pdf for all regimes and all…
A variational method is used to derive a self-consistent macro-particle model for relativistic electromagnetic kinetic plasma simulations. Extending earlier work [E. G. Evstatiev and B. A. Shadwick, J. Comput. Phys., vol. 245, pp. 376-398,…
In this paper we describe an adaptive softening length formalism for collisionless N-body and self-gravitating Smoothed Particle Hydrodynamics (SPH) calculations which conserves momentum and energy exactly. This means that spatially…
This paper focusses on the barely understood gap between the weakly nonlinear regime of structure formation and the onset of the virialized regime. While the former is accessed through perturbative calculations and the latter through…
Numerical simulations of the air in the atmosphere and water in the oceans are essential for numerical weather prediction. The state-of-the-art for performing these fluid simulations relies on an Eulerian viewpoint, in which the fluid…
The time irreversibility and fast relaxation of collapsing $N$-body gravitating systems (as opposed to the time reversibility of the equations of motion for individual stars or particles) are traditionally attributed to information loss due…
A novel method allowing to compute density, velocity and other fields in cosmological N--body simulations with unprecedentedly high spatial resolution is described. It is based on the tessellation of the three-dimensional manifold…
Particle-mesh simulations trade small-scale accuracy for speed compared to traditional, computationally expensive N-body codes in cosmological simulations. In this work, we show how a data-driven model could be used to learn an effective…
The evolution of closed gravitational systems is studied by means of $N$-body simulations. This, as well as being interesting in its own right, provides insight into the dynamical and statistical mechanical properties of gravitational…
The recently introduced relativistic Lagrangian darkon fluid model (EPJ C (2015) 75:9) is generalized to a self-gravitating, irrotational, pressure-less and stress free geodesic fluid, whose energy-momentum tensor is dust-like with…
$N$-body simulations study the dynamics of $N$ particles under the influence of mutual long-distant forces such as gravity. In practice, $N$-body codes will violate Newton's third law if they use either an approximate Poisson solver or…
This posting announces public availability of the GALAXY software package developed by the author over the past 40 years. It is a highly efficient code for the evolution of (almost) isolated, collisionless stellar systems, both disk-like…
A particle method for reproducing the phase space of collisionless stellar systems is described. The key idea originates in Liouville's theorem which states that the distribution function (DF) at time t can be derived from tracing necessary…
In this paper, we present a novel interface-driven adaptive variational procedure using a fully Eulerian description of fluid-structure interaction. The proposed fully-Eulerian procedure involves a fixed background unstructured mesh on…
This work deals with the numerical solution of the Vlasov equation. This equation gives a kinetic description of the evolution of a plasma, and is coupled with Poisson's equation for the computation of the self-consistent electric field.…