A two dimensional numerical model in the axisymmetric approximation that describes the flow structure in the magnetosphere of the white dwarf in the EX Hya system has been developed. Results of simulations show that the accretion in EX Hya proceeds via accretion columns, that are not closed and have curtain-like shapes. The thickness of the accretion curtains depends only weakly on the thickness of the accretion disk. This thickness developed in the simulations does not agree with observations. It is concluded that the main reason for the formation of thick accretion curtains in the used model is the assumption that the magnetic field penetrates fully into the plasma of the disk. An analysis based on simple estimates shows that a diamagnetic disk that fully or partially shields the magnetic field of the star may be a more attractive explanation for the observed features of the accretion in EX Hya.

The structures of plasma flows in close binary systems whose accretors have strong intrinsic magnetic fields are studied. A close binary system with the parameters of a typical polar is considered. The results of three dimensional numerical simulations of the material flow from the donor into the accretor Roche lobe are presented. Special attention is given to the flow structure in the vicinity of the inner Lagrangian point, where the accretion flow is formed. The interaction of the accretion flow material from the envelope of the donor with the magnetic field of the accretor results in the formation of a hierarchical structure of the magnetosphere, because less dense areas of the accretion flow are captured by the magnetic field of the white dwarf earlier than more dense regions. Taking into account this kind of magnetosphere structure can affect analysis results and interpretation of the observations.

It is now established that the hot Jupiters have extensive gaseous (ionospheric) envelopes, which expanding far beyond the Roche lobe. The envelopes are weakly bound to the planet and affected by strong influence of stellar wind fluctuations. Also, the hot Jupiters are located close to the parent star and therefore the magnetic field of stellar wind is an important factor, determining the structure their magnetosphere. At the same time, for a typical hot Jupiter, the velocity of stellar wind plasma, flowing around the atmosphere, is close to the Alfv\'en velocity. This should result in stellar wind parameters fluctuations (density, velocity, magnetic field), that can affect the conditions of formation of bow shock waves around a hot Jupiter, i. e. to switch flow from sub-Alfv\'en to super-Alfv\'en mode and back. In this paper, based on the results of three-dimensional numerical MHD modeling, it is confirmed that in the envelope of hot Jupiter, which is in Alfv\'en point vicinity of the stellar wind, both disappearance and appearance of the bow shock wave occures under the action of coronal mass ejection. The paper also shows that this effect can influence the observational manifestations of hot Jupiter, including luminosity in energetic part of the spectrum.

We show that the orbits of exoplanets of the "hot Jupiter" type, as a rule, are located close to the Alf\'{v}en point of the stellar wind of the parent star. At this, many hot Jupiters can be located in the sub-Alf\'{v}en zone in which the magnetic pressure of the stellar wind exceeds its dynamic pressure. Therefore, magnetic field of the wind must play an extremely important role for the flow of the stellar wind around the atmospheres of the hot Jupiters. This factor must be considered both in theoretical models and in the interpretation of observational data. The analysis shows that many typical hot Jupiters should have shock-less intrinsic magnetospheres, which, apparently, do not have counterparts in the Solar System. Such magnetospheres are characterized, primarily, by the absence of the bow shock, while the magnetic barrier (ionopause) is formed by the induced currents in the upper layers of the ionosphere. We confirmed this inference by the three-dimensional numerical simulation of the flow of the parent star stellar wind around the hot Jupiter HD 209458b in which we took into account both proper magnetic field of the planet and magnetic field of the wind.

We consider scenario of merger of two stellar mass black holes surrounded by an accretion disk. Due to emission of gravitational waves, the mass of the central object decreases and accretion disk experiences perturbation. Calculations show that the main consequence of this disturbance is formation of a shock wave propagating from the center of the disk to its periphery. Light curve is computed and duration of the flash is estimated under assumption that the flash terminates when the luminosity returns to the initial value. It is shown that, if the total mass of the merging binary is 55 M_sol (like in the event GW170814), the flash produced by the shock will increase bolometric luminosity of the disk by 4--6 orders of magnitude, up to 10^{45} erg/s (absolute stellar magnitude -23.8^m). With account of the distance to the source (540 Mpc) and for reasonable assumptions on the parameters of the accretion disk, it turns out that the apparent magnitude of the flash at the maximum of the spectral flux density should be 12.8^m -- 14.2^m, while duration of the flash --- few minutes. The main part of the shock radiation flux is emitted in the X-ray and gamma-ray ranges. In the spectral band of the EPIC instrument of the XMM-Newton observatory or the telescope eROSITA of the Spectrum-RG observatory (0.3--10 keV), luminosity will increase by 3--4 orders of magnitude, up to 10^{44} erg/s, corresponding to the apparent stellar magnitude about 17^m. Luminosity is at maximum in the observational band of the IBIS instrument of the INTEGRAL observatory (20 keV -- 10 MeV) and will be 10^{44} -- 10^{45} erg/s, corresponding to the apparent flux 10^{-4} photons per cm^2/s/keV at the wavelength ~100 keV. From the far UV to the longer wavelengths, the brightening is virtually absent.

We suggest a mechanism for generation of radio emission from polars. It is based on the cyclotron radiation of thermal electrons at the background of fluctuating magnetic field. The source of fluctuations is Alfv\'{e}n wave turbulence. Expressions for the radiation spectrum and degree of polarization are obtained. By the example of the polar AM~Her, the fluxes of radio emission from the accretion stream are calculated. Within the framework of proposed emission model, with realistic plasma characteristics, it is possible to reproduce observed radiation fluxes in the VLA bands.

The structure of plasma flows in close binary systems in which one of the components is a rapidly rotating magnetic white dwarf is studied. The main example considered is the AE Aquarii system; the spin period of the white dwarf is about a factor of 1000 shorter than the orbital period, and the magnetic field on the white dwarf surface is of order of 50 MG. The mass transfer in this system was analyzed via numerical solution of the system of MHD equations. These computations show that the magnetic field of the white dwarf does not significantly influence the velocity field of the material in its Roche lobe in the case of laminar flow regime, so that the field does not hinder the formation of a transient disk (ring) surrounding the magnetosphere. However, the efficiency of the energy and angular momentum exchange between the white dwarf and the surrounding material increases considerably with the development of turbulent motions in the matter, resulting in its acceleration at the magnetospheric boundary and further escape from the system at a high rate. The time scales of the transition of the system between the laminar and turbulent modes are close to those of the AE\,Aqr transition between its quiet and active phases.

We present a method that can be used to recover the spectrum of turbulence from observations of optically thin emission lines formed in astrophysical disks. Within this method we analyze how line intensity fluctuations depend on the angular resolution of the instrument, used for the observations. The method allows us to restore the slope of the power spectrum of velocity turbulent pulsations and estimate the upper boundary of the turbulence scale.

A three-dimensional numerical model for an accretion process investigation in the magnetosphere of a white dwarf in magnetic cataclysmic variables is developed. The model assumes that the white dwarf has a dipole magnetic field with its symmetry axis inclined to the rotation axis. The model is based on the equations of modified MHD, that describe the mean flow parameters in the wave MHD turbulence. Diffusion of the magnetic field and radiative heating and cooling are taken into account. The suitability of the model is confirmed by modeling the accretion in a typical intermediate polar. The computations show that a magnetosphere forms around the accretor, with the accretion occurring via columns. The accretion columns have a curtain-like shape, and arc-shaped zones of energy release form on the surface of the white dwarf in the magnetic poles area as a result of the matter infall.

A numerical hydrodynamical model for the evolution of spherically symmetric collapsing clouds, designed for the calculation of the thermal structure of these objects in both the prestellar and protostellar stages of their evolution, is presented. Distinctive features of the model include the possibility of independently describing the temperatures of the gas and dust, which is extremely important when calculating the thermal structure of prestellar and protostellar clouds, and the account of the radiation flux from the central protostar. This model is used to compare the theoretical density and temperature distributions with observations for nearby sites of star formation obtained with the Herschel Space Observatory. Application of the diffusion approximation with a flux limiter describes well the radial density and temperature distributions in protostellar clouds. However, significant differences between the model and observational density profiles were found for prestellar stages, suggesting the presence of appreciable deviations from equilibrium in the prestellar clouds. An approximate method for calculating the thermal structure of a cloud based on the adaptive $\tau$-approximation is presented. Application of the $\tau$-approximation yields good agreement with the diffusion approximation for the prestellar phase, but produces appreciable discrepancies for the protostellar phase, when the thermal structure of the accreting envelope is determined by the radiation of the protostar.

Context. In binary star systems, the winds from the two components impact each other, leading to strong shocks and regions of enhanced density and temperature. Potentially habitable circumbinary planets must continually be exposed to these interactions regions. Aims. We study, for the first time, the interactions between winds from low-mass stars in a binary system, to show the wind conditions seen by potentially habitable circumbinary planets. Methods. We use the advanced 3D numerical hydrodynamic code Nurgush to model the wind interactions of two identical winds from two solar mass stars with circular orbits and a binary separation of 0.5 AU. As input into this model, we use a 1D hydrodynamic simulation of the solar wind, run using the Versatile Advection Code. We derive the locations of stable and habitable orbits in this system to explore what wind conditions potentially habitable planets will be exposed to during their orbits. Results. Our wind interaction simulations result in the formation of two strong shock waves separated by a region of enhanced density and temperature. The wind-wind interaction region has a spiral shape due to Coriolis forces generated by the orbital motions of the two stars. The stable and habitable zone in this system extends from approximately 1.4 AU to 2.4 AU. (TRUNCATED)

The flow structure in magnetic cataclysmic variables was investigated taking into account the effects of strong magnetic field and fast rotation of the white dwarf. We modeled the AE Aqr system as a unique object that has the rotation period of the white dwarf is about 1000 times shorter than the orbital period of the binary system. Observations show that in spite of fast rotation of the white dwarf some part of the stream from the inner Lagrange point comes into the Roche lobe region. We analyzed possible mechanisms preventing material to outflow from the system.

We performed 3D MHD calculations of stream accretion in cataclysmic variable stars for which the white dwarf primary star possesses a strong and complex magnetic field. These calculations are motivated by observations of polars; cataclysmic variables containing white dwarfs with magnetic fields sufficiently strong to prevent the formation of an accretion disk. So an accretion stream flows from the L1 point and impacts directly onto one or more spots on the surface of the white dwarf. Observations indicate that the white dwarf, in some binaries, possesses a complex (non-dipolar) magnetic field. We perform simulations of 10 polars or equivalently one asynchronous polar at 10 different beat phases. Our models have an aligned dipole plus quadrupole magnetic field centered on the white dwarf primary. We find that for a sufficiently strong quadrupole component an accretion spot occurs near the magnetic equator for slightly less than half of our simulations while a polar accretion zone is active for most of the rest of the simulations. For one or two configurations; accretion at the dominant polar region and at an equatorial zone occurs simultaneously. These are the first 3D MHD calculations to confirm the existence of complex magnetic fields in magnetic CVs. We conclude, that it might be difficult observationally determine if the field is a pure dipole or if it is complex for polars, but there will be indications for some systems. Specifically, a complex magnetic field should be considered if the there is an accretion zone near the white dwarf's spin equator (orbital plane) or if there are two or more accretion regions that cannot be fit by a dipole magnetic field. For asynchronous polars, magnetic field constraints are expected to be substantially stronger, with clearer indicators of complex field geometry due to changes in accretion flow structure as a function of spin-orbit beat phase.

We have developed a three-dimensional numerical model and applied it to simulate plasma flows in semi-detached binary systems whose accretor possesses a strong intrinsic magnetic field. The model is based on the assumption that the plasma dynamics are determined by the slow mean flow, which forms a backdrop for the rapid propagation of MHD waves. The equations describing the slow motion of matter were obtained by averaging over rapidly propagating pulsations. The numerical model includes the diffusion of magnetic field by current dissipation in turbulent vortices, magnetic buoyancy, and wave MHD turbulence. A modified three-dimensional, parallel, numerical code was used to simulate the flow structure in close binary systems with various accretor magnetic fields, from $10^5$ to $10^8$ G. The conditions for the formation of the accretion disk and the criteria distinguishing the two types of flow corresponding to intermediate polars and polars are discussed.

The results of three-dimensional MHD numerical simulations are used to investigate the characteristic properties of the magnetic-field structures in the accretion disks of semi-detached binary systems. It is assumed that the intrinsic magnetic field of the accretor star is dipolar. Turbulent diffusion of the magnetic field in the disk is taken into account. The SS Cyg system is considered as an example. The results of the numerical simulations show the intense generation of a predominantly toroidal magnetic field in the accretion disk. Magnetic zones with well defined structures for the toroidal magnetic field form in the disk, which are separated by current sheets in which there is magnetic reconnection and current dissipation. Possible observational manifestations of such structures are discussed. It is shown that the interaction of a spiral precessional wave with the accretor's magnetosphere could lead to quasi-periodic oscillations of the accretion rate.

A physical model and two-dimensional numerical method for computing the evolution and spectra of protostellar clouds are described. The physical model is based on a system of magneto-gasdynamical equations, including ohmic and ambipolar diffusion, and a scheme for calculating the thermal and ionization structure of a cloud. The dust and gas temperatures are determined during the calculations of the thermal structure of the cloud. The results of computing the dynamical and thermal structure of the cloud are used to model the radiative transfer in continuum and in molecular lines. We presented the results for clouds in hydrostatic and thermal equilibrium. The evolution of a rotating magnetic protostellar cloud starting from a quasi-static state is also considered. Spectral maps for optically thick lines of linear molecules are analyzed. We have shown that the influence of the magnetic field and rotation can lead to a redistribution of angular momentum in the cloud and the formation of a characteristic rotational velocity structure. As a result, the distribution of the velocity centroid of the molecular lines can acquire an hourglass shape. We plan to use the developed program package together with a model for the chemical evolution to interpret and model observed starless and protostellar cores.

We investigate the formation and evolution of isothermal collapse nonuniformity for rotating magnetic interstellar clouds. The initial and boundary conditions correspond to the statement of the problem of homogeneous cloud contraction from a pressure equilibrium with the external medium. The initial uniform magnetic field is collinear with the angular velocity. Fast and slow magnetosonic rarefaction waves are shown to be formed and propagate from the boundary of the cloud toward its center in the early collapse stages. The front of the fast rarefaction wave divides the gas mass into two parts. The density, angular velocity, and magnetic field remain uniform in the inner region and have nonuniform profiles in the outer region. The rarefaction wave front surface can take both prolate and oblate shapes along the rotation axis, depending on the relationship between the initial angular velocity and magnetic field. We derive a criterion that separates the two regimes of rarefaction wave dynamics with the dominant role of electromagnetic and centrifugal forces. Based on analytical estimations and numerical calculations, we discuss possible scenarios for the evolution of collapse nonuniformity for rotating magnetic interstellar clouds.

Formation and evolution of general relativistic collapse inhomogeneity is investigated. It is shown that a rarefaction wave forms at initial collapse stages and propagates inside from the cloud boundary to its center. The focusing time of the rarefaction wave is found. In massive clouds, the rarefaction wave focusing time equals the free-fall time. The collapse of such clouds leads to black hole formation. In low-mass clouds, the rarefaction wave focusing time is less than the free-fall time. After focusing, the collapse of such clouds becomes fully inhomogeneous and can be sufficiently decelerated by the pressure gradient.

The explicit quasi-monotonic conservative TVD scheme and numerical method for the solution of the gravitational MHD equations are developed. The 2D numerical code for the simulation of multidimensional selfgravitating MHD flows on the Eulerian cylindrical grid is constructed. The results of test calculations show that the code has a good mathematical and computational properties and can be applied to the solution of a wide class of plasma physics and astrophysics problems. We have simulated, for example, a collapse of magnetized rotating protostellar clouds.