Source author record

Oscar Reula

Oscar Reula appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

8 published item(s)

preprint2022arXiv

Skylight: a new code for general-relativistic ray-tracing and radiative transfer in arbitrary spacetimes

To reproduce the observed spectra and light curves originated in the neighborhood of compact objects requires accurate relativistic ray-tracing codes. In this work, we present Skylight, a new numerical code for general-relativistic ray tracing and radiative transfer in arbitrary space-time geometries and coordinate systems. The code is capable of producing images, spectra, and light curves from astrophysical models of compact objects as seen by distant observers. We incorporate two different schemes, namely Monte Carlo radiative transfer, integrating geodesics from the astrophysical region to distant observers, and camera techniques with backwards integration from the observer to the emission region. The code is validated by successfully passing several test cases, among them: thin accretion disks and neutron star hot spot emission.

preprint2016arXiv

Covariant Hyperbolization of Force-free Electrodynamics

Force-Free Flectrodynamics (FFE) is a non-linear system of equations modeling the evolution of the electromagnetic field, in the presence of a magnetically dominated relativistic plasma. This configuration arises on several astrophysical scenarios, which represent exciting laboratories to understand physics in extreme regimes. We show that this system, when restricted to the correct constraint submanifold, is symmetric hyperbolic. In numerical applications is not feasible to keep the system in that submanifold, and so, it is necessary to analyze its structure first in the tangent space of that submanifold and then in a whole neighborhood of it. As already shown by Pfeiffer, a direct (or naive) formulation of this system (in the whole tangent space) results in a weakly hyperbolic system of evolution equations for which well-possednes for the initial value formulation does not follows. Using the generalized symmetric hyperbolic formalism due to Geroch, we introduce here a covariant hyperbolization for the FFE system. In fact, in analogy to the usual Maxwell case, a complete family of hyperbolizers is found, both for the restricted system on the constraint submanifold as well as for a suitably extended system defined in a whole neighborhood of it. A particular symmetrizer among the family is then used to write down the pertaining evolution equations, in a generic ($3+1$)-decomposition on a background spacetime. Interestingly, it turns out that for a particular choice of the lapse and shift functions of the foliation, our symmetrized system reduces to the one found by Pfeiffer. Finally we analyze the characteristic structure of the resulting evolution system.

preprint2015arXiv

Nonlinear electrodynamics as a symmetric hyperbolic system

Nonlinear theories generalizing Maxwell's electromagnetism and arising from a Lagrangian formalism have dispersion relations in which propagation planes factor into null planes corresponding to two effective metrics which depend on the point-wise values of the electromagnetic field. These effective Lorentzian metrics share the null (generically two) directions of the electromagnetic field. We show that, the theory is symmetric hyperbolic if and only if the cones these metrics give rise to have a non-empty intersection. Namely that there exist families of symmetrizers in the sense of Geroch which are positive definite for all covectors in the interior of the cones intersection. Thus, for these theories, the initial value problem is well-posed. We illustrate the power of this approach with several nonlinear models of physical interest such as Born-Infeld, Gauss-Bonnet and Euler-Heisenberg.

preprint2012arXiv

Turbulent flows for relativistic conformal fluids in 2+1 dimensions

We demonstrate that relativistic conformal hydrodynamics in 2+1 dimensions displays a turbulent behaviour which cascades energy to longer wavelengths on both flat and spherical manifolds. Our motivation for this study is to understand the implications for gravitational solutions through the AdS/CFT correspondence. The observed behaviour implies gravitational perturbations of the corresponding black brane/black hole spacetimes (for sufficiently large scales/temperatures) will display a similar cascade towards longer wavelengths.

preprint2011arXiv

Robustness of the Blandford-Znajek mechanism

The Blandford-Znajek mechanism has long been regarded as a key ingredient in models attempting to explain powerful jets in AGNs, quasars, blazzars etc. In such mechanism, energy is extracted from a rotating black hole and dissipated at a load at far distances. In the current work we examine the behaviour of the BZ mechanism with respect to different boundary conditions, revealing the mechanism robustness upon variation of these conditions. Consequently, this work closes a gap in our understanding of this important scenario.

preprint2008arXiv

Beyond ideal MHD: towards a more realistic modeling of relativistic astrophysical plasmas

Many astrophysical processes involving magnetic fields and quasi-stationary processes are well described when assuming the fluid as a perfect conductor. For these systems, the ideal-magnetohydrodynamics (MHD) description captures the dynamics effectively and a number of well-tested techniques exist for its numerical solution. Yet, there are several astrophysical processes involving magnetic fields which are highly dynamical and for which resistive effects can play an important role. The numerical modeling of such non-ideal MHD flows is significantly more challenging as the resistivity is expected to change of several orders of magnitude across the flow and the equations are then either of hyperbolic-parabolic nature or hyperbolic with stiff terms. We here present a novel approach for the solution of these relativistic resistive MHD equations exploiting the properties of implicit-explicit (IMEX) Runge Kutta methods. By examining a number of tests we illustrate the accuracy of our approach under a variety of conditions and highlight its robustness when compared with alternative methods, such as the Strang-splitting. Most importantly, we show that our approach allows one to treat, within a unified framework, both those regions of the flow which are fluid-pressure dominated (such as in the interior of compact objects) and those which are instead magnetic-pressure dominated (such as in their magnetospheres)

preprint2004arXiv

A model problem for the initial-boundary value formulation of Einstein's field equations

In many numerical implementations of the Cauchy formulation of Einstein's field equations one encounters artificial boundaries which raises the issue of specifying boundary conditions. Such conditions have to be chosen carefully. In particular, they should be compatible with the constraints, yield a well posed initial-boundary value formulation and incorporate some physically desirable properties like, for instance, minimizing reflections of gravitational radiation. Motivated by the problem in General Relativity, we analyze a model problem, consisting of a formulation of Maxwell's equations on a spatially compact region of spacetime with timelike boundaries. The form in which the equations are written is such that their structure is very similar to the Einstein-Christoffel symmetric hyperbolic formulations of Einstein's field equations. For this model problem, we specify a family of Sommerfeld-type constraint-preserving boundary conditions and show that the resulting initial-boundary value formulations are well posed. We expect that these results can be generalized to the Einstein-Christoffel formulations of General Relativity, at least in the case of linearizations about a stationary background.