-
Scattering images from autocorrelation functions of P-wave seismic velocity images: the case of Tenerife Island (Canary Islands, Spain)
Authors:
A. García-Yeguas,
A. Sánchez-Alzola,
L. De Siena,
J. Prudencio,
A. Díaz-Moreno,
J. M. Ibáñez
Abstract:
We present a P-wave scattering image of the volcanic structures under Tenerife Island using the autocorrelation functions of P-wave vertical velocity fluctuations. We have applied cluster analysis to total quality factor attenuation (Q inv t) and scattering quality factor attenuation (Q inv PSc) images to interpret the structures in terms of intrinsic and scattering attenuation variations on a 2D…
▽ More
We present a P-wave scattering image of the volcanic structures under Tenerife Island using the autocorrelation functions of P-wave vertical velocity fluctuations. We have applied cluster analysis to total quality factor attenuation (Q inv t) and scattering quality factor attenuation (Q inv PSc) images to interpret the structures in terms of intrinsic and scattering attenuation variations on a 2D plane, corresponding to a depth of 2000 m, and check the robustness of the scattering imaging. The results show that scattering patterns are similar to total attenuation patterns in the South of the island. There are two main areas where patterns differ: at Cañadas-Teide-Pico Viejo Complex high total attenuation and average-to-low scattering values are observed. We interpret the difference as induced by intrinsic attenuation. In the Santiago Ridge Zone (SRZ) region, high scattering values correspond to average total attenuation. In our interpretation, the anomaly is induced by an extended scatterer, geometrically related to the surficial traces of Garachico and El Chinyero historical eruptions and the area of highest seismic activity during the 2004-2008 seismic crises.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
A 3D joint interpretation of magnetotelluric and seismic tomographic models: the case of the volcanic island of Tenerife
Authors:
A. García-Yeguas,
J. Ledo,
P. Piña-Varas,
J. Prudencio,
P. Queralt,
A. Marcuello,
J. M. Ibañez,
B. Benjumea,
A. Sánchez-Alzola,
N. Pérez
Abstract:
In this work we have done a 3D joint interpretation of magnetotelluric and seismic tomography models. Previously we have described different techniques to infer the inner structure of the Earth. We have focused on volcanic regions, specifically on Tenerife Island volcano (Canary Islands, Spain). In this area, magnetotelluric and seismic tomography studies have been done separately. The novelty of…
▽ More
In this work we have done a 3D joint interpretation of magnetotelluric and seismic tomography models. Previously we have described different techniques to infer the inner structure of the Earth. We have focused on volcanic regions, specifically on Tenerife Island volcano (Canary Islands, Spain). In this area, magnetotelluric and seismic tomography studies have been done separately. The novelty of the present work is the combination of both techniques in Tenerife Island. For this aim we have applied Fuzzy Clusters Method at different depths obtaining several clusters or classes. From the results, a geothermal system has been inferred below Teide volcano, in the center of Tenerife Island. An edifice hydrothermally altered and full of fluids is situated below Teide, ending at 600 m below sea level. From this depth the resistivity and VP values increase downwards. We also observe a clay cap structure, a typical feature in geothermal systems related with low resistivity and low VP values.
△ Less
Submitted 4 February, 2024;
originally announced February 2024.
-
Total-variation methods for gravitational-wave denoising: performance tests on Advanced LIGO data
Authors:
Alejandro Torres-Forné,
Elena Cuoco,
Antonio Marquina,
José A. Font,
José M. Ibáñez
Abstract:
We assess total-variation methods to denoise gravitational-wave signals in real noise conditions, by injecting numerical-relativity waveforms from core-collapse supernovae and binary black hole mergers in data from the first observing run of Advanced LIGO. This work is an extension of our previous investigation where only Gaussian noise was used. Since the quality of the results depends on the reg…
▽ More
We assess total-variation methods to denoise gravitational-wave signals in real noise conditions, by injecting numerical-relativity waveforms from core-collapse supernovae and binary black hole mergers in data from the first observing run of Advanced LIGO. This work is an extension of our previous investigation where only Gaussian noise was used. Since the quality of the results depends on the regularization parameter of the model, we perform an heuristic search for the value that produces the best results. We discuss various approaches for the selection of this parameter, either based on the optimal, mean, or multiple values, and compare the results of the denoising upon these choices. Moreover, we also present a machine-learning-informed approach to obtain the Lagrange multiplier of the method through an automatic search. Our results provide further evidence that total-variation methods can be useful in the field of Gravitational-Wave Astronomy as a tool to remove noise.
△ Less
Submitted 19 June, 2018;
originally announced June 2018.
-
Neutron star collapse and gravitational waves with a non-convex equation of state
Authors:
M. A. Aloy,
J. M. Ibáñez,
N. Sanchis-Gual,
M. Obergaulinger,
J. A. Font,
S. Serna,
A. Marquina
Abstract:
The thermodynamical properties of the equation of state (EoS) of high-density matter (above nuclear saturation density) and the possible existence of exotic states such as phase transitions from nuclear/hadronic matter into quark-gluon plasma, or the appearance of hyperons, may critically influence the stability and dynamics of compact relativistic stars. From a theoretical point of view, establis…
▽ More
The thermodynamical properties of the equation of state (EoS) of high-density matter (above nuclear saturation density) and the possible existence of exotic states such as phase transitions from nuclear/hadronic matter into quark-gluon plasma, or the appearance of hyperons, may critically influence the stability and dynamics of compact relativistic stars. From a theoretical point of view, establishing the existence of those states requires the analysis of the `convexity' of the EoS. We show indications of the existence of regions in the dense-matter EoS where the thermodynamics may be non-convex as a result of a non-monotonic dependence of the sound speed with the rest-mass density. When this happens, non-conventional dynamics may develop. In this paper we investigate the effects of a phenomenological, non-convex EoS on the equilibrium structure of stable compact stars and on the dynamics of unstable neutron stars that collapse gravitationally to black holes, both for spherically symmetric and uniformly-rotating configurations. We show how the dynamics of the collapse with a non-convex EoS departs from the convex case, leaving distinctive imprints on the gravitational waveforms. The astrophysical significance of these results for microphysical EoSs is discussed.
△ Less
Submitted 24 January, 2019; v1 submitted 8 June, 2018;
originally announced June 2018.
-
Anomalous dynamics triggered by a non-convex equation of state in relativistic flows
Authors:
J. M. Ibáñez,
A. Marquina,
S. Serna,
M. A. Aloy
Abstract:
The non-monotonicity of the local speed of sound in dense matter at baryon number densities much higher than the nuclear saturation density ($n_0 \approx 0.16\,$fm$^{-3}$) suggests the possible existence of a non-convex thermodynamics which will lead to a non-convex dynamics. Here, we explore the rich and complex dynamics that an equation of state (EoS) with non-convex regions in the pressure-dens…
▽ More
The non-monotonicity of the local speed of sound in dense matter at baryon number densities much higher than the nuclear saturation density ($n_0 \approx 0.16\,$fm$^{-3}$) suggests the possible existence of a non-convex thermodynamics which will lead to a non-convex dynamics. Here, we explore the rich and complex dynamics that an equation of state (EoS) with non-convex regions in the pressure-density plane may develop as a result of genuinely relativistic effects, without a classical counterpart. To this end, we have introduced a phenomenological EoS, whose parameters can be restricted heeding to causality and thermodynamic stability constraints. This EoS shall be regarded as a toy-model with which we may mimic realistic (and far more complex) EoS of practical use in the realm of Relativistic Hydrodynamics.
△ Less
Submitted 7 January, 2018; v1 submitted 8 December, 2017;
originally announced December 2017.
-
Tracing accretion variability of high-mass YSOs via light echoes
Authors:
B. Stecklum,
S. Heese,
S. Wolf,
A. Caratti o Garatti,
J. M. Ibanez,
H. Linz
Abstract:
There is growing evidence for disk-mediated accretion being the dominant mode of star formation across nearly the whole stellar mass spectrum. The stochastic nature of this process has been realized which implies an inherent source variability. It can be traced more easily for low-mass YSOs (LMYSOs) since high-mass YSOs (HMYSOs) are still embedded even when reaching the ZAMS. While variable reflec…
▽ More
There is growing evidence for disk-mediated accretion being the dominant mode of star formation across nearly the whole stellar mass spectrum. The stochastic nature of this process has been realized which implies an inherent source variability. It can be traced more easily for low-mass YSOs (LMYSOs) since high-mass YSOs (HMYSOs) are still embedded even when reaching the ZAMS. While variable reflection nebulae around LMYSOs were among the earliest signs of star formation, little is known on the variability of scattered light from embedded clusters, the birthplaces of HMYSOs. Since the few most massive stars dominate this emission, their variability is literally reflected in scattered light. Moreover, because of their high luminosity, for a given ambient dust density and source distance, the associated nebulosities are much larger than those of LMYSOs. In this case, the light travel time becomes substantial. So the apparent brightness distribution constitutes a light echo, shaped by both the HMYSO variability history and the spatial distribution of the scattering medium. We report on early results of a NIR variability study of HMYSOs associated with Class II methanol masers which aims at revealing a possible correlation between maser flux density and infrared brightness. Additionally, relevant findings for the eruptive HMYSO S255IR-NIRS3 are presented.
△ Less
Submitted 4 December, 2017;
originally announced December 2017.
-
Disk-mediated accretion burst in a high-mass young stellar object
Authors:
A. Caratti o Garatti,
B. Stecklum,
R. Garcia Lopez,
J. Eislöffel,
T. P. Ray,
A. Sanna,
R. Cesaroni,
C. M. Walmsley,
R. D. Oudmaijer,
W. J. de Wit,
L. Moscadelli,
J. Greiner,
A. Krabbe,
C. Fischer,
R. Klein,
J. M. Ibañez
Abstract:
Solar-mass stars form via circumstellar disk accretion (disk-mediated accretion). Recent findings indicate that this process is likely episodic in the form of accretion bursts, possibly caused by disk fragmentation. Although it cannot be ruled out that high-mass young stellar objects (HMYSOs; $M>$8 M$_\odot$, $L_{bol}>$5$\times$10$^3$ L$_\odot$) arise from the coalescence of their low-mass brethre…
▽ More
Solar-mass stars form via circumstellar disk accretion (disk-mediated accretion). Recent findings indicate that this process is likely episodic in the form of accretion bursts, possibly caused by disk fragmentation. Although it cannot be ruled out that high-mass young stellar objects (HMYSOs; $M>$8 M$_\odot$, $L_{bol}>$5$\times$10$^3$ L$_\odot$) arise from the coalescence of their low-mass brethren, latest results suggest that they more likely form via disks. Accordingly, disk-mediated accretion bursts should occur. Here we report on the discovery of the first disk-mediated accretion burst from a $\sim$20 M$_\odot$ HMYSO. Our near-infrared images show the brightening of the central source and its outflow cavities. Near-infrared spectroscopy reveals emission lines typical of accretion bursts in low-mass protostars, but orders of magnitude more luminous. Moreover, the energy released and the inferred mass-accretion rate are also orders of magnitude larger. Our results identify disk accretion as the common mechanism of star formation across the entire stellar mass spectrum.
△ Less
Submitted 9 April, 2017;
originally announced April 2017.
-
Denoising of gravitational wave signals via dictionary learning algorithms
Authors:
Alejandro Torres-Forné,
Antonio Marquina,
José A. Font,
José M. Ibáñez
Abstract:
Gravitational wave astronomy has become a reality after the historical detections accomplished during the first observing run of the two advanced LIGO detectors. In the following years, the number of detections is expected to increase significantly with the full commissioning of the advanced LIGO, advanced Virgo and KAGRA detectors. The development of sophisticated data analysis techniques to impr…
▽ More
Gravitational wave astronomy has become a reality after the historical detections accomplished during the first observing run of the two advanced LIGO detectors. In the following years, the number of detections is expected to increase significantly with the full commissioning of the advanced LIGO, advanced Virgo and KAGRA detectors. The development of sophisticated data analysis techniques to improve the opportunities of detection for low signal-to-noise-ratio events is hence a most crucial effort. We present in this paper one such technique, dictionary-learning algorithms, which have been extensively developed in the last few years and successfully applied mostly in the context of image processing. However, to the best of our knowledge, such algorithms have not yet been employed to denoise gravitational wave signals. By building dictionaries from numerical relativity templates of both, binary black holes mergers and bursts of rotational core collapse, we show how machine-learning algorithms based on dictionaries can be also successfully applied for gravitational wave denoising. We use a subset of signals from both catalogs, embedded in non-white Gaussian noise, to assess our techniques with a large sample of tests and to find the best model parameters. The application of our method to the actual signal GW150914 shows promising results. Dictionary-learning algorithms could be a complementary addition to the gravitational wave data analysis toolkit. They may be used to extract signals from noise and to infer physical parameters if the data are in good enough agreement with the morphology of the dictionary atoms.
△ Less
Submitted 5 December, 2016;
originally announced December 2016.
-
Denoising of gravitational-wave signal GW150914 via total-variation methods
Authors:
Alejandro Torres-Forné,
Antonio Marquina,
José A. Font,
José M. Ibáñez
Abstract:
We apply a regularized Rudin-Osher-Fatemi total variation (TV) method to denoise the transient gravitational wave signal GW150914. We have previously applied TV techniques to denoise numerically generated grav- itational waves embedded in additive Gaussian noise, obtaining satisfactory results irrespective of the signal morphology or astrophysical origin. We find that the non-Gaussian, non-station…
▽ More
We apply a regularized Rudin-Osher-Fatemi total variation (TV) method to denoise the transient gravitational wave signal GW150914. We have previously applied TV techniques to denoise numerically generated grav- itational waves embedded in additive Gaussian noise, obtaining satisfactory results irrespective of the signal morphology or astrophysical origin. We find that the non-Gaussian, non-stationary noise from the gravitational wave event GW150914 can also be successfully removed with TV-denoising methods. The quality of the de- noised waveform is comparable to that obtained with the Bayesian approach used in the discovery paper [1]. TV-denoising techniques may thus offer an additional viable approach for waveform reconstruction.
△ Less
Submitted 22 February, 2016;
originally announced February 2016.
-
Total-variation-based methods for gravitational wave denoising
Authors:
Alejandro Torres,
Antonio Marquina,
José A. Font,
José M. Ibáñez
Abstract:
We describe new methods for denoising and detection of gravitational waves embedded in additive Gaussian noise. The methods are based on Total Variation denoising algorithms. These algorithms, which do not need any a priori information about the signals, have been originally developed and fully tested in the context of image processing. To illustrate the capabilities of our methods we apply them t…
▽ More
We describe new methods for denoising and detection of gravitational waves embedded in additive Gaussian noise. The methods are based on Total Variation denoising algorithms. These algorithms, which do not need any a priori information about the signals, have been originally developed and fully tested in the context of image processing. To illustrate the capabilities of our methods we apply them to two different types of numerically-simulated gravitational wave signals, namely bursts produced from the core collapse of rotating stars and waveforms from binary black hole mergers. We explore the parameter space of the methods to find the set of values best suited for denoising gravitational wave signals under different conditions such as waveform type and signal-to-noise ratio. Our results show that noise from gravitational wave signals can be successfully removed with our techniques, irrespective of the signal morphology or astrophysical origin. We also combine our methods with spectrograms and show how those can be used simultaneously with other common techniques in gravitational wave data analysis to improve the chances of detection.
△ Less
Submitted 28 September, 2014;
originally announced September 2014.
-
On the convexity of Relativistic Hydrodynamics
Authors:
José María Ibáñez,
Isabel Cordero-Carrión,
José María Martí,
Juan Antonio Miralles
Abstract:
The relativistic hydrodynamic system of equations for a perfect fluid obeying a causal equation of state is hyperbolic (Anile 1989 {\it Relativistic Fluids and Magneto-Fluids} (Cambridge: Cambridge University Press)). In this report, we derive the conditions for this system to be convex in terms of the fundamental derivative of the equation of state (Menikoff and Plohr 1989 {\it Rev. Mod. Phys.} {…
▽ More
The relativistic hydrodynamic system of equations for a perfect fluid obeying a causal equation of state is hyperbolic (Anile 1989 {\it Relativistic Fluids and Magneto-Fluids} (Cambridge: Cambridge University Press)). In this report, we derive the conditions for this system to be convex in terms of the fundamental derivative of the equation of state (Menikoff and Plohr 1989 {\it Rev. Mod. Phys.} {\bf 61} 75). The classical limit is recovered.
△ Less
Submitted 15 February, 2013;
originally announced February 2013.
-
On numerical relativistic hydrodynamics and barotropic equations of state
Authors:
J. M. Ibáñez,
I. Cordero-Carrión,
J. A. Miralles
Abstract:
The characteristic formulation of the relativistic hydrodynamic equations (Donat et al 1998 J. Comput. Phys. 146 58), which has been implemented in many relativistic hydro-codes that make use of Godunov-type methods, has to be slightly modified in the case of evolving barotropic flows. For a barotropic equation of state, a removable singularity appears in one of the eigenvectors. The singularity c…
▽ More
The characteristic formulation of the relativistic hydrodynamic equations (Donat et al 1998 J. Comput. Phys. 146 58), which has been implemented in many relativistic hydro-codes that make use of Godunov-type methods, has to be slightly modified in the case of evolving barotropic flows. For a barotropic equation of state, a removable singularity appears in one of the eigenvectors. The singularity can be avoided by means of a simple renormalization which makes the system of eigenvectors well defined and complete. An alternative strategy for the particular case of barotropic flows is discussed.
△ Less
Submitted 27 June, 2012; v1 submitted 26 June, 2012;
originally announced June 2012.
-
Maximal slicings in spherical symmetry: local existence and construction
Authors:
Isabel Cordero-Carrión,
José María Ibáñez,
Juan Antonio Morales-Lladosa
Abstract:
We show that any spherically symmetric spacetime locally admits a maximal spacelike slicing and we give a procedure allowing its construction. The construction procedure that we have designed is based on purely geometrical arguments and, in practice, leads to solve a decoupled system of first order quasi-linear partial differential equations. We have explicitly built up maximal foliations in Minko…
▽ More
We show that any spherically symmetric spacetime locally admits a maximal spacelike slicing and we give a procedure allowing its construction. The construction procedure that we have designed is based on purely geometrical arguments and, in practice, leads to solve a decoupled system of first order quasi-linear partial differential equations. We have explicitly built up maximal foliations in Minkowski and Friedmann spacetimes. Our approach admits further generalizations and efficient computational implementation. As by product, we suggest some applications of our work in the task of calibrating Numerical Relativity complex codes, usually written in Cartesian coordinates.
△ Less
Submitted 11 November, 2011;
originally announced November 2011.
-
Gravitational waves in dynamical spacetimes with matter content in the Fully Constrained Formulation
Authors:
Isabel Cordero-Carrión,
Pablo Cerdá-Durán,
José María Ibáñez
Abstract:
The Fully Constrained Formulation (FCF) of General Relativity is a novel framework introduced as an alternative to the hyperbolic formulations traditionally used in numerical relativity. The FCF equations form a hybrid elliptic-hyperbolic system of equations including explicitly the constraints. We present an implicit-explicit numerical algorithm to solve the hyperbolic part, whereas the elliptic…
▽ More
The Fully Constrained Formulation (FCF) of General Relativity is a novel framework introduced as an alternative to the hyperbolic formulations traditionally used in numerical relativity. The FCF equations form a hybrid elliptic-hyperbolic system of equations including explicitly the constraints. We present an implicit-explicit numerical algorithm to solve the hyperbolic part, whereas the elliptic sector shares the form and properties with the well known Conformally Flat Condition (CFC) approximation. We show the stability andconvergence properties of the numerical scheme with numerical simulations of vacuum solutions. We have performed the first numerical evolutions of the coupled system of hydrodynamics and Einstein equations within FCF. As a proof of principle of the viability of the formalism, we present 2D axisymmetric simulations of an oscillating neutron star. In order to simplify the analysis we have neglected the back-reaction of the gravitational waves into the dynamics, which is small (<2 %) for the system considered in this work. We use spherical coordinates grids which are well adapted for simulations of stars and allow for extended grids that marginally reach the wave zone. We have extracted the gravitational wave signature and compared to the Newtonian quadrupole and hexadecapole formulae. Both extraction methods show agreement within the numerical errors and the approximations used (~30 %).
△ Less
Submitted 3 April, 2012; v1 submitted 2 August, 2011;
originally announced August 2011.
-
A Roe-type Riemann solver based on the spectral decomposition of the equations of Relativistic Magnetohydrodynamics
Authors:
J. M. Ibáñez,
M. A. Aloy,
P. Mimica,
L. Antón,
J. A. Miralles,
J. M. Martí
Abstract:
In a recent paper (Antón et al. 2010) we have derived sets of right and left eigenvectors of the Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. We present a summary of the main steps followed in the above derivation and the numerical experiments carried out with the linearized (Roe-type) Riemann solver we ha…
▽ More
In a recent paper (Antón et al. 2010) we have derived sets of right and left eigenvectors of the Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. We present a summary of the main steps followed in the above derivation and the numerical experiments carried out with the linearized (Roe-type) Riemann solver we have developed, and some note on the (non-)convex character of the relativistic MHD equations.
△ Less
Submitted 11 December, 2010;
originally announced December 2010.
-
On the local existence of maximal slicings in spherically symmetric spacetimes
Authors:
Isabel Cordero-Carrión,
José María Ibáñez,
Juan Antonio Morales-Lladosa
Abstract:
In this talk we show that any spherically symmetric spacetime admits locally a maximal spacelike slicing. The above condition is reduced to solve a decoupled system of first order quasi-linear partial differential equations. The solution may be accomplished analytical or numerically. We provide a general procedure to construct such maximal slicings.
In this talk we show that any spherically symmetric spacetime admits locally a maximal spacelike slicing. The above condition is reduced to solve a decoupled system of first order quasi-linear partial differential equations. The solution may be accomplished analytical or numerically. We provide a general procedure to construct such maximal slicings.
△ Less
Submitted 4 March, 2010;
originally announced March 2010.
-
Dynamical spacetimes and gravitational radiation in a Fully Constrained Formulation
Authors:
Isabel Cordero-Carrión,
Pablo Cerdá-Durán,
José María Ibáñez
Abstract:
This contribution summarizes the recent work carried out to analyze the behavior of the hyperbolic sector of the Fully Constrained Formulation (FCF) derived in Bonazzola et al. 2004. The numerical experiments presented here allows one to be confident in the performances of the upgraded version of CoCoNuT's code by replacing the Conformally Flat Condition (CFC) approximation of the Einstein equat…
▽ More
This contribution summarizes the recent work carried out to analyze the behavior of the hyperbolic sector of the Fully Constrained Formulation (FCF) derived in Bonazzola et al. 2004. The numerical experiments presented here allows one to be confident in the performances of the upgraded version of CoCoNuT's code by replacing the Conformally Flat Condition (CFC) approximation of the Einstein equations by the FCF.
△ Less
Submitted 3 March, 2010;
originally announced March 2010.
-
Relativistic Magnetohydrodynamics: Renormalized eigenvectors and full wave decomposition Riemann solver
Authors:
L. Anton,
J. A. Miralles,
J. M. Marti,
J. M. Ibanez,
M. A. Aloy,
P. Mimica
Abstract:
We obtain renormalized sets of right and left eigenvectors of the flux vector Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. The renormalization procedure relies on the characterization of the degeneracy types in terms of the normal and tangential components of the magnetic field to the wavefront in the fl…
▽ More
We obtain renormalized sets of right and left eigenvectors of the flux vector Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. The renormalization procedure relies on the characterization of the degeneracy types in terms of the normal and tangential components of the magnetic field to the wavefront in the fluid rest frame. Proper expressions of the renormalized eigenvectors in conserved variables are obtained through the corresponding matrix transformations. Our work completes previous analysis that present different sets of right eigenvectors for non-degenerate and degenerate states, and can be seen as a relativistic generalization of earlier work performed in classical MHD. Based on the full wave decomposition (FWD) provided by the the renormalized set of eigenvectors in conserved variables, we have also developed a linearized (Roe-type) Riemann solver. Extensive testing against one- and two-dimensional standard numerical problems allows us to conclude that our solver is very robust. When compared with a family of simpler solvers that avoid the knowledge of the full characteristic structure of the equations in the computation of the numerical fluxes, our solver turns out to be less diffusive than HLL and HLLC, and comparable in accuracy to the HLLD solver. The amount of operations needed by the FWD solver makes it less efficient computationally than those of the HLL family in one-dimensional problems. However its relative efficiency increases in multidimensional simulations.
△ Less
Submitted 23 December, 2009;
originally announced December 2009.
-
Mathematical Issues in a Fully-Constrained Formulation of Einstein Equations
Authors:
I. Cordero-Carrión,
J. M. Ibáñez,
E. Gourgoulhon,
J. L. Jaramillo,
J. Novak
Abstract:
Bonazzola, Gourgoulhon, Grandclément, and Novak [Phys. Rev. D {\bf 70}, 104007 (2004)] proposed a new formulation for 3+1 numerical relativity. Einstein equations result, according to that formalism, in a coupled elliptic-hyperbolic system. We have carried out a preliminary analysis of the mathematical structure of that system, in particular focusing on the equations governing the evolution for…
▽ More
Bonazzola, Gourgoulhon, Grandclément, and Novak [Phys. Rev. D {\bf 70}, 104007 (2004)] proposed a new formulation for 3+1 numerical relativity. Einstein equations result, according to that formalism, in a coupled elliptic-hyperbolic system. We have carried out a preliminary analysis of the mathematical structure of that system, in particular focusing on the equations governing the evolution for the deviation of a conformal metric from a flat fiducial one. The choice of a Dirac's gauge for the spatial coordinates guarantees the mathematical characterization of that system as a (strongly) hyperbolic system of conservation laws. In the presence of boundaries, this characterization also depends on the boundary conditions for the shift vector in the elliptic subsystem. This interplay between the hyperbolic and elliptic parts of the complete evolution system is used to assess the prescription of inner boundary conditions for the hyperbolic part when using an excision approach to black hole spacetime evolutions.
△ Less
Submitted 8 April, 2008; v1 submitted 21 February, 2008;
originally announced February 2008.
-
Trapping Horizons as inner boundary conditions for black hole spacetimes
Authors:
J. L. Jaramillo,
E. Gourgoulhon,
I. Cordero-Carrion,
J. M. Ibanez
Abstract:
We present a set of inner boundary conditions for the numerical construction of dynamical black hole space-times, when employing a 3+1 constrained evolution scheme and an excision technique. These inner boundary conditions are heuristically motivated by the dynamical trapping horizon framework and are enforced in an elliptic subsystem of the full Einstein equation. In the stationary limit they r…
▽ More
We present a set of inner boundary conditions for the numerical construction of dynamical black hole space-times, when employing a 3+1 constrained evolution scheme and an excision technique. These inner boundary conditions are heuristically motivated by the dynamical trapping horizon framework and are enforced in an elliptic subsystem of the full Einstein equation. In the stationary limit they reduce to existing isolated horizon boundary conditions. A characteristic analysis completes the discussion of inner boundary conditions for the radiative modes.
△ Less
Submitted 18 September, 2007;
originally announced September 2007.
-
The exact solution of the Riemann problem in relativistic MHD with tangential magnetic fields
Authors:
Roberto Romero,
Jose M. Marti,
Jose A. Pons,
Jose M. Ibanez,
Juan A. Miralles
Abstract:
We have extended the procedure to find the exact solution of the Riemann problem in relativistic hydrodynamics to a particular case of relativistic magnetohydrodynamics in which the magnetic field of the initial states is tangential to the discontinuity and orthogonal to the flow velocity. The wave pattern produced after the break up of the initial discontinuity is analogous to the non--magnetic…
▽ More
We have extended the procedure to find the exact solution of the Riemann problem in relativistic hydrodynamics to a particular case of relativistic magnetohydrodynamics in which the magnetic field of the initial states is tangential to the discontinuity and orthogonal to the flow velocity. The wave pattern produced after the break up of the initial discontinuity is analogous to the non--magnetic case and we show that the problem can be understood as a purely relativistic hydrodynamical problem with a modified equation of state. The new degree of freedom introduced by the non-zero component of the magnetic field results in interesting effects consisting in the change of the wave patterns for given initial thermodynamical states, in a similar way to the effects arising from the introduction of tangential velocities. Secondly, when the magnetic field dominates the thermodynamical pressure and energy, the wave speeds approach the speed of light leading to fast shocks and fast and arbitrarily thin rarefaction waves. Our approach is the first non-trivial exact solution of a Riemann problem in relativistic magnetohydrodynamics and it can also be of great interest to test numerical codes against known analytical or exact solutions.
△ Less
Submitted 22 June, 2005;
originally announced June 2005.
-
Numerical 3+1 general relativistic magnetohydrodynamics: a local characteristic approach
Authors:
L. Anton,
O. Zanotti,
J. A. Miralles,
J. M. Marti,
J. M. Ibanez,
J. A. Font,
J. A. Pons
Abstract:
We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbol…
▽ More
We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbolic character in the solution procedure. All theoretical ingredients necessary to build up high-resolution shock-capturing schemes based on the solution of local Riemann problems (i.e. Godunov-type schemes) are described. In particular, we use a renormalized set of regular eigenvectors of the flux Jacobians of the relativistic magnetohydrodynamics equations. In addition, the paper describes a procedure based on the equivalence principle of general relativity that allows the use of Riemann solvers designed for special relativistic magnetohydrodynamics in GRMHD. Our formulation and numerical methodology are assessed by performing various test simulations recently considered by different authors. These include magnetized shock tubes, spherical accretion onto a Schwarzschild black hole, equatorial accretion onto a Kerr black hole, and magnetized thick accretion disks around a black hole prone to the magnetorotational instability.
△ Less
Submitted 2 June, 2005;
originally announced June 2005.
-
CFC+: Improved dynamics and gravitational waveforms from relativistic core collapse simulations
Authors:
P. Cerda-Duran,
G. Faye,
H. Dimmelmeier,
J. A. Font,
J. M. Ibanez,
E. Mueller,
G. Schaefer
Abstract:
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 (C…
▽ More
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 the coupled system of metric and fluid equations using the 3+1 formalism, specialized to a new framework for the gravitational field equations which we call CFC+. In this approach we add new degrees of freedom to the original CFC equations, which extend them by terms of second post-Newtonian order. The corrections for CFC+ are computed solving a system of elliptic linear equations. The new formalism is assessed with time evolutions of both rotating neutron stars in equilibrium and gravitational core collapse of rotating polytropes. Gravitational wave signals for a comprehensive sample of collapse models are extracted using either the quadrupole formula or directly from the metric. We discuss our results on the dynamics and the gravitational wave emission through a detailed comparison between CFC and CFC+ simulations. The main conclusion is that, for the neutron star spacetimes analyzed in the present work, no significant differences are found among CFC, CFC+, and full general relativity, which highlights the suitability of the former.
△ Less
Submitted 23 December, 2004;
originally announced December 2004.
-
Assessment of a high-resolution central scheme for the solution of the relativistic hydrodynamics equations
Authors:
Arturo Lucas-Serrano,
Jose A. Font,
Jose M. Ibanez,
Jose M. Marti
Abstract:
We assess the suitability of a recent high-resolution central scheme developed by Kurganov & Tadmor (2000) for the solution of the relativistic hydrodynamics equations. The novelty of this approach relies on the absence of Riemann solvers in the solution procedure. The computations we present are performed in one and two spatial dimensions in Minkowski spacetime. Standard numerical experiments s…
▽ More
We assess the suitability of a recent high-resolution central scheme developed by Kurganov & Tadmor (2000) for the solution of the relativistic hydrodynamics equations. The novelty of this approach relies on the absence of Riemann solvers in the solution procedure. The computations we present are performed in one and two spatial dimensions in Minkowski spacetime. Standard numerical experiments such as shock tubes and the relativistic flat-faced step test are performed. As an astrophysical application the article includes two-dimensional simulations of the propagation of relativistic jets using both Cartesian and cylindrical coordinates. The simulations reported clearly show the capabilities of the numerical scheme to yield satisfactory results, with an accuracy comparable to that obtained by the so-called high-resolution shock-capturing schemes based upon Riemann solvers (Godunov-type schemes), even well inside the ultrarelativistic regime. Such central scheme can be straightforwardly applied to hyperbolic systems of conservation laws for which the characteristic structure is not explicitly known, or in cases where the exact solution of the Riemann problem is prohibitively expensive to compute numerically. Finally, we present comparisons with results obtained using various Godunov-type schemes as well as with those obtained using other high-resolution central schemes which have recently been reported in the literature.
△ Less
Submitted 27 July, 2004;
originally announced July 2004.
-
"Mariage des Maillages": A new numerical approach for 3D relativistic core collapse simulations
Authors:
Harald Dimmelmeier,
Jerome Novak,
Jose A. Font,
Jose M. Ibanez,
Ewald Mueller
Abstract:
We present a new 3D general relativistic hydrodynamics code for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. It uses spectral methods for solving the metric equations, assuming the conformal flatness approximation for the three-metric. The matter equations are solved by high-resolution shock-capturing schemes. We…
▽ More
We present a new 3D general relativistic hydrodynamics code for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. It uses spectral methods for solving the metric equations, assuming the conformal flatness approximation for the three-metric. The matter equations are solved by high-resolution shock-capturing schemes. We demonstrate that the combination of a finite difference grid and a spectral grid can be successfully accomplished. This "Mariage des Maillages" (French for grid wedding) approach results in high accuracy of the metric solver and allows for fully 3D applications using computationally affordable resources, and ensures long term numerical stability of the evolution. We compare our new approach to two other, finite difference based, methods to solve the metric equations. A variety of tests in 2D and 3D is presented, involving highly perturbed neutron star spacetimes and (axisymmetric) stellar core collapse, demonstrating the ability to handle spacetimes with and without symmetries in strong gravity. These tests are also employed to assess gravitational waveform extraction, which is based on the quadrupole formula.
△ Less
Submitted 2 March, 2005; v1 submitted 8 July, 2004;
originally announced July 2004.
-
Jet stability and the generation of superluminal and stationary components
Authors:
I. Agudo,
J. L. Gomez,
J. M. Marti,
J. M Ibanez,
A. P. Marscher,
A. Alberdi,
M. A. Aloy,
P. E. Hardee
Abstract:
We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet r…
▽ More
We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component, instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness, and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component.
△ Less
Submitted 11 January, 2001;
originally announced January 2001.
-
Hyperbolic character of the angular moment equations of radiative transfer and numerical methods
Authors:
Jose A. Pons,
J. Ma. Ibanez,
Juan A. Miralles
Abstract:
We study the mathematical character of the angular moment equations of radiative transfer in spherical symmetry and conclude that the system is hyperbolic for general forms of the closure relation found in the literature. Hyperbolicity and causality preservation lead to mathematical conditions allowing to establish a useful characterization of the closure relations. We apply numerical methods sp…
▽ More
We study the mathematical character of the angular moment equations of radiative transfer in spherical symmetry and conclude that the system is hyperbolic for general forms of the closure relation found in the literature. Hyperbolicity and causality preservation lead to mathematical conditions allowing to establish a useful characterization of the closure relations. We apply numerical methods specifically designed to solve hyperbolic systems of conservation laws (the so-called Godunov-type methods), to calculate numerical solutions of the radiation transport equations in a static background. The feasibility of the method in any kind of regime, from diffusion to free-streaming, is demonstrated by a number of numerical tests and the effect of the choice of the closure relation on the results is discussed.
△ Less
Submitted 15 May, 2000;
originally announced May 2000.
-
Gravitational waves from the collapse and bounce of a stellar core in tensor-scalar gravity
Authors:
J. Novak,
J. M. Ibanez
Abstract:
Tensor-scalar theory of gravity allows the generation of gravitational waves from astrophysical sources, like Supernovæ, even in the spherical case. That motivated us to study the collapse of a degenerate stellar core, within tensor-scalar gravity, leading to the formation of a neutron star through a bounce and the formation of a shock. We discuss in this paper the effects of the scalar field on…
▽ More
Tensor-scalar theory of gravity allows the generation of gravitational waves from astrophysical sources, like Supernovæ, even in the spherical case. That motivated us to study the collapse of a degenerate stellar core, within tensor-scalar gravity, leading to the formation of a neutron star through a bounce and the formation of a shock. We discuss in this paper the effects of the scalar field on the evolution of the system, as well as the appearance of strong non-perturbative effects of this scalar field (the so-called ``spontaneous scalarization''). As a main result, we describe the resulting gravitational monopolar radiation (form and amplitude) and discuss the possibility of its detection by the gravitational detectors currently under construction, taking into account the existing constraints on the scalar field.
From the numerical point of view it is worthy to point out that we have developed a combined code which uses pseudo-spectral methods, for the evolution of the scalar field, and High Resolution Shock-Capturing schemes, for the evolution of the hydrodynamical system. Although this code has been used to integrate the field equations of that theory of gravity, in the spherically symmetric case, a by-product of the present work is to gain experience for an ulterior extension to multidimensional problems in Numerical Relativity of such numerical strategy.
△ Less
Submitted 16 November, 1999;
originally announced November 1999.
-
Relativistic Jets from Collapsars
Authors:
Miguel A. Aloy,
Ewald Mueller,
Jose M. Ibanez,
Jose M. Marti,
Andrew MacFadyen
Abstract:
Using a collapsar progenitor model of MacFadyen & Woosley we have simulated the propagation of an axisymmetric jet through a collapsing rotating massive star with the GENESIS multi-dimensional relativistic hydrodynamic code. The jet forms as a consequence of an assumed (constant or variable) energy deposition in the range $10^{50}$ erg s$^{-1}$ to $10^{51}$ erg s$^{-1}$ within a $30^{\circ}$ con…
▽ More
Using a collapsar progenitor model of MacFadyen & Woosley we have simulated the propagation of an axisymmetric jet through a collapsing rotating massive star with the GENESIS multi-dimensional relativistic hydrodynamic code. The jet forms as a consequence of an assumed (constant or variable) energy deposition in the range $10^{50}$ erg s$^{-1}$ to $10^{51}$ erg s$^{-1}$ within a $30^{\circ}$ cone around the rotation axis. The jet flow is strongly beamed ($\la$ few degrees), spatially inhomogeneous, and time dependent. The jet reaches the surface of the stellar progenitor ($R_{\ast} = 2.98 \times 10^{10} $cm) intact. At breakout the maximum Lorentz factor of the jet flow is 33. After breakout the jet accelerates into the circumstellar medium, whose density is assumed to decrease exponentially and then being constant $ρ_{\rm ext} = 10^{-5}$ gcm$^{-3}$. Outside the star the flow begins to expand also laterally ($v \sim c$), but the beam remains very well collimated. At a distance of $2.54 R_{\ast}$, where the simulation ends, the Lorentz factor has increased to 44.
△ Less
Submitted 17 January, 2000; v1 submitted 5 November, 1999;
originally announced November 1999.
-
Riemann Solvers in General Relativistic Hydrodynamics
Authors:
J. M. Ibanez,
M. A. Aloy,
J. A. Font,
J. M. Marti,
J. A. Miralles,
J. A. Pons
Abstract:
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.,…
▽ More
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 terms of which the spacetime metric is locally Minkowskian and permits an accurate description of numerical general relativistic hydrodynamics.
△ Less
Submitted 3 November, 1999;
originally announced November 1999.
-
Relativistic Jets from Collapsars
Authors:
Miguel A. Aloy,
Ewald Mueller,
Jose M. Ibanez,
Jose M. Marti,
Andrew MacFadyen
Abstract:
We have studied the relativistic beamed outflow proposed to occur in the collapsar model of gamma-ray bursts. A jet forms as a consequence of an assumed energy deposition of $\sim 10^{50}- 10^{51}$ erg/s within a $30^{\circ}$ cone around the rotation axis of the progenitor star. The generated jet flow is strongly beamed ($\la$ few degrees) and reaches the surface of the stellar progenitor (r…
▽ More
We have studied the relativistic beamed outflow proposed to occur in the collapsar model of gamma-ray bursts. A jet forms as a consequence of an assumed energy deposition of $\sim 10^{50}- 10^{51}$ erg/s within a $30^{\circ}$ cone around the rotation axis of the progenitor star. The generated jet flow is strongly beamed ($\la$ few degrees) and reaches the surface of the stellar progenitor (r $\approx 3 10^{10} $cm) intact. At break-out the maximum Lorentz factor of the jet flow is about 33. Simulations have been performed with the GENESIS multi-dimensional relativistic hydrodynamic code.
△ Less
Submitted 26 October, 1999;
originally announced October 1999.
-
Simulations of Relativistic Jets with GENESIS
Authors:
Miguel A. Aloy,
Jose M. Ibanez,
Jose M. Marti,
Jose L. Gomez,
Ewald Mueller
Abstract:
The multidimensional relativistic hydrodynamical code GENESIS has been used to obtain first results of {\it 3D} simulations of relativistic jets. We have studied the influence of a slight perturbation of the injection velocity field on the morphodynamics of otherwise axisymmetric relativistic jets.
The multidimensional relativistic hydrodynamical code GENESIS has been used to obtain first results of {\it 3D} simulations of relativistic jets. We have studied the influence of a slight perturbation of the injection velocity field on the morphodynamics of otherwise axisymmetric relativistic jets.
△ Less
Submitted 26 October, 1999;
originally announced October 1999.
-
Numerical simulations of relativistic wind accretion on to black holes using Godunov-type methods
Authors:
Jose A. Font,
Jose M. Ibanez,
Philippos Papadopoulos
Abstract:
We have studied numerically the so-called Bondi-Hoyle (wind) accretion on to a rotating (Kerr) black hole in general relativity. We have used the Kerr-Schild form of the Kerr metric, free of coordinate singularities at the black hole horizon. The `test-fluid' approximation has been adopted, assuming no dynamical evolution of the gravitational field. We have used a recent formulation of the gener…
▽ More
We have studied numerically the so-called Bondi-Hoyle (wind) accretion on to a rotating (Kerr) black hole in general relativity. We have used the Kerr-Schild form of the Kerr metric, free of coordinate singularities at the black hole horizon. The `test-fluid' approximation has been adopted, assuming no dynamical evolution of the gravitational field. We have used a recent formulation of the general relativistic hydrodynamic equations which casts them into a first-order hyperbolic system of conservation laws. Our studies have been performed using a Godunov-type scheme based on Marquina's flux-formula.
We find that regardless of the value of the black hole spin the final accretion pattern is always stable, leading to constant accretion rates of mass and momentum. The flow is characterized by a strong tail shock which is increasingly wrapped around the central black hole as the hole angular momentum increases. The rotation induced asymmetry in the pressure field implies that besides the well known drag, the black hole will experience also a lift normal to the flow direction.
△ Less
Submitted 23 October, 1999;
originally announced October 1999.
-
Gravitational waves from galaxy cluster distributions
Authors:
Vicent Quilis,
Jose M. Ibanez,
Diego Saez
Abstract:
Galaxy clusters are sources of gravitational radiation. The main aim of this paper is to give numerical estimates and theoretical description of the relevant features of the gravitational radiation coming from an appropriate spatial distributions of galaxy clusters. Since no analytical approaches are currently available to describe the strongly nonlinear regime, our numerical approach - combinin…
▽ More
Galaxy clusters are sources of gravitational radiation. The main aim of this paper is to give numerical estimates and theoretical description of the relevant features of the gravitational radiation coming from an appropriate spatial distributions of galaxy clusters. Since no analytical approaches are currently available to describe the strongly nonlinear regime, our numerical approach - combining numerical simulations with statistical arguments - seems to be an useful way of studying the main features of that radiation. Although far to be detectable with present technology, we advance some ideas about future observational strategies and its cosmological implications.
△ Less
Submitted 14 October, 1999;
originally announced October 1999.
-
High--Resolution 3D Simulations of Relativistic Jets
Authors:
M. A. Aloy,
J. M. Ibanez,
J. M. Marti J. L. Gomez,
E. Muller
Abstract:
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…
▽ More
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 relativistic beam give rise to a lumpy distribution of apparent speeds like that observed in M87. The beam is surrounded by a boundary layer of high specific internal energy. The properties of this layer are briefly discussed.
△ Less
Submitted 26 June, 1999;
originally announced June 1999.
-
An Efficient Implementation of Flux Formulae in Multidimensional Relativistic Hydrodynamical Codes
Authors:
M. A. Aloy,
J. A. Pons,
J. M. Ibanez
Abstract:
We derive and analyze a simplified formulation of the numerical viscosity terms appearing in the expression of the numerical fluxes associated to several High-Resolution Shock-Capturing schemes. After some algebraic pre-processing, we give explicit expressions for the numerical viscosity terms of two of the most widely used flux formulae, which implementation saves computational time in multidim…
▽ More
We derive and analyze a simplified formulation of the numerical viscosity terms appearing in the expression of the numerical fluxes associated to several High-Resolution Shock-Capturing schemes. After some algebraic pre-processing, we give explicit expressions for the numerical viscosity terms of two of the most widely used flux formulae, which implementation saves computational time in multidimensional simulations of relativistic flows. Additionally, such treatment explicitely cancells and factorizes a number of terms helping to amortiguate the growing of round-off errors. We have checked the performance of our formulation running a 3D relativistic hydrodynamical code to solve a standard test-bed problem and found that the improvement in efficiency is of high practical interest in numerical simulations of relativistic flows in Astrophysics.
△ Less
Submitted 15 April, 1999;
originally announced April 1999.
-
GENESIS: A high-resolution code for 3D relativistic hydrodynamics
Authors:
M. A. Aloy,
J. M. Ibanez,
J. M. Marti,
E. Muller
Abstract:
The main features of a three dimensional, high-resolution special relativistic hydro code based on relativistic Riemann solvers are described. The capabilities and performance of the code are discussed. In particular, we present the results of extensive test calculations which demonstrate that the code can accurately and efficiently handle strong shocks in three spatial dimensions. Results of th…
▽ More
The main features of a three dimensional, high-resolution special relativistic hydro code based on relativistic Riemann solvers are described. The capabilities and performance of the code are discussed. In particular, we present the results of extensive test calculations which demonstrate that the code can accurately and efficiently handle strong shocks in three spatial dimensions. Results of the performance of the code on single and multi-processor machines are given. Simulations (in double precision) with 7x10^6 computational cells require less than 1 Gb of RAM memory and ~7x10^-5 CPU seconds per zone and time step (on a SCI Cray-Origin 2000 with a R10000 processor). Currently, a version of the numerical code is under development, which is suited for massively parallel computers with distributed memory architecture (like, e.g., Cray T3E).
△ Less
Submitted 23 March, 1999;
originally announced March 1999.
-
Non-axisymmetric relativistic Bondi-Hoyle accretion onto a Kerr black hole
Authors:
Jose A. Font,
Jose M. Ibanez,
Philippos Papadopoulos
Abstract:
In our program of studying numerically the so-called Bondi-Hoyle accretion in the fully relativistic regime, we present here first results concerning the evolution of matter accreting supersonically onto a rotating (Kerr) black hole. These computations generalize previous results where the non-rotating (Schwarzschild) case was extensively considered. We parametrize our initial data by the asympt…
▽ More
In our program of studying numerically the so-called Bondi-Hoyle accretion in the fully relativistic regime, we present here first results concerning the evolution of matter accreting supersonically onto a rotating (Kerr) black hole. These computations generalize previous results where the non-rotating (Schwarzschild) case was extensively considered. We parametrize our initial data by the asymptotic conditions for the fluid and explore the dependence of the solution on the angular momentum of the black hole. Towards quantifying the robustness of our numerical results, we use two different geometrical foliations of the black hole spacetime, the standard form of the Kerr metric in Boyer-Lindquist coordinates as well as its Kerr-Schild form, which is free of coordinate singularities at the black hole horizon. We demonstrate some important advantages of using such horizon adapted coordinate systems.
Our numerical study indicates that regardless of the value of the black hole spin the final accretion pattern is always stable, leading to constant accretion rates of mass and momentum. The flow is characterized by a strong tail shock, which, unlike the Schwarzschild case, is increasingly wrapped around the central black hole as the hole angular momentum increases. The rotation induced asymmetry in the pressure field implies that besides the well known drag, the black hole will experience also a lift normal to the flow direction. This situation exhibits some analogies with the Magnus effect of classical fluid dynamics.
△ Less
Submitted 21 October, 1998;
originally announced October 1998.
-
Numerical evolution of matter in dynamical axisymmetric black hole spacetimes. I. Methods and tests
Authors:
S. Brandt,
J. A. Font,
J. M. Ibanez,
J. Masso,
E. Seidel
Abstract:
We have developed a numerical code to study the evolution of self-gravitating matter in dynamic black hole axisymmetric spacetimes in general relativity. The matter fields are evolved with a high-resolution shock-capturing scheme that uses the characteristic information of the general relativistic hydrodynamic equations to build up a linearized Riemann solver. The spacetime is evolved with an ax…
▽ More
We have developed a numerical code to study the evolution of self-gravitating matter in dynamic black hole axisymmetric spacetimes in general relativity. The matter fields are evolved with a high-resolution shock-capturing scheme that uses the characteristic information of the general relativistic hydrodynamic equations to build up a linearized Riemann solver. The spacetime is evolved with an axisymmetric ADM code designed to evolve a wormhole in full general relativity. We discuss the numerical and algorithmic issues related to the effective coupling of the hydrodynamical and spacetime pieces of the code, as well as the numerical methods and gauge conditions we use to evolve such spacetimes. The code has been put through a series of tests that verify that it functions correctly. Particularly, we develop and describe a new set of testbed calculations and techniques designed to handle dynamically sliced, self-gravitating matter flows on black holes, and subject the code to these tests. We make some studies of the spherical and axisymmetric accretion onto a dynamic black hole, the fully dynamical evolution of imploding shells of dust with a black hole, the evolution of matter in rotating spacetimes, the gravitational radiation induced by the presence of the matter fields and the behavior of apparent horizons through the evolution.
△ Less
Submitted 9 July, 1998;
originally announced July 1998.
-
General Relativistic Hydrodynamics with Special Relativistic Riemann Solvers
Authors:
Jose A. Pons,
Jose A. Font,
Jose M. Ibanez,
Jose M. Marti,
Juan A. Miralles
Abstract:
We present a general and practical procedure to solve the general relativistic hydrodynamic equations by 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 terms of which the spacetime metric is locally Minkowskian and permits accurate numerical calculations o…
▽ More
We present a general and practical procedure to solve the general relativistic hydrodynamic equations by 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 terms of which the spacetime metric is locally Minkowskian and permits accurate numerical calculations of general relativistic hydrodynamics problems using the numerical tools developed for the special relativistic case with negligible computational cost. The feasibility of the method has been confirmed by a number of numerical experiments.
△ Less
Submitted 21 July, 1998;
originally announced July 1998.
-
Gravitational waves from galaxy clusters: a new observable effect
Authors:
Vicent Quilis,
Jose Ma. Ibanez,
Diego Saez
Abstract:
A rich galaxy cluster showing strong resemblance with the observed ones is simulated. Cold dark matter spectrum, Gaussian statistics, flat universe, and two components -- baryonic gas plus dark matter particles -- are considered. We have calculated the gravitational-wave output during the epoch of the fully nonlinear and nonsymmetric cluster evolution. The amplitudes and frequencies of the resul…
▽ More
A rich galaxy cluster showing strong resemblance with the observed ones is simulated. Cold dark matter spectrum, Gaussian statistics, flat universe, and two components -- baryonic gas plus dark matter particles -- are considered. We have calculated the gravitational-wave output during the epoch of the fully nonlinear and nonsymmetric cluster evolution. The amplitudes and frequencies of the resulting gravitational waves are estimated. Since frequencies are very small --of the order of $10^{-17} Hz$ -- a complete pulse cannot be observed during an admissible integration time; nevertheless, it is proved that these waves can produce an interesting secular effect which appears to be observable with current technology.
△ Less
Submitted 22 May, 1998;
originally announced May 1998.
-
A "horizon adapted" approach to the study of relativistic accretion flows onto rotating black holes
Authors:
Jose A. Font,
Jose M. Ibanez,
Philippos Papadopoulos
Abstract:
We present a new geometrical approach to the study of accretion flows onto rotating (Kerr) black holes. Instead of Boyer-Lindquist coordinates, the standard choice in all existing numerical simulations in the literature, we employ the simplest example of a horizon adapted coordinate system, the Kerr-Schild coordinates. This choice eliminates boundary ambiguities and unphysical divergent behavior…
▽ More
We present a new geometrical approach to the study of accretion flows onto rotating (Kerr) black holes. Instead of Boyer-Lindquist coordinates, the standard choice in all existing numerical simulations in the literature, we employ the simplest example of a horizon adapted coordinate system, the Kerr-Schild coordinates. This choice eliminates boundary ambiguities and unphysical divergent behavior at the event horizon. Computations of Bondi-Hoyle accretion onto extreme Kerr black holes, performed here for the first time, demonstrate the key advantages of this procedure. We argue it offers the best approach to the numerical study of the, observationally, increasingly more accesible relativistic inner region around black holes.
△ Less
Submitted 21 August, 1998; v1 submitted 20 May, 1998;
originally announced May 1998.
-
Non-axisymmetric relativistic Bondi-Hoyle accretion onto a Schwarzschild black hole
Authors:
Jose A. Font,
Jose M. Ibanez
Abstract:
We present the results of an exhaustive numerical study of fully relativistic non-axisymmetric Bondi-Hoyle accretion onto a moving Schwarzschild black hole. We have solved the equations of general relativistic hydrodynamics with a high-resolution shock-capturing numerical scheme based on a linearized Riemann solver. The numerical code was previously used to study axisymmetric flow configurations…
▽ More
We present the results of an exhaustive numerical study of fully relativistic non-axisymmetric Bondi-Hoyle accretion onto a moving Schwarzschild black hole. We have solved the equations of general relativistic hydrodynamics with a high-resolution shock-capturing numerical scheme based on a linearized Riemann solver. The numerical code was previously used to study axisymmetric flow configurations past a Schwarzschild hole. We have analyzed and discussed the flow morphology for a sample of asymptotically high Mach number models. The results of this work reveal that initially asymptotic uniform flows always accrete onto the hole in a stationary way which closely resembles the previous axisymmetric patterns. This is in contrast with some Newtonian numerical studies where violent flip-flop instabilities were found. As discussed in the text, the reason can be found in the initial conditions used in the relativistic regime, as they can not exactly duplicate the previous Newtonian setups where the instability appeared. The dependence of the final solution with the inner boundary condition as well as with the grid resolution has also been studied. Finally, we have computed the accretion rates of mass and linear and angular momentum.
△ Less
Submitted 24 April, 1998;
originally announced April 1998.
-
On the role of shock waves in galaxy cluster evolution
Authors:
Vicent Quilis,
Jose Ma. Ibanez,
Diego Saez
Abstract:
Numerical simulations of galaxy clusters including two species -- baryonic gas and dark matter particles --are presented. Cold Dark Matter spectrum, Gaussian statistics and flat universe are assumed. The dark matter component is evolved numerically by means of a standard particle mesh method. The evolution of the baryonic component has been studied numerically by using a multidimensional (3D) hy…
▽ More
Numerical simulations of galaxy clusters including two species -- baryonic gas and dark matter particles --are presented. Cold Dark Matter spectrum, Gaussian statistics and flat universe are assumed. The dark matter component is evolved numerically by means of a standard particle mesh method. The evolution of the baryonic component has been studied numerically by using a multidimensional (3D) hydrodynamical code based on modern high resolution shock capturing techniques. These techniques are specially designed for treating accurately complex flows in which shocks appear and interact. With this picture, the role of shock waves in the formation and evolution of rich galaxy clusters is analyzed. Our results display two well differenced morphologies of the shocked baryonic matter: filamentary at early epochs and quasi-spherical at low redshifts.
△ Less
Submitted 6 March, 1998;
originally announced March 1998.
-
Legendre expansion of the neutrino-antineutrino annihilation kernel: Influence of high order terms
Authors:
J. A. Pons,
J. A. Miralles,
J. M. Ibanez
Abstract:
We calculate the Legendre expansion of the rate of the process $ν+ \barν \leftrightarrow e^+ + e^-$ up to 3rd order extending previous results of other authors which only consider the 0th and 1st order terms. Using different closure relations for the moment equations of the radiative transfer equation we discuss the physical implications of taking into account quadratic and cubic terms on the en…
▽ More
We calculate the Legendre expansion of the rate of the process $ν+ \barν \leftrightarrow e^+ + e^-$ up to 3rd order extending previous results of other authors which only consider the 0th and 1st order terms. Using different closure relations for the moment equations of the radiative transfer equation we discuss the physical implications of taking into account quadratic and cubic terms on the energy deposition outside the neutrinosphere in a simplified model. The main conclusion is that 2nd order is necessary in the semi-transparent region and gives good results if an appropriate closure relation is used.
△ Less
Submitted 26 February, 1998;
originally announced February 1998.
-
A multidimensional hydrodynamic code for structure evolution in cosmology
Authors:
Vicent Quilis,
Jose Ma. Ibanez,
Diego Saez
Abstract:
A cosmological multidimensional hydrodynamic code is described and tested. This code is based on modern high-resolution shock-capturing techniques. It can make use of a linear or a parabolic cell reconstruction as well as an approximate Riemann solver. The code has been specifically designed for cosmological applications. Two tests including shocks have been considered: the first one is a standa…
▽ More
A cosmological multidimensional hydrodynamic code is described and tested. This code is based on modern high-resolution shock-capturing techniques. It can make use of a linear or a parabolic cell reconstruction as well as an approximate Riemann solver. The code has been specifically designed for cosmological applications. Two tests including shocks have been considered: the first one is a standard shock tube and the second test involves a spherically symmetric shock. Various additional cosmological tests are also presented. In this way, the performance of the code is proved. The usefulness of the code is discussed; in particular, this powerful tool is expected to be useful in order to study the evolution of the hot gas component located inside nonsymmetric cosmological structures.
△ Less
Submitted 9 April, 1996;
originally announced April 1996.
-
A new spherically symmetric general relativistic hydrodynamical code
Authors:
Jose V. Romero,
Jose M. Ibanez,
Jose M. Marti,
Juan A. Miralles
Abstract:
In this paper we present a full general relativistic one-dimensional hydro-code which incorporates a modern high-resolution shock-capturing algorithm, with an approximate Riemann solver, for the correct modelling of formation and propagation of strong shocks. The efficiency of this code in treating strong shocks is demonstrated by some numerical experiments. The interest of this technique in sev…
▽ More
In this paper we present a full general relativistic one-dimensional hydro-code which incorporates a modern high-resolution shock-capturing algorithm, with an approximate Riemann solver, for the correct modelling of formation and propagation of strong shocks. The efficiency of this code in treating strong shocks is demonstrated by some numerical experiments. The interest of this technique in several astrophysical scenarios is discussed.
△ Less
Submitted 26 September, 1995; v1 submitted 25 September, 1995;
originally announced September 1995.
-
Neutrino pair synchrotron radiation from relativistic electrons in strong magnetic fields
Authors:
A. Vidaurre,
A. Perez,
H. Sivak,
J. Bernabeu,
J. M. Ibanez
Abstract:
The emissivity for the neutrino pair synchrotron radiation in strong magnetic fields has been calculated both analytically and numerically for high densities and moderate temperatures, as can be found in neutron stars. Under these conditions, the electrons are relativistic and degenerate. We give here our results in terms of an universal function of a single variable. For two different regimes o…
▽ More
The emissivity for the neutrino pair synchrotron radiation in strong magnetic fields has been calculated both analytically and numerically for high densities and moderate temperatures, as can be found in neutron stars. Under these conditions, the electrons are relativistic and degenerate. We give here our results in terms of an universal function of a single variable. For two different regimes of the electron gas we present a simplified calculation and compare our results to those of Kaminker et al. Agreement is found for the classical region, where many Landau levels contribute to the emissivity , but some differences arise in the quantum regime. One finds that the emissivity for neutrino pair synchrotron radiation is competitive, and can be dominant, with other neutrino processes for magnetic fields of the order $B \sim 10^{14} - 10^{15} G $.This indicates the relevance of this process for some astrophysical scenarios, such as neutron stars and supernovae.
△ Less
Submitted 10 July, 1995; v1 submitted 7 July, 1995;
originally announced July 1995.
-
Galaxy clusters and microwave background anisotropy
Authors:
V. Quilis,
J. M. Ibanez,
D. Saez
Abstract:
Previous estimates of the microwave background anisotropies produced by freely falling spherical clusters are discussed. These estimates are based on the Swiss-Cheese and Tolman-Bondi models. It is proved that these models give only upper limits to the anisotropies produced by the observed galaxy clusters. By using spherically symmetric codes including pressureless matter and a hot baryonic gas,…
▽ More
Previous estimates of the microwave background anisotropies produced by freely falling spherical clusters are discussed. These estimates are based on the Swiss-Cheese and Tolman-Bondi models. It is proved that these models give only upper limits to the anisotropies produced by the observed galaxy clusters. By using spherically symmetric codes including pressureless matter and a hot baryonic gas, new upper limits are obtained. The contributions of the hot gas and the pressureless component to the total anisotropy are compared. The effects produced by the pressure are proved to be negligible; hence, estimations of the cluster anisotropies based on N-body simulations are hereafter justified. After the phenomenon of violent relaxation, any realistic rich cluster can only produce small anisotropies with amplitudes of order $10^{-7}$. During the rapid process of violent relaxation, the anisotropies produced by nonlinear clusters are expected to range in the interval $(10^{-6},10^{-5})$. The angular scales of these anisotropies are discussed.
△ Less
Submitted 8 July, 1995; v1 submitted 6 July, 1995;
originally announced July 1995.