Related papers: A New Time-Dependent Finite Difference Method for …
The outflows from gamma ray bursts, active galactic nuclei and relativistic jets in general interact with the surrounding media through collisionless shocks. With three dimensional relativistic particle-in-cell simulations we investigate…
In recent years, Rectified flow (RF) has gained considerable popularity largely due to its generation efficiency and state-of-the-art performance. In this paper, we investigate the degree to which RF automatically adapts to the intrinsic…
A one-dimensional particle-in-cell (PIC) simulation tracks a fast magnetosonic shock over time scales comparable to an inverse ion gyrofrequency. The magnetic pressure is comparable to the thermal pressure upstream. The shock propagates…
Turbulent relative dispersion is studied theoretically with a focus on the evolution of probability distribution of the relative separation of two passive particles. A finite separation speed and a finite correlation of relative velocity,…
Drag laws for particles in fluids are often expressed in terms of the undisturbed fluid velocity, defined as the fluid velocity a particle sees before the disturbance develops in the fluid. In two-way coupled point-particle simulations the…
We describe a new hybrid framework to model non-thermal spectral signatures from highly energetic particles embedded in a large-scale classical or relativistic MHD flow. Our method makes use of \textit{Lagrangian} particles moving through…
In this paper, we propose stochastic structure-preserving schemes to compute the effective diffusivity for particles moving in random flows. We first introduce the motion of particles using the Lagrangian formulation, which is modeled by…
We consider the acceleration of charged particles near ultra-relativistic shocks, with Lorentz factor Gamma_s >> 1. We present simulations of the acceleration process and compare these with results from semi-analytical calculations. We show…
We present a new numerical dissipation algorithm, which can be efficiently used in combination with centered finite-difference methods. We start from a formulation of centered finite-volume methods for Numerical Relativity, in which…
We present results from a fully relativistic Monte Carlo simulation of diffusive shock acceleration (DSA) in unmodified shocks. The computer code uses a single algorithmic sequence to smoothly span the range from nonrelativistic speeds to…
This paper proposes a novel and simple linear model to capture line losses for use in linearized DC models, such as optimal power flow (DC-OPF) and security-constrained economic dispatch (SCED). The \textit{Line Loss Outer Approximation}…
The problem of anomalous diffusion in momentum space is considered for plasma-like systems on the basis of a new collision integral, which is appropriate for consideration of the probability transition function (PTF) with long tails in…
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,…
Ultrahigh peak power femtosecond laser pulses create extreme states of matter that are currently being probed with great interest. Plasma optics have been proposed for shaping and amplifying high-power pulses, but they are subject to huge…
A very intuitive description of nucleus-nucleus collision phenomena is provided by the relativistic fluid dynamics. We consider a 1+1 dimensional relativistic imperfect fluid flow to approximate the high energy heavy ion collision. The…
Symplectic schemes are powerful methods for numerically integrating Hamiltonian systems, and their long-term accuracy and fidelity have been proved both theoretically and numerically. However direct applications of standard symplectic…
We present a computational algorithm for computing short range forces between particles. The algorithm has two distinguishing features. First, it is optimized for multi-processor computers, and will use as many processors as are available.…
Collisionless shocks are frequently analyzed using the magnetohydrodynamics (MHD) formalism, even though MHD assumes a small mean free path. Yet, isotropy of pressure, fruit of binary collisions and assumed in MHD, may not apply in…
In this article, on the basis of the Langevin equation applied to velocity fluctuations, we numerically model the Partial Variance of Increments, which is a useful tool to measure time and spatial correlations in space plasmas. We consider…
We have examined cutoffs and pile-ups due to various processes in the spectra of particles produced by shock acceleration, and found that, even in the absence of energy losses, the shape of the spectrum of accelerated particles at energies…