Vol 43, No 2 (2017)

The instability of anisotropic disk systems with elongated stellar orbits has been investigated. N-body generalized polytropic models of stellar disks have been constructed. They are shown to be unstable with respect to the bar formation at any degree of anisotropy. This result differs from the results of the studies of such models by other authors. The bar pattern speed and amplitude have been found. The initial distribution of precession rates and the adiabatic invariants of stellar orbits have been calculated. A bar is shown to be formed in such systems due to the radial orbit instability.

Using the GADGET-2 code modified by us, we have computed hydrodynamic models of a protoplanetary disk perturbed by a low-mass companion. We have considered the cases of circular and eccentric orbits coplanar with the disk and inclined relative to its midplane. During our simulations we computed the column density of test particles on the line of sight between the central star and observer. On this basis we computed the column density of circumstellar dust by assuming the dust and gas to be well mixed with a mass ratio of 1: 100. To study the influence of the disk orientation relative to the observer on the interstellar extinction, we performed our computations for four inclinations of the line of sight to the disk plane and eight azimuthal directions. The column densities in the circumstellar disk of the central star and the circumbinary disk were computed separately. Our computations have shown that periodic column density oscillations can arise in both inner and circumbinary disks. The amplitude and shape of these oscillations depend on the system’s parameters (the orbital eccentricity and inclination, the component mass ratio) and its orientation in space. The results of our simulations can be used to explain the cyclic brightness variations of young UX Ori stars.

The “mass loading” of the solar wind by cometary ions produced by the photoionization of neutral molecules outflowing from the cometary nucleus plays a major role in the interaction of the solar wind with cometary atmospheres. In particular, this process leads to a decrease in the solar wind velocity with a transition from supersonic velocities to subsonic ones through the bow shock. The so-called single-fluid approximation, in which the interacting plasma flows are considered as a single fluid, is commonly used in modeling such an interaction. However, it is occasionally necessary to know the distribution of parameters for the components of the interacting plasma flows. For example, when the flow of the cometary dust component in the interplanetary magnetic field is considered, the dust particle charge, which depends significantly on the composition of the surrounding plasma, needs to be known. In this paper, within the framework of a three-dimensional magnetohydrodynamic model of the solar wind flow around cometary ionospheres, we have managed to separately obtain the density distributions of solar wind protons and cometary ions between the bow shock and the cometary ionopause (in the shock layer). The influence of the interplanetary magnetic field on the position of the point of intersection between the densities with the formation of a region near the ionopause where the proton density is essentially negligible compared to the density of cometary ions is investigated. Such a region was experimentally detected by the Vega-2 spacecraft when investigating Comet Halley in March 1986. The results of the model considered below are compared with some experimental data obtained by the Giotto spacecraft under the conditions of flow around Comets Halley and Grigg–Skjellerup in 1986 and 1992, respectively. Unfortunately, our results of calculations on Comet Churyumov–Gerasimenko are only predictive in character, because the trajectory of the Rosetta spacecraft, which manoeuvred near its surface for several months, is complex.