Related papers: A New Time-Dependent Finite Difference Method for …
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…
In the paper, we propose a new effective mathematical formulation and resulting universal numerical algorithm capable of tackling various HF models in the framework of a unified approach. The presented numerical scheme is not limited to any…
In this paper, we study the problem of computing the effective diffusivity for particles moving in chaotic flows. Instead of solving a convection-diffusion type cell problem in the Eulerian formulation (arising from homogenization theory…
We develop a relativistic lattice Boltzmann model capable of describing relativistic fluid dynamics at ultra-high velocities, with Lorentz factors up to $\gamma \sim 10$. To this purpose, we first build a new lattice kinetic scheme by…
A universal particle velocity based algorithm for simulating hydraulic fracture with leak-off, previously demonstrated for the PKN and KGD models, is extended to obtain solutions for a penny-shaped crack. The numerical scheme is capable of…
To simulate the dynamics of fluid with polydisperse particles on macroscale level, one has to solve hydrodynamic equations with several relaxation terms, representing momentum transfer from fluid to particles and vice versa. For small…
Relativistic shocks accelerate particles by the first order Fermi mechanism. These particles then emit synchrotron emission in the post shock gas. We have developed a numerical code which integrates the relativistic Euler equations for…
Accurate mesh-free simulation of fluid flows involving complex boundaries requires that the boundaries be captured accurately in terms of particles. In the context of incompressible/weakly-compressible fluid flow, the SPH method is more…
We present results from a Monte Carlo simulation of a parallel collisionless shock undergoing particle acceleration. Our simulation, which contains parameterized scattering and a particular thermal leakage injection model, calculates the…
We review recent progress on collisionless relativistic shocks. Kinetic instability theory is briefed including its predictions and limitations. The main focus is on numerical experiments in (i) pair and (ii) electron-nucleon plasmas. The…
A general, fast, and effective approach is developed for numerical calculation of kinetic plasma dispersion relations. The plasma dispersion function is approximated by $J$-pole expansion. Subsequently, the dispersion relation is…
We consider the simulation of isentropic flow in pipelines and pipe networks. Standard operating conditions in pipe networks suggest an emphasis to simulate low Mach and high friction regimes -- however, the system is stiff in these regimes…
An "exact" method for scalar one-dimensional hyperbolic conservation laws is presented. The approach is based on the evolution of shock particles, separated by local similarity solutions. The numerical solution is defined everywhere, and is…
Generalization of the self-similar solution for ultrarelativistic shock waves (Blandford & McKee, 1976) is obtained in presence of losses localized on the shock front or distributed in the downstream medium. It is shown that there are two…
This paper introduces a framework for simulating finite dimensional representations of (jump) diffusion sample paths over finite intervals, without discretisation error (exactly), in such a way that the sample path can be restored at any…
Based on the particle-in-cell (PIC) plasma simulation method, the speed-limited PIC (SLPIC) method delivers faster kinetic plasma simulation in cases where the particle distributions evolve slowly compared with the maximum stable PIC…
We have developed a new, very efficient numerical scheme to solve the CR diffusion convection equation that can be applied to the study of the nonlinear time evolution of CR modified shocks for arbitrary spatial diffusion properties. The…
Immersed boundary-lattice Boltzmann method (IB-LBM) has been widely used for simulation of particle-laden flows recently. However, it was limited to small-scale simulations with no more than O(103) particles. Here, we expand IB-LBM for…
The paper presents a solution for axisymmetric propagation of a penny-shaped crack driven by a thinning fluid. The solution to the accuracy of four significant digits, at least, is obtained on the basis of the modified formulation of…
Physical processes ranging from the Lamb shift to the energy loss dE/dx of a charged particle traversing a plasma entail processes that occur over a wide range of energy or length scales. Different physical mechanisms dominate at one or the…