BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Olindo Zanotti

Olindo Zanotti contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEgr-qcmath.NAastro-ph.IMphysics.comp-phastro-ph.SRNumerical Analysisphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - Emerging
27works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2015arXiv

A Posteriori Subcell Limiting of the Discontinuous Galerkin Finite Element Method for Hyperbolic Conservation Laws

The purpose of this work is to propose a novel a posteriori finite volume subcell limiter technique for the Discontinuous Galerkin finite element method for nonlinear systems of hyperbolic conservation laws in multiple space dimensions that works well for arbitrary high order of accuracy in space and time and that does not destroy the natural subcell resolution properties of the DG method. High order time discretization is achieved via a one-step ADER approach that uses a local space-time discontinuous Galerkin predictor method. Our new limiting strategy is based on the so-called MOOD paradigm, which aposteriori verifies the validity of a discrete candidate solution against physical and numerical detection criteria. Within the DG scheme on the main grid, the discrete solution is represented by piecewise polynomials of degree N. For those troubled cells that need limiting, our new limiter approach recomputes the discrete solution by scattering the DG polynomials at the previous time step onto a set of N_s=2N+1 finite volume subcells per space dimension. A robust but accurate ADER-WENO finite volume scheme then updates the subcell averages of the conservative variables within the detected troubled cells. The choice of N_s=2N+1 subcells is optimal since it allows to match the maximum admissible time step of the finite volume scheme on the subgrid with the maximum admissible time step of the DG scheme on the main grid. We illustrate the performance of the new scheme via the simulation of numerous test cases in two and three space dimensions, using DG schemes of up to tenth order of accuracy in space and time (N=9). The method is also able to run on massively parallel large scale supercomputing infrastructure, which is shown via one 3D test problem that uses 10 billion space-time degrees of freedom per time step.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement

We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method. Due to the one-step nature of the underlying scheme, the resulting algorithm is particularly well suited for an AMR strategy on space-time adaptive meshes, i.e.with time-accurate local time stepping. The AMR property has been implemented 'cell-by-cell', with a standard tree-type algorithm, while the scheme has been parallelized via the Message Passing Interface (MPI) paradigm. The new scheme has been tested over a wide range of examples for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations of compressible gas dynamics and the equations of magnetohydrodynamics (MHD). High order in space and time have been confirmed via a numerical convergence study and a detailed analysis of the computational speed-up with respect to highly refined uniform meshes is also presented. We also show test problems where the presented high order AMR scheme behaves clearly better than traditional second order AMR methods. The proposed scheme that combines for the first time high order ADER methods with space--time adaptive grids in two and three space dimensions is likely to become a useful tool in several fields of computational physics, applied mathematics and mechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

General relativistic radiation hydrodynamics of accretion flows. I: Bondi-Hoyle accretion

We present a new code for performing general-relativistic radiation-hydrodynamics simulations of accretion flows onto black holes. The radiation field is treated in the optically-thick approximation, with the opacity contributed by Thomson scattering and thermal bremsstrahlung. Our analysis is concentrated on a detailed numerical investigation of hot two-dimensional, Bondi-Hoyle accretion flows with various Mach numbers. We find significant differences with respect to purely hydrodynamical evolutions. In particular, once the system relaxes to a radiation-pressure dominated regime, the accretion rates become about two orders of magnitude smaller than in the purely hydrodynamical case, remaining however super-Eddington as are the luminosities. Furthermore, when increasing the Mach number of the inflowing gas, the accretion rates become smaller because of the smaller cross section of the black hole, but the luminosities increase as a result a stronger emission in the shocked regions. Overall, our approach provides the first self-consistent calculation of the Bondi-Hoyle luminosity, most of which is emitted within r~100 M from the black hole, with typical values L/L_Edd ~ 1-7, and corresponding energy efficiencies eta_BH ~ 0.09-0.5. The possibility of computing luminosities self-consistently has also allowed us to compare with the bremsstrahlung luminosity often used in modelling the electromagnetic counterparts to supermassive black-hole binaries, to find that in the optically-thick regime these more crude estimates are about 20 times larger than our radiation-hydrodynamics results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: viscous heat-conducting fluids and elastic solids

This paper is concerned with the numerical solution of the unified first order hyperbolic formulation of continuum mechanics recently proposed by Peshkov & Romenski, denoted as HPR model. In that framework, the viscous stresses are computed from the so-called distortion tensor A, which is one of the primary state variables. A very important key feature of the model is its ability to describe at the same time the behavior of inviscid and viscous compressible Newtonian and non-Newtonian fluids with heat conduction, as well as the behavior of elastic and visco-plastic solids. This is achieved via a stiff source term that accounts for strain relaxation in the evolution equations of A. Also heat conduction is included via a first order hyperbolic evolution equation of the thermal impulse, from which the heat flux is computed. The governing PDE system is hyperbolic and fully consistent with the principles of thermodynamics. It is also fundamentally different from first order Maxwell-Cattaneo-type relaxation models based on extended irreversible thermodynamics. The connection between the HPR model and the classical hyperbolic-parabolic Navier-Stokes-Fourier theory is established via a formal asymptotic analysis in the stiff relaxation limit. From a numerical point of view, the governing partial differential equations are very challenging, since they form a large nonlinear hyperbolic PDE system that includes stiff source terms and non-conservative products. We apply the successful family of one-step ADER-WENO finite volume and ADER discontinuous Galerkin finite element schemes in the stiff relaxation limit, and compare the numerical results with exact or numerical reference solutions obtained for the Euler and Navier-Stokes equations. To show the universality of the model, the paper is rounded-off with an application to wave propagation in elastic solids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Model for an optically thick torus in local thermodynamic equilibrium around a black hole

We propose a simple model for an optically thick radiative torus in local thermodynamic equilibrium around a Kerr black hole. The hydrodynamics structure, which is not affected by the radiation field, is the same as for the so--called polish doughnuts. Under the assumption of isentropic fluid and polytropic equation of state, a simple stationary and axisymmetric solution to the relativistic radiation hydrodynamics equations is possible, for which the temperature of the torus scales like the specific enthalpy. The astrophysical relevance of the model is briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Oscillations of radiation pressure supported tori near black holes

We study the dynamics of radiation pressure supported tori around Schwarzschild black holes, focusing on their oscillatory response to an external perturbation. Using KORAL, a general relativistic radiation hydrodynamics code capable of modeling all radiative regimes from the optically thick to the optically thin, we monitor a sample of models at different initial temperatures and opacities, evolving them in two spatial dimensions for $\sim 165$ orbital periods. The dynamics of models with high opacity is very similar to that of purely hydrodynamics models, and it is characterized by regular oscillations which are visible also in the light curves. As the opacity is decreased, the tori quickly and violently migrate towards the gas-pressure dominated regime, collapsing towards the equatorial plane. When the spectra of the $L_2$ norm of the mass density are considered, high frequency inertial-acoustic modes of oscillations are detected (with the fundamental mode at a frequency $68 M_{\rm BH}^{-1}\,\rm Hz$), in close analogy to the phenomenology of purely hydrodynamic models. An additional mode of oscillation, at a frequency $129 M_{\rm BH}^{-1}\,\rm Hz$, is also found, which can be unambiguously attributed to the radiation. The spectra extracted from the light curves are typically more noisy, indicating that in a real observation such modes would not be easily detected.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Particle acceleration in the polar cap region of an oscillating neutron star

We revisit particle acceleration in the polar cap region of a neutron star by taking into account both general relativistic effects and the presence of toroidal oscillations at the star surface. In particular, we address the question of whether toroidal oscillations at the stellar surface can affect the acceleration properties in the polar cap. We solve numerically the relativistic electrodynamics equations in the stationary regime, focusing on the computation of the Lorentz factor of a space-charge-limited electron flow accelerated in the polar cap region of a rotating and oscillating pulsar. To this extent, we adopt the correct expression of the general relativistic Goldreich-Julian charge density in the presence of toroidal oscillations. Depending on the ratio of the actual charge density of the pulsar magnetosphere to the Goldreich-Julian charge density, we distinguish two different regimes of the Lorentz factor of the particle flow, namely an oscillatory regime produced for sub-GJ current density configurations, which does not produce an efficient acceleration, and a true accelerating regime for super-GJ current density configurations. We find that star oscillations may be responsible for a significant asymmetry in the pulse profile that depends on the orientation of the oscillations with respect to the pulsar magnetic field. In particular, significant enhancements of the Lorentz factor are produced by stellar oscillations in the super-GJ current density regime.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solving the relativistic magnetohydrodynamics equations with ADER discontinuous Galerkin methods, a posteriori subcell limiting and adaptive mesh refinement

We present a new numerical tool for solving the special relativistic ideal MHD equations that is based on the combination of the following three key features: (i) a one-step ADER discontinuous Galerkin (DG) scheme that allows for an arbitrary order of accuracy in both space and time, (ii) an a posteriori subcell finite volume limiter that is activated to avoid spurious oscillations at discontinuities without destroying the natural subcell resolution capabilities of the DG finite element framework and finally (iii) a space-time adaptive mesh refinement (AMR) framework with time-accurate local time-stepping. The divergence-free character of the magnetic field is instead taken into account through the so-called "divergence-cleaning" approach. The convergence of the new scheme is verified up to 5th order in space and time and the results for a set of significant numerical tests including shock tube problems, the RMHD rotor and blast wave problems, as well as the Orszag-Tang vortex system are shown. We also consider a simple case of the relativistic Kelvin-Helmholtz instability with a magnetic field, emphasizing the potential of the new method for studying turbulent RMHD flows. We discuss the advantages of our new approach when the equations of relativistic MHD need to be solved with high accuracy within various astrophysical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Space-time adaptive ADER discontinuous Galerkin finite element schemes with a posteriori sub-cell finite volume limiting

In this paper we present a novel arbitrary high order accurate discontinuous Galerkin (DG) finite element method on space-time adaptive Cartesian meshes (AMR) for hyperbolic conservation laws in multiple space dimensions, using a high order \aposteriori sub-cell ADER-WENO finite volume \emph{limiter}. Notoriously, the original DG method produces strong oscillations in the presence of discontinuous solutions and several types of limiters have been introduced over the years to cope with this problem. Following the innovative idea recently proposed in \cite{Dumbser2014}, the discrete solution within the troubled cells is \textit{recomputed} by scattering the DG polynomial at the previous time step onto a suitable number of sub-cells along each direction. Relying on the robustness of classical finite volume WENO schemes, the sub-cell averages are recomputed and then gathered back into the DG polynomials over the main grid. In this paper this approach is implemented for the first time within a space-time adaptive AMR framework in two and three space dimensions, after assuring the proper averaging and projection between sub-cells that belong to different levels of refinement. The combination of the sub-cell resolution with the advantages of AMR allows for an unprecedented ability in resolving even the finest details in the dynamics of the fluid. The spectacular resolution properties of the new scheme have been shown through a wide number of test cases performed in two and in three space dimensions, both for the Euler equations of compressible gas dynamics and for the magnetohydrodynamics (MHD) equations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A high order special relativistic hydrodynamic and magnetohydrodynamic code with space-time adaptive mesh refinement

We present a high order one-step ADER-WENO finite volume scheme with space-time adaptive mesh refinement (AMR) for the solution of the special relativistic hydrodynamic and magnetohydrodynamic equations. By adopting a local discontinuous Galerkin predictor method, a high order one-step time discretization is obtained, with no need for Runge--Kutta sub-steps. This turns out to be particularly advantageous in combination with space-time adaptive mesh refinement, which has been implemented following a "cell-by-cell" approach. As in existing second order AMR methods, also the present higher order AMR algorithm features time-accurate local time stepping (LTS), where grids on different spatial refinement levels are allowed to use different time steps. We also compare two different Riemann solvers for the computation of the numerical fluxes at the cell interfaces. The new scheme has been validated over a sample of numerical test problems in one, two and three spatial dimensions, exploring its ability in resolving the propagation of relativistic hydrodynamical and magnetohydrodynamical waves in different physical regimes. The astrophysical relevance of the new code for the study of the Richtmyer--Meshkov instability is briefly discussed in view of future applications.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An ADER-WENO Finite Volume AMR code for Astrophysics

A high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions is presented. A high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method, while a high order spatial accuracy is obtained through a WENO reconstruction. Thanks to the one-step nature of the underlying scheme, the resulting algorithm can be efficiently imported within an AMR framework on space-time adaptive meshes. We provide convincing evidence that the presented high order AMR scheme behaves better than traditional second order AMR methods. Tests are shown of the new scheme for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations and the equations of ideal magnetohydrodynamics. The proposed scheme is likely to become a useful tool in several astrophysical scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Explaining the subpulse drift velocity of pulsar magnetosphere within the space-charge limited flow model

We try to explain the subpulse drift phenomena adopting the space-charge limited flow (SCLF) model and comparing the plasma drift velocity in the inner region of pulsar magnetospheres with the observed velocity of drifting subpulses. We apply the approach described in a recent paper of van Leeuwen & Timokhin (2012), where it was shown that the standard estimation of the subpulse drift velocity through the total value of the scalar potential drop in the inner gap gives inaccurate results, while the exact expression relating the drift velocity to the gradient of the scalar potential should be used instead. After considering a selected sample of sources taken from the catalog of Weltevrede, Edwards & Stappers (2006) with coherently drifting subpulses and reasonably known observing geometry, we show that their subpulse drift velocities would correspond to the drift of the plasma located very close or above the pair formation front. Moreover, a detailed analysis of PSR B0826-34 and PSR B0818-41 reveals that the variation of the subpulse separation with the pulse longitude can be successfully explained by the dependence of the plasma drift velocity on the angular coordinates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High Order Cell-Centered Lagrangian-Type Finite Volume Schemes with Time-Accurate Local Time Stepping on Unstructured Triangular Meshes

We present a novel cell-centered direct Arbitrary-Lagrangian-Eulerian (ALE) finite volume scheme on unstructured triangular meshes that is high order accurate in space and time and that also allows for time-accurate local time stepping (LTS). The new scheme uses the following basic ingredients: a high order WENO reconstruction in space on unstructured meshes, an element-local high-order accurate space-time Galerkin predictor that performs the time evolution of the reconstructed polynomials within each element, the computation of numerical ALE fluxes at the moving element interfaces through approximate Riemann solvers, and a one-step finite volume scheme for the time update which is directly based on the integral form of the conservation equations in space-time. The inclusion of the LTS algorithm requires a number of crucial extensions, such as a proper scheduling criterion for the time update of each element and for each node; a virtual projection of the elements contained in the reconstruction stencils of the element that has to perform the WENO reconstruction; and the proper computation of the fluxes through the space-time boundary surfaces that will inevitably contain hanging nodes in time due to the LTS algorithm. We have validated our new unstructured Lagrangian LTS approach over a wide sample of test cases solving the Euler equations of compressible gasdynamics in two space dimensions, including shock tube problems, cylindrical explosion problems, as well as specific tests typically adopted in Lagrangian calculations, such as the Kidder and the Saltzman problem. When compared to the traditional global time stepping (GTS) method, the newly proposed LTS algorithm allows to reduce the number of element updates in a given simulation by a factor that may depend on the complexity of the dynamics, but which can be as large as 4.7.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High order numerical simulations of the Richtmyer Meshkov instability in a relativistic fluid

We study the Richtmyer--Meshkov (RM) instability of a relativistic perfect fluid by means of high order numerical simulations with adaptive mesh refinement (AMR). The numerical scheme adopts a finite volume Weighted Essentially Non-Oscillatory (WENO) reconstruction to increase accuracy in space, a local space-time discontinuous Galerkin predictor method to obtain high order of accuracy in time and a high order one-step time update scheme together with a &#34;cell-by-cell&#34; space-time AMR strategy with time-accurate local time stepping. In this way, third order accurate (both in space and in time) numerical simulations of the RM instability are performed, spanning a wide parameter space. We present results both for the case in which a light fluid penetrates into a higher density one (Atwood number $A>0$), and for the case in which a heavy fluid penetrates into a lower density one (Atwood number $A<0$). We find that, for large Lorentz factors γ_s of the incident shock wave, the relativistic RM instability is substantially weakened and ultimately suppressed. More specifically, the growth rate of the RM instability in the linear phase has a local maximum which occurs at a critical value of γ_s ~ [1.2,2]. Moreover, we have also revealed a genuine relativistic effect, absent in Newtonian hydrodynamics, which arises in three dimensional configurations with a non-zero velocity component tangent to the incident shock front. In this case, the RM instability is strongly affected, typically resulting in less efficient mixing of the fluid.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Arbitrary-Lagrangian-Eulerian One-Step WENO Schemes for Stiff Hyperbolic Balance Laws

In this article we present a new family of high order accurate Arbitrary Lagrangian-Eulerian one-step WENO finite volume schemes for the solution of stiff hyperbolic balance laws. High order accuracy in space is obtained with a standard WENO reconstruction algorithm and high order in time is obtained using the local space-time discontinuous Galerkin method recently proposed in Dumbser, Enaux, and Toro (2008). In the Lagrangian framework considered here, the local space-time DG predictor is based on a weak formulation of the governing PDE on a moving space-time element. For the space-time basis and test functions we use Lagrange interpolation polynomials defined by tensor-product Gauss-Legendre quadrature points. The moving space-time elements are mapped to a reference element using an isoparametric approach, i.e. the space-time mapping is defined by the same basis functions as the weak solution of the PDE. We show some computational examples in one space-dimension for non-stiff and for stiff balance laws, in particular for the Euler equations of compressible gas dynamics, for the resistive relativistic MHD equations, and for the relativistic radiation hydrodynamics equations. Numerical convergence results are presented for the stiff case up to sixth order of accuracy in space and time and for the non-stiff case up to eighth order of accuracy in space and time.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High Order Space-Time Adaptive WENO Finite Volume Schemes for Non-Conservative Hyperbolic Systems

We present a class of high order finite volume schemes for the solution of non-conservative hyperbolic systems that combines the one-step ADER-WENO finite volume approach with space-time adaptive mesh refinement (AMR). The resulting algorithm, which is particularly well suited for the treatment of material interfaces in compressible multi-phase flows, is based on: (i) high order of accuracy in space obtained through WENO reconstruction, (ii) a high order one-step time discretization via a local space-time discontinuous Galerkin predictor method, and (iii) the use of a path conservative scheme for handling the non-conservative terms of the equations. The AMR property with time accurate local time stepping, which has been treated according to a &#39;cell-by-cell&#39; strategy, strongly relies on the high order one-step time discretization, which naturally allows a high order accurate and consistent computation of the jump terms at interfaces between elements using different time steps. The new scheme has been successfully validated on some test problems for the Baer-Nunziato model of compressible multiphase flows.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accurate Simulations of Binary Black-Hole Mergers in Force-Free Electrodynamics

We provide additional information on our recent study of the electromagnetic emission produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. As anticipated in a recent letter, our results show that although a dual-jet structure is present, the associated luminosity is ~ 100 times smaller than the total one, which is predominantly quadrupolar. We here discuss the details of our implementation of the equations in which the force-free condition is not implemented at a discrete level, but rather obtained via a damping scheme which drives the solution to satisfy the correct condition. We show that this is important for a correct and accurate description of the current sheets that can develop in the course of the simulation. We also study in greater detail the three-dimensional charge distribution produced as a consequence of the inspiral and show that during the inspiral it possesses a complex but ordered structure which traces the motion of the two black holes. Finally, we provide quantitative estimates of the scaling of the electromagnetic emission with frequency, with the diffused part having a dependence that is the same as the gravitational-wave one and that scales as L ~ Omega^{10/3-8/3}, while the collimated one scales as L ~ Omega^{5/3-6/3}, thus with a steeper dependence than previously estimated. We discuss the impact of these results on the potential detectability of dual jets from supermassive black holes and the steps necessary for more accurate estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

General relativistic radiation hydrodynamics of accretion flows: II. Treating stiff source terms and exploring physical limitations

We present the implementation of an implicit-explicit (IMEX) Runge-Kutta numerical scheme for general relativistic hydrodynamics coupled to an optically thick radiation field in two existing GR-hydrodynamics codes. We argue that the necessity of such an improvement arises naturally in astrophysically relevant regimes where the optical thickness is high as the equations become stiff. By performing several 1D tests we verify the codes&#39; new ability to deal with this stiffness and show consistency. Then, still in 1D, we compute a luminosity versus accretion rate diagram for the setup of spherical accretion onto a Schwarzschild black hole and find good agreement with previous work. Lastly, we revisit the supersonic Bondi Hoyle Lyttleton (BHL) accretion in 2D where we can now present simulations of realistic temperatures, down to T~10^6 K. Here we find that radiation pressure plays an important role, but also that these highly dynamical set-ups push our approximate treatment towards the limit of physical applicability. The main features of radiation hydrodynamics BHL flows manifest as (i) an effective adiabatic index approaching gamma_effective ~ 4/3; (ii) accretion rates two orders of magnitude lower than without radiation pressure; (iii) luminosity estimates around the Eddington limit, hence with an overall radiative efficiency as small as eta ~ 10^{-2}; (iv) strong departures from thermal equilibrium in shocked regions; (v) no appearance of the flip-flop instability. We conclude that the current optically thick approximation to the radiation transfer does give physically substantial improvements over the pure hydro also in set-ups departing from equilibrium, and, once accompanied by an optically thin treatment, is likely to provide a fundamental tool for investigating accretion flows in a large variety of astrophysical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the detectability of dual jets from binary black holes

We revisit the suggestion that dual jets can be produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. By performing independent calculations and by computing the electromagnetic emission in a way which is consistent with estimates using the Poynting flux, we show that a dual-jet structure is present but energetically subdominant with respect to a non-collimated and predominantly quadrupolar emission, which is similar to the one computed when the binary is in electrovacuum. While our findings set serious restrictions on the detectability of dual jets from coalescing binaries, they also increase the chances of detecting an EM counterpart from these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Electromagnetic counterparts from counter-rotating relativistic kicked discs

We show the results of two dimensional general relativistic inviscid and isothermal hydrodynamical simulations comparing the behavior of co-rotating (with respect to the black hole rotation) and counter-rotating circumbinary quasi-Keplerian discs in the post merger phase of a supermassive binary black hole system. While confirming the spiral shock generation within the disc due to the combined effects of mass loss and recoil velocity of the black hole, we find that the maximum luminosity of counter-rotating discs is a factor ~(2-12) higher than in the co-rotating case, depending on the spin of the black hole. On the other hand, the luminosity peak happens ~10 days later with respect to the co-rotating case, for a binary with a total mass M~10^6 M_\odot. Although the global dynamics of counter-rotating discs in the post merger phase of a merging event is very similar to that for co-rotating discs, an important difference has been found. In fact, increasing the spin of the central black hole produces more luminous co-rotating discs while less luminous counter-rotating ones.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Explaining radio emission of magnetars via rotating and oscillating magnetospheres of neutron stars

We investigate the conditions for radio emission in rotating and oscillating magnetars, by focusing on the main physical processes determining the position of their death-lines in the P-\dot{P} diagram, i.e. of those lines that separate the regions where the neutron star may be radio-loud or radio-quiet. After using the general relativistic expression for the electromagnetic scalar potential in the magnetar magnetosphere, we find that larger compactness parameters of the star as well as larger inclination angles between the rotation axis and the magnetic moment produce death-lines well above the majority of known magnetars. This is consistent with the observational evidence of no regular radio emission from the magnetars in the frequency range typical for the ordinary pulsars. On the contrary, when oscillations of the magnetar are taken into account, the death-lines shift downward and the conditions necessary for the generation of radio emission in the magnetosphere are met. Present observations showing a close connection between the burst activity of magnetars and the generation of the radio emission in the magnetar magnetosphere are naturally accounted for within our interpretation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical simulations of high Lundquist number relativistic magnetic reconnection

We present the results of two-dimensional and three-dimensional magnetohydrodynamical numerical simulations of relativistic magnetic reconnection, with particular emphasis on the dynamics of the plasma in a Petschek-type configuration with high Lundquist numbers, S\sim 10^5-10^8. The numerical scheme adopted, allowing for unprecedented accuracy for this type of calculations, is based on high order finite volume and discontinuous Galerkin methods as recently proposed by \citet{Dumbser2009}. The possibility of producing high Lorentz factors is discussed, showing that Lorentz factors close to \sim 4 can be produced for a plasma parameter σ_m=20. Moreover, we find that the Sweet-Parker layers are unstable, generating secondary magnetic islands, but only for S > S_c = 10^8, much larger than what is reported in the Newtonian regime. Finally, the effects of a mildly anisotropic Ohm law are considered in a configuration with a guide magnetic field. Such effects produce only slightly faster reconnection rates and Lorentz factors of about 1% larger with respect to the perfectly isotropic Ohm law.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical simulations of relativistic magnetic reconnection with Galerkin methods

We present the results of two-dimensional magnetohydrodynamical numerical simulations of relativistic magnetic reconnection, with particular emphasis on the dynamics of Petschek-type configurations with high Lundquist numbers, S ~ 10^5-10^8. The numerical scheme adopted, allowing for unprecedented accuracy for this type of calculations, is based on high order finite volume and discontinuous Galerkin methods as recently proposed by Dumbser & Zanotti (2009). The possibility of producing high Lorentz factors is discussed, by studying the effects produced on the dynamics by different magnetization and resistivity regimes. We show that Lorentz factors close to ~4 can be produced for a plasma magnetization parameter sigma=20. Moreover, we find that the Sweet-Parker layers are unstable, generating secondary magnetic islands, but only for S>S_c~10^8, much larger than what is reported in the Newtonian regime.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the development of QPOs in Bondi-Hoyle accretion flows

The numerical investigation of Bondi-Hoyle accretion onto a moving black hole has a long history, both in Newtonian and in general-relativistic physics. By performing new two-dimensional and general-relativistic simulations onto a rotating black hole, we point out a novel feature, namely, that quasi-periodic oscillations (QPOs) are naturally produced in the shock cone that develops in the downstream part of the flow. Because the shock cone in the downstream part of the flow acts as a cavity trapping pressure perturbations, modes with frequencies in the integer ratios 2:1 and 3:1 are easily produced. The frequencies of these modes depend on the black-hole spin and on the properties of the flow, and scale linearly with the inverse of the black-hole mass. Our results may be relevant for explaining the detection of QPOs in Sagittarius A*, once such detection is confirmed by further observations. Finally, we report on the development of the flip-flop instability, which can affect the shock cone under suitable conditions; such an instability has been discussed before in Newtonian simulations but was never found in a relativistic regime.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Oscillations of relativistic axisymmetric tori and implications for modelling kHz-QPOs in neutron-star X-ray binaries

We perform a global linear perturbative analysis, and investigate the oscillation properties of relativistic, non-selfgravitating tori orbiting around neutron stars in the slow rotation limit approximation. Extending the work done in Schwarzschild and Kerr backgrounds, we consider the axisymmetric oscillations of vertically integrated tori in the Hartle-Thorne spacetime. The equilibrium models are constructed by selecting a number of different non-Keplerian distributions of specific angular momentum, allowing for disc sizes $L \sim 0.5 - 600$ gravitational radii. Our results, obtained after solving a global eigenvalue problem to compute the xisymmetric $p$-modes, indicate that such oscillation modes could account with most observed lower ($ν_L$) and upper ($ν_U$) high frequency quasi-periodic oscillations for Sco X-1, and for some Z sources and Atoll sources with $ν_L\ gtrsim 500$ Hz. However, when $ν_L \lesssim 500$ Hz, $p$-modes oscillations do not account for the linear relation $ν_U = A ν_L + B$, $B \neq 0$ between the upper and the lower high frequency quasi-periodic oscillations that are observed in neutron star low-mass X-ray binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

EM counterparts of recoiling black holes: general relativistic simulations of non-Keplerian discs

We investigate the dynamics of a circumbinary disc that responds to the loss of mass and to the recoil velocity of the black hole produced by the merger of a binary system of supermassive black holes. We perform the first two-dimensional general relativistic hydrodynamics simulations of \textit{extended} non-Keplerian discs and employ a new technique to construct a &#34;shock detector&#34;, thus determining the precise location of the shocks produced in the accreting disc by the recoiling black hole. In this way we can study how the properties of the system, such as the spin, mass and recoil velocity of the black hole, affect the mass accretion rate and are imprinted on the electromagnetic emission from these sources. We argue that the estimates of the bremsstrahlung luminosity computed without properly taking into account the radiation transfer yield cooling times that are unrealistically short. At the same time we show, through an approximation based on the relativistic isothermal evolution, that the luminosity produced can reach a peak value above $L \simeq 10^{43} \ {\rm erg/s} $ at about $\sim 30\,{\rm d}$ after the merger of a binary with total mass $M\simeq 10^6 M_\odot$ and persist for several days at values which are a factor of a few smaller. If confirmed by more sophisticated calculations such a signal could indeed lead to an electromagnetic counterpart of the merger of binary black-hole system.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Very High Order $\PNM$ Schemes on Unstructured Meshes for the Resistive Relativistic MHD Equations

In this paper we propose the first better than second order accurate method in space and time for the numerical solution of the resistive relativistic magnetohydrodynamics (RRMHD) equations on unstructured meshes in multiple space dimensions. The nonlinear system under consideration is purely hyperbolic and contains a source term, the one for the evolution of the electric field, that becomes stiff for low values of the resistivity. For the spatial discretization we propose to use high order $\PNM$ schemes as introduced in \cite{Dumbser2008} for hyperbolic conservation laws and a high order accurate unsplit time discretization is achieved using the element-local space-time discontinuous Galerkin approach proposed in \cite{DumbserEnauxToro} for one-dimensional balance laws with stiff source terms. The divergence free character of the magnetic field is accounted for through the divergence cleaning procedure of Dedner et al. \cite{Dedneretal}. To validate our high order method we first solve some numerical test cases for which exact analytical reference solutions are known and we also show numerical convergence studies in the stiff limit of the RRMHD equations using $\PNM$ schemes from third to fifth order of accuracy in space and time. We also present some applications with shock waves such as a classical shock tube problem with different values for the conductivity as well as a relativistic MHD rotor problem and the relativistic equivalent of the Orszag-Tang vortex problem. We have verified that the proposed method can handle equally well the resistive regime and the stiff limit of ideal relativistic MHD. For these reasons it provides a powerful tool for relativistic astrophysical simulations involving the appearance of magnetic reconnection.

0 citations0 reviewsNo code linkedOpen access