Refining a relativistic, hydrodynamic solver: Admitting ultra-relativistic flows
arXiv:astro-ph/0606012 · doi:10.1016/j.jcp.2009.05.012
Abstract
We have undertaken the simulation of hydrodynamic flows with bulk Lorentz factors in the range 10^2--10^6. We discuss the application of an existing relativistic, hydrodynamic primitive-variable recovery algorithm to a study of pulsar winds, and, in particular, the refinement made to admit such ultra-relativistic flows. We show that an iterative quartic root finder breaks down for Lorentz factors above 10^2 and employ an analytic root finder as a solution. We find that the former, which is known to be robust for Lorentz factors up to at least 50, offers a 24% speed advantage. We demonstrate the existence of a simple diagnostic allowing for a hybrid primitives recovery algorithm that includes an automatic, real-time toggle between the iterative and analytical methods. We further determine the accuracy of the iterative and hybrid algorithms for a comprehensive selection of input parameters and demonstrate the latter's capability to elucidate the internal structure of ultra-relativistic plasmas. In particular, we discuss simulations showing that the interaction of a light, ultra-relativistic pulsar wind with a slow, dense ambient medium can give rise to asymmetry reminiscent of the Guitar nebula leading to the formation of a relativistic backflow harboring a series of internal shockwaves. The shockwaves provide thermalized energy that is available for the continued inflation of the PWN bubble. In turn, the bubble enhances the asymmetry, thereby providing positive feedback to the backflow.
33 pages, 11 figures
References in corpus (25)
- How Massive Single Stars End their Life
- The Evolution and Structure of Pulsar Wind Nebulae
- HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
- The Velocity Distribution of Isolated Radio Pulsars
- Relativistic Jets in Collapsars
- Primitive Variable Solvers for Conservative General Relativistic Magnetohydrodynamics
- Neutron Star Population Dynamics.II: 3D Space Velocities of Young Pulsars
- An efficient shock-capturing central-type scheme for multidimensional relativistic flows. II. Magnetohydrodynamics
- Pulsar Wind Nebulae in the Chandra Era
- Equation of State in Relativistic Magnetohydrodynamics: variable versus constant adiabatic index
- Axially symmetric relativistic MHD simulations of Pulsar Wind Nebulae in Supernova Remnants - On the origin of torus and jet-like features
- An efficient shock-capturing central-type scheme for multidimensional relativistic flows. I. Hydrodynamics
- Relativistic MHD simulations of pulsar bow-shock nebulae
- An evolutionary model for pulsar-driven supernova remnants
- Axisymmetric general relativistic hydrodynamics: Long-term evolution of neutron stars and stellar collapse to neutron stars and black holes
- Equation of State in Numerical Relativistic Hydrodynamics
- Bow Shocks from Neutron Stars: Scaling Laws and HST Observations of the Guitar Nebula
- Interaction of High-Velocity Pulsars with Supernova Remnant Shells
- 3D Hydrodynamic Simulations of Relativistic Extragalactic Jets
- Three-dimensional hydrodynamic simulations of asymmetric pulsar wind bow shocks
- The Origin of Complex Behavior of Linearly Polarized Components in Parsec-Scale Jets
- X-ray Emission from the Guitar Nebula
- Hydrodynamic models
- A Numerical Study of Relativistic Fluid Collapse
- Modeling Pulsar Wind Nebulae