Numerical modeling of transient structures in the disks of spiral galaxies

Nonstationary gas-dynamical processes occuring in the disks of spiral galaxies due to the external gravitational field have been studied using numerical simulations. A series of hydrodynamical discontinuities, including strong shock fronts and contact discontinuities, arises as a result of the nonlinear, supersonic interaction between emerging spiral formations and the flow of matter in the two-armed global morphology.     

Mathematical modeling of the structure of an accretion disk in a stellar binary system

Numerical simulations of gas-dynamical processes taking place in the accretion disk of a stellar binary system are presented. The initial state of the disk is an equilibrium gaseous configuration. Mechanisms for the development of spiral waves and associated variations in the angular momentum of the gas are considered. The influence of the ratio of the binary-component masses and the initial disk configuration are investigated. It is concluded that the existence of a steady-state disk is impossible without a flow of gas from the donor star.     

Comparison between two-and three-dimensional modeling of the structure of an accretion disk in a binary system

We compare two-and three-dimensional modeling of gas-dynamical processes in the accretion disk of a binary system. The origin of spiral waves and the loss of the angular momentum related to them are considered. It is concluded that a steady state of the disk cannot exist without taking into account t he gas inflow from the donor star.     

The structure of cool accretion disks in semidetached binaries

We present a qualitative analysis of possible changes in the structure of accretion disks that occur in the transition from hot to cool disks. We suggest that an additional spiral-density wave can exist in the inner parts of the disk, where gas-dynamical perturbations are negligible. We consider the formation of this wave and its parameters. The results of a three-dimensional gas-dynamical simulation of a cool accretion disk are presented; these results confirm the possibility of the formation of a new, “precessional,” spiral wave in the inner regions of a cool accretion disk. Possible observational manifestations of such a wave are discussed.     

A model for superoutbursts in SU UMa-type binaries

We suggest a new mechanism for the superoutbursts in SU UMa binaries, in which the increase in the accretion rate resulting in a superoutburst is associated with the formation of a spiral “precessional” wave in the inner parts of the disk, where gas-dynamical perturbations are negligible. The existence of such waves was suggested by us previously. The results of three-dimensional gas-dynamical simulations have shown that a considerable increase in the accretion rate (by up to an order of magnitude) is associated with the formation of the precessional wave. The features of the precessional spiral wave can explain both the energy release in the superoutburst and all its observational manifestations. One distinguishing feature of superoutbursts in SU UMa-type stars is the formation of a “superhump” in the light curve. Our model reproduces well both the formation of a superhump and its observational features, including its period, which is up to 3–7% longer than the orbital period, and the detectability of the superhump independent of the orbital inclination of the binary.     

Synthetic Doppler tomograms of gas flows in the binary system IP Peg

We have synthesized Doppler tomograms of gas flows in the binary system IP Peg using the results of three-dimensional gas-dynamical computations. Gas-dynamical modeling in combination with Doppler tomography enables identification of the key elements of flows in Doppler maps without solution of an ill-posed inverse problem. A comparison of the synthetic tomograms with observations shows that, in the quiescent state of the system, the most luminous components are (1) the shock wave induced by interaction between the circumbinary envelope and the stream from the Lagrange point L ₁ (the “hot line”) and (2) the gas condensation at the apogee of the quasi-elliptical disk. Both the single spiral shock wave arm in the gas-dynamical solution and the stream from L ₁ contribute little to the luminosity. In the active state of the system, when the stream from L ₁ does not play an appreciable role and the disk dominates, both areas of enhanced luminosity in the observational tomograms are associated with the two arms of the spiral shock wave in the disk.     

The structure of the common envelope in a close binary system

As a result of the interaction between an elliptical accretion disk and gas flowing into it from the circumbinary envelope in a close binary in the course of its orbital motion, the matter of the disk and the circum-disk halo is periodically ejected from the vicinity of the Lagrange point L₃, and a common envelope is formed in the system. Three-dimensional numerical gas-dynamical modeling is used to study the structure and dynamics of the envelope and determine its basic parameters. The evolution of the envelope is followed on timescales of the order of several orbital periods. The matter flow ejected through the vicinity of L₃ displays a spiral shape. The maximum size of the forming spiral structure is restricted by the self-intersection point, and is of the order of four to five times the component separation. We consider the dynamics of the regions directly adjacent to the spiral structure: an inner, rarified and outer, fragmented region, which further makes a transition to an expanding diffuse ring.     

The structure of matter flows in semi-detached binaries after the termination of mass transfer

We present the results of three-dimensional gas-dynamical simulations of matter flows in semi-detached binaries after termination of the mass transfer between the components of the system. The structure of the residual accretion disk is studied. When the mass transfer has ended, the quasi-elliptical disk becomes circular and its structure changes: tidal interactions result in the formation of a second arm in the spiral shock wave. In addition, a condensation (blob) moving through the disk with variable velocity is formed. The blob is maintained by interactions with the arms of the spiral shock and exists essentially over the entire lifetime of the disk. We also show that, for a viscosity corresponding to α～0.01 (typical for observed accretion disks), the lifetime of the residual disk is about 50 orbital periods.     

Supernova explosion mechanism taking into account large-scale convection and neutrino transport

Two types of supernovae are considered: thermonuclear supernovae, whose explosions are due to thermonuclear energy, and core-collapse supernovae, whose explosions are due to the gravitational energy of collapsing stars released in the form of neutrinos. Numerical models of supernovae are discussed. Themain problem in devising supernova explosion mechanisms is producing the energy required to disperse the envelope. In theoretical models, it is necessary to solve multi-dimensional problems involving complex physics (3D gas dynamics, neutrino transport, large-scale convective instability, and other important physical processes). In recent years, the development of large-scale convection during supernova explosions has been reconsidered. Self-consistent problems problems in three-dimensional, gas-dynamical instability have been considered. Two-dimensional gas-dynamical calculations taking into account neutrino absorption in the envelope have been performed. The spherically symmetric collapse and neutrino transport were calculated including all reactions, leading to a new understanding of possible paths for the development of supernova theory. The main emphasis is placed on the neutrino transport and the basis for promising multidimensional models taking into account large-scale convective instability.     

Polygonal structure of spiral galaxies

We have carried out numerical simulations of hydrodynamical processes occurring in the disks of spiral galaxies. The initial state of the disk is an equilibrium stellar-gaseous configuration. The spherical component is described by a standard analytical model for the gravitational potential. The behavior of the modeled disk in the presence of an external perturbation is analyzed. The results of numerical simulations of stellar-gaseous galactic disks aimed at studying the formation of polygonal structures in spiral galaxies are presented. The possible influence of spur-like formations on the appearance of polygonal structure is studied.     

The Influence of Superflares of Host Stars on the Dynamics of the Envelopes of Hot Jupiters

Results of a study of the influence of solar-type host stars superflares on the gas dynamics of the extended envelopes of giant exoplanets are presented. During flare events, the radiation intensity of the host star in the extreme ultraviolet and soft X-ray can increase by several orders of magnitude for a short time, leading to strong local heating of the exoplanet atmosphere on the side facing the star, with the formation of shocks in the atmosphere. Computations of the gas-dynamical response of the atmosphere of the hot Jupiter HD 209458b to characteristic superflares of solar-like stars were carried out earlier in [1] using a one-dimensional aeronomical model correctly taking into account heating and chemical processes in the atmosphere. To investigate the outflow of atmospheric gas, the results obtained with this onedimensional model were used as simple boundary conditions for computations of the three-dimensional flow structure after a flare. The results of these three-dimensional gas-dynamical computations show that the mass ejection of the flare increases the size of the envelope over several hours, which could be detected with existing observing facilities. It is shown that the mass-loss rates for the most powerful superflare considered could be enhanced by an order of magnitude over several tens of hours after the flare.

     

Doppler mapping of the SS Cyg system during outburst

Spectral observations of the SS Cyg system in its active state are used to construct Hβ and Hγ Doppler tomograms. These are compared with analogous tomograms for the quiescent state and synthetic tomograms derived from the results of three-dimensional gas-dynamical modeling. The parameters of the accretion disk during the outburst are estimated. An explanation for the observed flow pattern is proposed, based on a numerical model with an elliptical accretion disk.     

Interaction of the dark-matter cusp with the baryonic component in disk galaxies

The influence of the formation and evolution of a (disk) galaxy on the matter distribution in the dark-matter halo is considered. Calculations of the evolution of an isolated dark-matter halo were carried out with and without including a baryonic component. N-body simulations (for the dark-matter halo) and gas-dynamical numerical simulations (for the baryonic gas) were used for this analysis. Star formation, feedback, and heating and cooling of the interstellar medium were taken into account in the gas-dynamical calculations. The results of these numerical simulations with high spatial resolution indicate that 1) including the star formation resolves the so-called cusp problem (according to CDMcosmological models, the density distribution in the central regions of the dark-matter halo should have a distinct peak (cusp), which is not shown by observations); 2) the interaction of the dark matter with dynamical substructures of the stellar-gas galactic disk (spiralwaves, a bar) affects the shape of the dark-matter halo. In particular, the calculated dark-matter distribution in the plane of the disk is more symmetric when the baryonic component is taken into account.     

Observational Manifestations of the Precessional Wave in the Accretion Disk of a Cataclysmic Variable Star

Effects due to the interaction of the steam from the inner Lagrangian point with the accretion disk in a cataclysmic variable star are considered. The results of three-dimensional gas-dynamical numerical simulations confirm that the disk thickness in the vicinity of the interaction with the stream is minimum when the component-mass ratio is 0.6. As a consequence, some of the matter from the stream does not collide with the outer edge of the accretion disk, and continues its motion unperturbed toward the accretor. This part of the stream subsequent interacts (collides) with a thickening of the accretion disk due to the presence of a precessional wave in the disk, leading to the appearance of an additional zone of heating at the disk surface. This additional zone of enhanced luminosity (hot spot) is a direct observational manifestation of the precessional wave in the accretion disk.     

Origin of prolonged X-ray flares on active late-type stars

Soft X-ray data for prolonged flares in subgiants in RS CVn binary systems and some other active late-type stars (AB Dor, Algol) are analyzed. During these nonstationary events, a large amount of hot plasma with temperatures exceeding 10⁸ K exists for many hours. Numerical simulations of gas-dynamical processes in the X-ray source—giant loops—can yield reliable estimates of the plasma parameters and flare-source size. This confirms that such phenomena exist while considerable energy is supplied to the top part of a giant loop or system of loops. Refined estimates of the flare energy (up to 10³⁷ erg) and scales contradict the widely accepted idea that prolonged X-ray flares are associated with the evolution of local magnetic fields. The energy of the current component of the large-scale magnetic field arising during the ejection of magnetic field by plasma jets or stellar wind is estimated. Two cases are considered: a global stellar field and fields connecting regions with oppositely directed unipolar magnetic fields. The inferred energy of the current component of the magnetic field associated with distortion of the initial MHD configuration is close to the total flare energy, suggesting that large-scale magnetic fields play an important role in prolonged flares. The flare process encompasses some portion of a streamer belt and may propagate along the entire magnetic equator of the star during the most powerful prolonged events.     

A study of the outburst development in the classical symbiotic star Z And within the colliding-winds model

Two-dimensional gas-dynamical modeling of the mass-flow structures in binary systems is used to study the outburst development in the classical symbiotic star Z And. The stage-by-stage rise of the light during the outburst can be explained using a model with colliding winds. We suggest a scenario for the development of the outburst and study the influence of possible changes in the flow structure on the brightness of the system. The model variations of the luminosity due to the formation of a system of shocks are in good agreement with the observed brightness variations.     

Turbulence in the Accretion Disks in Semi-Detached Binary Systems with Various Component-Mass Ratios

Results of three-dimensional gas-dynamical numerical simulations of the structure of matter flows in semi-detached binary systems with various component-mass ratios are considered. The main elements of the flows in the models studied are described. The characteristics of density waves inside the accretion disk for various component-mass ratios are considered. The influence of the precessional density wave on the development of instability in the accretion disks and the characteristics of developing turbulence are analyzed. Values of the Shakura–Syunyaev coefficient α for the simulated systems are calculated.     

Fine vertical structure of the galactic gaseous disk

The three-dimensional evolution of an ensemble of N particles (N = 8 × 10⁵) in the external gravitational field of a galaxy perturbed by a spiral density wave is considered. The particles simulate clouds of interstellar gas, and inelastic two-body collisions between them are taken into account. The three-dimensional structure of the gaseous galactic layer and the vertical profile of the spiral arms are computed. It is shown that: (1) the local thickness of the gaseous galactic disk has a minimum where the volume gas density has a maximum (the maximum density of the interstellar medium is shifted downstream relative to the galactic shock front), (2) the configuration of the vertical profile of the spiral arms changes radically when the corotation region is crossed. Our first result explains the negative correlation between the thickness of the gas layer and the density derived from neutral-hydrogen observations. The second result can be used in the next generation of neutral-hydrogen observations to localize the corotation radius in the Galaxy.     

Structure of the circumbinary envelopes of young binary stars with elliptical orbits

The results of two-dimensional gas-dynamical numerical simulations of the structure of matter flows in the envelopes of a number of T Tauri binary systems with elliptical orbits are considered. The main flow elements in inner regions of protoplanetary disks of these stars are described. The influence of shocks on the size of the gap—a rarified region in the inner parts of the protoplanetary disk—is analyzed. A method is proposed for estimating the size of this gap from the numerical simulations, and the gap sizes for the studied stars are determined and compared with observational results. The flow dynamics in the gap is considered, and the periodic variations of the gap size on time scales of several orbital periods are analyzed. Possible observational manifestations of the studied flows are discussed.     

Impact of the magnetic field on the structure of accretion disks in semi-detached binaries

We discuss characteristic features of the magnetic gas-dynamical structure of the flows in a semi-detached binary system obtained from three-dimensional simulations, assuming that the intrinsic magnetic field of the accreting star is dipolar. The turbulent diffusion of the magnetic field is taken into account. The SS Cyg system is considered as an example. Including the magnetic field can alter the basic parameters of the accretion disk, such as the accretion rate and the characteristic density. The magnetic field in the disk is primarily toroidal.