We performed axisymmetric, grid-based, ideal magnetohydrodynamic (MHD) simulations of oscillating cusp-filling tori orbiting a non-rotating neutron star. A pseudo-Newtonian potential was used to construct the constant angular momentum tori in equilibrium. The inner edge of the torus is terminated by a "cusp" in the effective potential. The initial motion of the model tori was perturbed with uniform sub-sonic vertical and diagonal velocity fields. As the configuration evolved in time, we measured the mass accretion rate on the neutron star surface and obtained the power spectrum. The prominent mode of oscillation in the cusp torus is the radial epicyclic mode. It would appear that vertical oscillations are suppressed by the presence of the cusp. From our analysis it follows that the mass accretion rate carries a modulation imprint of the oscillating torus, and hence so does the boundary layer luminosity.

We perform numerical simulations in the close vicinity of a slowly rotating young stellar object. Using our own resistive MHD Zeus347 code in 2D axisymmetry, magnetospheric interaction experiences four robust stages in evolution. Quasi-stationary axial and conical streams of outflowing matter can last many orbital periods as results of the resistivity-facilitated magnetic reconnections. The shape of the magnetic field depends on resistivity in the magnetosphere.

Stellar magnetosphere and accretion disk interact, and a result should be outflow launched from the innermost vicinity of a protostellar object. We simulated physical conditions in this region by resistive MHD simulations. Outflows resembling the observed ones do not happen in the closest vicinity, except for quasi-stationary funnel flows onto the star, but could occur at few tens of stellar radii above the star. Numerical simulations we performed using our resistive version of ZEUS-3D code, ZEUS347.

Axisymmetric resistive MHD simulations for radially self-similar initial conditions are performed, using the NIRVANA code. The magnetic diffusivity could occur in outflows above an accretion disk, being transferred from the underlying disk into the disk corona by MHD turbulence (anomalous turbulent diffusivity), or as a result of ambipolar diffusion in partially ionized flows. We introduce, in addition to the classical magnetic Reynolds number Rm, which measures the importance of resistive effects in the induction equation, a new number Rb, which measures the importance of the resistive effects in the energy equation. We find two distinct regimes of solutions in our simulations. One is the low-resistivity regime, in which results do not differ much from ideal-MHD solutions. In the high-resistivity regime, results seem to show some periodicity in time-evolution, and depart significantly from the ideal-MHD case. Whether this departure is caused by numerical or physical reasons is of considerable interest for numerical simulations and theory of astrophysical outflows and is currently investigated.

Numerical simulations with self-similar initial and boundary conditions provide a link between theoretical and numerical investigations of jet dynamics. We perform axisymmetric resistive magnetohydrodynamic (MHD) simulations for a generalised solution of the Blandford & Payne type, and compare them with the corresponding analytical and numerical ideal-MHD solutions. We disentangle the effects of the numerical and physical diffusivity. The latter could occur in outflows above an accretion disk, being transferred from the underlying disk into the disk corona by MHD turbulence (anomalous turbulent diffusivity), or as a result of ambipolar diffusion in partially ionized flows. We conclude that while the classical magnetic Reynolds number $R_{\rm m}$ measures the importance of resistive effects in the induction equation, a new introduced number, $\rbeta=(\beta/2)R_{\rm m}$ with $\beta$ the plasma beta, measures the importance of the resistive effects in the energy equation. Thus, in magnetised jets with $\beta<2$, when $\rbeta \la 1$ resistive effects are non-negligible and affect mostly the energy equation. The presented simulations indeed show that for a range of magnetic diffusivities corresponding to $\rbeta \ga 1$ the flow remains close to the ideal-MHD self-similar solution.

Turbulent accretion disks are the most probable origin of jets from young stellar objects. Using an extended version of ZEUS-3D code in the axisymmetry option we investigate the evolution of a disk wind into a collimated jet. Resistive jets propagate slower with increasing magnetic diffusivity, and the degree of collimation decreases. The resistive disk included in the simulations ejects persistent, well collimated jet.

We investigate the evolution of a disk wind into a collimated jet under the influence of magnetic diffusivity. Using the ZEUS-3D code in the axisymmetry option we solve the time-dependent resistive MHD equations for a model setup of a central star surrounded by an accretion disk. The disk is taken as a time-independent boundary condition for the mass flow rate and the magnetic flux distribution. We find that the diffusive jets propagate slower into the ambient medium, most probably due to the lower mass flow rate in axial direction. Close to the star we find that a quasi stationary state evolves after several hundreds (weak diffusion) or thousands (strong diffusion) of disk rotations. Magnetic diffusivity affects the protostellar jet structure de-collimating it. We explain these effects in the framework of the Lorentz force.

Accretion powered pulsars exhibit a variety of lightcurves. In this paper we propose to classify the observed lightcurves as single or double pulsed.We analyze the lightcurves of 86 accretion powered pulsars and assign them to these classes. We present three datasets: first in which the classification can be easily done, second for which the classification is more difficult and not certain, and third for which we were unable to classify the pulsar because of lack of published data. We analyze a simple model in which the angles between the magnetic and rotation axis $\beta$, and between the rotation axis and the line of sight $\theta$ are random, and show that it is inconsistent with the data. We also present a model in whichthe angle between the magnetic axis and rotation axis is restricted and compare it with the data. This leads to an upper limit on the angle $\beta < 40^\circ$. We conclude that there must be a mechanism that leads to alignment of the magnetic and spin axis in X-ray pulsars.

We investigate the evolution of a disk wind into a collimated jet under the influence of magnetic diffusivity, assuming that the turbulent pattern in the disk will also enter the disk corona and the jet. Using the ZEUS-3D code in the axisymmetry option we solve the time-dependent resistive MHD equations for a model setup of a central star surrounded by an accretion disk. We find that the diffusive jets propagate slower into the ambient medium. Close to the star we find that a quasi stationary state evolves after several hundred (weak diffusion) or thousand (strong diffusion) disk rotations. Magnetic diffusivity affects the protostellar jet structure as follows. The jet poloidal magnetic field becomes de-collimated. The jet velocity increases with increasing diffusivity, while the degree of collimation for the hydrodynamic flow remains more or less the same. We suggest that the mass flux is a proper tracer for the degree of jet collimation and find indications of a critical value for the magnetic diffusivity above which the jet collimation is only weak.

Flow of matter onto strongly magnetized neutron stars in X-ray binaries proceeds through accretion funnels that roughly follow geometry of the magnetic field. X-rays originate near surface of the neutron star, and it may happen that the accretion flow passes through the line of sight as the star rotates. We consider the effects of such accretion flow eclipses on the X-ray light curves of accretion powered pulsars, and present a set of X-ray light-curves measured by BATSE for A0535+262 for which this phenomenon is very likely to take place.