Related papers: A Particle Module for the PLUTO Code: III -- Dust
The simulation of geometrically resolved rigid particles in a fluid relies on coupling algorithms to transfer momentum both ways between the particles and the fluid. In this article, the fluid flow is modeled with a parallel Lattice…
We develop an inertial coupling method for modeling the dynamics of point-like 'blob' particles immersed in an incompressible fluid, generalizing previous work for compressible fluids. The coupling consistently includes excess (positive or…
This paper exposes a novel exploratory formalism, which end goal is the numerical simulation of the dynamics of a cloud of particles weakly or strongly coupled with a turbulent fluid. Giventhe large panel of expertise of the list of…
A method to avoid the explicit time integration of small dust grains in the two fluid gas/dust smoothed particle hydrodynamics (SPH) approach is proposed. By assuming a very simple exponential decay model for the relative velocity between…
Approaches used in modern numerical simulations of the dynamics of dust and gas in circumstellar disks are tested. The gas and dust are treated like interpenetrating continuous media that can exchange momentum. A stiff coupling between the…
We present a new Eulerian framework for the computation of turbulent compressible multiphase channel flows, specifically to assess turbulence modulation by dispersed particulate matter in dilute concentrations but with significant mass…
We describe a simple method for simulating the dynamics of small grains in a dusty gas, relevant to micron-sized grains in the interstellar medium and grains of centimetre size and smaller in protoplanetary discs. The method involves…
The key role that dust plays in the interstellar medium has motivated the development of numerical codes designed to study the coupled evolution of dust and gas in systems such as turbulent molecular clouds and protoplanetary discs. Drift…
Piernik is a multi-fluid grid magnetohydrodynamic (MHD) code based on the Relaxing Total Variation Diminishing (RTVD) conservative scheme. The original code has been extended by addition of dust described within the particle approximation.…
We propose an efficient algorithm for the evaluation of the potential and its gradient of gravitational/electrostatic $N$-body systems, which we call particle mesh multipole method (PMMM or PM$^3$). PMMM can be understood both as an…
Context: Planet-disk interactions play a crucial role in the understanding of planet formation and disk evolution. There are multiple numerical tools available to simulate these interactions, including the often-used FARGO code and its…
Dust-gas mixtures are the simplest example of a two fluid mixture. We show that when simulating such mixtures with particles or with particles coupled to grids a problem arises due to the need to resolve a very small length scale when the…
Exascale supercomputing unleashes the potential for simulations of astrophysical systems with unprecedented resolution. Taking full advantage of this computing power requires the development of new algorithms and numerical methods that are…
We have developed a parallel Particle-Particle, Particle-Mesh (P3M) simulation code for the Cray T3E parallel supercomputer that is well suited to studying the time evolution of systems of particles interacting via gravity and gas forces in…
In a companion paper we have shown how the equations describing gas and dust as two fluids coupled by a drag term can be reformulated to describe the system as a single fluid mixture. Here we present a numerical implementation of the…
We develop a numerical a framework to study phoretic particle dynamics in two dimensions. The particles are modeled as chemically active rigid circles, which can emit or absorb a solute into surrounding fluid. The interaction between…
Traditional cosmological hydrodynamical simulations usually assume equal-numbered but unequal-mass dark matter and baryonic particles, which can lead to spurious collisional heating due to energy equipartition. To avoid such a numerical…
Particle-in-cell merging algorithms aim to resample dynamically the six-dimensional phase space occupied by particles without distorting substantially the physical description of the system. Whereas various approaches have been proposed in…
This is the first in a series of papers on the effects of dust on the formation, propagation, and structure of nonlinear MHD waves and MHD shocks in weakly-ionized plasmas. We model the plasma as a system of 9 interacting fluids, consisting…
The charge on an aggregate immersed in a plasma environment distributes itself over the aggregate's surface; this can be approximated theoretically by assuming a multipole distribution. The dipole-dipole (or higher order) charge…