Search

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., Comment: 6 pages, 3 figures, Conference Modern stellar astronomy-2019, SAO RAS, October 7-11, 2019 (Submitted: Research in Astronomy and Astrophysics)

A new dispersion equation is obtained for a non-equilibrium medium with an exponential relaxation model of a vibrationally excited gas. We have researched the dependencies of the pump source and the heat removal on the medium thermodynamic parameters. The boundaries of sound waves stability regions in a non-equilibrium gas have been determined. The nonlinear stage of sound waves instability development in a vibrationally excited gas has been investigated within CSPH-TVD and MUSCL numerical schemes using parallel technologies OpenMP-CUDA. We have obtained a good agreement of numerical simulation results with the linear perturbations dynamics at the initial stage of the sound waves growth caused by instability. At the nonlinear stage, the sound waves amplitude reaches the maximum value that leads to the formation of shock waves system., Comment: 11 pages, 5 figures

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., Comment: 12 pages, 7 figures

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., Comment: 15 pages, 15 figures, accepted for publication in A&A

The problem of choice of boundary conditions are discussed for the case of numerical integration of the shallow water equations on a substantially irregular relief. In modeling of unsteady surface water flows has a dynamic boundary partitioning liquid and dry bottom. The situation is complicated by the emergence of sub- and supercritical flow regimes for the problems of seasonal floodplain flooding, flash floods, tsunami landfalls. Analysis of the use of various methods of setting conditions for the physical quantities of liquid when the settlement of the boundary shows the advantages of using the waterfall type conditions in the presence of strong inhomogeneities landforms. When there is a waterfall on the border of the computational domain and heterogeneity of the relief in the vicinity of the boundary portion may occur, which is formed by the region of critical flow with the formation of a hydraulic jump, which greatly weakens the effect of the waterfall on the flow pattern upstream., Comment: 16 p., in Russian, 6 figures