DIPLODOCUS II: Implementation of transport equations and test cases relevant to micro-scale physics of jetted astrophysical sources

Abstract

DIPLODOCUS (Distribution-In-PLateaux methODOlogy for the CompUtation of transport equationS) is a framework being developed for the mesoscopic modelling of astrophysical systems via the transport of particle distribution functions through the seven dimensions of phase space, including continuous forces and discrete interactions between particles. Following Paper I, which details the mathematical background, this second paper provides an overview of the numerical implementation in the form of the code package Diplodocus$.$jl, written in Julia, including the description of a novel Monte-Carlo sampling technique for the pre-computation of anisotropic collision integrals. In addition to the discussion of numerical implementation, a selection of test cases are presented to examine the package's capabilities. These test cases focus on micro-scale physical effects: binary collisions, emissive interactions, and external forces that are relevant to the modelling of jetted astrophysical sources, such as Active Galactic Nuclei and X-Ray Binaries.