0000000000799144
AUTHOR
J. M. Ibanez
A Divergence-Free High-Resolution Code for MHD
We describe a 2.5D numerical code to solve the equations of ideal magnetohydrodynamics (MHD). The numerical code, based on high-resolution shock-capturing (HRSC) techniques, solves the equations written in conservation form and computes the numerical fluxes using a linearized Riemann solver. A special procedure is used to force the conservation of magnetic flux along the time.
Riemann Solvers in General Relativistic Hydrodynamics
Our contribution concerns with the numerical solution of the 3D general relativistic hydrodynamical system of equations within the framework of the 3+1 formalism. We summarize the theoretical ingredients which are necessary in order to build up a numerical scheme based on the solution of local Riemann problems. Hence, the full spectral decomposition of the Jacobian matrices of the system, i.e., the eigenvalues and the right and left eigenvectors, is explicitly shown. An alternative approach consists in using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows. Our proposal relies on a local change of coordinates in te…
Hydrodynamical and Emission Simulations of Relativistic Jets: Stability and Generation of Superluminal and Stationary Components
We present 2D hydrodynamical and emission simulations of the jet stabilityafter the introduction of strong perturbations on a relativistic jet. These simulations show that the interaction of a single strong perturbation with the underlying jet results in the formation of multiple conical shocks with very specific observational properties.
A Numerical Study of Relativistic Jets
Numerical simulations of supersonic jets are able to explain the structures observed in many VLA images of radio sources. The improvements achieved in classical simulations (see Hardee, these proceedings) are in contrast with the almost complete lack of relativistic simulations the reason being that numerical difficulties arise from the highly relativistic flows typical of extragalactic jets. For our study, we have developed a two-dimensional code which is based on (i) an explicit conservative differencing of the special relativistic hydrodynamics (SRH) equations and (ii) the use of an approximate Riemann solver (see Marti et al. 1995a,b and references therein).
Relativistic MHD simulations of extragalactic jets
We have performed a comprehensive parameter study of the morphology and dynamics of axisymmetric, magnetized, relativistic jets by means of numerical simulations. The simulations have been performed with an upgraded version of the GENESIS code which is based on a second-order accurate finite volume method involving an approximate Riemann solver suitable for relativistic ideal magnetohydrodynamic flows, and a method of lines. Starting from pure hydrodynamic models we consider the effect of a magnetic field of increasing strength (up to β ≡ |b|2/2p ≈ 3.3 times the equipartition value) and different topology (purely toroidal or poloidal). We computed several series of models investigating the …
High--Resolution 3D Simulations of Relativistic Jets
We have performed high-resolution 3D simulations of relativistic jets with beam flow Lorentz factors up to 7, a spatial resolution of 8 cells per beam radius, and for up to 75 normalized time units to study the morphology and dynamics of 3D relativistic jets. Our simulations show that the coherent fast backflows found in axisymmetric models are not present in 3D models. We further find that when the jet is exposed to non-axisymmetric perturbations, (i) it does not display the strong perturbations found for 3D classical hydrodynamic and MHD jets (at least during the period of time covered by our simulations), and (ii) it does propagate according to the 1D estimate. Small 3D effects in the re…
Relativistic simulations of superluminal sources
Abstract We present numerical simulations of the radio emission from hydrodynamical relativistic jets. The quiescent-state jet emission consists of quasi-periodic knots of high emission, associated with internal recollimation shocks. Superluminal components can be reproduced by introducing a square-wave perturbation in the injection velocity of the jet. Strong interactions of the resulting moving shock and the standing recollimations result in a “drag” and increase in emission of the latter.
Upwind Relativistic MHD Code for Astrophysical Applications
We describe the status of devolpment of a 2.5D numerical code to solve the equations of ideal relativistic magnetohydrodynamics. The numerical code, based on high-resolution shock-capturing techniques, solves the equations written in conservation form and computes the numerical fluxes using a linearized Riemann solver. A special procedure is used to force the conservation of magnetic flux along the evolution.
Analysis of the Characteristics in the Meudon Constrained Evolution Scheme
A first analysis of the characteristics associated with the evolving modes in the constraint evolution scheme proposed by the Meudon group in 2004 is presented. The system is written in a first-order hyperbolic form and a so-called generalized Dirac gauge is considered. Applications to inner boundary conditions in an excised approach to black hole evolutions are discussed.
Simulations of Precessing Jets
We report on the results of a three-dimensional, relativistic, hydrodynamical simulation of a precessing jet through which a compact blob of matter is set to propagate. We conclude that the morphology of superluminal sources is the result of a complex combination of phase motions, viewing angle selection effects, and non-linear interactions between perturbations and the underlying jet and/or external medium.
CFC+: Improved dynamics and gravitational waveforms from relativistic core collapse simulations
Core collapse supernovae are a promising source of detectable gravitational waves. Most of the existing (multidimensional) numerical simulations of core collapse in general relativity have been done using approximations of the Einstein field equations. As recently shown by Dimmelmeier et al (2002a,b), one of the most interesting such approximation is the so-called conformal flatness condition (CFC) of Isenberg, Wilson and Mathews. Building on this previous work we present here new results from numerical simulations of relativistic rotational core collapse in axisymmetry, aiming at improving the dynamics and the gravitational waveforms. The computer code used for these simulations evolves th…
On the convexity of relativistic ideal magnetohydrodynamics
We analyze the influence of the magnetic field in the convexity properties of the relativistic magnetohydrodynamics system of equations. To this purpose we use the approach of Lax, based on the analysis of the linearly degenerate/genuinely non-linear nature of the characteristic fields. Degenerate and non-degenerate states are discussed separately and the non-relativistic, unmagnetized limits are properly recovered. The characteristic fields corresponding to the material and Alfv\'en waves are linearly degenerate and, then, not affected by the convexity issue. The analysis of the characteristic fields associated with the magnetosonic waves reveals, however, a dependence of the convexity con…