Related papers: Particle-Gas Dynamics with Athena: Method and Conv…
We describe the algorithm, implementation and numerical tests of a multifluid dust module in the Athena++ magnetohydrodynamic (MHD) code. The module can accommodate an arbitrary number of dust species interacting with the gas via…
The implementation of a new particle module describing the physics of dust grains coupled to the gas via drag forces is the subject of this work. The proposed particle-gas hybrid scheme has been designed to work in Cartesian as well as in…
We describe Pegasus, a new hybrid-kinetic particle-in-cell code tailored for the study of astrophysical plasma dynamics. The code incorporates an energy-conserving particle integrator into a stable, second-order--accurate, three-stage…
We present a new magnetohydrodynamic-particle-in-cell (MHD-PIC) code integrated into the Athena++ framework. It treats energetic particles as in conventional PIC codes while the rest of thermal plasmas are treated as background fluid…
(Abridged) We present local 2D and 3D hybrid numerical simulations of particles and gas in the midplane of protoplanetary disks (PPDs) using the Athena code. The particles are coupled to gas aerodynamically, with particle-to-gas feedback…
We investigate the utility of deep learning for modeling the clustering of particles that are aerodynamically coupled to turbulent fluids. Using a Lagrangian particle module within the Athena++ hydrodynamics code, we simulate the dynamics…
Numerical simulation of numerous mm/cm-sized particles embedded in a gaseous disk has become an important tool in the study of planet formation and in understanding the dust distribution in observed protoplanetary disks. However, the mutual…
We perform a systematic study of the dynamics of dust particles in protoplanetary disks with embedded planets using global 2-D and 3-D inviscid hydrodynamic simulations. Lagrangian particles have been implemented into magnetohydrodynamic…
Forming macroscopic solid bodies in circumstellar discs requires local dust concentration levels significantly higher than the mean. Interactions of the dust particles with the gas must serve to augment local particle densities, and…
We present local simulations that verify the linear streaming instability that arises from aerodynamic coupling between solids and gas in protoplanetary disks. This robust instability creates enhancements in the particle density in order to…
We describe a module for the Athena code that solves the gray equations of radiation hydrodynamics (RHD), based on the first two moments of the radiative transfer equation. We use a combination of explicit Godunov methods to advance the gas…
We present a stable and convergent method for studying a system of gas and dust, coupled through viscous drag in both non-stiff and stiff regimes. To account for the effects of dust drag in the update of the fluid quantities, we employ a…
We present an implementation of Smoothed Particle Hydrodynamics (SPH) in an adaptive-mesh PPPM algorithm. The code evolves a mixture of purely gravitational particles and gas particles. The code retains the desirable properties of previous…
A method for implementing cylindrical coordinates in the Athena magnetohydrodynamics (MHD) code is described. The extension follows the approach of Athena's original developers and has been designed to alter the existing…
Hybrid particle-field methods are computationally efficient approaches for modelling soft matter systems. So far applications of these methodologies have been limited to constant volume conditions. Here, we reformulate particle-field…
A new code for astrophysical magnetohydrodynamics (MHD) is described. The code has been designed to be easily extensible for use with static and adaptive mesh refinement. It combines higher-order Godunov methods with the constrained…
The streaming instability is a leading mechanism for concentrating solids and initiating planetesimal formation in protoplanetary disks. Although numerous studies have explored its linear growth, nonlinear evolution, and implications for…
We present the formulation, algorithm and numerical tests of the magnetohydrodynamic-particle-in-cell (MHD-PIC) method with particles treated under the guiding center approximation, which we term the MHD-gPIC method, and it is implemented…
The quasi-neutral hybrid model with kinetic ions and fluid electrons is a promising approach for bridging the inherent multi-scale nature of many problems in space and laboratory plasmas. Here, a novel, implicit, particle-in-cell based…
We investigate streaming instability using two-fluid approximation (neutral gas and dust) in a quasi-global, unstratified protoplanetary disc, with the help of PIERNIK code. We compare amplification rate of the eigen-mode in numerical…