Cosmic ray (CR) protons are an important component in many astrophysical systems. Processes like CR injection, cooling, adiabatic changes as well as active CR transport through the medium strongly modify the CR momentum distribution and have to be taken into account in hydrodynamical simulations. We present an efficient novel numerical scheme to accurately compute the evolution of the particle distribution function by solving the Fokker-Planck equation with a low number of spectral bins (10 - 20), which is required to include a full spectrum for every computational fluid element. The distribution function is represented by piecewise power laws and is not forced to be continuous, which enables an optimal representation of the spectrum. The Fokker-Planck equation is solved with a two-moment approach evolving the CR number and energy density. The low numerical diffusion of the scheme reduces the numerical errors by orders of magnitude in comparison to classical schemes with piecewise constant spectral representations. With this method not only the spectral evolution of CRs can be computed accurately in magnetohydrodynamic simulations but also their dynamical impact as well as CR ionisation. This allows for more accurate models for astrophysical plasmas, like the interstellar medium, and direct comparisons with observations.

We investigate the impact of cosmic rays (CRs) on galactic outflows from a multi-phase interstellar medium with solar neighbourhood conditions. The three-dimensional magneto-hydrodynamical simulations include CRs as a relativistic fluid in the advection-diffusion approximation. The thermal and chemical state of the ISM is computed with a non-equilibrium chemical network. We find that CRs (injected with 10 \% of the supernova energy) efficiently support the launching of outflows and strongly affect their phase structure. Outflows leaving the midplane are denser ($\rho \sim 10^{-26}\,\mathrm{g\,cm}^{-3}$), colder ($\sim 10^4\,\mathrm{K}$), and slower ($\sim 30\,\mathrm{km\,s}^{-1}$) if CRs are considered in addition to thermal SNe. The CR supported outflows are also smoother, in particular at larger heights ($> 1\,\mathrm{kpc}$ above the midplane) without the direct impact of SN explosions. Approximately $5\% - 25\%$ of the injected CR energy is lost via hadronic cooling. Smaller diffusion coefficients lead to slightly larger hadronic losses but allow for steeper CR pressure gradients, stronger outflows and larger accelerations. Up to a height of $z \sim1\,\mathrm{kpc}$ there are large volumes in approximate pressure equilibrium between thermal and CR component. At larger altitudes the CR pressure is $10-100$ times as large as the thermal counterpart. More than $\sim 1\,\mathrm{kpc}$ away from the midplane, CRs provide the dominant gas acceleration mechanism.

We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is computed using a chemical network that follows the abundances of H+, H, H2, CO, C+, and free electrons and includes (self-)shielding of the gas and dust. We find that CRs perceptibly thicken the disk with the heights of 90% (70%) enclosed mass reaching ~1.5 kpc (~0.2 kpc). The simulations indicate that CRs alone can launch and sustain strong outflows of atomic and ionized gas with mass loading factors of order unity, even in solar neighborhood conditions and with a CR energy injection per SN of 10^50 erg, 10% of the fiducial thermal energy of an SN. The CR-driven outflows have moderate launching velocities close to the midplane (~100 km/s) and are denser (\rho~1e-24 - 1e-26 g/cm^3), smoother, and colder than the (thermal) SN-driven winds. The simulations support the importance of CRs for setting the vertical structure of the disk as well as the driving of winds.

We present a novel implementation of cosmic rays (CR) in the magneto-hydrodynamic code FLASH. CRs are described as separate fluids with different energies. CR advection, energy dependent anisotropic diffusion with respect to the magnetic field and adiabatic losses to follow the evolution of spectra are taken into account. We present a first study of the transport and immediate (~150 kyr) dynamical impact of CRs on the turbulent magnetised interstellar medium around supernova remnants on scales up to 80 pc. CR diffusion quickly leads to an efficient acceleration of low-density gas (mainly perpendicular to the magnetic field) with accelerations up to two orders of magnitude above the thermal values. Peaked (at 1 GeV) CR injection spectra have a stronger impact on the dynamics than power-law spectra. For realistic magnetic field configurations low energy CRs (with smaller diffusion coefficients) distribute anisotropically with large spatial variations of a factor of ten and more. Adiabatic losses can change the local spectra perceptibly but do not have an integral effect on the dynamics at the spatial and temporal scales considered here. We discuss the potential global impact of CRs and anisotropic transport on the dynamical structure of the ISM and also detail the limitations of the model.

Star-forming dwarf galaxies can be seen as the local proxies of the high-redshift building blocks of more massive galaxies according to the current paradigm of the hierarchical galaxy formation. They are low-mass objects, and therefore their rotation speed is very low. Several galaxies are observed to show quite strong magnetic fields. These cases of strong ordered magnetic fields seem to correlate with a high, but not extremely high, star formation rate. We investigate whether these magnetic fields could be generated by the cosmic-ray-driven dynamo. The environment of a dwarf galaxy is unfavourable for the large-scale dynamo action because of the very slow rotation that is required to create the regular component of the magnetic field. We built a 3D global model of a dwarf galaxy that consists of two gravitational components: the stars and the dark-matter halo described by the purely phenomenological profile proposed previously. We solved a system of magnetohydrodynamic (MHD) equations that include an additional cosmic-ray component described by the fluid approximation. We found that the cosmic-ray-driven dynamo can amplify the magnetic field with an exponential growth rate. The $e$-folding time is correlated with the initial rotation speed. The final mean value of the azimuthal flux for our models is of the order of few $\mu$G and the system reaches its equipartition level. The results indicate that the cosmic-ray-driven dynamo is a process that can explain the magnetic fields in dwarf galaxies.

We present simulations of the magnetized interstellar medium (ISM) in models of massive star forming (40 Msun / yr) disk galaxies with high gas surface densities (~100 Msun / pc^2) similar to observed star forming high-redshift disks. We assume that type II supernovae deposit 10 per cent of their energy into the ISM as cosmic rays and neglect the additional deposition of thermal energy or momentum. With a typical Galactic diffusion coefficient for CRs (3e28 cm^2 / s) we demonstrate that this process alone can trigger the local formation of a strong low density galactic wind maintaining vertically open field lines. Driven by the additional pressure gradient of the relativistic fluid the wind speed can exceed 1000 km/s, much higher than the escape velocity of the galaxy. The global mass loading, i.e. the ratio of the gas mass leaving the galactic disk in a wind to the star formation rate becomes of order unity once the system has settled into an equilibrium. We conclude that relativistic particles accelerated in supernova remnants alone provide a natural and efficient mechanism to trigger winds similar to observed mass-loaded galactic winds in high-redshift galaxies. These winds also help explaining the low efficiencies for the conversion of gas into stars in galaxies as well as the early enrichment of the intergalactic medium with metals. This mechanism can be at least of similar importance than the traditionally considered momentum feedback from massive stars and thermal and kinetic feedback from supernova explosions.

We investigate streaming instability using two-fluid approximation (neutral gas and dust) in a quasi-global, unstratified protoplanetary disc, with the help of PIERNIK code. We compare amplification rate of the eigen-mode in numerical simulations, with the corresponding growth resulting from the linear stability analysis of full system of Euler's equation including aerodynamic drag. Following Youdin & Goodman (2005) we show that (1) rapid dust clumping occurs due to the difference in azimuthal velocities of gas and dust, coupled by the drag force, (2) initial density perturbations are amplified by several orders of magnitude. We demonstrate that the multi-fluid extension of the simple and efficient Relaxing TVD scheme, implemented in PIERNIK, leads to results, which are compatible with those obtained with other methods.

We study the possible magnetization of cosmic voids by void galaxies. Recently, observations revealed isolated starforming galaxies within the voids. Furthermore, a major fraction of a voids volume is expected to be filled with magnetic fields of a minimum strength of about $10^{-15}$ G on Mpc scales. We estimate the transport of magnetic energy by cosmic rays (CR) from the void galaxies into the voids. We assume that CRs and winds are able to leave small isolated void galaxies shortly after they assembled, and then propagate within the voids. For a typical void, we estimate the magnetic field strength and volume filling factor depending on its void galaxy population and possible contributions of strong active galactic nuclei (AGN) which border the voids. We argue that the lower limit on the void magnetic field can be recovered, if a small fraction of the magnetic energy contained in the void galaxies or void bordering AGNs is distributed within the voids.

We present three-dimensional global numerical simulations of the cosmic-ray (CR) driven dynamo in barred galaxies. We study the evolution of the interstellar medium of the barred galaxy in the presence of non-axisymmetric component of the potential, i.e., the bar. The magnetohydrodynamical dynamo is driven by CRs, which are continuously supplied to the disk by supernova (SN) remnants. No magnetic field is present at the beginning of simulations but one-tenth of SN explosions is a source of a small-scale randomly oriented dipolar magnetic field. In all models we assume that 10% of 10^51 erg SN kinetic energy output is converted into CR energy. To compare our results directly with the observed properties of galaxies, we construct realistic maps of polarized radio emission. The main result is that the CR-driven dynamo can amplify weak magnetic fields up to a few {\mu}G within a few Gyr in barred galaxies. The obtained e-folding time is equal to 300 Myr and the magnetic field reaches equipartition at time t ~ 4.0 Gyr. Initially, the completely random magnetic field evolves into large-scale structures. An even (quadrupole-type) configuration of the magnetic field with respect to the galactic plane can be observed. Additionally, the modeled magnetic field configuration resembles maps of the polarized intensity observed in barred galaxies. Polarization vectors are distributed along the bar and between spiral arms. Moreover, the drift of magnetic arms with respect to the spiral pattern in the gas density distribution is observed during the entire simulation time.

We present recent developments of global galactic-scale numerical models of the Cosmic Ray (CR) driven dynamo, which was originally proposed by Parker (1992). We conduct a series of direct CR+MHD numerical simulations of the dynamics of the interstellar medium (ISM), composed of gas, magnetic fields and CR components. We take into account CRs accelerated in randomly distributed supernova (SN) remnants, and assume that SNe deposit small-scale, randomly oriented, dipolar magnetic fields into the ISM. The amplification timescale of the large-scale magnetic field resulting from the CR-driven dynamo is comparable to the galactic rotation period. The process efficiently converts small-scale magnetic fields of SN-remnants into galactic-scale magnetic fields. The resulting magnetic field structure resembles the X-shaped magnetic fields observed in edge-on galaxies.

Radio observations show that magnetic fields are present in dwarf irregular galaxies (dIrr) and its strength is comparable to that found in spiral galaxies. Slow rotation, weak shear and shallow gravitational potential are the main features of a typical dIrr galaxy. These conditions of the interstellar medium in a dIrr galaxy seem to unfavourable for amplification of the magnetic field through the dynamo process. Cosmic-ray driven dynamo is one of the galactic dynamo model, which has been successfully tested in case of the spiral galaxies. We investigate this dynamo model in the ISM of a dIrr galaxy. We study its efficiency under the influence of slow rotation, weak shear and shallow gravitational potential. Additionally, the exploding supernovae are parametrised by the frequency of star formation and its modulation, to reproduce bursts and quiescent phases. We found that even slow galactic rotation with a low shearing rate amplifies the magnetic field, and that rapid rotation with a low value of the shear enhances the efficiency of the dynamo. Our simulations have shown that a high amount of magnetic energy leaves the simulation box becoming an efficient source of intergalactic magnetic fields.

We study the global evolution of the magnetic field and interstellar medium (ISM) of the barred and ringed galaxies in the presence of non-axisymmetric components of the potential, i.e. the bar and/or the oval perturbations. The magnetohydrodynamical dynamo is driven by cosmic rays (CR), which are continuously supplied to the disk by supernova (SN) remnants. Additionally, weak, dipolar and randomly oriented magnetic field is injected to the galactic disk during SN explosions. To compare our results directly with the observed properties of galaxies we construct realistic maps of high-frequency polarized radio emission. The main result is that CR driven dynamo can amplify weak magnetic fields up to few $\mu$G within few Gyr in barred and ringed galaxies. What is more, the modelled magnetic field configuration resembles maps of the polarized intensity observed in barred and ringed galaxies.

Irregular galaxies are usually smaller and less massive than their spiral, S0, and elliptical counterparts. Radio observations indicate that a magnetic field is present in irregular galaxies whose value is similar to that in spiral galaxies. However, the conditions in the interstellar medium of an irregular galaxy are unfavorable for amplification of the magnetic field because of the slow rotation and low shearing rate. We investigate the cosmic-ray driven dynamo in the interstellar medium of an irregular galaxy. We study its efficiency under the conditions of slow rotation and weak shear. The star formation is also taken into account in our model and is parametrized by the frequency of explosions and modulations of activity. The numerical model includes a magnetohydrodynamical dynamo driven by cosmic rays that is injected into the interstellar medium by randomly exploding supernovae. In the model, we also include essential elements such as vertical gravity of the disk, differential rotation approximated by the shearing box, and resistivity leading to magnetic reconnection. We find that even slow galactic rotation with a low shearing rate amplifies the magnetic field, and that rapid rotation with a low value of the shear enhances the efficiency of the dynamo. Our simulations have shown that a high amount of magnetic energy leaves the simulation box becoming an efficient source of intergalactic magnetic fields.

We present a new particle module of the magnetohydrodynamic (MHD) Piernik code. The original multi-fluid grid code based on the Relaxing Total Variation Diminishing (RTVD) scheme has been extended by addition of dust described within the particle approximation. The dust is now described as a system of interacting particles. The particles can interact with gas, which is described as a fluid. In this poster we introduce the scheme used to solve equations of motion for the particles and present the first results coming from the module. The results of test problems are also compared with the results coming from fluid simulations made with Piernik-MHD code. The comparison shows the most important differences between fluid and particle approximations used to describe dynamical evolution of dust under astrophysical conditions.

Piernik is a multi-fluid grid magnetohydrodynamic (MHD) code based on the Relaxing Total Variation Diminishing (RTVD) conservative scheme. The original code has been extended by addition of dust described within the particle approximation. The dust is now described as a system of interacting particles. The particles can interact with gas, which is described as a fluid. The comparison between the test problem results and the results coming from fluid simulations made with Piernik code shows the most important differences between fluid and particle approximations used to describe dynamical evolution of dust under astrophysical conditions.

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.

We report first results of first global galactic-scale CR-MHD simulations of cosmic-ray-driven dynamo. We investigate the dynamics of magnetized interstellar medium (ISM), which is dynamically coupled with the cosmic-ray (CR) gas. We assume that exploding stars deposit small-scale, randomly oriented, dipolar magnetic fields into the differentially rotating ISM, together with a portion of cosmic rays, accelerated in supernova shocks. We conduct numerical simulations with the aid of a new parallel MHD code PIERNIK. We find that the initial magnetization of galactic disks by exploding magnetized stars forms a favourable conditions for the cosmic-ray-driven dynamo. We demonstrate that dipolar magnetic fields supplied on small SN-remnant scales, can be amplified exponentially, by the CR-driven dynamo, to the present equipartition values, and transformed simultaneously to large galactic-scales. The resulting magnetic field structure in an evolved galaxy appears spiral in the face-on view and reveals the so called X-shaped structure in the edge-on view.

We investigate the cosmic ray driven dynamo in the interstellar medium of irregular galaxy. The observations (Chyzy et al. 2000, 2003) show that the magnetic field in irregular galaxies is present and its value reaches the same level as in spiral galaxies. However the conditions in the medium of irregular galaxy are very unfavorable for amplification the magnetic field due to slow rotation and low shearing rate. In this work we present numerical model of the interstellar medium in irregular galaxies. The model includes magnetohydrodynamical dynamo driven by cosmic rays in the interstellar medium provided by random supernova explosions. We describe models characterized by different shear and rotation. We find that even slow galactic rotation with low shearing rate gives amplification of the magnetic field. Simulations have shown that high amount of the magnetic energy flow out off the simulation region becoming an efficient source of intergalactic magnetic fields.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin & Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we present an extension of PIERNIK, which is designed for simulations of diffusive propagation of the Cosmic-Ray (CR) component in the magnetized ISM.

We conduct global galactic--scale magnetohydrodynamical (MHD) simulations of the cosmic--ray driven dynamo. We assume that exploding stars deposit small--scale, randomly oriented, dipolar magnetic fields into the differentially rotating ISM, together with a portion of cosmic rays, accelerated in supernova shocks. Our simulations are performed with the aid of a new parallel MHD code PIERNIK. We demonstrate that dipolar magnetic fields supplied on small SN--remnant scales, can be amplified exponentially by the CR--driven dynamo to the present equipartition values, and transformed simultaneously to large galactic--scales by an inverse cascade promoted by resistive processes.

We present new developments on the Cosmic--Ray driven, galactic dynamo, modeled by means of direct, resistive CR--MHD simulations, performed with ZEUS and PIERNIK codes. The dynamo action, leading to the amplification of large--scale galactic magnetic fields on galactic rotation timescales, appears as a result of galactic differential rotation, buoyancy of the cosmic ray component and resistive dissipation of small--scale turbulent magnetic fields. Our new results include demonstration of the global--galactic dynamo action driven by Cosmic Rays supplied in supernova remnants. An essential outcome of the new series of global galactic dynamo models is the equipartition of the gas turbulent energy with magnetic field energy and cosmic ray energy, in saturated states of the dynamo on large galactic scales.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin & Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we present Ohmic resistivity extension of the original Relaxing TVD MHD scheme, and show examples of magnetic reconnection in cases of uniform and current-dependent resistivity prescriptions.

We present a parameter study of the magnetohydrodynamical dynamo driven by cosmic rays in the interstellar medium (ISM) focusing on the efficiency of magnetic field amplification and the issue of energy equipartition between magnetic, kinetic and cosmic ray (CR) energies. We perform numerical CR-MHD simulations of the ISM using the extended version of ZEUS-3D code in the shearing box approximation and taking into account the presence of Ohmic resistivity, tidal forces and vertical disk gravity. CRs are supplied in randomly distributed supernova (SN) remnants and are described by the diffusion-advection equation, which incorporates an anisotropic diffusion tensor. The azimuthal magnetic flux and total magnetic energy are amplified depending on a particular choice of model parameters. We find that the most favorable conditions for magnetic field amplification correspond to magnetic diffusivity of the order of $3\times 10^{25} \cm^2\s^{-1}$, SN rates close to those observed in the Milky Way, periodic SN activity corresponding to spiral arms, and highly anisotropic and field-aligned CR diffusion. The rate of magnetic field amplification is relatively insensitive to the magnitude of SN rates in a rage of spanning 10% up to 100% of realistic values. The timescale of magnetic field amplification in the most favorable conditions is 150 Myr, at galactocentric radius equal to 5 kpc. The final magnetic field energies fluctuate near equipartition with the gas kinetic energy. In all models CR energy exceeds the equipartition values by a least an order of magnitude, in contrary to the expected equipartition. We suggest that the excess of cosmic rays can be attributed to the fact that the shearing-box does not permit cosmic rays to leave the system along the horizontal magnetic field.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin and Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we introduce the multifluid extension of Relaxing TVD scheme and present a test case of dust migration in a two-fluid disk composed of gas and dust. We demonstrate that due to the difference in azimuthal velocities of gas and dust and the drag force acting on both components dust drifts towards maxima of gas pressure distribution.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme. The original scheme has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we shortly introduce basic elements of the Relaxing TVD MHD algorithm, following Trac & Pen (2003) and Pen et al. (2003), and then focus on the conservative implementation of the shearing box model, constructed with the aid of the Masset's (2000) method. We present results of a test example of a formation of a gravitationally bounded object (planet) in a self-gravitating and differentially rotating fluid.

In the present paper we construct maps of polarized synchrotron radio emission of a whole galaxy, based on local models of the cosmic ray (CR) driven dynamo. We perform numerical simulations of the dynamo in local Cartesian domains, with shear-periodic boundary conditions, placed at the different galactocentric radii. Those local solutions are concatenated together to construct the synchrotron images of the whole galaxy. The main aim of the paper is to compare the model results with the observed radio continuum emission from nearly edge-on spiral galaxy. On the basis of the modeled evolution of the magnetic field structure, the polarization maps can be calculated at different time-steps and at any orientation of the modeled galaxy. For the first time a self-consistent cosmic-ray electron distribution is used to integrate synchrotron emissivity along the line of sight. Finally, our maps are convolved with the given radiotelescope beam. We show that it is possible to reconstruct the extended magnetic halo structures of the edge-on galaxies (so called X-shaped structures).

We conduct a series of magnetohydrodynamical (MHD) simulations of magnetized interstellar medium (ISM) disturbed by exploding stars. Each star deposits a randomly oriented, dipolar magnetic field into ISM. The simulations are performed in a Cartesian box, in a reference frame that is corotating with the galactic disk. The medium is stratified by vertical galactic gravity. The resulting turbulent state of ISM magnetized by the stellar explosions is processed with the aid of Fourier analysis. The results leads to the conclusion that the input of magnetic energy from exploding stars is additionally multiplied by differential rotation. The resulting magnetic field appears to grow up in small-scale component, while the total magnetic flux remains limited. Our results indicate that magnetic field originating from exploding stars can be a source of initial magnetic fields for a subsequent dynamo process.

We analyze the electromotive force (EMF) terms and basic assumptions of the linear and nonlinear dynamo theories in our three-dimensional (3D) numerical model of the Parker instability with cosmic rays and shear in a galactic disk. We also apply the well known prescriptions of the EMF obtained by the nonlinear dynamo theory (Blackman & Field 2002 and Kleeorin et al. 2003) to check if the EMF reconstructed from their prescriptions corresponds to the EMF obtained directly from our numerical models. We show that our modeled EMF is fully nonlinear and it is not possible to apply any of the considered nonlinear dynamo approximations due to the fact that the conditions for the scale separation are not fulfilled.

Qualitatively new aspects of the (linear and non-linear) stability of sheared relativistic (slab) jets are analyzed. The linear problem has been solved for a wide range of jet models well inside the ultrarelativistic domain (flow Lorentz factors up to 20; specific internal energies $\approx 60c^2$). As a distinct feature of our work, we have combined the analytical linear approach with high-resolution relativistic hydrodynamical simulations, which has allowed us i) to identify, in the linear regime, resonant modes specific to the relativistic shear layer ii) to confirm the result of the linear analysis with numerical simulations and, iii) more interestingly, to follow the instability development through the non-linear regime. We find that very high-order reflection modes with dominant growth rates can modify the global, long-term stability of the relativistic flow. We discuss the dependence of these resonant modes on the jet flow Lorentz factor and specific internal energy, and on the shear layer thickness. The results could have potential applications in the field of extragalactic relativistic jets.

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.

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.

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.

We present the first numerical model of the magnetohydrodynamical cosmic-ray (CR) driven dynamo of the type proposed by Parker (1992). The driving force of the amplification process comes from CRs injected into the galactic disk in randomly distributed spherical regions representing supernova remnants. The underlying disk is differentially rotating. An explicit resistivity is responsible for the dissipation of the small-scale magnetic field component. We obtain amplification of the large-scale magnetic on a timescale 250 Myr.

We present a numerical algorithm for the incorporation of the active cosmic ray transport, into the ZEUS-3D magnetohydrodynamical code. The cosmic ray transport is described by the diffusion-advection equation. The applied form of the diffusion tensor allows for anisotropic diffusion of cosmic rays along and across the magnetic field direction, which is controlled by two parameters: the parallel and perpendicular diffusion coefficients. The implemented numerical algorithm is tested by comparison of the diffusive transport of cosmic rays to analytical solutions of the diffusion equation. Our method is numerically stable for a wide range of diffusion coefficients, including the realistic values inferred from the observational data for the Milky Way of about $ 6 \times 10^{28} \cm^2 \s^{-1}$. The presented algorithm is applied for exemplary simulations of the the Parker instability triggered by cosmic rays injected by a single SN remnant.

Parker instability leads to the formation of tangential discontinuities in a magnetic field and subsequent magnetic reconnection due to a numerical and/or an explicit resistivity. In this paper we investigate the role of the uniform, localized and numerical resistivity on the diffusion of magnetic field lines during the growth phase of Parker instability modes. We propose a new method to quantify the diffusion of magnetic field lines which is attributed to the presence of resistivity in ideal and non-ideal MHD codes. The method relies (1) on integration of magnetic lines in between periodic boundaries, (2) on measurements of the dispersion of magnetic field lines with the left and the right periodic boundaries and (3) on a statistical analysis of shifts of a large set of magnetic lines. The proposed method makes it possible to detect topological evolution of magnetic field. We perform a series of resistive MHD simulations of the Parker instability in uniformly rotating galactic disks. We follow the topological evolution of the magnetic field evolving due to the Parker instability and relate it to the ratio of total to uniform magnetic field in galactic disks. We find that after the onset of the Parker instability, the magnetic field becomes first tangled and later on it evolves toward a uniform state due to the presence of resistivity. A similar effect of a varying contribution of a turbulent magnetic field is observed in arms and inter-arm regions of galaxies.

We investigate favourable circumstances for fast magnetic reconnection in astrophysical plasmas based on recent results by Rogers et al. (2001). Given that a critical magnetic field structure with antiparallel field lines exists, our analysis demonstrates that a sufficient condition for fast reconnection is that the ratio of the thermal pressure to the magnetic field pressure $\beta$ should be significantly larger than $ 2 m_e/m_p$ (twice the ratio of electron mass to proton mass). Using several examples (like the different components of the interstellar medium, the intergalactic medium, active galactic nuclei and jets) we show that in almost all astrophysical plasmas, magnetic reconnection proceeds fast i.e. independent of the resistivity, with a few percent of the Alfv{\'e}n speed. Only for special cases like neutron stars and white dwarfs is $\beta$ smaller than $2 m_e/ m_p$.

We present a scenario for efficient magnetization of very young galaxies about 0.5 Gigayears after the Big-Bang by a cosmic ray-driven dynamo. These objects experience a phase of strong star formation during this first $10^9$ years. We transfer the knowledge of the connection between star formation and the production rate of cosmic rays by supernova remnants to such high redshift objects. Since the supernova rate is a direct measure for the production rate of cosmic rays we conclude that very young galaxies must be strong sources of cosmic rays. The key argument of our model is the finding that magnetic fields and cosmic rays are dynamically coupled, i.e. a strong cosmic ray source contains strong magnetic fields since the relativistic particles drive an efficient dynamo in a galaxy via their buoyancy. We construct a phenomenological model of a dynamo driven by buoyancy of cosmic rays and show that if azimuthal shearing is strong enough the dynamo amplification timescale is close to the buoyancy timescale of the order of several $10^7 \div 10^8$ yr. We predict that young galaxies are strongly magnetized and may contribute significantly to the gamma-ray-background.

A set of simulations concerning the influence of internal energy on the stability of relativistic jets is presented. Results show that perturbations saturate when the amplitude of the velocity perturbation approaches the speed of light limit. Also, contrary to what predicted by linear stability theory, jets with higher specific internal energy appear to be more stable.

We present a new fast dynamo model for galactic magnetic fields, which is based on the Parker-shearing instability and magnetic reconnection, in the spirit of the model proposed by Parker (1992). We introduce a new scenario of flux tube interactions and estimate the dynamo transport coefficient basing on simple geometrical arguments. The obtained expressions are equivalent to the formally derived helicity $\alpha_d$ and diffusivity $\eta_d$ in the first paper of this series. The model we propose predicts that the $\alpha$-effect in galactic discs has opposite sign with respect to that resulting directly from the sign of the Coriolis force. We estimate the rate of magnetic heating due to the reconnection of magnetic flux tubes, which plays an important role in our dynamo model. The corresponding luminosities of the diffuse X-ray emission are consistent with the ROSAT observations of nearby galaxies. The present considerations synthesize the ideas of Parker with our own results presented in the preceding papers (Hanasz & Lesch 1993, 1997; Hanasz 1997).

In this paper we investigate the idea of Hanasz & Lesch 1993 that the galactic dynamo effect is due to the Parker instability of magnetic flux tubes. In addition to the former approach, we take into account more general physical conditions in this paper, by incorporating cosmic rays and differential forces due to the axisymmetric differential rotation and the density waves as well. We present the theory of slender magnetic flux tube dynamics in the thin flux tube approximation and the Lagrange description. This is the application of the formalism obtained for solar magnetic flux tubes by Spruit (1981), to the galactic conditions. We perform a linear stability analysis for the Parker-shearing instability of magnetic flux tubes in galactic discs and then calculate the dynamo coefficients. We present a number of new effects which are very essential for cosmological and contemporary evolution of galactic magnetic fields. First of all we demonstrate that a very strong dynamo $\alpha$-effect is possible in the limit of weak magnetic fields in presence of cosmic rays. Second, we show that the differential force resulting from axisymmetric differential rotation and the linear density waves causes that the $\alpha$-effect is essentially magnified in galactic arms and switched off in the interarm regions. Moreover, we predict a non-uniform magnetic field in spiral arms and well aligned one in interarm regions. These properties are well confirmed by recent observational results by Beck & Hoernes (1996)

The observed wiggles and knots in astrophysical jets as well as the curvilinear motion of radio emitting features are frequently interpreted as signatures of the Kelvin-Helmholtz (KH) instability (eg. Hardee 1987). We investigate the KH instability of a hydrodynamic jet composed of a relativistic gas, surrounded by a nonrelativistic external medium and moving with a relativistic bulk speed. We show basic nonlinear effects, which become important for a finite amplitude KH modes. Since the KH instability in supersonic jets involves acoustic waves over-reflected on jet boundaries, the basic nonlinear effect relies on the steepening of the acoustic wave fronts, leading to the formation of shocks. It turns our that the shocks appear predominantly in the external nonrelativistic gas, while the internal acoustic waves remain linear for a much longer time. In addition, the external medium "hardens" as soon as the boundary oscillation velocity becomes comparable to the external sound speed. On the other hand, the amplification of internal waves due to the over-reflection is limited by a nonlinearity of the Lorentz $\gamma$ factor. This implies that the sidereal oscillations of the jet boundary, resulting from the K-H instability, are limited to very small amplitudes comparable to a fraction of the jet radius.

In this paper we continue investigations of the Parker-shearing instability performing numerical simulations of the magnetic flux tube dynamics in the thin flux-tube approximation. We show that evolution of flux tubes resulting from numerical simulations is very similar to that of linear solutions if the vertical displacements are smaller than the vertical scale height $H$ of the galactic disc. If the vertical displacements are comparable to $H$, the vertical growth of perturbations is faster in the nonlinear range than in the linear one and we observe a rapid inflation of the flux tube at its top, which leads to a singularity in numerical simulations, if only the cosmic rays are taken into account. Then we perform simulations for the case of nonuniform external medium, which show that the dominating wavelength of the Parker instability is the same as the wavelength of modulations of external medium. As a consequence of this fact, in the case of dominating cosmic ray pressure, the dynamo $\alpha$ effect related to these short wavelength modulations is much more efficient than that related to the linearly most unstable long wavelengths modes of the Parker instability. Under the influence of differential forces resulting from differential rotation and the density waves, the $\alpha$-effect is essentially magnified in the spiral arms and diminished in the interarm regions, what confirms our previous results obtained in the linear approximation.

The recent development of tomographic analysis of radio data shows the presence of sheaths around several extragalactic jets. We have studied the Kelvin-Helmholtz instability of such stratified jets, describing them as inner relativistic beams embedded in slower plasma envelopes. We find that the para- meters of the envelope generally determine the Kelvin-Helmholtz properties of the jet which thus appears isolated from any ambient external medium.

We investigate the Kelvin-Helmholtz instability of stratified jets. The internal component (core) is made of a relativistic gas moving with a relativistic bulk speed. The second component (sheath or envelope) flows between the core and external gas with a nonrelativistic speed. Such a two-component jet describes a variety of possible astrophysical jet configurations like e.g. (1) a relativistic electron-positron beam penetrating a classical electron-proton disc wind or (2) a beam-cocoon structure. We perform a linear stability analysis of such a configuration in the hydrodynamic, plane-parallel, vortex-sheet approximation. The obtained solutions of the dispersion relation show very apparent differences with respect to the single-jet solutions. Due to the reflection of sound waves at the boundary between sheet and external gas, the growth rate as a function of wavenumber presents a specific oscillation pattern. Overdense sheets can slow down the growth rate and contribute to stabilize the configuration. Moreover, we obtain the result that even for relatively small sheet widths the properties of sheet start to dominate the jet dynamics. Such effects could have important astrophysical implications, for instance on the origin of the dichotomy between radio-loud and radio-quiet objects.