Particle module of Piernik MHD code

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.Comment: 4 pages, submitted to the Proceedings of the 17'th Young Scientists' Conference on Astronomy and Space Physics (April 26 - May 1, 2010 Kyiv, Ukraine

Comparison between particle and fluid approximations to dust dynamics

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.Comment: 2 pages, submitted to the Proceedings of the International Conference of Young Astronomers (September 7-13, 2009 Krak\'ow, Poland

Cosmic rays can drive strong outflows from gas-rich high-redshift disk 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.Comment: 5 pages, 5 figures, accepted in ApJL; corrected titl

Global simulations of galactic dynamo driven by cosmic-rays and exploding magnetized stars

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

Cosmic ray driven dynamo in barred and ringed galaxies

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

Global galactic dynamo driven by cosmic-rays and exploding magnetized stars

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.Comment: 11 pages, 4 figure