-
Efficient WENO schemes for nonuniform grids
Authors:
M. C. Martí,
P. Mulet,
D. F. Yáñez,
D. Zorío
Abstract:
A set of arbitrarily high-order WENO schemes for reconstructions on nonuniform grids is presented. These non-linear interpolation methods use simple smoothness indicators with a linear cost with respect to the order, making them easy to implement and computationally efficient. The theoretical analysis to verify the accuracy and the essentially non-oscillatory properties are presented together with…
▽ More
A set of arbitrarily high-order WENO schemes for reconstructions on nonuniform grids is presented. These non-linear interpolation methods use simple smoothness indicators with a linear cost with respect to the order, making them easy to implement and computationally efficient. The theoretical analysis to verify the accuracy and the essentially non-oscillatory properties are presented together with some numerical experiments involving algebraic problems in order to validate them. Also, these general schemes are applied for the solution of conservation laws and hyperbolic systems in the context of finite volume methods.
△ Less
Submitted 15 May, 2024;
originally announced May 2024.
-
WENO scheme on characteristics for the equilibrium dispersive model of chromatography with generalized Langmuir isotherms
Authors:
R. Donat,
M. C. Martí,
P. Mulet
Abstract:
Column chromatography is a laboratory and industrial technique used to separate different substances mixed in a solution. Mathematically, it can be modelled using non-linear partial differential equations whose main ingredients are the adsorption isotherms, which are non-linear functions modelling the affinity between the different substances in the solution and the solid stationary phase filling…
▽ More
Column chromatography is a laboratory and industrial technique used to separate different substances mixed in a solution. Mathematically, it can be modelled using non-linear partial differential equations whose main ingredients are the adsorption isotherms, which are non-linear functions modelling the affinity between the different substances in the solution and the solid stationary phase filling the column. The goal of this work is twofold. Firstly, we aim to extend the techniques of Donat, Guerrero and Mulet (Appl. Numer. Math. 123 (2018) 22-42) to other adsorption isotherms. In particular, we propose a family of generalized Langmuir-type isotherms and prove that the correspondence between the concentrations of solutes in the liquid phase (the primitive variables) and the conserved variables is well defined and admits a global smooth inverse that can be computed numerically. Secondly, to establish the well-posedness of the mathematical model, we study the eigenstructure of the Jacobian of the mentioned correspondence and use this characteristic information to get oscillation-free sharp interfaces on the numerical approximate solutions. To do so, we determine the structure of the Jacobian matrix of the system and use it to deduce its eigenstructure. We combine the use of characteristic-based numerical fluxes with a second-order implicit-explicit scheme proposed in the cited reference and perform some numerical experiments with Tóth's adsorption isotherms to demonstrate that the characteristic-based schemes produce accurate numerical solutions with no oscillations, even when steep gradients appear in the solutions.
△ Less
Submitted 15 May, 2024;
originally announced May 2024.
-
On approximate implicit Taylor methods for ordinary differential equations
Authors:
Antonio Baeza,
Raimund Bürger,
María del Carmen Martí,
Pep Mulet,
David Zorío
Abstract:
An efficient approximate version of implicit Taylor methods for initial-value problems of systems of ordinary differential equations (ODEs) is introduced. The approach, based on an approximate formulation of Taylor methods, produces a method that requires less evaluations of the function that defines the ODE and its derivatives than the usual version. On the other hand, an efficient numerical solu…
▽ More
An efficient approximate version of implicit Taylor methods for initial-value problems of systems of ordinary differential equations (ODEs) is introduced. The approach, based on an approximate formulation of Taylor methods, produces a method that requires less evaluations of the function that defines the ODE and its derivatives than the usual version. On the other hand, an efficient numerical solution of the equation that arises from the discretization by means of Newton's method is introduced for an implicit scheme of any order. Numerical experiments illustrate that the resulting algorithm is simpler to implement and has better performance than its exact counterpart.
△ Less
Submitted 2 February, 2024;
originally announced February 2024.
-
Magnetic and thermal acceleration in extragalactic jets: An application to NGC 315
Authors:
L. Ricci,
M. Perucho,
J. López-Miralles,
J. M. Martí,
B. Boccardi
Abstract:
Aims. Relativistic jets launched from active galactic nuclei accelerate up to highly relativistic velocities within a few parsecs to tens of parsecs. The precise way in which this process takes place is still under study. While magnetic acceleration is known to be able to accelerate relativistic outflows, little attention has been paid to the role of thermal acceleration. The latter has been assum…
▽ More
Aims. Relativistic jets launched from active galactic nuclei accelerate up to highly relativistic velocities within a few parsecs to tens of parsecs. The precise way in which this process takes place is still under study. While magnetic acceleration is known to be able to accelerate relativistic outflows, little attention has been paid to the role of thermal acceleration. The latter has been assumed to act only on compact regions, very close to the central engine, and to become negligible on parsec scales. However, this holds under the assumption of small internal energies as compared to the magnetic ones, and whether this is true or what happens when we drop this assumption is currently uncertain. Methods. We use a 2D relativistic magnetohydrodynamical code to explore jet acceleration from sub-parsec to parsec scales. As initial conditions for our models, we use observational constraints on jet properties derived by means of very long baseline interferometry observations for a Fanaroff Riley I radio galaxy, NGC\,315. We investigate the parameter space established for this source and perform a number of simulations of magnetically, thermally or kinetically dominated jets at injection, and compare our results with the observed ones. Results. Our simulated jets show that when thermal energy is comparable to or exceeds magnetic energy, thermal acceleration becomes significant at parsec scales. This result has important consequences, potentially extending the acceleration region far beyond the collimation scales, as thermal acceleration can effectively operate within a conically expanding jet. In all the models, we observe acceleration to be driven by expansion, as expected. A number of our models allow us to reproduce the acceleration and opening angles observed in NGC\,315. Finally, our results indicate that disk-launched winds might play an important role in the jet propagation.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
Filamentary structures as the origin of blazar jet radio variability
Authors:
Antonio Fuentes,
José L. Gómez,
José M. Martí,
Manel Perucho,
Guang-Yao Zhao,
Rocco Lico,
Andrei P. Lobanov,
Gabriele Bruni,
Yuri Y. Kovalev,
Andrew Chael,
Kazunori Akiyama,
Katherine L. Bouman,
He Sun,
Ilje Cho,
Efthalia Traianou,
Teresa Toscano,
Rohan Dahale,
Marianna Foschi,
Leonid I. Gurvits,
Svetlana Jorstad,
Jae-Young Kim,
Alan P. Marscher,
Yosuke Mizuno,
Eduardo Ros,
Tuomas Savolainen
Abstract:
Supermassive black holes at the centre of active galactic nuclei power some of the most luminous objects in the Universe. Typically, very long baseline interferometric (VLBI) observations of blazars have revealed only funnel-like morphologies with little information of the ejected plasma internal structure, or lacked the sufficient dynamic range to reconstruct the extended jet emission. Here we sh…
▽ More
Supermassive black holes at the centre of active galactic nuclei power some of the most luminous objects in the Universe. Typically, very long baseline interferometric (VLBI) observations of blazars have revealed only funnel-like morphologies with little information of the ejected plasma internal structure, or lacked the sufficient dynamic range to reconstruct the extended jet emission. Here we show microarcsecond-scale angular resolution images of the blazar 3C 279 obtained at 22 GHz with the space VLBI mission RadioAstron, which allowed us to resolve the jet transversely and reveal several filaments produced by plasma instabilities in a kinetically dominated flow. Our high angular resolution and dynamic range image suggests that emission features traveling down the jet may manifest as a result of differential Doppler-boosting within the filaments, as opposed to the standard shock-in-jet model invoked to explain blazar jet radio variability. Moreover, we infer that the filaments in 3C 279 are possibly threaded by a helical magnetic field rotating clockwise, as seen in the direction of the flow motion, with an intrinsic helix pitch angle of ~45 degrees in a jet with a Lorentz factor of ~13 at the time of observation.
△ Less
Submitted 3 November, 2023;
originally announced November 2023.
-
Polarised light modular microscopy (Pol-ModMicro) for identifying hemozoin crystals in Plasmodium
Authors:
Fraser Eadie,
Matthew P Gibbins,
Graham Gibson,
Matthias Marti,
Akhil Kallepalli
Abstract:
Plasmodium spp. are the protozoan parasites responsible for malaria. Plasmodium spp. synthesise a biocrystal, hemozoin, which can be observed under cross-polarised light. These birefringent crystals can be seen due to different refractive indices of the hemozoin crystal and red blood cells. Here, we present a polarised light modular microscopy (Pol-ModMicro) solution, complete with illumination so…
▽ More
Plasmodium spp. are the protozoan parasites responsible for malaria. Plasmodium spp. synthesise a biocrystal, hemozoin, which can be observed under cross-polarised light. These birefringent crystals can be seen due to different refractive indices of the hemozoin crystal and red blood cells. Here, we present a polarised light modular microscopy (Pol-ModMicro) solution, complete with illumination sources and a robust imaging system to capture birefringence and identify the parasites. We achieve this by combining ModLight light sources with bespoke components designed for the OpenFlexure microscope to image blood smears. Further, a simple and robust algorithm capable of enhancing birefringence is presented. This solution provides image quality that is comparable to a substantially more expensive proprietary microscope.
△ Less
Submitted 26 April, 2023;
originally announced April 2023.
-
On the application of Jacobian-free Riemann solvers for relativistic radiation magnetohydrodynamics under M1 closure
Authors:
Jose López-Miralles,
Jose María Martí,
Manel Perucho
Abstract:
Radiative transfer plays a major role in high-energy astrophysics. In multiple scenarios and in a broad range of energy scales, the coupling between matter and radiation is essential to understand the interplay between theory, observations and numerical simulations. In this paper, we present a novel scheme for solving the equations of radiation relativistic magnetohydrodynamics within the parallel…
▽ More
Radiative transfer plays a major role in high-energy astrophysics. In multiple scenarios and in a broad range of energy scales, the coupling between matter and radiation is essential to understand the interplay between theory, observations and numerical simulations. In this paper, we present a novel scheme for solving the equations of radiation relativistic magnetohydrodynamics within the parallel code Lóstrego. These equations, which are formulated taking successive moments of the Boltzmann radiative transfer equation, are solved under the gray-body approximation and the M1 closure using an IMEX time integration scheme. The main novelty of our scheme is that we introduce for the first time in the context of radiation magnetohydrodynamics a family of Jacobian-free Riemann solvers based on internal approximations to the Polynomial Viscosity Matrix, which were demonstrated to be robust and accurate for non-radiative applications. The robustness and the limitations of the new algorithms are tested by solving a collection of one-dimensional and multi-dimensional test problems, both in the free-streaming and in the diffusion radiation transport limits. Due to its stable performance, the applicability of the scheme presented in this paper to real astrophysical scenarios in high-energy astrophysics is promising. In future simulations, we expect to be able to explore the dynamical relevance of photon-matter interactions in the context of relativistic jets and accretion discs, from microquasars and AGN to gamma-ray bursts.
△ Less
Submitted 1 December, 2022;
originally announced December 2022.
-
Perspectives for self-driving labs in synthetic biology
Authors:
Hector Garcia Martin,
Tijana Radivojevic,
Jeremy Zucker,
Kristofer Bouchard,
Jess Sustarich,
Sean Peisert,
Dan Arnold,
Nathan Hillson,
Gyorgy Babnigg,
Jose Manuel Marti,
Christopher J. Mungall,
Gregg T. Beckham,
Lucas Waldburger,
James Carothers,
ShivShankar Sundaram,
Deb Agarwal,
Blake A. Simmons,
Tyler Backman,
Deepanwita Banerjee,
Deepti Tanjore,
Lavanya Ramakrishnan,
Anup Singh
Abstract:
Self-driving labs (SDLs) combine fully automated experiments with artificial intelligence (AI) that decides the next set of experiments. Taken to their ultimate expression, SDLs could usher a new paradigm of scientific research, where the world is probed, interpreted, and explained by machines for human benefit. While there are functioning SDLs in the fields of chemistry and materials science, we…
▽ More
Self-driving labs (SDLs) combine fully automated experiments with artificial intelligence (AI) that decides the next set of experiments. Taken to their ultimate expression, SDLs could usher a new paradigm of scientific research, where the world is probed, interpreted, and explained by machines for human benefit. While there are functioning SDLs in the fields of chemistry and materials science, we contend that synthetic biology provides a unique opportunity since the genome provides a single target for affecting the incredibly wide repertoire of biological cell behavior. However, the level of investment required for the creation of biological SDLs is only warranted if directed towards solving difficult and enabling biological questions. Here, we discuss challenges and opportunities in creating SDLs for synthetic biology.
△ Less
Submitted 1 November, 2022; v1 submitted 14 October, 2022;
originally announced October 2022.
-
3D RMHD simulations of jet-wind interactions in High Mass X-ray Binaries
Authors:
J. López-Miralles,
M. Perucho,
J. M. Martí,
S. Migliari,
V. Bosch-Ramon
Abstract:
The interaction of jets in High-Mass X-ray Binaries (HMXBs) with the strong winds driven by the hot companion star in the vicinity of the compact object is fundamental to understand the jet dynamics, non-thermal emission and long-term stability. However, the role of the jet magnetic field in this process is unclear. We study the dynamical role of weak and moderate-to-strong toroidal magnetic field…
▽ More
The interaction of jets in High-Mass X-ray Binaries (HMXBs) with the strong winds driven by the hot companion star in the vicinity of the compact object is fundamental to understand the jet dynamics, non-thermal emission and long-term stability. However, the role of the jet magnetic field in this process is unclear. We study the dynamical role of weak and moderate-to-strong toroidal magnetic fields during the first hundreds of seconds of jet propagation, focusing on the magnetized flow dynamics and the mechanisms of energy conversion. We have developed the code Lóstrego v1.0, a new 3D RMHD code to simulate astrophysical plasmas in Cartesian coordinates. Using this tool, we performed the first 3D RMHD numerical simulations of relativistic magnetized jets propagating through the clumpy stellar wind in a HMXB. The overall morphology and dynamics of weakly magnetized jet models is similar to previous hydrodynamical simulations, where the jet head generates a strong shock in the ambient medium and the initial over-pressure with respect to the stellar wind drives one or more recollimation shocks. In the time scales of our simulations, these jets are ballistic and seem to be more stable against internal instabilities than jets with the same power in the absence of fields. However, moderate-to-strong toroidal magnetic fields favour the development of current-driven instabilities and the disruption of the jet within the binary. A detailed analysis of the energy distribution in the relativistic outflow and the ambient medium reveals that both magnetic and internal energies can contribute to the effective acceleration of the jet. We certify that the jet feedback into the ambient medium is highly dependent on the jet energy distribution at injection, where hotter, more dilute and/or more magnetized jets are more efficient, as anticipated by feedback studies in the case of jets in active galaxies.
△ Less
Submitted 22 February, 2022;
originally announced February 2022.
-
A degenerating convection-diffusion system modelling froth flotation with drainage
Authors:
Raimund Bürger,
Stefan Diehl,
M. Carmen Martí,
Yolanda Vásquez
Abstract:
Froth flotation is a common unit operation used in mineral processing. It serves to separate valuable mineral particles from worthless gangue particles in finely ground ores. The valuable mineral particles are hydrophobic and attach to bubbles of air injected into the pulp. This creates bubble-particle aggregates that rise to the top of the flotation column where they accumulate to a froth or foam…
▽ More
Froth flotation is a common unit operation used in mineral processing. It serves to separate valuable mineral particles from worthless gangue particles in finely ground ores. The valuable mineral particles are hydrophobic and attach to bubbles of air injected into the pulp. This creates bubble-particle aggregates that rise to the top of the flotation column where they accumulate to a froth or foam layer that is removed through a launder for further processing. At the same time, the hydrophilic gangue particles settle and are removed continuously. The drainage of liquid due to capillarity is essential for the formation of a stable froth layer. This effect is included into a previously formulated hyperbolic system of partial differential equations that models the volume fractions of floating aggregates and settling hydrophilic solids [R. Bürger, S. Diehl and M.C. Martí, {\it IMA J. Appl. Math.} {\bf 84} (2019) 930--973]. The construction of desired steady-state solutions with a froth layer is detailed and feasibility conditions on the feed volume fractions and the volumetric flows of feed, underflow and wash water are visualized in so-called operating charts. A monotone numerical scheme is derived and employed to simulate the dynamic behaviour of a flotation column. It is also proven that, under a suitable Courant-Friedrichs-Lewy (CFL) condition, the approximate volume fractions are bounded between zero and one when the initial data are.
△ Less
Submitted 9 February, 2022;
originally announced February 2022.
-
Reading M87's DNA: A Double Helix revealing a large scale Helical Magnetic Field
Authors:
Alice Pasetto,
Carlos Carrasco-Gonzalez,
Jose L. Gomez,
Jose M. Marti,
Manel Perucho,
Shane P. O'Sullivan,
Craig Anderson,
Daniel Jacobo Diaz-Gonzalez,
Antonio Fuentes,
John Wardle
Abstract:
We present unprecedented high fidelity radio images of the M87 jet. We analyzed Jansky Very Large Array (VLA) broadband, full polarization, radio data from 4 to 18 GHz. The observations were taken with the most extended configuration (A configuration), which allow the study of the emission of the jet up to kpc scales with a linear resolution $\sim$10 pc. The high sensitivity and resolution of our…
▽ More
We present unprecedented high fidelity radio images of the M87 jet. We analyzed Jansky Very Large Array (VLA) broadband, full polarization, radio data from 4 to 18 GHz. The observations were taken with the most extended configuration (A configuration), which allow the study of the emission of the jet up to kpc scales with a linear resolution $\sim$10 pc. The high sensitivity and resolution of our data allow to resolve the jet width. We confirm a double-helix morphology of the jet material between $\sim$300 pc and $\sim$1 kpc. We found a gradient of the polarization degree with a minimum at the projected axis and maxima at the jet edges, and a gradient in the Faraday depth with opposite signs at the jet edges. We also found that the behavior of the polarization properties along the wide range of frequencies is consistent with internal Faraday depolarization. All these characteristics strongly support the presence of a helical magnetic field in the M87 jet up to 1 kpc from the central black hole although the jet is most likely particle dominated at these large scales. Therefore, we propose a plausible scenario in which the helical configuration of the magnetic field has been maintained to large scales thanks to the presence of Kelvin-Helmholtz instabilities.
△ Less
Submitted 13 December, 2021;
originally announced December 2021.
-
Long-term FRII jet evolution in dense environments
Authors:
Manel Perucho,
José María Martí,
Vicent Quilis
Abstract:
We present long-term numerical three-dimensional simulations of a relativistic outflow propagating through a galactic ambient medium and environment, up to distances $\sim 100$~kpc. Our aim is to study the role of dense media in the global dynamics of the radio source. We use a relativistic gas equation of state, and a basic description of thermal cooling terms. In previous work, we showed that a…
▽ More
We present long-term numerical three-dimensional simulations of a relativistic outflow propagating through a galactic ambient medium and environment, up to distances $\sim 100$~kpc. Our aim is to study the role of dense media in the global dynamics of the radio source. We use a relativistic gas equation of state, and a basic description of thermal cooling terms. In previous work, we showed that a linear perturbation could enhance the jet propagation during the early phases of evolution, by introducing obliquity to the jet reverse shock. Here, we show that this effect is reduced in denser media. We find that the dentist-drill effect acts earlier, due to slower jet propagation and an increased growth of the helical instability. The global morphology of the jet is less elongated, with more prominent lobes. The fundamental physical parameters of the jet generated structure derived from our simulations fall within the estimated values derived for FRII jets in the 3C sample. In agreement with previous axisymmetric and three dimensional simulations in lower density media, we conclude that shock heating of the interstellar and intergalactic media is very efficient in the case of powerful, relativistic jets.
△ Less
Submitted 7 December, 2021; v1 submitted 6 December, 2021;
originally announced December 2021.
-
On the deceleration of Fanaroff-Riley Class I jets: mass loading of magnetized jets by stellar winds
Authors:
Andreu Anglés-Castillo,
Manel Perucho,
José María Martí,
Robert A. Laing
Abstract:
In this paper we present steady-state RMHD simulations that include a mass-load term to study the process of jet deceleration. The mass-load mimics the injection of a proton-electron plasma from stellar winds within the host galaxy into initially pair plasma jets, with mean stellar mass-losses ranging from $10^{-14}$ to $10^{-9}\,{M_\odot\,yr^{-1}}$. The spatial jet evolution covers…
▽ More
In this paper we present steady-state RMHD simulations that include a mass-load term to study the process of jet deceleration. The mass-load mimics the injection of a proton-electron plasma from stellar winds within the host galaxy into initially pair plasma jets, with mean stellar mass-losses ranging from $10^{-14}$ to $10^{-9}\,{M_\odot\,yr^{-1}}$. The spatial jet evolution covers $\sim 500\,{\rm pc}$ from jet injection in the grid at 10~pc from the jet nozzle. Our simulations use a relativistic gas equation of state and a pressure profile for the ambient medium. We compare these simulations with previous dynamical simulations of relativistic, non-magnetised jets. Our results show that toroidal magnetic fields can prevent fast jet expansion and the subsequent embedding of further stars via magnetic tension. In this sense, magnetic fields avoid a runaway deceleration process. Furthermore, when the mass-load is large enough to increase the jet density and produce fast, differential jet expansion, the conversion of magnetic energy flux into kinetic energy flux (i.e., magnetic acceleration), helps to delay the deceleration process with respect to non-magnetised jets. We conclude that the typical stellar population in elliptical galaxies cannot explain jet deceleration in classical FRI radio galaxies. However, we observe a significant change in the jet composition, thermodynamical parameters and energy dissipation along its evolution, even for moderate values of the mass-load.
△ Less
Submitted 20 October, 2020;
originally announced October 2020.
-
Cosmic magnetic fields with MASCLET: an application to galaxy clusters
Authors:
Vicent Quilis,
José Ma. Martí,
Susana Planelles
Abstract:
We describe and test a new version of the adaptive mesh refinement (AMR) cosmological code MASCLET. The new version of the code includes all the ingredients of its previous version plus a description of the evolution of the magnetic field under the approximation of the ideal magneto-hydrodynamics (MHD). To preserve the divergence-free condition of MHD, the original divergence cleaning algorithm of…
▽ More
We describe and test a new version of the adaptive mesh refinement (AMR) cosmological code MASCLET. The new version of the code includes all the ingredients of its previous version plus a description of the evolution of the magnetic field under the approximation of the ideal magneto-hydrodynamics (MHD). To preserve the divergence-free condition of MHD, the original divergence cleaning algorithm of Dedner et al. (2002) is implemented. We present a set of well-known 1D and 2D tests, such as several shock-tube problems, the fast rotor and the Orszag-Tang vortex. The performance of the code in all the tests is excellent with estimated median relative errors of $\nabla \cdot {\bf B}$ in the 2D tests smaller than $5 \times 10^{-5}$ for the fast rotor test, and $5 \times 10^{-3}$ for the Orszag-Tang vortex. As an astrophysical application of the code, we present a simulation of a cosmological box of 40 comoving Mpc side length in which a primordial uniform comoving magnetic field of strength 0.1 nG is seeded. The simulation shows how the magnetic field is channelled along the filaments of gas and is concentrated and amplified within galaxy clusters. Comparison with the values expected from pure compression reveals an additional amplification of the magnetic field caused by turbulence in the central region of the cluster. Values of the order of $\sim 1 μ$G are obtained in clusters at $z\sim 0$ with median relative errors of $\nabla \cdot {\bf B}$ below 0.4\%. The implications of a proper description of the dynamics of the magnetic field and their possible observational counterparts in future facilities are discussed.
△ Less
Submitted 27 March, 2020;
originally announced March 2020.
-
Real-Time Semantic Stereo Matching
Authors:
Pier Luigi Dovesi,
Matteo Poggi,
Lorenzo Andraghetti,
Miquel Martí,
Hedvig Kjellström,
Alessandro Pieropan,
Stefano Mattoccia
Abstract:
Scene understanding is paramount in robotics, self-navigation, augmented reality, and many other fields. To fully accomplish this task, an autonomous agent has to infer the 3D structure of the sensed scene (to know where it looks at) and its content (to know what it sees). To tackle the two tasks, deep neural networks trained to infer semantic segmentation and depth from stereo images are often th…
▽ More
Scene understanding is paramount in robotics, self-navigation, augmented reality, and many other fields. To fully accomplish this task, an autonomous agent has to infer the 3D structure of the sensed scene (to know where it looks at) and its content (to know what it sees). To tackle the two tasks, deep neural networks trained to infer semantic segmentation and depth from stereo images are often the preferred choices. Specifically, Semantic Stereo Matching can be tackled by either standalone models trained for the two tasks independently or joint end-to-end architectures. Nonetheless, as proposed so far, both solutions are inefficient because requiring two forward passes in the former case or due to the complexity of a single network in the latter, although jointly tackling both tasks is usually beneficial in terms of accuracy. In this paper, we propose a single compact and lightweight architecture for real-time semantic stereo matching. Our framework relies on coarse-to-fine estimations in a multi-stage fashion, allowing: i) very fast inference even on embedded devices, with marginal drops in accuracy, compared to state-of-the-art networks, ii) trade accuracy for speed, according to the specific application requirements. Experimental results on high-end GPUs as well as on an embedded Jetson TX2 confirm the superiority of semantic stereo matching compared to standalone tasks and highlight the versatility of our framework on any hardware and for any application.
△ Less
Submitted 24 February, 2020; v1 submitted 1 October, 2019;
originally announced October 2019.
-
GANtruth - an unpaired image-to-image translation method for driving scenarios
Authors:
Sebastian Bujwid,
Miquel Martí,
Hossein Azizpour,
Alessandro Pieropan
Abstract:
Synthetic image translation has significant potentials in autonomous transportation systems. That is due to the expense of data collection and annotation as well as the unmanageable diversity of real-words situations. The main issue with unpaired image-to-image translation is the ill-posed nature of the problem. In this work, we propose a novel method for constraining the output space of unpaired…
▽ More
Synthetic image translation has significant potentials in autonomous transportation systems. That is due to the expense of data collection and annotation as well as the unmanageable diversity of real-words situations. The main issue with unpaired image-to-image translation is the ill-posed nature of the problem. In this work, we propose a novel method for constraining the output space of unpaired image-to-image translation. We make the assumption that the environment of the source domain is known (e.g. synthetically generated), and we propose to explicitly enforce preservation of the ground-truth labels on the translated images.
We experiment on preserving ground-truth information such as semantic segmentation, disparity, and instance segmentation. We show significant evidence that our method achieves improved performance over the state-of-the-art model of UNIT for translating images from SYNTHIA to Cityscapes. The generated images are perceived as more realistic in human surveys and outperforms UNIT when used in a domain adaptation scenario for semantic segmentation.
△ Less
Submitted 26 November, 2018;
originally announced December 2018.
-
Applying Faster R-CNN for Object Detection on Malaria Images
Authors:
Jane Hung,
Deepali Ravel,
Stefanie C. P. Lopes,
Gabriel Rangel,
Odailton Amaral Nery,
Benoit Malleret,
Francois Nosten,
Marcus V. G. Lacerda,
Marcelo U. Ferreira,
Laurent Rénia,
Manoj T. Duraisingh,
Fabio T. M. Costa,
Matthias Marti,
Anne E. Carpenter
Abstract:
Deep learning based models have had great success in object detection, but the state of the art models have not yet been widely applied to biological image data. We apply for the first time an object detection model previously used on natural images to identify cells and recognize their stages in brightfield microscopy images of malaria-infected blood. Many micro-organisms like malaria parasites a…
▽ More
Deep learning based models have had great success in object detection, but the state of the art models have not yet been widely applied to biological image data. We apply for the first time an object detection model previously used on natural images to identify cells and recognize their stages in brightfield microscopy images of malaria-infected blood. Many micro-organisms like malaria parasites are still studied by expert manual inspection and hand counting. This type of object detection task is challenging due to factors like variations in cell shape, density, and color, and uncertainty of some cell classes. In addition, annotated data useful for training is scarce, and the class distribution is inherently highly imbalanced due to the dominance of uninfected red blood cells. We use Faster Region-based Convolutional Neural Network (Faster R-CNN), one of the top performing object detection models in recent years, pre-trained on ImageNet but fine tuned with our data, and compare it to a baseline, which is based on a traditional approach consisting of cell segmentation, extraction of several single-cell features, and classification using random forests. To conduct our initial study, we collect and label a dataset of 1300 fields of view consisting of around 100,000 individual cells. We demonstrate that Faster R-CNN outperforms our baseline and put the results in context of human performance.
△ Less
Submitted 11 March, 2019; v1 submitted 25 April, 2018;
originally announced April 2018.
-
Total and linearly polarized synchrotron emission from overpressured magnetized relativistic jets
Authors:
Antonio Fuentes,
José L. Gómez,
José M. Martí,
Manel Perucho
Abstract:
We present relativistic magnetohydrodynamic (RMHD) simulations of stationary overpressured magnetized relativistic jets which are characterized by their dominant type of energy, namely internal, kinetic, or magnetic. Each model is threaded by a helical magnetic field with a pitch angle of $45^\circ$ and features a series of recollimation shocks produced by the initial pressure mismatch, whose stre…
▽ More
We present relativistic magnetohydrodynamic (RMHD) simulations of stationary overpressured magnetized relativistic jets which are characterized by their dominant type of energy, namely internal, kinetic, or magnetic. Each model is threaded by a helical magnetic field with a pitch angle of $45^\circ$ and features a series of recollimation shocks produced by the initial pressure mismatch, whose strength and number varies as a function of the dominant type of energy. We perform a study of the polarization signatures from these models by integrating the radiative transfer equations for synchrotron radiation using as inputs the RMHD solutions. These simulations show a top-down emission asymmetry produced by the helical magnetic field and a progressive confinement of the emission into a jet spine as the magnetization increases and the internal energy of the non-thermal population is considered to be a constant fraction of the thermal one. Bright stationary components associated with the recollimation shocks appear presenting a relative intensity modulated by the Doppler boosting ratio between the pre-shock and post-shock states. Small viewing angles show a roughly bimodal distribution in the polarization angle due to the helical structure of the magnetic field, which is also responsible for the highly stratified degree of linear polarization across the jet width. In addition, small variations of the order of $26^\circ$ are observed in the polarization angle of the stationary components, which can be used to identify recollimation shocks in astrophysical jets.
△ Less
Submitted 5 April, 2018;
originally announced April 2018.
-
A multitask deep learning model for real-time deployment in embedded systems
Authors:
Miquel Martí,
Atsuto Maki
Abstract:
We propose an approach to Multitask Learning (MTL) to make deep learning models faster and lighter for applications in which multiple tasks need to be solved simultaneously, which is particularly useful in embedded, real-time systems. We develop a multitask model for both Object Detection and Semantic Segmentation and analyze the challenges that appear during its training. Our multitask network is…
▽ More
We propose an approach to Multitask Learning (MTL) to make deep learning models faster and lighter for applications in which multiple tasks need to be solved simultaneously, which is particularly useful in embedded, real-time systems. We develop a multitask model for both Object Detection and Semantic Segmentation and analyze the challenges that appear during its training. Our multitask network is 1.6x faster, lighter and uses less memory than deploying the single-task models in parallel. We conclude that MTL has the potential to give superior performance in exchange of a more complex training process that introduces challenges not present in single-task models.
△ Less
Submitted 31 October, 2017;
originally announced November 2017.
-
Magnetized, Relativistic Jets
Authors:
Manel Perucho,
José María Martí,
José Luís Gómez,
Antonio Fuentes
Abstract:
Extragalactic relativistic jets are composed by charged particles and magnetic fields, as inferred from the synchrotron emission that we receive from them. The Larmor radii of the particles propagating along the magnetic field are much smaller than the scales of the problem, providing the necessary coherence to the system to treat is as a flow. We can thus study them using relativistic magnetohydr…
▽ More
Extragalactic relativistic jets are composed by charged particles and magnetic fields, as inferred from the synchrotron emission that we receive from them. The Larmor radii of the particles propagating along the magnetic field are much smaller than the scales of the problem, providing the necessary coherence to the system to treat is as a flow. We can thus study them using relativistic magnetohydrodynamics. As a first step, we have studied the structure of steady-state configurations of jets by using numerical simulations. We have used a helical field configuration and have changed different relevant parameters that control the way in which the energy flux is distributed in jets (namely, the proportion of the energy flux carried by internal, kinetic or magnetic energy). Our results show significant differences among the different kinds of jets. Finally, we also report on results based on synthetic maps of our simulated jets.
△ Less
Submitted 4 August, 2017;
originally announced August 2017.
-
Radio mode feedback: Does relativity matter?
Authors:
Manel Perucho,
José María Martí,
Vicent Quilis,
Marina Borja-Lloret
Abstract:
Radio mode feedback, associated with the propagation of powerful outflows in active galaxies, is a crucial ingredient in galaxy evolution. Extragalactic jets are well collimated and relativistic, both in terms of thermodynamics and kinematics. They generate strong shocks in the ambient medium, associated with observed hotspots, and carve cavities that are filled with the shocked jet flow. In this…
▽ More
Radio mode feedback, associated with the propagation of powerful outflows in active galaxies, is a crucial ingredient in galaxy evolution. Extragalactic jets are well collimated and relativistic, both in terms of thermodynamics and kinematics. They generate strong shocks in the ambient medium, associated with observed hotspots, and carve cavities that are filled with the shocked jet flow. In this Letter, we compare the pressure evolution in the hotspot and the cavity generated by relativistic and classical jets. Our results show that the classical approach underestimates the cavity pressure by a factor larger or equal to 2 for a given shocked volume during the whole active phase. The tension between both approaches can only be alleviated by unrealistic jet flow densities or gigantic jet areas in the classical case. As a consequence, the efficiency of a relativistic jet heating the ambient is typically around 20% larger compared with a classical jet, and the heated volume is 2 to 10 times larger during the time evolution. This conflict translates into two substantially disparate manners, both spatially and temporal, of heating the ambient medium. These differences are expected to have relevant implications on the star formation rates of the host galaxies and their evolution.
△ Less
Submitted 18 July, 2017;
originally announced July 2017.
-
Okutama-Action: An Aerial View Video Dataset for Concurrent Human Action Detection
Authors:
Mohammadamin Barekatain,
Miquel Martí,
Hsueh-Fu Shih,
Samuel Murray,
Kotaro Nakayama,
Yutaka Matsuo,
Helmut Prendinger
Abstract:
Despite significant progress in the development of human action detection datasets and algorithms, no current dataset is representative of real-world aerial view scenarios. We present Okutama-Action, a new video dataset for aerial view concurrent human action detection. It consists of 43 minute-long fully-annotated sequences with 12 action classes. Okutama-Action features many challenges missing i…
▽ More
Despite significant progress in the development of human action detection datasets and algorithms, no current dataset is representative of real-world aerial view scenarios. We present Okutama-Action, a new video dataset for aerial view concurrent human action detection. It consists of 43 minute-long fully-annotated sequences with 12 action classes. Okutama-Action features many challenges missing in current datasets, including dynamic transition of actions, significant changes in scale and aspect ratio, abrupt camera movement, as well as multi-labeled actors. As a result, our dataset is more challenging than existing ones, and will help push the field forward to enable real-world applications.
△ Less
Submitted 15 June, 2017; v1 submitted 9 June, 2017;
originally announced June 2017.
-
The internal structure of overpressured, magnetized, relativistic jets
Authors:
José M. Martí,
Manel Perucho,
José L. Gómez
Abstract:
This work presents the first characterization of the internal structure of overpressured steady superfast magnetosonic relativistic jets in connection with their dominant type of energy. To this aim, relativistic magnetohydrodynamic simulations of different jet models threaded by a helical magnetic field have been analyzed covering a wide region in the magnetosonic Mach number - specific internal…
▽ More
This work presents the first characterization of the internal structure of overpressured steady superfast magnetosonic relativistic jets in connection with their dominant type of energy. To this aim, relativistic magnetohydrodynamic simulations of different jet models threaded by a helical magnetic field have been analyzed covering a wide region in the magnetosonic Mach number - specific internal energy plane. The merit of this plane is that models dominated by different types of energy (internal energy: hot jets; rest-mass energy: kinetically dominated jets; magnetic energy: Poynting-flux dominated jets) occupy well separated regions. The analyzed models also cover a wide range of magnetizations. Models dominated by the internal energy (i.e., hot models, or Poynting-flux dominated jets with magnetizations larger than but close to 1) have a rich internal structure characterized by a series of recollimation shocks and present the largest variations in the flow Lorentz factor (and internal energy density). Conversely, in kinetically dominated models there is not much internal nor magnetic energy to be converted into kinetic one and the jets are featureless, with small variations in the flow Lorentz factor. The presence of a significant toroidal magnetic field threading the jet produces large gradients in the transversal profile of the internal energy density. Poynting-flux dominated models with high magnetization (>10) are prone to be unstable against magnetic pinch modes, which sets limits to the expected magnetization in parsec-scale AGN jets {and/or constrains their magnetic field configuration}.
△ Less
Submitted 12 September, 2016; v1 submitted 2 September, 2016;
originally announced September 2016.
-
Modular Models for Intuitionistic Justification Logic
Authors:
Michel Marti,
Thomas Studer
Abstract:
We present the logic iJT4, which is an explicit version of intuitionistic S4 and establish soundness and completeness with respect to modular models.
We present the logic iJT4, which is an explicit version of intuitionistic S4 and establish soundness and completeness with respect to modular models.
△ Less
Submitted 25 April, 2016;
originally announced April 2016.
-
The structure of steady, relativistic, magnetised jets with rotation
Authors:
José-María Martí
Abstract:
We present equilibrium models of relativistic magnetised, infinite, axisymmetric jets with rotation propagating through an homogeneous, unmagnetised ambient medium at rest. The jet models are characterised by six functions defining the radial profiles of density, pressure, and the toroidal and axial components of velocity and magnetic field. Fixing the ambient pressure and the jet rest-mass densit…
▽ More
We present equilibrium models of relativistic magnetised, infinite, axisymmetric jets with rotation propagating through an homogeneous, unmagnetised ambient medium at rest. The jet models are characterised by six functions defining the radial profiles of density, pressure, and the toroidal and axial components of velocity and magnetic field. Fixing the ambient pressure and the jet rest-mass density and axial components of the flow velocity and magnetic field, we analyze the influence of the toroidal magnetic field and several rotation laws on the structure of the equilibrium models. Our approach excludes by construction the analysis of the self-consistently magnetically launched jet models or the force-free equilibrium solutions. Several forbidden regions in the magnetic pitch angle/magnetization plane are found where models of the class considered in our study could not be settled. These forbidden regions are associated with the existence of maximum axial and toroidal magnetic field components compatible with the prescribed equilibrium condition at the jet surface, and/or an excess of centrifugal force producing gaps with negative pressures in the jet. The present study can be easily extended to jet models with different transversal profiles and magnetic field configurations.
In the last part of the paper, we test the ability of our RMHD code to maintain steady equilibrium models of axisymmetric RMHD jets in one and two spatial dimensions. The one dimensional numerical simulations serve also as a consistency proof of the fidelity of the analytical steady solutions discussed in the first part of the paper. The present study allows us to build initial equilibrium jet models with selected properties for dynamical (and emission) simulations of magnetised relativistic jets with rotation.
△ Less
Submitted 22 June, 2015;
originally announced June 2015.
-
The Spanish Square Kilometre Array White Book
Authors:
J. A. Acosta-Pulido,
I. Agudo,
A. Alberdi,
J. Alcolea,
E. J. Alfaro,
A. Alonso-Herrero,
G. Anglada,
P. Arnalte-Mur,
Y. Ascasibar,
B. Ascaso,
R. Azulay,
R. Bachiller,
A. Baez-Rubio,
E. Battaner,
J. Blasco,
C. B. Brook,
V. Bujarrabal,
G. Busquet,
M. D. Caballero-Garcia,
C. Carrasco-Gonzalez,
J. Casares,
A. J. Castro-Tirado,
L. Colina,
F. Colomer,
I. de Gregorio-Monsalvo
, et al. (94 additional authors not shown)
Abstract:
The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on…
▽ More
The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.
△ Less
Submitted 17 June, 2015; v1 submitted 10 June, 2015;
originally announced June 2015.
-
On the correction of conserved variables for numerical RMHD with staggered constrained transport
Authors:
José M. Martí
Abstract:
Despite the success of the combination of conservative schemes and staggered constrained transport algorithms in the last fifteen years, the accurate description of highly magnetized, relativistic flows with strong shocks represents still a challenge in numerical RMHD. The present paper focusses in the accuracy and robustness of several correction algorithms for the conserved variables, which has…
▽ More
Despite the success of the combination of conservative schemes and staggered constrained transport algorithms in the last fifteen years, the accurate description of highly magnetized, relativistic flows with strong shocks represents still a challenge in numerical RMHD. The present paper focusses in the accuracy and robustness of several correction algorithms for the conserved variables, which has become a crucial ingredient in the numerical simulation of problems where the magnetic pressure dominates over the thermal pressure by more than two orders of magnitude.
Two versions of non-relativistic and fully relativistic corrections have been tested and compared using a magnetized cylindrical explosion with high magnetization ($ \ge 10^4$) as test. In the non-relativistic corrections, the total energy is corrected for the difference in the classical magnetic energy term between the average of the staggered fields and the conservative ones, before (CA1) and after (CA1') recovering the primitive variables. These corrections are unable to pass the test at any numerical resolution. The two relativistic approaches (CA2 and CA2'), correcting also the magnetic terms depending on the flow speed in both the momentum and the total energy, reveal as much more robust. These algorithms pass the test succesfully and with very small deviations of the energy conservation ($\le 10^{-4}$), and very low values of the total momentum ($\le 10^{-8}$). In particular, the algorithm CA2' (that corrects the conserved variables after recovering the primitive variables) passes the test at all resolutions.
The numerical code used to run all the test cases is briefly described.
△ Less
Submitted 12 January, 2015;
originally announced January 2015.
-
Large-scale jets from active galactic nuclei as a source of ICM heating: cavities and shocks
Authors:
Manel Perucho,
José María Martí,
Vicent Quilis,
Elena Ricciardelli
Abstract:
The evolution of powerful extragalactic jets is not only interesting by itself, but also for its impact on the evolution of the host galaxy and its surroundings. We have performed long-term axisymmetric numerical simulations of relativistic jets with different powers to study their evolution through an environment with a pressure and density gradient. Our results show key differences in the evolut…
▽ More
The evolution of powerful extragalactic jets is not only interesting by itself, but also for its impact on the evolution of the host galaxy and its surroundings. We have performed long-term axisymmetric numerical simulations of relativistic jets with different powers to study their evolution through an environment with a pressure and density gradient. Our results show key differences in the evolution of jets with different powers in terms of the spatial and temporal scales of energy deposition. According to our results, the observed morphology in X-ray cavities requires that an important fraction of the jet's energetic budget is in the form of internal energy. Thus, light, lepton-dominated jets are favoured. In all cases, heating is mainly produced by shocks. Cavity overpressure is sustained by an important population of thermal particles. Our simulations reproduce the cool-core structure in projected, luminosity-weighted temperature. We have performed an additional simulation of a slow, massive jet and discuss the differences with its relativistic counterparts. Important qualitative and quantitative differences are found between the non-relativistic and the relativistic jets. Our conclusions point towards a dual-mode of AGN kinetic feedback, depending on the jet power.
△ Less
Submitted 11 September, 2014;
originally announced September 2014.
-
On the deceleration of FRI jets: mass loading by stellar-winds
Authors:
Manel Perucho,
José María Martí,
Robert A. Laing,
Philip E. Hardee
Abstract:
Jets in low-luminosity radio galaxies are known to decelerate from relativistic speeds on parsec scales to mildly or sub-relativistic speeds on kiloparsec scales. Several mechanisms have been proposed to explain this effect, including strong reconfinement shocks and the growth of instabilities (both leading to boundary-layer entrainment) and mass loading from stellar winds or molecular clouds. We…
▽ More
Jets in low-luminosity radio galaxies are known to decelerate from relativistic speeds on parsec scales to mildly or sub-relativistic speeds on kiloparsec scales. Several mechanisms have been proposed to explain this effect, including strong reconfinement shocks and the growth of instabilities (both leading to boundary-layer entrainment) and mass loading from stellar winds or molecular clouds. We have performed a series of axisymmetric simulations of the early evolution of jets in a realistic ambient medium to probe the effects of mass loading from stellar winds using the code Ratpenat. We study the evolution of Fanaroff-Riley Class I (FRI) jets, with kinetic powers L_j \sim 1.e41-1.e44 erg/s, within the first 2 kpc of their evolution, where deceleration by stellar mass loading should be most effective. Mass entrainment rates consistent with present models of stellar mass loss in elliptical galaxies produce deceleration and effective decollimation of weak FRI jets within the first kiloparsec. However, powerful FRI jets are not decelerated significantly. In those cases where the mass loading is important, the jets show larger opening angles and decollimate at smaller distances, but the overall structure and dynamics of the bow-shock are similar to those of unloaded jets with the same power and thrust. According to our results, the flaring observed on kpc scales is initiated by mass loading in the weaker FRI jets and by reconfinement shocks or the growth of instabilities in the more powerful jets. The final mechanism of decollimation and deceleration is always the development of disruptive pinching modes.
△ Less
Submitted 4 April, 2014;
originally announced April 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.
-
Efficient modeling of laser plasma interactions in high energy density scenarios
Authors:
F. Fiuza,
M. Marti,
R. A. Fonseca,
L. O. Silva,
J. Tonge,
J. May,
W. B. Mori
Abstract:
We describe how a new framework for coupling a full-PIC algorithm with a reduced PIC algorithm has been implemented into the code OSIRIS. We show that OSIRIS with this new hybrid-PIC algorithm can efficiently and accurately model high energy density scenarios such as ion acceleration in laser-solid interactions and fast ignition of fusion targets. We model for the first time the full density range…
▽ More
We describe how a new framework for coupling a full-PIC algorithm with a reduced PIC algorithm has been implemented into the code OSIRIS. We show that OSIRIS with this new hybrid-PIC algorithm can efficiently and accurately model high energy density scenarios such as ion acceleration in laser-solid interactions and fast ignition of fusion targets. We model for the first time the full density range of a fast ignition target in a fully self-consistent hybrid-PIC simulation, illustrating the possibility of stopping the laser generated electron flux at the core region with relatively high efficiencies. Computational speedups greater than 1000 times are demonstrated, opening the way for full-scale multi-dimensional modeling of high energy density scenarios and for the guiding of future experiments.
△ Less
Submitted 14 May, 2012;
originally announced May 2012.
-
Intracluster Medium reheating by relativistic jets
Authors:
Manel Perucho,
Vicent Quilis,
José María Martí
Abstract:
Galactic jets are powerful energy sources reheating the intra-cluster medium in galaxy clusters. Their crucial role in the cosmic puzzle, motivated by observations, has been established by a great number of numerical simulations missing the relativistic nature of these jets. We present the first relativistic simulations of the very long term evolution of realistic galactic jets. Unexpectedly, our…
▽ More
Galactic jets are powerful energy sources reheating the intra-cluster medium in galaxy clusters. Their crucial role in the cosmic puzzle, motivated by observations, has been established by a great number of numerical simulations missing the relativistic nature of these jets. We present the first relativistic simulations of the very long term evolution of realistic galactic jets. Unexpectedly, our results show no buoyant bubbles, but large cocoon regions compatible with the observed X-ray cavities. The reheating is more efficient and faster than in previous scenarios, and it is produced by the shock wave driven by the jet, that survives for several hundreds of Myrs. Therefore, the X-ray cavities in clusters produced by powerful relativistic jets would remain confined by weak shocks for extremely long periods, whose detection could be an observational challenge.
△ Less
Submitted 13 October, 2011; v1 submitted 27 July, 2011;
originally announced July 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.
-
Stability of three-dimensional relativistic jets: implications for jet collimation
Authors:
M. Perucho,
J. M. Marti,
J. M. Cela,
M. Hanasz,
R. de la Cruz,
F. Rubio
Abstract:
The stable propagation of jets in FRII sources is remarkable if one takes into account that large-scale jets are subjected to potentially highly disruptive three-dimensional (3D) Kelvin-Helmholtz instabilities. Numerical simulations can address this problem and help clarify the causes of this remarkable stability. Following previous studies of the stability of relativistic flows in two dimensions…
▽ More
The stable propagation of jets in FRII sources is remarkable if one takes into account that large-scale jets are subjected to potentially highly disruptive three-dimensional (3D) Kelvin-Helmholtz instabilities. Numerical simulations can address this problem and help clarify the causes of this remarkable stability. Following previous studies of the stability of relativistic flows in two dimensions (2D), it is our aim to test and extend the conclusions of such works to three dimensions. We present numerical simulations for the study of the stability properties of 3D, sheared, relativistic flows. This work uses a fully parallelized code Ratpenat that solves equations of relativistic hydrodynamics in 3D. The results of the present simulations confirm those in 2D. We conclude that the growth of resonant modes in sheared relativistic flows could be important in explaining the long-term collimation of extragalactic jets.
△ Less
Submitted 24 May, 2010;
originally announced May 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.
-
Statistical kinetic treatment of relativistic binary collisions
Authors:
F. Peano,
M. Marti,
L. O. Silva,
G. Coppa
Abstract:
In particle-based algorithms, the effect of binary collisions is commonly described in a statistical way, using Monte Carlo techniques. It is shown that, in the relativistic regime, stringent constraints should be considered on the sampling of particle pairs for collision, which are critical to ensure physically meaningful results, and that nonrelativistic sampling criteria (e.g., uniform random…
▽ More
In particle-based algorithms, the effect of binary collisions is commonly described in a statistical way, using Monte Carlo techniques. It is shown that, in the relativistic regime, stringent constraints should be considered on the sampling of particle pairs for collision, which are critical to ensure physically meaningful results, and that nonrelativistic sampling criteria (e.g., uniform random pairing) yield qualitatively wrong results, including equilibrium distributions that differ from the theoretical Jüttner distribution. A general procedure for relativistically consistent algorithms is provided, and verified with three-dimensional Monte Carlo simulations, thus opening the way to the numerical exploration of the statistical properties of collisional relativistic systems.
△ Less
Submitted 10 February, 2009;
originally announced February 2009.
-
Spectral evolution of superluminal components in parsec-scale jets
Authors:
P. Mimica,
M. A. Aloy,
I. Agudo,
J. M. Marti,
J. L. Gomez,
J. A. Miralles
Abstract:
(Abridged) We present numerical simulations of the spectral evolution and emission of radio components in relativistic jets. We have developed an algorithm (SPEV) for the transport of a population of non-thermal electrons including radiative losses. For large values of the ratio of gas pressure to magnetic field energy density, \ab \sim 6\times 10^4, quiescent jet models show substantial spectra…
▽ More
(Abridged) We present numerical simulations of the spectral evolution and emission of radio components in relativistic jets. We have developed an algorithm (SPEV) for the transport of a population of non-thermal electrons including radiative losses. For large values of the ratio of gas pressure to magnetic field energy density, \ab \sim 6\times 10^4, quiescent jet models show substantial spectral evolution, with observational consequences only above radio frequencies. Larger values of the magnetic field (\ab \sim 6\times 10^2), such that synchrotron losses are moderately important at radio frequencies, present a larger ratio of shocked-to-unshocked regions brightness than the models without radiative losses, despite the fact that they correspond to the same underlying hydrodynamic structure. We also show that jets with a positive photon spectral index result if the lower limit γ_min of the non-thermal particle energy distribution is large enough. A temporary increase of the Lorentz factor at the jet inlet produces a traveling perturbation that appears in the synthetic maps as a superluminal component. We show that trailing components can be originated not only in pressure matched jets, but also in over-pressured ones, where the existence of recollimation shocks does not allow for a direct identification of such features as Kelvin-Helmholtz modes, and its observational imprint depends on the observing frequency. If the magnetic field is large (\ab \sim 6\times 10^2), the spectral index in the rarefaction trailing the traveling perturbation does not change much with respect to the same model without any hydrodynamic perturbation. If the synchrotron losses are considered the spectral index displays a smaller value than in the corresponding region of the quiescent jet model.
△ Less
Submitted 20 February, 2009; v1 submitted 7 November, 2008;
originally announced November 2008.
-
The Influence of Helical Magnetic Fields in the Dynamics and Emission of Relativistic Jets
Authors:
M. Roca-Sogorb,
M. Perucho,
J. L. Gomez,
J. M. Marti,
L. Anton,
M. A. Aloy,
I. Agudo
Abstract:
We present numerical relativistic magnetohydrodynamic and emission simulations aimed to study the role played by the magnetic field in the dynamics and emission of relativistic jets in Active Galactic Nuclei. We focus our analysis on the study of the emission from recollimation shocks since they may provide an interpretation for the stationary components seen at parsec-scales in multiple sources…
▽ More
We present numerical relativistic magnetohydrodynamic and emission simulations aimed to study the role played by the magnetic field in the dynamics and emission of relativistic jets in Active Galactic Nuclei. We focus our analysis on the study of the emission from recollimation shocks since they may provide an interpretation for the stationary components seen at parsec-scales in multiple sources. We show that the relative brightness of the knots associated with the recollimation shocks decreases with increasing jet magnetization, suggesting that jets presenting stationary components may have a relatively weak magnetization, with magnetic fields of the order of equipartition or below.
△ Less
Submitted 26 February, 2008;
originally announced February 2008.
-
A numerical simulation of the evolution and fate of a FRI jet. The case of 3C 31
Authors:
M. Perucho,
J. M. Marti
Abstract:
The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. In this paper, we test the present theoretical and observational models via a relativistic numerical simulation of the jets in the radio galaxy 3C 31. We use the parameters derived from the modelling presented by \cite{lb02a,lb02b} as input parameters for the simulation of the evolution of the source, thu…
▽ More
The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. In this paper, we test the present theoretical and observational models via a relativistic numerical simulation of the jets in the radio galaxy 3C 31. We use the parameters derived from the modelling presented by \cite{lb02a,lb02b} as input parameters for the simulation of the evolution of the source, thus assuming that they have not varied over the lifetime of the source. We simulate about 10 % of the total lifetime of the jets in 3C 31. Realistic density and pressure gradients for the atmosphere are used. The simulation includes an equation of state for a two-component relativistic gas that allows a separate treatment of leptonic and baryonic matter. We compare our results with the modelling of the observational data of the source. Our results show that the bow shock evolves self-similarly at a quasi-constant speed, with slight deceleration by the end of the simulation, in agreement with recent X-ray observations that show the presence of bow shocks in FRI sources. The jet expands until it becomes underpressured with respect to the ambient medium, and then recollimates. Subsequent oscillations around pressure equilibrium and generation of standing shocks lead to the mass loading and disruption of the jet flow. We derive an estimate for the minimum age of the source of $t>1. 10^8 \rm{yrs}$, which may imply continuous activity of 3C 31 since the triggering of its activity. The simulation shows that weak CSS sources may be the young counterparts of FRIs. We conclude that the observed properties of the jets in 3C 31 are basically recovered by the standing shock scenario.
△ Less
Submitted 12 September, 2007;
originally announced September 2007.
-
Analytical and numerical studies of fluid instabilities in relativistic jets
Authors:
M. Perucho,
A. P. Lobanov,
J. M. Marti
Abstract:
Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed…
▽ More
Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed in the parsec-scale jet of the well known quasar 3C 273. We compare the model with the radio structure observed in 3C 273 on parsec scales using very long baseline interferometry (VLBI) and constrain the basic properties of the flow. Our results show that Kelvin-Helmholtz instabilities are the most plausible mechanism to generate the observed structures.
△ Less
Submitted 8 February, 2006;
originally announced February 2006.
-
Very High Mach Number Electrostatic Shocks in Collisionless Plasmas
Authors:
G. Sorasio,
M. Marti,
R. Fonseca,
L. O. Silva
Abstract:
The kinetic theory of collisionless electrostatic shocks resulting from the collision of plasma slabs with different temperatures and densities is presented. The theoretical results are confirmed by self-consistent particle-in-cell simulations, revealing the formation and stable propagation of electrostatic shocks with very high Mach numbers ($M \gg 10$), well above the predictions of the classi…
▽ More
The kinetic theory of collisionless electrostatic shocks resulting from the collision of plasma slabs with different temperatures and densities is presented. The theoretical results are confirmed by self-consistent particle-in-cell simulations, revealing the formation and stable propagation of electrostatic shocks with very high Mach numbers ($M \gg 10$), well above the predictions of the classical theories for electrostatic shocks.
△ Less
Submitted 23 December, 2005;
originally announced December 2005.
-
Nonlinear stability of relativistic sheared planar jets
Authors:
M. Perucho,
J. M. Marti,
M. Hanasz
Abstract:
The linear and non-linear stability of sheared, relativistic planar jets is studied by means of linear stability analysis and numerical hydrodynamical simulations. Our results extend the previous Kelvin-Hemlholtz stability studies for relativistic, planar jets in the vortex sheet approximation performed by Perucho et al. (2004a,b) by including a shear layer between the jet and the external mediu…
▽ More
The linear and non-linear stability of sheared, relativistic planar jets is studied by means of linear stability analysis and numerical hydrodynamical simulations. Our results extend the previous Kelvin-Hemlholtz stability studies for relativistic, planar jets in the vortex sheet approximation performed by Perucho et al. (2004a,b) by including a shear layer between the jet and the external medium and more general perturbations. The models considered span a wide range of Lorentz factors ($2.5-20$) and internal energies ($0.08 c^2-60 c^2$) and are classified into three classes according to the main characteristics of their long-term, non-linear evolution. We observe a clear separation of these three groups in a relativistic Mach-number Lorentz-factor plane. Jets with a low Lorentz factor and small relativistic Mach number are disrupted after saturation. Those with a large Lorentz factor and large relativistic Mach number are the stablest, due to the appearance of short wavelength resonant modes which generate local mixing and heating in the shear layer around a fast, unmixed core, giving a plausible solution for the problem of the long-term stability of relativistic jets. A third group is present between them, including jets with intermediate values of Lorentz factor and relativistic Mach number, which are disrupted by a slow process of mixing favored by an efficient and continuous conversion of kinetic into internal energy. In the long term, all the models develop a distinct transversal structure (shear/transition layers) as a consequence of KH perturbation growth, depending on the class they belong to. The properties of these shear layers are analyzed in connection with the parameters of the original jet models.
△ Less
Submitted 14 October, 2005;
originally announced October 2005.
-
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.
-
Real vs. simulated relativistic jets
Authors:
J. L. Gomez,
J. M. Marti,
I. Agudo,
A. P. Marscher,
S. G. Jorstad,
M. A. Aloy
Abstract:
Intensive VLBI monitoring programs of jets in AGN are showing the existence of intricate emission patterns, such as upstream motions or slow moving and quasi-stationary componentes trailing superluminal features. Relativistic hydrodynamic and emission simulations of jets are in very good agreement with these observations, proving as a powerful tool for the understanding of the physical processes…
▽ More
Intensive VLBI monitoring programs of jets in AGN are showing the existence of intricate emission patterns, such as upstream motions or slow moving and quasi-stationary componentes trailing superluminal features. Relativistic hydrodynamic and emission simulations of jets are in very good agreement with these observations, proving as a powerful tool for the understanding of the physical processes taking place in the jets of AGN, microquasars and GRBs. These simulations show that the variability of the jet emission is the result of a complex combination of phase motions, viewing angle selection effects, and non-linear interactions between perturbations and the underlying jet and/or ambient medium. Both observations and simulations suggest that shock-in-jet models may be an overly simplistic idealization when interpreting the emission structure observed in actual jets.
△ Less
Submitted 3 January, 2005;
originally announced January 2005.
-
Simulations of the relativistic parsec-scale jet in 3C273
Authors:
M. Perucho,
A. P. Lobanov,
J. M. Marti
Abstract:
We present a hydrodynamical 3D simulation of the relativistic jet in 3C273, in comparison to previous linear perturbation analysis of Kelvin-Helmholtz instability developing in the jet. Our aim is to assess advantages and limitations of both analytical and numerical approaches and to identify spatial and temporal scales on which the linear regime of Kelvin-Helmholtz instability can be applied in…
▽ More
We present a hydrodynamical 3D simulation of the relativistic jet in 3C273, in comparison to previous linear perturbation analysis of Kelvin-Helmholtz instability developing in the jet. Our aim is to assess advantages and limitations of both analytical and numerical approaches and to identify spatial and temporal scales on which the linear regime of Kelvin-Helmholtz instability can be applied in studies of morphology and kinematics of parsec-scale jets.
△ Less
Submitted 17 December, 2004;
originally announced December 2004.
-
A dynamical model for the CSO-MSO-FRII evolution: hints from hot spot properties
Authors:
M. Perucho,
J. M. Marti
Abstract:
Compact Symmetric Objects (CSOs) are considered the young counterparts of large doubles according to advance speeds measured or inferred from spectral ageing. Here we present a simple power law model for the CSO/FRII evolution based on the study of sources with well defined hot-spots. The luminosity of the hot spots is estimated under minimum energy conditions. The advance of the source is consi…
▽ More
Compact Symmetric Objects (CSOs) are considered the young counterparts of large doubles according to advance speeds measured or inferred from spectral ageing. Here we present a simple power law model for the CSO/FRII evolution based on the study of sources with well defined hot-spots. The luminosity of the hot spots is estimated under minimum energy conditions. The advance of the source is considered to proceed in ram pressure equilibrium with the ambient medium. Finally, we also assume that the jets feeding the hot spots are relativistic and have a time dependent power. Comparison with observational data allows to interpret the CSO-FRII evolution in terms of decreasing jet power with time.
△ Less
Submitted 28 July, 2004;
originally announced July 2004.
-
Stability of hydrodynamical relativistic planar jets. II. Long-term nonlinear evolution
Authors:
M. Perucho,
J. M. Marti,
M. Hanasz
Abstract:
In this paper we continue our study of the Kelvin-Helmholtz (KH) instability in relativistic planar jets following the long-term evolution of the numerical simulations which were introduced in Paper I. The models have been classified into four classes (I to IV) with regard to their evolution in the nonlinear phase, characterized by the process of jet/ambient mixing and momentum transfer. Models…
▽ More
In this paper we continue our study of the Kelvin-Helmholtz (KH) instability in relativistic planar jets following the long-term evolution of the numerical simulations which were introduced in Paper I. The models have been classified into four classes (I to IV) with regard to their evolution in the nonlinear phase, characterized by the process of jet/ambient mixing and momentum transfer. Models undergoing qualitatively different non-linear evolution are clearly grouped in well-separated regions in a jet Lorentz factor/jet-to-ambient enthalpy diagram. Jets with a low Lorentz factor and small enthalpy ratio are disrupted by a strong shock after saturation. Those with a large Lorentz factor and enthalpy ratio are unstable although the process of mixing and momentum exchange proceeds to a longer time scale due to a steady conversion of kinetic to internal energy in the jet. In these cases, the high value of the initial Lorentz seems to prevent transversal velocity from growing far enough to generate the strong shock that breaks the slower jets. Finally, jets with either high Lorentz factors and small enthalpy ratios or low Lorentz factors and large enthalpy ratios appear as the most stable.
△ Less
Submitted 28 July, 2004;
originally announced July 2004.
-
Stability of hydrodynamical relativistic planar jets. I. Linear evolution and saturation of Kelvin-Helmholtz modes
Authors:
M. Perucho,
M. Hanasz,
J. M. Marti,
H. Sol
Abstract:
The effects of relativistic dynamics and thermodynamics in the development of Kelvin-Helmholtz instabilities in planar, relativistic jets along the early phases (namely linear and saturation phases) of evolution has been studied by a combination of linear stability analysis and high-resolution numerical simulations for the most unstable first reflection modes in the temporal approach. Three diff…
▽ More
The effects of relativistic dynamics and thermodynamics in the development of Kelvin-Helmholtz instabilities in planar, relativistic jets along the early phases (namely linear and saturation phases) of evolution has been studied by a combination of linear stability analysis and high-resolution numerical simulations for the most unstable first reflection modes in the temporal approach. Three different values of the jet Lorentz factor (5, 10 and 20) and a few different values of specific internal energy of the jet matter (from 0.08 to $60.0 c^2$) have been considered. Figures illustrating the evolution of the perturbations are also shown.
△ Less
Submitted 27 July, 2004;
originally announced July 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.