Related papers: An Efficient Numerical Scheme for Simulating Parti…
A numerical scheme that incorporates a thermal leakage injection model into a combined gas dynamics and cosmic ray (CR, hereafter) diffusion-convection code has been developed. The particle injection is followed numerically by filtering the…
Numerical solutions of the cosmic-ray (CR) magneto-hydrodynamic equations are dogged by a powerful numerical instability, which arises from the constraint that CRs can only stream down their gradient. The standard cure is to regularize by…
We have calculated the evolution of cosmic ray (CR) modified astrophysical shocks for a wide range of shock Mach numbers and shock speeds through numerical simulations of diffusive shock acceleration (DSA) in 1D quasi- parallel plane…
We have developed a novel computer code designed to follow the evolution of cosmic-ray modified shocks, including the full momentum dependence of the particles for a realistic diffusion coefficient model. In this form the problem is…
We consider cosmic ray (CR) modified shocks with both streaming and diffusion in the two-fluid description. Previously, numerical codes were unable to incorporate streaming in this demanding regime, and have never been compared against…
We have developed a cosmic ray (CR) shock code in one dimensional spherical geometry with which the particle distribution, the gas flow and their nonlinear interaction can be followed numerically in a frame comoving with an expanding shock.…
One of the key questions to understanding the efficiency of diffusive shock acceleration of the cosmic rays (CRs) is the injection process from thermal particles. A self-consistent injection model based on the interactions of the…
We present a more accurate numerical scheme for the calculation of diffusive shock acceleration of cosmic rays using Stochastic Differential Equations. The accuracy of this scheme is demonstrated using a simple analytical flow profile that…
Cosmological hydrodynamic simulations of large scale structure in the universe have shown that accretion shocks and merger shocks form due to flow motions associated with the gravitational collapse of nonlinear structures. Estimated speed…
A popular numerical method to model the dynamics of a 'full spectrum' of cosmic rays (CRs), also applicable to radiation/neutrino hydrodynamics (RHD), is to discretize the spectrum at each location/cell as a piecewise power law in 'bins' of…
We discuss new methods to integrate the cosmic ray (CR) evolution equations coupled to magneto-hydrodynamics (MHD) on an unstructured moving mesh, as realised in the massively parallel AREPO code for cosmological simulations. We account for…
We present new results characterizing cosmological shocks within adaptive mesh refinement N-Body/hydrodynamic simulations that are used to predict non-thermal components of large-scale structure. This represents the first study of shocks…
We present a new code aimed at the simulation of diffusive shock acceleration (DSA), and discuss various test cases which demonstrate its ability to study DSA in its full time-dependent and non-linear developments. We present the numerical…
Kinetic equations model the position-velocity distribution of particles subject to transport and collision effects. Under a diffusive scaling, these combined effects converge to a diffusion equation for the position density in the limit of…
Using an improved version of the previously introduced CRASH (Cosmic Ray Acceleration SHock) code, we have calculated the time evolution of cosmic-ray (CR) modified quasi-parallel plane shocks for Bohm-like diffusion, including…
In this paper, we develop and analyze a stochastic algorithm for solving space-time fractional diffusion models, which are widely used to describe anomalous diffusion dynamics. These models pose substantial numerical challenges due to the…
This work proposes a novel shape optimization framework for geometric inverse problems governed by the advection--diffusion equation, based on the coupled complex boundary method (CCBM). Building on recent developments [Afr22, Rab23, Rab25,…
This article presents a new finite element method for convection-diffusion equations by enhancing the continuous finite element space with a flux space for flux approximations that preserve the important mass conservation locally on each…
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…
First of all, this paper presents some improvements of DSMC method in the form of new schemes and approaches, that, for a wide class of problems, increase performance and reduce the demands on computer resources. The most important…