Excitation of Turbulence in Protoplanetary Disks around Binary Stars

Accretion disks around young binary stars are subject to strong forces exerted by the system components. Gas–dynamical interactions excite strong non-linear perturbations in the disk, which can give rise to turbulence. This study considers a mechanism for the excitation of turbulence based on the instability of non-linear waves in a rotating flow. It is shown that the spectrum of non-linear perturbations excited in the inner part of the disk leads to turbularization of the flow. Estimates of the Shakura–Syunyaev index, α ～ 0.01?0.05, are obtained for use in numerical models of accretion disks.     

Non-steady-state accretion disks in X-ray novae: Outburst models for Nova Monocerotis 1975 and Nova Muscae 1991

We have fit outbursts of two X-ray novae (Nova Monocerotis 1975=A0620-00 and Nova Muscae GS 1991=1124-683) using a non-steady-state accretion-disk model. The model is based on a new solution for a diffusion-type equation for non-steady-state accretion and describes the evolution of a viscous α disk in a binary system after the peak of the outburst, when the matter in the disk is totally ionized. The accretion rate in the disk decreases according to a power law. We derive formulas for the accretion rate and effective temperature of the disk. The model has three free input parameters: the mass of the central object M, the turbulence parameter α, and the normalization parameter δt. The results of the modeling are compared with the observed X-ray and optical B and V light curves. The estimates for the turbulence parameter α are similar: 0.2–0.4 for A 0620-00 and 0.45–0.65 for GS 1124-683, suggesting a similar nature for the viscosity in the accretion disks around the compact objects in these sources. We have also derived the distances to these systems as functions of the masses of their compact objects.     

The development of large-scale instability in Keplerian stellar accretion disks

Our previous studies of large-scale vortical flows arising in shear flows of stellar accretion disks with Keplerian azimuthal velocity distributions as a result of the development of small perturbations are continued. The development of large-scale instability in an accretion disk is investigated via mathematical modeling. One result obtained is the change of the disk flow structure due to the formation of large vortices. In the limiting case, sufficiently long evolution leads to the formation of several asymmetric spiral structures of the flow of disk matter. The presence of large-scale structures leads to angular-momentum redistribution in the disk.     

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.     

The impact of large-scale turbulence on the redistribution of angular momentum in stellar accretion disks

We consider the origin and development of large-scale turbulence in a shear flow in a stellar accretion disk. The ratio of the kinetic energy of vortices originating in the turbulent flow and the total initial kinetic energy of the rotating disk is essentially constant. The large-scale structures that form are able to redistribute the angular momentum without any appreciable heating of the matter.     

The development of large-scale instability in stellar accretion disks and its influence on the redistribution of angular momentum

The paper continues our studies of large-scale instability arising during shearmotions in stellar accretion disks due to the development of small perturbations. The evolution of a local perturbation introduced into the outer part of a stationary accretion disk is modeled mathematically. The possible formation of large-scale structures that propagate throughout the disk, leading to an appreciable redistribution of angular momentum, is demonstrated.     

Amplification of fluctuations and current dynamo due to helical and chiral effects in geophysical and plasma-like media

Dynamo effect is considered in more general case than well known one when not only mean motion but also conducting components possess nonzero mean helicity. Dispersion equation for helical motions is studied, including inhomogeneous case. Criteria for development of instability are found. Interaction between large-scale internal wave and small-scale helical turbulence in plane Couette flow of fluid with statically stable uniform density gradient in gravitational field is studied basing on a set of equations for scalar and quasi-scalar wave fields and for helicity of turbulence. This system of equation describes acceleration of wave growth due to its backward action on turbulence, i.e. wave-turbulent instability. Chiral media are further considered and effects topologically close to helical ones are analyzed. It is shown that in such media ranges of scales exist, in which fluctuations of electric field can be anomalously amplified.     

Three-dimensional gas-dynamical modeling of changes in the flow structure during the transition from the quiescent to the active state in symbiotic stars

The results of three-dimensional modeling of the flow structure in the classical symbiotic system Z Andromedae are presented. Outbursts in systems of this type occur when the accretion rate exceeds the upper limit of the steady-burning range. Therefore, in order to realize the transition from a quiescent to an active state, it is necessary to find a mechanism capable of sufficiently increasing the accretion rate on the time scales typical for outburst development. Our calculations provide support for a mechanism for the transition from quiescence to outburst in classical symbiotic systems suggested earlier based on two-dimensional calculations (Bisikalo et al., 2002). Our results show that an accretion disk forms in the system for a wind velocity of 20 km s^?1. The accretion rate for the solution with the disk is ～22.5–25% of the mass-loss rate of the donor, which is ～4.5?5 × 10^?8M_⊙ yr^?1 for Z And. This value is in agreement with the steady-burning range for the white-dwarf masses usually accepted for this system. When the wind velocity increases from 20 to 30 km s^?1, the accretion disk is destroyed and the disk material falls onto the accretor surface. This process is followed by an approximately twofold jump in the accretion rate. The resulting growth in the accretion rate is sufficient so as to exceed the upper limit of the steady-burning range, thus bringing the system into an active state. The time during which the accretion rate is above the steady-burning value is in very good agreement with observations. Our analysis leads us to conclude that small variations in the donor wind velocity can lead to the transition from disk accretion to wind accretion and, as a consequence, to the transition from a quiescent to an active state in classical symbiotic stars.     

Morphology and evolution of bars in a wandering gravel-bed river; lower Fraser river, British Columbia, Canada

A hierarchical typology for the channels and bars within aggradational wandering gravel-bed rivers is developed from an examination of a 50 km reach of lower Fraser River, British Columbia, Canada. Unit bars, built by stacking of gravelly bedload sheets, are the key dynamic element of the sediment transfer system, linking sediment transport during individual freshets to the creation, development and remoulding of compound bar platforms that have either a lateral or medial style. Primary and secondary unit bars are identified, respectively, as those that deliver sediment to compound bars from the principal channel and those that redistribute sediment across the compound bar via seasonal anabranches and smaller channels. The record of bar accretion evident in ground-penetrating radar sequences is consistent with the long-term development of bar complexes derived from historical aerial photographs. For two compound bars, inter-annual changes associated with individual sediment transport episodes are measured using detailed topographic surveys and longer-term changes are quantified using sediment budgets derived for individual bars from periodic channel surveys. Annual sediment turnover on the bars is comparable with the bed material transfer rate along the channel, indicating that relatively little bed material bypasses the bars. Bar construction and change are accomplished mainly by lateral accretion as the river has limited capacity to raise bed load onto higher surfaces. Styles of accretion and erosion and, therefore, the major bar form morphologies on Fraser River are familiar and consistent with those in gravelly braided channels but the wandering style does exhibit some distinctive features. For example, 65-year histories reveal the potential for long sequences of uninterrupted accretion in relatively stable wandering rivers that are unlikely in braided rivers.     

Are catch-up reefs an artefact of coring?

Drill cores through modern coral reefs commonly show a time lag in reef initiation followed by a phase of rapid accretion to sea level from submerged foundations – the so-called ‘catch-up response’. But because of the difficulty of drilling in these environments, core distribution is usually restricted to accessible areas that may not fully represent reef history, especially if the reef initiated in patches or developed with a prograde or retrograde geometry. As a consequence, core data have the potential to give a misleading impression of reef development, particularly with respect to the timing of initiation and response to sea-level rise. Here, we use computer models to simulate keep-up reef development and, from them, quantify variations in the timing of reef initiation and accretion rate using mock cores taken through the completed simulations. The results demonstrate that cores consistently underestimate the timing of reef initiation and overestimate the reef accretion rate so that, statistically, a core through a keep-up reef will most likely produce a catch-up pattern – an initiation lag followed by a phase of rapid accretion to sea level. This implies that catch-up signatures may be an artefact of coring and that keep-up reefs are significantly more common than previous core studies claim.     

Magneto-rotational instability in the accreting envelope of a protostar and the formation of the large-scale magnetic field

We investigate the role of the magnetic field in the collapse of a gas-dust cloud into a massive gravitating object. Observations of one such object (G31.41) indicate that the magnetic field has an hourglass shape oriented along the rotation axis of the matter, due to the freezing-in of the magnetic-field lines in the accreting matter. It is believed that accretion in stellar disks is associated with the transport of angular momentum from the center to the periphery, which could be initiated by large-scale vortex structures arising in the presence of unstable rotational flows of matter. The numerical simulations have established that the equilibrium configuration of a gas-dust disk rotating in a spherically symmetrical gravitational potential is subject to the development of strong instability in the presence of a weak magnetic field. It is shown that the development of instability leads to a transport of angular momentum to the disk periphery by large-scale vortex structures, together with the accretion of matter onto the gravitating object. The magnetic-field lines near the equator take on a chaotic character, but an hourglass configuration is observed near the rotation axis, in agreement with observations.     

Ring structures in Orion KL

We present the results of studies of the superfine structure of H₂O maser sources in the Orion Nebula. Powerful, low-velocity, compact maser sources are distributed in eight active zones. Highly organized structures in the form of chains of compact components were revealed in two of these, in the molecular cloud OMC-1. The component sizes are ～0.1 AU and their brightness temperatures are T _b =10¹²?10¹⁶ K. The structures correspond to tangential sections of concentric rings viewed edge-on. The ring emission is concentrated in the azimuthal plane, decreasing the probability of their discovery. The formation of protostars is accompanied by the development of accretion disks and bipolar flows, with associated H₂O maser emission. The accretion disks are in the stage of fragmentation into protoplanetary rings. In a Keplerian approximation, the protostars have low masses, possibly evidence for instability of the systems. Supermaser emission of the rings is probably triggered by precession of the accretion disk. The molecular cloud’s radial velocity is V _LSR=7.74 km/s and its optical depth is τ≈5. The emission from components with velocities within the maser window is additionally amplified. The components’ emission is linearly polarized via anisotropic pumping.     

Formation of accretion disks in close-binary systems with magnetic fields

We have developed a three-dimensional numerical model and applied it to simulate plasma flows in semi-detached binary systems whose accretor possesses a strong intrinsic magnetic field. The model is based on the assumption that the plasma dynamics are determined by the slow mean flow, which forms a backdrop for the rapid propagation of MHD waves. The equations describing the slow motion of matter were obtained by averaging over rapidly propagating pulsations. The numerical model includes the diffusion of magnetic field by current dissipation in turbulent vortices, magnetic buoyancy, and wave MHD turbulence. A modified three-dimensional, parallel, numerical code was used to simulate the flow structure in close binary systems with various accretor magnetic fields, from 10⁵ to 10⁸ G. The conditions for the formation of the accretion disk and the criteria distinguishing the two types of flow corresponding to intermediate polars and polars are discussed.     

Features of the accretion in the EX Hydrae system: Results of numerical simulation

A two-dimensional numerical model in the axisymmetric approximation that describes the flow structure in the magnetosphere of the white dwarf in the EX Hya system has been developed. Results of simulations show that the accretion in EX Hya proceeds via accretion columns, which are not closed and have curtain-like shapes. The thickness of the accretion curtains depends only weakly on the thickness of the accretion disk. This thickness developed in the simulations does not agree with observations. It is concluded that the main reason for the formation of thick accretion curtains in the model is the assumption that the magnetic field penetrates fully into the plasma of the disk. An analysis based on simple estimates shows that a diamagnetic disk that fully or partially shields the magnetic field of the star may be a more attractive explanation for the observed features of the accretion in EX Hya.     

Excitation of electromagnetic lower hybrid instability in the region of bouncing electron beams

The region of bouncing electron beams in the earth’s magnetosphere can be unstable against a non-resonant electromagnetic lower hybrid instability. The instability is purely growing in the rest frame of the plasma, and can be excited either by the temperature anisotropy or the drift velocity of the bouncing electron beams. The growth rates of the instability decrease with the increase of cold electron density. Consequently the growth rate is maximum at the equator where the cold electron density is minimum. The intense turbulence generated by this instability could broaden the bouncing electron beams thereby explaining the observed wider cone width of the beams at the equator. The instability could generate magnetic pulsations in the frequency range of orderPc ₁?Pc₃ with typical wavelength ≈ (3–10) km in the, magnetosphere during magnetic storms or substorms.     

A new look at anomalous X-ray Pulsars

We explore the possibility of explaining Anomalous X-ray Pulsars (AXPs) and Soft Gammaray Repeaters (SGRs) in a scenario with fall-back magnetic accretion onto a young isolated neutron star. The X-ray emission of the pulsar in this case originates due to the accretion of matter onto the surface of the neutron star from a magnetic slab surrounding its magnetosphere. The spin-down rate of the neutron star expected in this picture is close to the observed value. We show that such neutron stars are relatively young and are going through the transition from the propeller state to the accretor state. The pulsar’s activity in gamma-rays is connected with its relative youth, and is enabled by energy stored in a non-equilibrium layer located in the crust of the low-mass neutron star. This energy can be released due to the mixing of matter in the neutron star crust with super heavy nuclei approaching its surface and becoming unstable. The fission of nuclei in the low-density region initiates chain reactions leading to a nuclear explosion. Outbursts are probably triggered by instability developing in the region where the matter accreted by the neutron star accumulates in the magnetic polar regions.     

Neutronization of matter in a stellar core and convection during gravitational collapse

The roles of neutrinos and convective instability in collapsing supernovae are considered. Spherically symmetrical computations of the collapse using the Boltzmann equation for the neutrinos lead to the formation of the condition of convective instability, \({\left( {\frac{{\partial P}}{{\partial s}}} \right)_{\rho {Y_l}}}\frac{{ds}}{{dr}} + {\left( {\frac{{\partial P}}{{\partial {Y_L}}}} \right)_{\rho s}}\frac{{d{Y_L}}}{{dr}}\) < 0, in a narrow region of matter accretion above the neutrinosphere. If instability arises in this region, the three-dimensional solution will represent a correction to the spherically symmetrical solution for the gravitational collapse. The mean neutrino energies change only negligibly in the narrow region of accretion. Nuclear statistical equilibrium is usually assumed in the hot proto-neutron stellar core, to simplify the computations of the collapse. Neutronization with the participation of free neutrons is most efficient. However, the decay of nuclei into nucleons is hindered during the collapse, because the density grows too rapidly compared to the growth in the temperature, and an appreciable fraction of the energy is carried away by neutrinos. The entropy of the matter per nucleon is modest at the stellar center. All the energy is in degenerate electrons during the collapse. If the large energy of these degenerate electrons is taken into account, neutrons are efficiently formed, even in cool matter with reduced Y_e (the difference between the numbers of electrons and positrons per nucleon). This process brings about an increase in the optical depth to neutrinos, the appearance of free neutrons, and an increase in the entropy per nucleon at the center. The convectively unstable region at the center increases. The development of large-scale convection is illustrated using a multi-dimensional gas-dynamical model for the evolution of a stationary, unstable state (without taking into account neutrino transport). The time for the development of convective instability (several milliseconds) does not exceed the time for the existence of the unstable region at the center (10ms). The realization of this type of instability is fundamentally different from a spherically symmetrical model. The flux of neutrinos changes and the mean energy of the neutrinos is increased, which has important implications for the detection of neutrinos from supernovae. For these same reasons, the energy absorped in the supernova envelope also changes in the transition to such a multi-dimensional model.     

The influence of vegetation on turbulence and flow velocities in European salt-marshes

Flow hindrance by salt‐marsh vegetation is manifested in the structure of the tidal current; it has a significant impact on sediment transport and it has been related to increased sediment accretion. The flow characteristics in three different vegetation types (Spartina maritima, Sp. anglica and Salicornia sp./Suaeda maritima) were measured on three salt‐marshes in Portugal and England. These in situ measurements differ from laboratory flume experiments with ‘clean’ vegetation by the complexity of natural canopies. Skimming flow develops above the Spartina canopy when the vegetation is fully submerged. In this situation, a low turbulence zone with nearly constant velocity in the denser canopy is separated from the skimming flow above by an interface characterized by high Reynolds stresses. In the low turbulence zone, a positive relationship exists between turbulence intensity and shoot density, which is due to wake turbulence generated locally in the canopy. The rate of particle settling should be increased in that zone. The lower limit of skimming flow is best predicted by the height within the canopy that includes 85% of the biomass. For emergent Spartina canopies and the short Salicornia/Suaeda marsh, the maximal velocity‐gradient is shifted upwards compared to a standard boundary layer over bare sediment and the turbulence is attenuated near the bed, but to a lesser extent than for fully submerged Spartina canopies. A turbulence reduction near the bed was observed in all measured profiles; that should enhance sediment deposition and protects the bed against subsequent erosion.     

太阳风中磁流体湍流的特征和本质

在空间探测的早期,太阳风湍流现象就已经被空间飞船的观测发现了。但是太阳风湍流不同于过去学术界熟悉的流体湍流。虽然国际学术界已进行了长期的研究,但太阳风湍流的本质仍不能得到明确的认识,许多观测现象都得不到解释。涂传诒等^［1～5］研究了这些现象,促进了理论上的新进展。这些研究指出,太阳风起伏中存在着非线性湍流相互作用,但不是完全发展的湍流,而是正处于由阿尔芬波向着完全发展的湍流过渡的状态。这一新的概念把阿尔芬波传播理论与磁流体湍流理论结合起来,导致了理论上的新进展,从而揭示了太阳风湍流的本质,阐明了一系列过去不能解释的观测现象。着重介绍了这些研究的成果和意义。     

Periodic instabilities in protostellar disks with azimuthal magnetic fields

The stability of magnetohydrodynamic oscillations in a protostellar disk with a toroidal magnetic field is analyzed. It is shown that, apart from the aperiodic magnetorotational instability, two other types of periodic instabilities of non-axisymmetric perturbations can exist. The simultaneous presence of azimuthal and vertical components of the wave vector are necessary for these to exist. One instability is due to the inductive winding-up of the azimuthal magnetic field of the wave, and the other arises when the field amplitude is increased by a comoving Hall wave, transferring magnetic field into a region of enhanced field intensity. The bandwidths of the unstable wave numbers are analyzed as a function of the Hall current, the β parameter of a plasma, and the angle between the direction of wave propagation and the plane of the disk. Regions in the accretion disks typical of T Tauri stars are indentified where these instabilities could be most active.