Related papers: An adaptively refined phase-space element method f…
We describe a technique for solving the combined collisionless Boltzmann and Poisson equations in a discretised, or lattice, phase space. The time and the positions and velocities of `particles' take on integer values, and the forces are…
We develop a diffuse solid method that is versatile and accurate for modeling wetting and multiphase flows in highly complex geometries. In this scheme, we harness N + 1-component phase field models to investigate interface shapes and flow…
Modeling dispersed solid phases in fluids still represents a computational challenge when considering a small-scale coupling in wide systems, such as the atmosphere or industrial processes at high Reynolds numbers. A numerical method is…
In collisionless and weakly collisional plasmas, the particle distribution function is a rich tapestry of the underlying physics. However, actually leveraging the particle distribution function to understand the dynamics of a weakly…
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…
We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust,…
We introduce a deterministic discrete-particle simulation approach, the Linearly-Transformed Particle-In-Cell (LTPIC) method, that employs linear deformations of the particles to reduce the noise traditionally associated with particle…
Tidal dissipation plays an important role in the dynamical evolution of moons, planets, stars and compact remnants. The interesting complexity originates from the interplay between the internal structure and external tidal forcing. Recent…
A kinetic equation for the collisional evolution of stable, bound, self gravitating and slowly relaxing systems is established, which is valid when the number of constituents is very large. It accounts for the detailed dynamics and self…
In this paper the development of a physically consistent phase-field theory of solidification shrinkage is presented. The coarse-grained hydrodynamic equations are derived directly from the N-body Hamiltonian equations in the framework of…
The Vlasov-Poisson system describes interacting systems of collisionless particles. For solutions with small initial data in three dimensions it is known that the spatial density of particles decays like $t^{-3}$ at late times. In this…
Active Flux is a modified Finite Volume method that evolves additional Degrees of Freedom for each cell that are located on the interface by a non-conservative method to compute high-order approximations to the numerical fluxes through the…
Most approaches in Lagrangian fluid dynamics simulations proceed from the definition of particle volumes, from which discrete versions of the spatial differential operators are derived. Recently, Gallou\"et and M\'erigot [1] simultaneously…
We present a framework for the simulation of rigid and deformable bodies in the presence of contact and friction. Our method is based on a non-smooth Newton iteration that solves the underlying nonlinear complementarity problems (NCPs)…
We present an algorithm for quickly generating multiple realizations of N-body simulations to be used, for example, for cosmological parameter estimation from surveys of large-scale structure. Our algorithm uses a new method to resample the…
We consider a perturbative approach to the Vlasov-Poisson system for cosmic structure formation that does not rely on any truncation of the momentum-cumulant hierarchy. The generally non-trivial linear solution is computed by solving a…
This work outlines a new multi-physics-compatible immersed rigid body method for Eulerian finite-volume simulations. To achieve this, rigid bodies are represented as a diffuse scalar field and an interface seeding method is employed to…
Real-time simulation of elastic structures is essential in many applications, from computer-guided surgical interventions to interactive design in mechanical engineering. The Finite Element Method is often used as the numerical method of…
This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the…
We have developed an adaptive multigrid code for solving the Poisson equation in gravitational simulations. Finer rectangular subgrids are adaptively created in locations where the density exceeds a local level-dependent threshold. We…