The results of MHD numerical simulations of the formation and development of magnetized jets are presented. Similarity criteria for comparisons of the results of laboratory laser experiments and numerical simulations of astrophysical jets are discussed. The results of laboratory simulations of jets generated in experiments at the Neodim laser installation are presented.

Many pulsars propagate through the interstellar medium (ISM) with supersonic velocities, and their pulsar winds interact with the interstellar medium (ISM), forming bow shocks and magnetotails (PWN). We model the propagation of pulsars through the inhomogeneous ISM using non-relativistic axisymmetric magneto-hydrodynamic (MHD) simulations. We take into account the wind from the star, and the azimuthal and poloidal components of the magnetic field, and investigate the PWN at different levels of magnetization (the ratio of magnetic to matter energy-densities) in the wind. We consider the interaction of PWN with small-scale and large-scale imhomogeneities in the ISM at different values of magnetization. We conclude that the inhomogeneities in the ISM can change the shapes of the bow shocks and magnetotails at different values of the magnetization.We compare the results of our simulations with the images of the Guitar Nebula and other PWN that show irregularities in the shapes of their bow shocks and magnetotails. We conclude that these irregularities may be caused by the interaction of PWN with the inhomogeneities in the ISM.

Axisymmetric MHD simulations are used to investigate the Bondi-Hoyle accretion onto an isolated magnetized neutron star moving supersonically (with Mach number of 3) through the interstellar medium. The star is assumed to have a dipole magnetic field aligned with its motion and a magnetospheric radius R_m less then the accretion radius R_BH, so that the gravitational focusing is important. We find that the accretion rate to a magnetized star is smaller than that to a non-magnetized star for the parameters considered. Close to the star the accreting matter falls to the star's surface along the magnetic poles with a larger mass flow to the leeward pole of the star. In the case of a relatively large stellar magnetic field, the star's magnetic field is stretched in the direction of the matter flow outside of R_m (towards the windward side of the star). For weaker magnetic fields we observed oscillations of the closed magnetosphere frontward and backward. These are accompanied by strong oscillations of the mass accretion rate which varies by factors ~ 3. Old slowly rotating neutron stars with no radio emission may be visible in the X-ray band due to accretion of interstellar matter. In general, the star's velocity, magnetic moment, and angular velocity vectors may all be in different directions so that the accretion luminosity will be modulated at the star's rotation rate.

We investigate spherical accretion to a rotating magnetized star in the "propeller" regime using axisymmetric resistive magnetohydrodynamic simulations. The regime is predicted to occur if the magnetospheric radius is larger than the corotation radius and smaller than the light cylinder radius. The simulations show that accreting matter is expelled from the equatorial region of the magnetosphere and that it moves away from the star in a supersonic, disk-shaped outflow. At larger radial distances the outflow slows down and becomes subsonic. The equatorial matter outflow is initially driven by the centrifugal force, but at larger distances the pressure gradient force becomes significant. We find the fraction of the Bondi accretion rate which accretes to the surface of the star.

We use axisymmetric magnetohydrodynamic simulations to investigate the spinning-down of magnetars rotating in the propeller regime and moving supersonically through the interstellar medium. The simulations indicate that magnetars spin-down rapidly due to this interaction, faster than for the case of a non-moving star. From many simulation runs we have derived an approximate scaling laws for the angular momentum loss rate, \dot{L} \propto \~\eta_m^{0.3}\mu^{0.6}\rho^{0.8}{\cal M}^{-0.4} \Omega_*^{1.5}, where \rho is the density of the interstellar medium, \cal M is Mach number, \mu is the star's magnetic moment, \Omega_* is its angular velocity, and \eta_m is magnetic diffusivity. A magnetar with a surface magnetic field of 10^{13} - 10^{15} G is found to spin-down to a period P > 10^5-10^6 s in \sim 10^4 - 10^5 years. There is however uncertainty about the value of the magnetic diffusivity so that the time-scale may be longer. We discuss this model in respect of Soft Gamma Repeaters (SGRs) and the isolated neutron star candidate RXJ1856.5-3754.

This work investigates Bondi accretion to a rotating magnetized star in the "propeller" regime using axisymmetric resistive, magnetohydrodynamic simulations. In this regime accreting matter tends to be expelled from the equatorial region of the magnetosphere where the centrifugal force on matter rotating with the star exceeds the gravitational force. The regime is predicted to occur if the magnetospheric radius larger than the corotation radius and less than the light cylinder radius. The simulations show that accreting matters is expelled from the equatorial region of the magnetosphere and that it moves away from the star in a supersonic, disk-shaped outflow. At larger radial distances the outflow slows down and becomes subsonic. The equatorial matter outflow is initially driven by the centrifugal force, but at larger distances the pressure gradient becomes significant. We find that the star is spun-down mainly by the magnetic torques at its surface with the rate of loss of angular momentum $\dot{L}$ proportional to $-\Omega_*^{1.3}\mu^{0.8}$, where $\Omega_*$ is the star's rotation rate and $\mu$ is its magnetic moment. Further, we find that $\dot{L}$ is approximately independent of the magnetic diffusivity of the plasma $\eta_m$. The fraction of the Bondi accretion rate which accretes to the surface of the star is found to be $\propto\Omega_*^{-1.0}\mu^{-1.7}\eta_m^{0.4}$. Predictions of this work are important for the observability of isolated old neutron stars and for wind fed pulsars in X-ray binaries.

This work investigates the propagation of magnetized, isolated old neutron stars through the interstellar medium. We performed axisymmetric non relativistic magnetohydrodynamic (MHD) simulations of the propagation of a non-rotating star with dipole magnetic field aligned with its velocity through the interstellar medium (ISM). Effects of rotation will be discussed in a subsequent work. We consider two cases: (1) where theaccretion radius $R_{acc}$ is comparable to the magnetic standoff distance or Alfv\'en radius $R_A$ and gravitational focusing is important; and (2) where $R_{acc}<< R_A$ and the magnetized star interacts with the ISM as a ``magnetic plow'', without significant gravitational focusing. For the first case simulations were done at a low Mach number ${\cal M}=3$ for a range of values of the magnetic field $B_*$. For the second case, simulations were done for higher Mach numbers, ${\cal M}=10,$ 30, and 50. In both cases, the magnetosphere of the star represents an obstacle for the flow, and a shock wavestands in front of the star. Magnetic field lines are stretched downwind from the star and form a hollow elongated magnetotail. Reconnection of the magnetic field is observed in the tail which may lead to acceleration of particles. Similar powers are estimated to be released in the bow shock wave and in the magnetotail. The estimated powers are, however, below present detection limits. Results of our simulations maybe applied to other strongly magnetized stars, for example, white dwarfs and magnetic Ap stars. Future more sensitive observations may reveal long magnetotails of magnetized stars moving through the ISM.

Quasi-spherical supersonic (Bondi-type) accretion to a star with a dipole magnetic field is investigated using resistive magnetohydrodynamic simulations. A systematic study is made of accretion to a non-rotating star, while sample results for a rotating star are also presented. A new stationary subsonic accretion flow is found with a steady rate of accretion to the magnetized star smaller than the Bondi accretion rate. Dependences of the accretion rate and the flow pattern on the magnetic momentum of the star and the magnetic diffusivity are presented. For slow star's rotation the accretion flow is similar to that in non-rotating case, but in the case of fast rotation the structure of the subsonic accretion flow is fundamentally different and includes a region of ``propeller'' outflow. The methods and results described here are of general interest and can be applied to systems where matter accretes with low angular momentum.