Observational studies of collimation in jets in active galactic nuclei (AGN) are a key to understanding their formation and acceleration processes. We have performed an automated search for jet shape transitions in a sample of 367 AGN using VLBA data at 15 GHz and 1.4 GHz. This search has found ten out of 29 nearby jets at redshifts z<0.07 with a transition from a parabolic to conical shape, while the full analyzed sample is dominated by distant AGN with a typical z about 1. The ten AGN are UGC00773, NGC1052, 3C111, 3C120, TXS0815-094, Mrk180, PKS1514+00, NGC6251, 3C371, and BL Lac. We conclude that the geometry transition may be a common effect in AGN jets. It can be observed only when sufficient linear resolution is obtained. Supplementing these results with previously reported shape breaks in the nearby AGN 1H0323+342 and M87, we estimate that the break occurs at 10^5-10^6 gravitational radii from the nucleus. We suggest that the jet shape transition happens when the bulk plasma kinetic energy flux becomes equal to the Poynting energy flux, while the ambient medium pressure is assumed to be governed by Bondi accretion. In general, the break point may not coincide with the Bondi radius. The observational data supports our model predictions on the jet acceleration and properties of the break point.

At present, there is no direct information about the evolution of inclination angle $\chi$ between magnetic and rotational axes in radio pulsars. As to theoretical models of pulsar evolution, they predict both the alignment, i.e. evolution of inclination angle $\chi$ to $0^{\circ}$, and its counter-alignment, i.e. evolution to $90^{\circ}$. In this paper, we demonstrate that the statistics of interpulse pulsars can give us the key test to solve the alignment/counter-alignment problem as the number of orthogonal interpulse pulsars ($\chi \approx 90^{\circ}$) drastically depends on the evolutionary trajectory.

In this first paper from forthcoming series of works devoted to radio image of relativistic jets from active galactic nuclei the role of internal structure of a flow is discussed. We determine the radial profiles of all physical values for reasonable Michel magnetization parameter $\sigma_{\rm M}$ and ambient pressure $P_{\rm ext}$. Maps of Doppler boosting factor $\delta$ and observed directions of linear polarization of synchrotron emission are also constructed.

Thermal effects on the properties of four electromagnetic waves propagating in the pulsar magnetosphere are analyzed. It is shown that thermal effects change only quantitatively the dispersion properties of superluminal ordinary O-mode freely escaping the pulsar magnetosphere; properties of the extraordinary X-mode remain unchanged. As for two subluminal waves propagating along magnetic field lines, for them thermal effects result in essential absorption. However, this attenuation occurs at considerable distances from the neutron star, so there is no doubt in their existence.

We propose a new method of estimating a mass of a super massive black hole residing in the center of an active galaxy. The active galaxy M87 offers a convenient test case for the method due to the existence of a large amount of observational data on the jet and ambient environment properties in the central area of the object. We suggest that the observed transition of a jet boundary shape from a parabolic to a conical form is associated with the flow transiting from the magnetically dominated regime to the energy equipartition between plasma bulk motion and magnetic field. By coupling the unique set of observations available for the jet kinematics, environment and boundary profile with our MHD modelling under assumption on the presence of a dynamically important magnetic field in the M87 jet, we estimate the central black hole mass and spin. The method leads us to believe that the M87 super massive black hole has a mass somewhat larger than typically accepted so far.

A laboratory simulation of astrophysical processes is one of the intensively developed areas of plasma physics. A new series of experiments has been launched recently on the Plasma Focus type facility in NRC Kurchatov Institute. The main goal is to study the mechanisms of the jet stabilization, due to which it can propagate at distances much greater than their transverse dimensions. The experiments with stationary gas filling revealed regimes in which a narrowly collimated plasma jet was formed, the head of which was no wider than several centimeters at jet propagation distances of up to 100 cm. The PF-1000 (IFPiLM, Warsaw, Poland) and KPF-4 (SFTI, Sukhum, Abkhazia) experiments are aimed at creating profiled initial gas distributions to control the conditions of plasma jet propagation in the ambient plasma. Estimations of the dimensionless parameters, i.e. the Mach, Reynolds, and Peclet numbers which were achieved during the experiments, showed that the PF-facilities can be used for the YSO jets modelling. The future experiments, which can allow one to understand the nature of the stable plasma ejections observed in many astrophysical sources, are discussed.

We analyze the torque applied to a rotating magnetized sphere in a vacuum. It is shown that for the correct determination of one of the torque's component the angular momentum of the electromagnetic field within the body should be taken into account.

Although fifty years have passed since the discovery of radio pulsars, there is still no satisfactory understanding of how these amazing objects operate. While there has been significant progress in understanding the basic properties of radio pulsars, there is as yet no consensus on key issues, such as the nature of coherent radio emission or the conversion mechanism of the electromagnetic energy of the pulsar wind into particle energy. In this review, we present the main theoretical results on the magnetosphere of neutron stars. We formulate a number of apparently simple questions, which nevertheless remain unanswered since the very beginning of the field and which must be resolved before any further progress can be made.

We investigate the internal structure of the current sheet in the pulsar wind within force-free and two-fluid MHD approximations. Within the force-free approximation we obtain general asymptotic solution of the Grad-Shafranov equation for quasi-spherical pulsar wind up to the second order in small parameter $\varepsilon = (\Omega r/c)^{-1}$. The solution allows an arbitrary latitudinal structure of the radial magnetic field, including that obtained in the numerical simulations of oblique rotators. It is also shown that the shape of the current sheet does not depend on the latitudinal structure. For the internal region of the current sheet outside the fast magnetosonic surface where the force-free approximation is not valid we use two-fluid MHD approximation. Carrying out calculations in the comoving reference frame we succeed in determining intrinsic electric and magnetic fields of a sheet. It allows us to analyze time-dependent effects which were not investigated up to now. In particular, we estimate the efficiency of the particle acceleration inside the sheet. Finally, after investigating the motion of individual particles in the time-dependent current sheet, we find self-consistently the width of the sheet and its time evolution.

The effect of a mass-loading of the magnetohydrodynamic (MHD) flow in a relativistic jet from active galactic nuclei (AGN) due to $\gamma\gamma\rightarrow e^+e^-$ conversion is considered analytically. We argue that the effects of charge average separation due to specific initial pairs' motion lead to partial magnetic and electric fields screening or enhancement. The effect of the fields screening has not been considered earlier. The pairs with the center-of-mass moving faster or slower than the bulk plasma flow create a surface charge and a current which either screen or enhance both electric and magnetic fields in a pair creation domain. This impacts the bulk flow motion, which either accelerates or decelerates. The pairs with the center-of-mass moving with exactly the drift velocity does not induce the fields disturbance. In this case the flow decelerates due to pure mass-loading. For these different cases the Lorentz factor of the loaded outflow is calculated as a function of loading pairs number density. The effect may be important on sub-parsec to parsec scales due to conversion of TeV jet radiation on soft infra-red to ultra-violet external isotropic photon field. This leads to an outer shell acceleration. The conversion of MeV jet radiation on larger scales may account for the flow deceleration due to pure mass-loading. The proposed mechanism may be a source of internal shocks and instabilities in the pair creation region.

Our previous paper outlined the general aspects of the theory of radio light curve and polarization formation for pulsars. We predicted the one-to-one correspondence between the tilt of the linear polarization position angle of the and the circular polarization. However, some of the radio pulsars indicate a clear deviation from that correlation. In this paper we apply the theory of the radio wave propagation in the pulsar magnetosphere for the analysis of individual effects leading to these deviations. We show that within our theory the circular polarization of a given mode can switch its sign, without the need to introduce a new radiation mode or other effects. Moreover, we show that the generation of different emission modes on different altitudes can explain pulsars, that presumably have the X-O-X light-curve pattern, different from what we predict. General properties of radio emission within our propagation theory are also discussed. In particular, we calculate the intensity patterns for different radiation altitudes and present light curves for different observer viewing angles. In this context we also study the light curves and polarization profiles for pulsars with interpulses. Further, we explain the characteristic width of the position angle curves by introducing the concept of a wide emitting region. Another important feature of radio polarization profiles is the shift of the position angle from the center, which in some cases demonstrates a weak dependence on the observation frequency. Here we demonstrate that propagation effects do not necessarily imply a significant frequency-dependent change of the position angle curve.

At present, there are theoretical models of radio pulsar evolution which predict both the alignment, i.e., evolution of inclination angle $\chi$ between magnetic and rotational axes to $0^{\circ}$, and its counter-alignment, i.e., evolution to $90^{\circ}$. At the same time, both models well describe the pulsar distribution on $P$-$\dot P$ diagram. For this reason, up to now it was impossible to determine the braking mechanisms since it was rather difficult to estimate inclination angle evolution on the basis of observation. In this paper we demonstrate that statistics of interpulse pulsars can give us the key to solve alignment/counter-alignment problem as the number of interpulse pulsars (both, having $\chi \sim 0^{\circ}$ and $\chi \sim 90^{\circ}$) drastically depends on evolution of inclination angle.

Surprisingly, the chronology of nearly 50 years of the pulsar magnetosphere and pulsar wind research is quite similar to the history of our civilization. Using this analogy, I have tried to outline the main results obtained in this field. In addition to my talk, the possibility of particle acceleration due to different processes in the pulsar magnetosphere is discussed in more detail.

Deceleration of relativistic jets from active galactic nuclei (AGNs) detected recently by MOJAVE team is discussed in connection with the interaction of the jet material with the external photon field. Appropriate energy density of the isotropic photon field which is necessary to decelerate jets is determined. It is shown that the disturbances of the electric potential and magnetic surfaces play important role in general dynamics of particle deceleration.

Almost 50 years after radio pulsars were discovered in 1967, our understanding of these objects remains incomplete. On the one hand, within a few years it became clear that neutron star rotation gives rise to the extremely stable sequence of radio pulses, that the kinetic energy of rotation provides the reservoir of energy, and that electromagnetic fields are the braking mechanism. On the other hand, no consensus regarding the mechanism of coherent radio emission or the conversion of electromagnetic energy to particle energy yet exists. In this review, we report on three aspects of pulsar structure that have seen recent progress: the self-consistent theory of the magnetosphere of an oblique magnetic rotator; the location, geometry, and optics of radio emission; and evolution of the angle between spin and magnetic axes. These allow us to take the next step in understanding the physical nature of the pulsar activity.

We investigate the motion of the primary beam outside the light cylinder in the pulsar wind. Inside the light cylinder both primary and secondary plasma move along dipole magnetic field lines where their energies can be arbitrary. But at larger distances the theory predicts quasi-radial motion with the velocity exactly corresponding to the drift velocity which cannot be the same for primary and secondary plasma. Hence, the deceleration of the primary beam is to take place simultaneously resulting in the acceleration of the secondary plasma. We investigate this process in the three-fluid MHD approximation and demonstrate that for most pulsars the energy of the beam remains practically unchanged. Only for young radio pulsars (Crab, Vela) essential deceleration up to the energy of the secondary plasma takes place outside the fast magnetosonic surface $r_{\rm F} \sim (10$-$100) R_{\rm L}$, the energy of secondary plasma itself increasing insufficiently.

The analysis of the frequency dependence of the observed shift of the cores of relativistic jets in active galactic nuclei (AGN) allows us to evaluate the number density of the outflowing plasma n and, hence, the multiplicity parameter, which is a ratio of particle number density to Goldreich-Julian number density. We have obtained the median value for the multiplicity parameter to be 3x10^13 and the median value for the Michel magnetization parameter is 8 from an analysis of 97 sources. Since the magnetization parameter can be interpreted as the maximum possible Lorentz factor of the bulk motion which can be obtained for relativistic magnetohydrodynamic (MHD) flow, this estimate is in agreement with the observed superluminal motion of bright features in AGN jets. Moreover, knowing these key parameters, one can determine the transverse structure of the flow. We show that the poloidal magnetic field and particle number density are much larger in the center of the jet than near the at boundary. The MHD model can also explain the typical observed level of jet acceleration. Finally, casual connectivity of strongly collimated jets is discussed.

The sizes of pulsar radio pulses in the plane of the sky are determined. This is important not only in relation to the possibility of directly resolving the radio pulses spatially, but also for verifying and placing constraints on existing theories of the radio emission. The shape of the pulse radio image and its variation as a function of the pulse phase for the extraordinary (X) and ordinary (O) modes are determined for the simple case of a dipolar magnetic field. Images are obtained for pulsars with various angles between their magnetic and rotational axes, and also for various model parameters, such as the height of the emission, the size of the emission region as a function of the height, and the Lorentz factor of the secondary plasma. It is shown that the size of the resulting image is appreciably smaller at the center of the pulse than at its edges.

We analize the quasi-spherical accretion in the presence of axisymmetric vortex turbulence. It is shown that in this case the turbulence changes mainly the effective gravity potential but not the effective pressure.

Previously unexplored accretion regimes associated with the rotation of accreting matter, namely the perturbations of a quasi-spherical subsonic settling flow and Bondi-Hoyle accretion in the presence of axial rotation, are considered within the framework of ideal hydrodynamics. For subsonic settling accretion, the perturbations are shown to grow rapidly as the gravitating center is approached, so that the flow in the inner regions can no longer be considered as quasi-spherical. For Bondi-Hoyle accretion, a vacuum cylindrical cavity is shown to be formed at large distances from the gravitating center near the flow axis, with the flow velocity outside this cavity being virtually independent of the distance to the rotation axis.

The magnitude of the anomalous torque acting on a rotating magnetized ball in vacuum is specified. Its value is shown to depend on the magnetic field structure inside the body.

It was as early as the 1980s that A V Gurevich and his group proposed a theory to explain the magnetosphere of radio pulsars and the mechanism by which they produce coherent radio emission. The theory has been sharply criticized and is currently rarely mentioned when discussing the observational properties of radio pulsars, even though all the criticisms were in their time disproven in a most thorough and detailed manner. Recent results show even more conclusively that the theory has no internal inconsistencies. New observational data also demonstrate the validity of the basic conclusions of the theory. Based on the latest results on the effects of wave propagation in the magnetosphere of a neuron star, we show that the developed theory does indeed allow quantitative predictions of the evolution of neutron stars and the properties of the observed radio emission.

Using the Grad-Shafranov equation, we consider a new analytical model of the black hole magnetosphere based on the assumption that the magnetic field is radial near the horizon and uniform (cylindrical) in the jet region. Within this model, we have managed to show that the angular velocity of particles near the rotation axis of the black hole can be smaller than one half of the black hole angular velocity. This result is consistent with the latest numerical simulations.

We study the influence of the propagation effects on the mean profiles of radio pulsars using the Kravtsov-Orlov method of the wave propagation in the inhomogeneous media. This approach allows us firstly to include into consideration the transition from geometrical optics to vacuum propagation, the cyclotron absorption, and the wave refraction simultaneously. In addition, arbitrary non-dipole magnetic field configuration, drift motion of plasma particles, and their realistic energy distribution are taken into account. The one-to-one correspondence between the signs of circular polarization and position angle (p.a.) derivative along the profile for both ordinary and extraordinary waves is predicted. Using the numerical integration we now can model the main profiles of radio pulsars. It is shown that standard S-shape form of the p.a. swing can be realized for small enough pair production multiplicity and large enough bulk plasma Lorentz factor only. It is also shown that the value of p.a. maximum derivative, that is often used for determination the angle between magnetic dipole and rotation axis, depends on the plasma parameters and could differ from the rotation vector model (RVM) prediction.

The paper deals with the one possible mechanism of the pulsar radio emission, i.e., with the collective curvature radiation of the relativistic particle stream moving along the curved magnetospheric magnetic field lines. It is shown that the electromagnetic wave containing one cylindrical harmonic exp{is{\phi}} can not be radiated by the curvature radiation mechanism, that corresponds to radiation of a charged particle moving along curved magnetic field lines. The point is that the particle in vacuum radiates the triplex of harmonics (s, s \pm 1), so for the collective curvature radiation the wave polarization is very important and cannot be fixed a priori. For this reason the polarization of real unstable waves must be determined directly from the solution of wave equations for the media. Its electromagnetic properties should be described by the dielectric permittivity tensor \^{\epsilon}({\omega},k,r), that contains the information on the reaction on all possible types of radiation.

The observed shift of the core of the relativistic AGN jets as a function of frequency allows us to evaluate the number density of outflowing plasma and, hence, the magnetization parameter. New analysis of more than 20 sources shows that for most of jets the magnetization parameter ~ 10 - 100. Since the magnetization parameter is the maximum possible bulk Lorentz-factor of the relativistic flow, this estimate is consistent with the observed superluminal motion.

The key properties of the wave propagation theory in the magntosphere of radio pulsars based on the Kravtsov-Orlov equation are presented. It is shown that for radio pulsars with known circular polarization and the swing of the linear polarization position angle one can determine which mode, ordinary or extraordinary one, forms mainly the mean profile of the radio emission. The comparison of the observational data with the theory predictions demonstrates their good agreement.

In this review, analytical results obtained for a wide class of stationary axisymmetric flows in the vicinity of compact astrophysical objects are analyzed, with an emphasis on quantitative predictions for specific sources. Recent years have witnessed a great increase in understanding the formation and properties of astrophysical jets. This is due not only to new observations but also to advances in analytical theory which has produced fairly simple relations, and to what can undoubtedly be called a breakthrough in numerical simulation which has enabled confirmation of theoretical predictions. Of course, we are still very far from fully understanding the physical processes occurring in compact sources. Nevertheless, the progress made raises hopes for near-future test observations that can give insight into the physical processes occurring in active astrophysical objects.

We demonstrate that the 1D cylindrical version of the Grad-Shafranov equation is more rich than classical self-similar ones, and more suitable for the astrophysical jets we observe. In particular, it allows us to describe the central (and, hence, the most energetic) part of the flow. Both relativistic and non-relativistic versions are discussed. It is shown that taking into account the finite pressure of the external media one can determine the magnetic flux within the central core. We found as well that for non-relativistic flows which are magnetically dominated near the origin the solution can be constructed only in the presense of the oblique shock near the base of a jet where the additional heating is to take place.

The goal of this presentation is in paying attention to the 1D cylindrical version of the Grad-Shafranov (GS) equation. In our opinion, this approach is more rich than classical self-similar ones, and more suitable for astrophysical jets we observe. In particular, it allows us describing the central (and, hence, the most energetic) part of the flow.

We explore the possibilities for detecting pulsars that have ceased to radiate in the radio band. We consider two models: the model with hindered particle escape from the pulsar surface (first suggested by Ruderman and Sutherland 1975) and the model with free particle escape (Arons 1981; Mestel 1999). In the model with hindered particle escape, the number of particles that leave the pulsar magnetosphere is small and their radiation cannot be detected with currently available instruments. At the same time, for the free particle escape model, both the number of particles and the radiation intensity are high enough for such pulsars to be detectable with the presently available receivers such as GLAST and AGILE spacecrafts. It is also possible that extinct radio pulsars can be among the unidentified EGRET sources.

We study the statistical distribution of extinct radio pulsars -- isolated neutron stars, that just crossed the death line. An important element that distinguishes our study from other works is a consistent allowance for the evolution of the angle between the magnetic and spin axes. We find the distribution of extinct radio pulsars in spin period for two models: the model with hindered particle escape from the neutron star surface and the model with free particle escape. The total number of extinct radio pulsars is shown to be much smaller than for the model in which the evolution of an inclination angle is disregarded. The decrease in total number of extinct radio pulsars arises from the fact, that account for evolution of the inclination angle allows the transition of a neutron star to the propeller stage appear at much smaller periods (P~5-10 s) than it is commonly assumed.

The aim of this article is to draw attention to the importance of the electric current loss in the energy output of radio pulsars. We remind that even the losses attributed to the magneto-dipole radiation of a pulsar in vacuum can be written as a result of an Ampere force action of the electric currens flowing over the neutron star surface (Michel, 1991, Beskin et al., 1993). It is this force that is responsible for the transfer of angular momentum of a neutron star to an outgoing magneto-dipole wave. If a pulsar is surrounded by plasma, and there is no longitudinal current in its magnetosphere, there is no energy loss (Beskin et al., 1993, Mestel et al., 1999). It is the longitudinal current closing within the pulsar polar cap that exerts the retardation torque acting on the neutron star. This torque can be determined if the structure of longitudinal current is known. Here we remind of the solution by Beskin, Gurevitch & Istomin (1993) and discuss the validity of such an assumption. The behavior of the recently observed "part-time job" pulsar B1931+24 can be naturally explained within the model of current loss while the magneto-dipole model faces difficulties.

The aim of this article is to draw attention to the importance of the electric current loss in the energy output of radio pulsars. We remind that even the losses attributed to the magneto-dipole radiation of a pulsar in vacuum can be written as a result of Ampere force action of the electric currents flowing over the neutron star surface (Michel 1991, Beskin, Gurevich & Istomin 1993). It is this force that is responsible for the transfer of angular momentum of a neutron star to an outgoing magneto-dipole wave. If a pulsar is surrounded by plasma, and there is no longitudinal current in its magnetosphere, there is no energy loss (Beskin, Gurevich & Istomin 1993, Mestel, Panagi & Shibata 1999). It is the longitudinal current closing within the pulsar polar cap that exerts the retardation torque acting on the neutron star. This torque can be determined if the structure of longitudinal current is known. Here we remind of the solution by Beskin, Gurevich & Istomin (1993) and discuss the validity of such an assumption. Finally, it is shown that the behaviour of the recently observed "part-time job" pulsar B1931+24 can be naturally explained within the model of current loss while the magneto-dipole model faces difficulties.

We revisit the case of magneto-rotational evolution of neutron stars (NSs) with accounting for changes of the angle \chi between spin and magnetic axes. This element of the evolution of NSs is very important for age estimates and population modeling, but usually it is neglected. In the framework of two models of energy losses we demonstrate that unless NSs are born with the inclination angle \chi_0 very close (within 1 degree) to the position of maximal losses, pulsars with short initial periods quickly reach regions of small energy losses (aligned rotators in the case of magneto-dipole losses and orthogonal rotators -- in the case of current losses) without significant spin-down. This means that either most of known NSs should be born with relatively long periods close to presently observed values, or the initial inclination angles should be large (close to 90) for the magnetodipole model and very small (close to 0) for the current losses model, or that both models cannot be applied for the whole evolutionary track of a typical NS. In particular, magnetar candidates should be born with periods close to the observed ones within the scope of both models of energy losses. We discuss how these considerations can influence population synthesis models and age estimates for different types of isolated NSs. Focusing on the model of current losses we illustrate our conclusion with evolutionary track reconstructions on the P--\sin \chi plane. We conclude, that most probably both existing models of energy losses -- magneto-dipole and longitudinal current -- require serious modifications, on the other hand, estimates obtained under the standard assumption of \sin \chi=1=const does not have solid theoretical ground and should be taken with care.

The problem of the efficiency of particle acceleration for parabolic poloidal magnetic field is considered within the approach of steady axisymmetric MHD flow. For the large Michel magnetization parameter it is possible to linearize the stream equation near the force-free solution and to solve the problem self-consistently as was done by Beskin, Kuznetsova and Rafikov (1998) for monopole magnetic field. It is shown that on the fast magnetosonic surface the particle Lorentz factor does not exceed the standard value. On the other hand, in the supersonic region the Lorentz factor grows linearly with the distance from the rotational axis, reaching almost its maximal value. This corresponds to almost the full conversion of the Poynting energy flux into the particle kinetic one.

We analyze the statistical distribution of neutron stars at the stage of a supersonic propeller. An important point of our analysis is allowance for the evolution of the angle of inclination of the magnetic axis to the spin axis of the neutron star for the boundary of the transition to the supersonic propeller stage for two models: the model with hindered particle escape from the stellar surface and the model with free particle escape. As a result, we have shown that a consistent allowance for the evolution of the inclination angle in the region of extinct radio pulsars for the two models leads to an increase in the total number of neutron stars at the supersonic propeller stage. This increase stems from he fact that when allowing for the evolution of the inclination angle $\chi$ for neutron stars in the region of extinct radio pulsars and, hence, for the boundary of the transition to the propeller stage, this transition is possible at shorter spin periods (P~5-10 s) than assumed in the standard model.

We study the two-dimensional structure of thin transonic accretion discs in the vicinity of black holes. We use the hydrodynamical version of the Grad-Shafranov equation and focus on the region inside the marginally stable orbit. We show that all components of the dynamical force in the disc become significant near the sonic surface and especially important in the supersonic region. Under certain conditions, the disc structure is shown to be far from radial, and we review the affected disc properties, in particular the role of the critical condition at the sonic surface: it determines neither the accretion rate nor the angular momentum in the accretion disc. Finally, we present a simple model explaining quasi-periodical oscillations that have been observed in the infrared and X-ray radiation of the Galactic Centre.

My lecture is devoted to the analytical results available for a large class of axisymmetric stationary flows in the vicinity of compact astrophysical objects. First, the most general case is formulated corresponding to the axisymmetric stationary MHD flow in the Kerr metric. Then, I discuss the hydrodynamical version of the Grad-Shafranov equation. Although not so well-known as the full MHD one, it allows us to clarify the nontrivial structure of the Grad-Shafranov approach as well as to discuss the simplest version of the 3+1-split language -- the most convenient one for the description of ideal flows in the vicinity of rotating black holes. Finally, I consider several examples that demonstrate how this approach can be used to obtain the quantitative description of the real transonic flows in the vicinity of rotating and moving black holes.

We study the two-dimensional structure of thin transonic accretion discs in the vicinity of a non-spinning black hole within the framework of hydrodynamical version of the Grad-Shafranov equation. Our analysis focuses on the region inside the marginally stable orbit (MSO), r < r_{ms}. We show that all components of the dynamical force in the disc become significant near the sonic surface and (especially) in the supersonic region. Under certain conditions, the disc structure can be far from radial, and we review the affected disc properties, in particular the role of the critical condition at the sonic surface. Finally, we present a simple model aimed at explaining quasi-periodical oscillations that have been observed in the infra-red and X-ray radiation of the Galactic Centre.

The effects of radiation drag force on the structure of relativistic electron-positron and electron-proton outflows are considered within the one-fluid approximation for quasi-monopole cold outflow. It is shown that for a Poynting-dominated case the drag force does not change the particle energy inside the fast magnetosonic surface. In this region the action of the drag results in a diminishing of the Poynting flux, not the particle flux. Outside the fast magnetosonic surface, for intermediate photon density the drag force may result in additional acceleration of the plasma. This acceleration is a result of the disturbance of magnetic surfaces under the action of the drag. At even larger distances particles are not frozen into the magnetic field and the drag force decelerates them efficiently. In the case of extreme photon densities, the disturbance of magnetic surfaces becomes large and the drag force changes the total energy flux significantly, the particles becoming nonrelativistic. We find that for Active Galactic Nuclei the photon density is too low to disturb the parameters of an ideal MHD outflow. The drag action may result in additional acceleration of outgoing plasma only for central engines with very high luminosities. For cosmological gamma-ray bursts the drag force can strongly affect the process of formation of a Poynting-dominated outflow.

The goal of my lecture is to present the introduction into the hydrodynamical version of the Grad-Shafranov equation. Although not so well-known as the full MHD one, it allows us to clarify the nontrivial structure of the Grad-Shafranov approach as well as to discuss the simplest version of the 3+1-split language -- the most convenient one for the description of the ideal flows in the vicinity of a rotating black hole.

The internal structure of the thin transonic disk accreting onto a nonrotating black hole inside the last stable orbit ($r < 3r_{\rm g}$) is considered within the hydrodynamical version of the Grad-Shafranov equation. It is shown that in the vicinity of the sonic surface takes place a sharp diminishing of the disk thickness. As a result, in the vertical balance equation the dynamical forces $\rho[({\bf v}\nabla){\bf v}]_{\theta}$ become important, these on the sonic surface being of the same order as the pressure gradient $\nabla_{\theta} P$. In the supersonic region the thickness of the disk is determined by the form of ballistic trajectories rather than by the pressure gradient.

A magnetohydrodynamic model is constructed for a cylindrical jet immersed in an external uniform magnetic field. It is shown that, as in the force-free case, the total electric current within the jet can be zero. The particle energetics and the magnetic field structure are determined in a self-consistent way; all jet parameters depend on the physical conditions in the external medium. In particular, we show that a region with subsonic flow can exist in the central jet regions. In actual relativistic jets, most of the energy is transferred by the electromagnetic field only when the magnetization parameter is sufficiently large, $\sigma>10^6$. We also show that, in general, the well-known solution with a central core, $B_z = B_0/(1+\varpi^2/\varpi_c^2)$, can not be realized in the presence of an external medium.

The most general form of the nonrelativistic Grad-Shafranov equation describing anisotropic pressure effects is formulated within the double adiabatic approximation. It gives a possibility to analyze quantitatively how the anisotropic pressure affects the 2D structure of the ideal magnetohydrodynamical flows.

The Grad--Shafranov approach to the problem of the structure of the black hole magnetosphere is discussed. For the double transonic flow, the number of boundary conditions in the pair creation region is shown to be sufficient to determine not only the longitudinal electric current, but also the angular velocity of a flow as a solution to a problem. As a result, the energy loss is determined by the physical parameters at the particle creation region rather than the "boundary conditions" at the event horizon.

The two-fluid effects on the radial outflow of relativistic electron-positron plasma are considered. It is shown that for large enough Michel magnetization parameter (1969) and multiplication parameter, the one-fluid MHD approximation remains correct in the whole region where |E| < |B|. In the case in which the longitudinal electric current is smaller than the Goldreich-Julian one, the acceleration of particles is determined near the light surface where |E| = |B|. It is shown that, as in the previously considered cylindrical geometry (Beskin, Gurevich & Istomin 1983), almost all electromagnetic energy is transformed into the energy of particles in the narrow boundary layer. The relative thickness of this layer is shown to be of the order of the inverse of multiplication parameter.

We consider a model of axisymmetric neutron star magnetosphere. In our approach, the current density in the region of open field lines is constant and the returning current flows in a narrow layer along the separatrix. In this case, the stream equation describing the magnetic field structure is linear both in the open and closed regions, the main problem lying in matching the solutions along the separatrix (Okamoto 1974; Lyubarskii 1990). We demonstrate that it is the stability condition on the separatrix that allows us to obtain a unique solution of the problem. In particular, the zero point of magnetic field is shown to locate near the light cylinder. Moreover, the hypothesis of the existence of the nonlinear Ohm's Law (Beskin, Gurevich & Istomin 1983) connecting the potential drop in the pair creation region and the longitudinal electric current flowing in the magnetosphere is confirmed.