In this paper we introduce a new method for analysing Milky Way phase-space which allows us to reveal the imprint left by the Milky Way bar and spiral arms on the stars with full phase-space data in Gaia Data Release 2. The unprecedented quality and extended spatial coverage of these data enable us to discover six prominent stellar density structures in the disc to a distance of 5 kpc from the Sun. Four of these structures correspond to the spiral arms detected previously in the gas and young stars (Scutum-Centaurus, Sagittarius, Local and Perseus). The remaining two are associated with the main resonances of the Milky Way bar where corotation is placed at around 6.2 kpc and the outer Lindblad resonance beyond the Solar radius, at around 9 kpc. For the first time we provide evidence of the imprint left by spiral arms and resonances in the stellar densities not relying on a specific tracer, through enhancing the signatures left by these asymmetries. Our method offers new avenues for studying how the stellar populations in our Galaxy are shaped.

Using observational data on the kinematical properties and density distributions of the subsystems of the Milky Way galaxy, we construct a set of multi-component equilibrium models of its disk. The dynamics of the disk is studied numerically using collisionless-gaseous numerical simulations. After approximately one Gyr, a prominent central bar is formed with a semi-axis of about three kiloparsecs. Outside the central regions, a multi-armed spiral pattern develops, which can be characterized by the superposition of m=2, 3 and 4-armed spiral patterns. The spiral structure and the bar exist for at least 3 Gyr in our simulations. The presence of the bar in the disk of the MilkyWay galaxy imposes rather strict limitations on the density distributions in the subsystems of the MilkyWay galaxy.We find that bar does not form if the radial scale length of the surface density distribution of the disk is more than 2.6 kpc. Analogously, the formation of bar is suppressed in the Milky Way disk in models with a massive and compact stellar bulge. If future GAIA observations confirm the existence of the three-four armed spiral pattern in the disk of the MilkyWay, this will prove the long-term existence of spiral patterns in galactic disks.

Parallel implementation features of self-gravitating gas dynamics modeling on multiple GPUs are considered applying the GPU-Direct technology. The parallel algorithm for solving of the self-gravitating gas dynamics problem based on hybrid OpenMP-CUDA parallel programming model has been described in detail. The gas-dynamic forces are calculated by the modified SPH-method (Smoothed Particle Hydrodynamics) while the N-body problem gravitational interaction is obtained by the direct method (so-called Particle-Particle algorithm). The key factor in the SPH-method performance is creation of the neighbor lists of the particles which contribute into the gas-dynamic forces calculation. Our implementation is based on hierarchical grid sorting method using a cascading algorithm for parallel computations of partial sums at CUDA block. The parallelization efficiency of the algorithm for various GPUs of the Nvidia Tesla line (K20, K40, K80) is studied in the framework of galactic' gaseous halos collisions models by the SPH-method.

This paper focuses on the parallel implementation of a direct $N$-body method~(particle-particle algorithm) and the application of multiple GPUs for galactic dynamics simulations. Application of a hybrid OpenMP-CUDA technology is considered for models with a number of particles $N \sim 10^5 \div 10^7$. By means of $N$-body simulations of gravitationally unstable stellar galactic we have investigated the algorithms parallelization efficiency for various Nvidia Tesla graphics processors~(K20, K40, K80). Particular attention was paid to the parallel performance of simulations and accuracy of the numerical solution by comparing single and double floating-point precisions~(SP and DP). We showed that the double-precision simulations are slower by a factor of~$1.7$ than the single-precision runs performed on Nvidia Tesla K-Series processors. We also claim that application of the single-precision operations leads to incorrect result in the evolution of the non-axisymmetric gravitating $N$-body systems. In particular, it leads to significant quantitative and even qualitative distortions in the galactic disk evolution. For instance, after $10^4$ integration time steps for the single-precision numbers the total energy, momentum, and angular momentum of a system with $N = 2^{20}$ conserve with accuracy of $10^{-3}$, $10^{-2}$ and $10^{-3}$ respectively, in comparison to the double-precision simulations these values are $10^{-5}$, $10^{-15}$ and $10^{-13}$, respectively. Our estimations evidence in favour of usage of the second-order accuracy schemes with double-precision numbers since it is more efficient than in the fourth-order schemes with single-precision numbers.

We study numerically large-scale magnetic field evolution and its enhancement in gaseous disks of spiral galaxies. We consider a set of models with the various spiral pattern parameters and the initial magnetic field strength with taking into account gas self-gravity and cooling/heating. In agreement with previous studies, we find out that galactic magnetic field is mostly aligned with gaseous structures, however small-scale gaseous structures (spurs and clumps) are more chaotic than the magnetic field structure. In spiral arms magnetic field strongly coexists with the gas distribution, in the inter-arm region we see filamentary magnetic field structure. Simulations reveal the presence of the small-scale irregularities of the magnetic field as well as the reversal of magnetic field at the outer edge of the large-scale spurs. We provide evidences that the magnetic field in the spiral arms has a stronger mean-field component, and there is a clear inverse correlation between gas density and plasma-beta parameter, compared to the rest of the disk with a more turbulent component of the field and an absence of correlation between gas density and plasma-beta. We show the mean field growth up to 3-10$\mu G$ in the cold gas during several rotation periods (500-800 Myr), whereas ratio between azimuthal and radial field is equal to 4/1. Mean field strength increases by a factor of 1.5-2.5 for models with various spiral pattern parameters. Random magnetic field component can reach up to 25 % from the total strength. By making an analysis of the time-depended evolution of radial Poynting flux we point out that the magnetic field strength is enhanced stronger at the galactic outskirts which is due to the radial transfer of magnetic energy by the spiral arms pushing the magnetic field outward. Our results also support the presence of sufficient conditions for development of MRI at distances >11 kpc.

Using $N$-body simulations ($N\sim 10^6 - 10^7$), we examine how a non-axisymmetric dark halo affects the dynamical evolution of the structure in collisionless (stellar) discs. We demonstrate how the model parameters such as mass of the halo, initial conditions in the disc and the halo axes ratio affect morphology and kinematics of the stellar discs. We show that a non-axisymmetric halo can generate a large-scale spiral density pattern in the embedded stellar disc. The pattern is observed in the disc for many periods of its revolution, even if the disc is gravitationally over-stable. The growth of the spiral arms is not accompanied by significant dynamical heating of the disc, irrelevant to its initial parameters. We also investigate transformation of the dark halo's shape driven by the long-lived spiral pattern in the disc . We show that the analysis of the velocity field in the stellar disc and in the spiral pattern gives us a possibility to figure out the spatial orientation of the triaxial-shaped dark halo and to measure the triaxiality.

We present results of numerical modeling made for the galactic stellar and stellar-gas disk embedded in the spherical halo and bulge. The stellar disk is simulated by N-body system, the equations of hydrodynamics are solved by TVD-method. We used TREEcode-algorithm for calculation of a self-gravity in stellar and gaseous components. The possibility of bars birth in a hot stellar disk because of gravitational instability of a cold gas component is investigated. The conditions of occurrence lopsided-galaxies from a axisymmetric disk as a result of gravitational instability are explored. The self-consistent models of double bars are constructed and the dynamical stability of these structures is discussed.

The results of numerical simulations of the dynamics of galactic disks, which are submerged into nonaxisymmetric dark massive halo are discussed. Galaxy disks dynamics in the nonaxisymmetric (triaxial) dark halo were investigated in detail by the high resolution numerical hydrodynamical methods (TVD \& SPH) and N-body models. The long-lived two arms spiral structure generates for a wide range of parameters. The spiral structure is global and number of turns can be 2-3 in depends of model parameters. Morphology and kinematics of spiral structures were investigated in depends of the halo and the disk parameters. The spiral structure rotates slowly and the angular velocity varies is quasiperiodic.

In Part I of this paper we described an equilibrium model of a jet in the gravitational field corresponding to the rigid-rotation region of the galactic disk. We used linear stability analysis to find the waveguide-resonance instability of internal gravity waves due to the superreflection of these waves from the jet boundary. In this part of the paper, we perform nonlinear numerical 2D and 3D simulations of the development of this instability. We show that the shocks produced by this instability in the ambient medium of the jet are localized inside a cone with a large opening angle and are capable of producing features that are morphologically similar to those observed in galaxies with active nuclei (NGC 5252 for example).

We discuss causes of the formation of the observed kinematics and morphology of cones of ionized matter in the neighborhood of the nuclei of Seyfert galaxies. The results of linear stability analysis of an optically thin conic jet where radiation cooling and gravity play an important part are reported. The allowance for radiation cooling is shown to result in strong damping of all acoustic modes and to have insignificant effect on unstable surface Kelvin--Helmholtz modes. In the case of waveguide--resonance internal gravity modes radiative cooling suppresses completely the instability of waves propagating away from the ejection source and, vice versa, reduces substantially the growth time scale of unstable sourceward propagating modes. The results obtained can be used to study ionization cones in Seyfert galaxies with radio jets. In particular, our analysis shows that surface Kelvin--Helmholtz modes and volume harmonics are capable of producing regular features observed in optical emission-line images of such galaxies.

The Seyfert galaxies with Z-shaped emission filaments in the Narrow Line Region (NLR) are considered. We assume that observable Z-shaped structures and velocity pattern of NLR may be explained as tridimensional helical waves in the ionization cone.

Results of 2D spectral observations for 4 Sy galaxies with ionization cones are presented. Images in the [OIII] line, velocity fields of the ionized gas, stellar velocity and velocity dispersion fields were obtained at the SAO RAS 6m telescope. Non-circular gas motions and Z-shaped emission filaments could be explained as helical waves located in the ionization cone. This waves are generated by the hydrodynamical instability due to the velocity break between galactic ISM and outflowing matter from the AGN engine. The axis of cone lies close to the direction of a nuclear mini bar.

Dynamics of linear perturbations in a differentially rotating accretion disk with non-homogeneous vertical structure is investigated. It has been found that turbulent viscosity results in instability of both pinching oscillations, and bending modes. Not only the low-frequency fundamental modes, but also the high-frequency reflective harmonics appear to be unstable. The question of the limits of applicability of the thin disk model (MTD) is also investigated. The insignificant distinctions in the dispersion properties of MTD and three-dimensional model appear for wave numbers k<(1-3)/h (h is the half-thickness of a disk). In the long-wavelenght limit, the relative difference between eigenfrequencies of the unstable acoustic mode in the 3D-model and the MTD is smaller than 5%.

The model of a geometrically thin gaseous disk in the external gravitational potential is considered. The dinamics of small nonaxisymmetric perturbations in the plane of the accretion disk with dissipative effects is investigated. It is showed, that conditions of development and parameters of unstable oscillation modes in the opticaly thick accretion disk are strongly depended on the models of viscosity and opacity.