Related papers: An implicit scheme for solving the anisotropic dif…
Anisotropic diffusion is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and the interplay of cosmic rays with the Galactic magnetic field. This diffusion term contributes to the highly stiff nature of…
In the interstellar medium of galaxies and the intracluster gas of galaxy clusters, the charged particles making up cosmic rays are moving almost exclusively along (but not across) magnetic field lines. The resulting anisotropic transport…
We present a new multi-fluid, multi-temperature plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (non-linear, implicit) scheme for collisional low-temperature plasma. The particle transport is described using the…
The anisotropic diffusion equation is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and its interplay with the ambient magnetic field. This diffusion term contributes to the highly stiff nature of the…
Simple finite differencing of the anisotropic diffusion equation, where diffusion is only along a given direction, does not ensure that the numerically calculated heat fluxes are in the correct direction. This can lead to negative…
The macroscopic behaviour of cosmic rays in turbulent magnetic fields is discussed. An implementation of anisotropic diffusion of cosmic rays with respect to the magnetic field in a non-conservative, high-order, finite-difference…
We develop a numerical scheme for solving a fully special relativistic resistive radiation magnetohydrodynamics. Our code guarantees conservations of total mass, momentum and energy. Radiation energy density and radiation flux are…
In weakly collisional, strongly magnetised plasmas such as the intracluster medium (ICM), hot accretion flows and the solar corona, the transport of heat and momentum occurs primarily along magnetic field lines. In this paper we present a…
In this paper, we present an interior penalty discontinuous Galerkin finite element scheme for solving diffusion problems with strong anisotropy arising in magnetized plasmas for fusion applications. We demonstrate the accuracy produced by…
The acceleration of energetic particle transport in high amplitude magnetosonic and Alfvenic turbulence is considered using the method of Monte Carlo particle simulations, involving integration of particle equations of motion. We derive the…
We present a robust and accurate numerical method for the anisotropic diffusion equation in curvilinear coordinates. This study extends the recent work [Muir et al., Computer Physics Communications, 2025] for solving the anisotropic…
We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial…
We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following [Du Toit et al., Comp. Phys.…
In many astrophysical plasmas, the Coulomb collision is insufficient to maintain an isotropic temperature, and the system is driven to the anisotropic regime. In this case, magnetohydrodynamic (MHD) models with anisotropic pressure are…
The aim of this work is to introduce a numerical method to cope with the multiscale nature of confined plasma physics. These investigations are focused on fluid plasma description under large magnetic field. The difficulties in this context…
Cosmic rays (CRs) are an important source of feedback in a variety of astrophysical contexts. Magneto-hydrodynamical (MHD) simulations treating CRs as a fluid have shown that how their feedback operates is strongly dependent on their…
Energetic particle transport in a finite amplitude magnetosonic and Alfvenic turbulence is considered using Monte Carlo particle simulations, which involve an integration of particle equation of motion. We show that in a low-Betha plasma…
Implicit particle-in-cell codes offer advantages over their explicit counterparts in that they suffer weaker stability constraints on the need to resolve the higher frequency modes of the system. This feature may prove particularly valuable…
We propose an explicit-implicit scheme for numerically solving Special Relativistic Radiation Hydrodynamic (RRHD) equations, which ensures a conservation of total energy and momentum (matter and radiation). In our scheme, 0th and 1st moment…
In magnetized plasma, the magnetic field confines particles around field lines. The ratio between the intensity of the parallel and perpendicular viscosity or heat conduction may reach the order of $10^{12}$. When the magnetic fields have…