We present the results of 3D-hydrodynamical simulations of accretion flow in the eclipsing dwarf nova V1239 Her in quiescence. The model includes the optical star filling its Roche lobe, a gas stream emanating from the inner Lagrangian point of the binary system, and the accretion disc structure. A cold hydrogen gas stream is initially emitted towards a point-like gravitational centre. A stationary accretion disc is formed in about 15 orbital periods after the beginning of accretion. The model takes into account partial ionization of hydrogen and uses realistic cooling function for hydrogen. The light curve of the system is calculated as the volume emission of optically thin layers along the line of sight up to the optical depth \tau =2/3 calculated using Planck-averaged opacities. The calculated eclipse light curves show good agreement with observations, with the changing shape of pre-eclipse and post-eclipse light curves being explained entirely due to ~ 50% variations in the mass accretion rate through the gas stream.

The results of three-dimensional numerical simulations of mass transfer in semidetached binaries after mass transfer termination is presented and the structure of residual accretion disk was investigated. It was shown that after the mass transfer termination the quasi-elliptic accretion disk becomes circular. It was also shown that the cessation of mass transfer results in changing of the structure of the accretion disk: the second arm of spiral shock is appeared as well as a new dense formation (blob), the latter moving through the disk with variable velocity. The blob doesn't smear out under the action of dissipation but is sustained by interaction with arms of spiral shock for practically all the lifetime of the disk. We also analyze the dependence of disk's lifetime on the value of viscosity. For the value of parameter $\alpha$ that is typical for observing accretion disk ($\alpha\sim0.01$) the lifetime of residual disk is found to be equal to 50 orbital periods.

We present synthetic Doppler maps of gaseous flows in semidetached binaries based on the results of 3D gas dynamical simulations. Using of gas dynamical calculations alongside with Doppler tomography technique permits to identify main features of the flow on the Doppler maps without solution of ill-posed inverse problem. Comparison of synthetic tomograms with observations makes possible both to refine the gas dynamical model and to interpret the observational data.

We present the results of three-dimensional numerical simulations of mass transfer in semi-detached binary with a donor star whose rotation vector precesses around the orbital rotation axis of the binary in the observer's coordinate frame. The calculations support our previous model of flow without a `hot spot'. Characteristic features of the flow in this model, such as the formation of an circumbinary envelope, the absence of a `hot spot' at the edge of the accretion disk, and the formation of a shock wave along the edge of the stream, are also present in the solution for a binary with precessing donor star. The parameters of accretion disk and of the structure of the near-disk regions recur with the precessional period of the rotation axis of donor star.

Results of 3D gas dynamical simulation of mass transfer in binaries are presented for systems with various values of viscosity. Analysis of obtained solutions shows that in the systems with low value of viscosity the flow structure is qualitatively similar to one for systems with high viscosity. Presented calculations confirm that there is no shock interaction between the stream from L1 and the forming accretion disk (`hot spot') at any value of viscosity.

We present the results of systematic comparison of 2D and 3D numerical models of mass transfer in semidetached binaries. It is shown that only for case of \gamma \sim 1 (near-isothermal case) the obtained 2D and 3D solutions are qualitatively similar. For higher value of \gamma the 3D flow structure is drastically changed and for \gamma=1.2 the accretion disk is not formed in the system. Numerical results show that for case of \gamma=1.2 the 2D and 3D solutions are different. We discuss the spiral-shaped shock waves obtained in numerical models. It is shown that these shocks are not intrinsic spirals and are caused by collisions of the gaseous flows in the system.

We present the results of three-dimensional numerical simulations of mass transfer in semi-detached binary systems in which the mass-losing star is rotating. The cases of aligned and misaligned non-synchronous rotation of the donor star are considered; the resulting flow patterns are compared to the synchronous case. The main properties of the flow, such as the formation of an circumbinary envelope, the absence of a "hot spot" on the edge of the accretion disk, and the formation of a shock wave along the flow edge, are qualitatively similar to those obtained earlier. For the case of misaligned, non-synchronous rotation, the behavior of the disk and surrounding matter in established flow regime reflects changes in the boundary conditions at the surface of the donor star; in other words, a "driven disk" model is realized in the calculations.

We consider the dependence of the internal structure of a neutron star in a close binary system on the semi-major axis of the binary orbit, focusing on the case when the Roche lobes of the components are nearly filled. We adopt a polytropic equation of state. The temporal evolution of the semi-major axis and its dependence on the mass ratio of the binary components and the polytropic index are determined. The calculation are carried out right up to the moment of contact, when quasi-stationary model becomes invalid. We analyze differences in the shapes of the pulses of gravitational radiation emitted by a pair of point masses and by a binary neutron star, taking into account its internal structure and tidal deformations.

Stationary magnetohydrodynamic (MHD) outflows from a rotating, conducting Keplerian accretion disk threaded by B-field are investigated numerically by time-dependent, axisymmetric (2.5D) simulations using a Godunov-type code. A large class of stationary magneto-centrifugally driven winds are found where matter is accelerated from a thermal speed at the disk to much larger velocity, greater than the fast magnetosonic speed and larger than the escape speed. The flows are approximately spherical outflows with only small collimation within the simulation region. Numerical results are shown to coincide with the theoretical predictions of ideal, axisymmetric MHD to high accuracy. Investigation of the influence of outer boundary conditions, particularly that on the toroidal component of magnetic field shows that the commonly used ``free'' boundary condition leads to artificial magnetic forces which can act to give spurious collimation. New boundary conditions are proposed which do not generate artificial forces. Artificial results may also arise for cases where the Mach cones on the outer boundaries are partially directed into the simulation region.

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.

Hydrodynamic simulations have been used to study accretion disks consisting of counterrotating components with an intervening shear layer(s). Configurations of this type can arise from the accretion of newly supplied counterrotating matter onto an existing corotating disk. The grid-dependent numerical viscosity of our hydro code is used to simulate the influence of a turbulent viscosity of the disk. Firstly, we consider the case where the gas well above the disk midplane rotates with angular rate +\Omega(r) and that well below has the same properties but rotates with rate -\Omega(r). We find that there is angular momentum annihilation in a narrow equatorial boundary layer in which matter accretes supersonically with a velocity which approaches the free-fall velocity and the average accretion speed of the disk can be enormously larger than that for a conventional \alpha-disk rotating in one direction. Secondly, we consider the case of a corotating accretion disk for r<r_t and a counterrotating disk for r>r_t. In this case we observed, that matter from the annihilation layer lost its stability and propagated inward pushing matter of inner regions of the disk to accrete. Thirdly, we investigated the case where counterrotating matter inflowing from large radial distances encounters an existing corotating disk. Friction between the inflowing matter and the existing disk is found to lead to fast boundary layer accretion along the disk surfaces and to enhanced accretion in the main disk. These models are pertinent to the formation of counterrotating disks in galaxies and possibly in Active Galactic Nuclei and in X-ray pulsars in binary systems.

The basic parameters determining the flow pattern for a nonviscous, non-heat-conducting gas in a semidetached binary system without a magnetic field are identified. Three-dimensional gas-dynamical modeling of the mass transfer enables investigation of the influence of these parameters on the structure of gas flows. The parameter on which the flow pattern depends most strongly is the adiabatic index $\gamma$. The effect of other parameters is small, and only leads to unimportant quantitative changes in the solutions obtained. The main properties of flows typical of semidetached binaries without magnetic fields are summarized.

The results of 3D hydrodynamic simulation of mass transfer in semidetached binaries of different types (cataclysmic variables and low-mass X-ray binaries) are presented. We find that taking into account of a circumbinary envelope leads to significant changes in the stream-disc morphology. In particular, the obtained steady-state self-consistent solutions show an absence of impact between gas stream from the inner Lagrangian point L1 and forming accretion disc. The stream deviates under the action of gas of circumbinary envelope, and does not cause the shock perturbation of the disc boundary (traditional `hotspot'). At the same time, the gas of circumbinary envelope interacts with the stream and causes the formation of an extended shock wave, located on the stream edge. We discuss the implication of this model without `hotspot' (but with a shock wave located outside the disc) for interpretation of observations. The comparison of synthetic light curves with observations proves the validity of the discussed hydrodynamic model without `hotspot'. We also consider the influence of a circumbinary envelope on the mass transfer rate in semidetached binaries. The obtained features of flow structure in the vicinity of L1 show that the gas of circumbinary envelope plays an important role in the flow dynamics, and that it leads to significant (in order of magnitude) increasing of the mass transfer rate. The comparison of gaseous flows structure obtained in 2D and 3D approaches is presented. We discuss the common features of the flow structures and the possible reasons of revealed differences.

The results of three-dimensional numerical simulations of the flow of matter in non-magnetic semidetached binary systems are presented. Self-consistent solutions indicate the absence of a shock interaction between the stream of matter flowing from the inner Lagrange point and the accretion disk (a "hot spot"). At the same time, the interaction between the stream and the common envelope of the system forms an extended shock wave along the edge of the stream, whose observational properties are roughly equivalent to those of a hot spot in the disk. Comparison of synethesized and observed light curves confirm the adequacy of the proposed model for the flow of matter without formation of a hot spot. The model makes it possible for the first time to reproduce the variety of the observed light curves of cataclysmic variables in the framework of a single model, without invoking additional assumptions.

Results of numerical simulations of the impact of a common envelope on the matter flow pattern near the outflowing component in a semidetached binary system are presented. Three-dimensional modeling of the matter transfer gas dynamics in a low-mass X-ray binary X1822-371 enable investigation of the structure of flows in the vicinity of the inner Lagrange point L1. Taking into account the common envelope of the system substantially changes the flow pattern near the Roche surface of the outflowing component. In a stationary regime, accretion of common envelope gas is observed over a significant fraction of the donor star's surface, which inhibits the flow of gas along the Roche surface to L1. The change in the flow pattern is particularly significant near L1, where the stream of common envelope gas strips matter off the stellar surface. This, in turn, significantly increases (by an order of magnitude) the gas flow from the donor surface in comparison with the estimates of standard models.

Results of numerical simulations of matter flows in a semidetached binary system similar to the low-mass X-ray binary X1822--371 are presented. Three-dimensional modeling of the mass transfer gas dynamics makes it possible to investigate gas streams in the system and to study the influence of a common envelope. The presence of the common envelope leads to the absence of shock interaction between the stream of matter flowing from the inner Lagrange point and the gas in the accretion disk. The stream is deflected by the gas in the common envelope and approaches the disk tangentially, so that it does not cause any shock perturbation ("hot spot") on the disk. At the same time, the interaction of the stream with the common envelope leads to the formation of an extended shock wave along the edge of the stream. The observational manifestation of this shock is estimated to be equivalent to that of a hot spot in the disk. The calculated accretion disk parameters are presented.