Related papers: A numerical approach to the non-uniqueness problem…
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…
Cosmic rays (CRs) play a pivotal role in various astrophysical systems, delivering feedback over a broad range of scales. However, modeling CR transport remains challenging due to its inherently multi-scale nature and complex microphysics.…
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…
When a magnetised, thermal plasma containing a shock wave efficiently accelerates Cosmic Rays (CR) a reaction is exerted on the fluid structure. The simplest description of this process involves three conservation equations for the…
The acceleration of relativistic particles due to repeated scattering across a shock wave remains the most attractive model for the production of energetic cosmic rays. This process has been analyzed extensively during the past two decades…
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.…
The stability of the cosmic ray modified shock (CRMS) is studied by means of numerical simulations. Owing to the nonlinear feedback of cosmic-ray (CR) acceleration, a downstream state of the modified shock can no longer be uniquely…
The inability of the Navier-Stokes-Fourier equations to capture rarefaction effects motivates us to adopt the extended hydrodynamic equations. In the present work, a hydrodynamic model, which consists of the conservation laws closed with…
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…
It is well known that cosmic rays (CRs) contribute significantly to the pressure of the interstellar medium in our own Galaxy, suggesting that they may play an important role in regulating star formation during the formation and evolution…
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…
Many recent numerical studies have argued that cosmic rays (CRs) from supernovae (SNe) or active galactic nuclei (AGN) could play a crucial role in galaxy formation, in particular by establishing a CR-pressure dominated circum-galactic…
A numerical method for integrating the equations describing a dynamically coupled system made of a fluid and cosmic-rays is developed. In smooth flows the effect of CR pressure is accounted for by modification of the characteristic…
We study numerically the standard one pressure model of two fluid flows with energy equations. This system is not solved in time derivative. It has been transformed into an equivalent system solved in time derivative. We show that the…
The speed of sound in two-phase pipe flow systems is often several orders of magnitude greater than the travelling speed of hydraulic information (volume fractions.) Dynamically simulating such flows requires resolution of acoustic and…
Convection schemes are a large source of error in global weather and climate models, and modern resolutions are often too fine to parameterise convection but are still too coarse to fully resolve it. Recently, numerical solutions of…
We analyze the cosmic-ray magnetohydrodynamic (CR MHD) equations to improve the numerical simulations. We propose to solve them in the fully conservation form, which is equivalent to the conventional CR MHD equations. In the fully…
The problem of two stiff fluids (energy density = pressure) moving radially in spherical symmetry is treated. The metric ansatz is chosen spherically symmetric, conformally static with a multiplicative separation of variables. The first…
We propose a new numerical approach to compute nonclassical solutions to hyperbolic conservation laws. The class of finite difference schemes presented here is fully conservative and keep nonclassical shock waves as sharp interfaces,…
We present COSMOCR, a numerical code for the investigation of cosmic ray related studies in computational cosmology. The code follows the diffusive shock acceleration, the mechanical and radiative energy losses and the spatial transport of…