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.

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.

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.

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.

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.

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.

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.

We consider the transverse structure and stability properties of relativistic jets formed in the course of the collapse of a massive progenitor. Our numerical simulations show the presence of a strong shear in the bulk velocity of such jets. This shear can be responsible for a very rapid shear--driven instability that arises for any velocity profile. This conclusion has been confirmed both by numerical simulations and theoretical analysis. The instability leads to rapid fluctuations of the main hydrodynamical parameters (density, pressure, Lorentz factor, etc.). However, the perturbations of the density are effectively decoupled from those of the pressure because the beam of the jet is radiation--dominated. The characteristic growth time of instability is much shorter than the life time of the jet and, therefore, may lead to a complete turbulent beam. In the course of the non-linear evolution, these fluctuations may yield to internal shocks which can be randomly distributed in the jet. In the case that internal shocks in a ultrarelativistic outflow are responsible for the observed phenomenology of gamma-ray bursts, the proposed instability can well account for the short-term variability of gamma-ray light curves down to milliseconds.

We study the asymptotic behaviour of scaling solutions with a dissipative fluid and we show that, contrary to recent claims, the existence of stable accelerating attractor solution which solves the `energy' coincidence problem depends crucially on the chosen equations of state for the thermodynamical variables. We discuss two types of equations of state, one which contradicts this claim, and one which supports it.

Two possible interpretations of FRW cosmologies (perfect fluid or dissipative fluid)are considered as consecutive phases of the system. Necessary conditions are found, for the transition from perfect fluid to dissipative regime to occur, bringing out the conspicuous role played by a particular state of the system (the ''critical point '').

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.

The evolution of an homogeneous and isotropic dissipative fluid is analyzed using dynamical systems techniques. The dissipation is driven by bulk viscous pressure and the truncated Israel-Stewart theory is used. Although almost all solutions inflate, we show that only few of them can be considered as physical solutions since the dominant energy condition is not satisfied.

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.

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 $\rho_{\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.

Tensor-scalar theory of gravity allows the generation of gravitational waves from astrophysical sources, like Supernov\ae{}, 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.

We present the first radio emission simulations from high resolution three dimensional relativistic hydrodynamic jets, which allow for a study of the observational implications of the interaction between the jet and external medium. This interaction gives rise to a stratification of the jet where a fast spine is surrounded by a slow high energy shear layer. The stratification, and in particular the large specific internal energy and slow flow in the shear layer largely determines the emission from the jet. If the magnetic field in the shear layer becomes helical (e.g., resulting from an initial toroidal field and an aligned field component generated by shear) the emission shows a cross section asymmetry, in which either the top or the bottom of the jet dominates the emission. This, as well as limb or spine brightening, is a function of the viewing angle and flow velocity, and the top/bottom jet emission predominance can be reversed if the jet changes direction with respect to the observer, or presents a change in velocity. The asymmetry is more prominent in the polarized flux, because of field cancellation (or amplification) along the line of sight. Recent observations of jet cross section emission asymmetries in the blazar 1055+018 can be explained assuming the existence of a shear layer with a helical magnetic field.

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.

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.

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.

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.

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.

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.

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.

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).

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.

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.

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.

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.

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.

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.

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.

We calculate the Legendre expansion of the rate of the process $\nu + \bar{\nu} \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.

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.

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.

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.

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.