No Arabic abstract
Radial transport of particles, elements and fluid driven by internal stresses in three-dimensional (3D) astrophysical accretion disks is an important phenomenon, potentially relevant for the outward dust transport in protoplanetary disks, origin of the refractory particles in comets, isotopic equilibration in the Earth-Moon system, etc. To gain better insight into these processes, we explore the dependence of meridional circulation in 3D disks with shear viscosity on their thermal stratification, and demonstrate strong effect of the latter on the radial flow. Previous locally isothermal studies have normally found a pattern of the radial outflow near the midplane, switching to inflow higher up. Here we show, both analytically and numerically, that a flow, which is inward at all altitudes, is possible in disks with entropy and temperature steeply increasing with height. Such thermodynamic conditions may be typical in the optically thin, viscously heated accretion disks. Disks in which these conditions do not hold should feature radial outflow near the midplane, as long as their internal stress is provided by the shear viscosity. Our results can also be used for designing hydrodynamical disk simulations with a prescribed pattern of the meridional circulation.
In standard accretion discs, outward angular momentum transfer by viscous forces is compensated by the inward motion of the accreting matter. However, the vertical structure of real accretion discs leads to meridional circulation with comparable amplitudes of poloidal velocities. Using thin-disc approximation, we consider different regimes of disc accretion with different vertical viscosity scalings. We show that, while gas-pressure-dominated discs can easily have a midplane outflow, standard thin radiation-pressure-dominated disc is normally moving inwards at all the heights. However, quasi-spherical scaling for pressure ($ppropto varpi^{-5/2}$) leads to a midplane outflow for a very broad range of parameters. It particular, this may lead to a reversed, outward heat advection in geometrically thick discs when the temperature decreases rapidly enough with height. While the overall direction of heat advection depends on the unknown details of vertical structure and viscosity mechanisms, existence of the midplane counterflow in quasi-spherical flows is a robust result weakly dependent on the parameters and the assumptions of the model. Future models of thick radiatively inefficient flows should take meridional circulation into account.
We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies citep{mun19}, we consider circular binaries with various mass ratios ($0.1leq q_{rm{b}}leq1$) and study accretion from ``infinite, steady-supply disks and from finite-sized, viscously spreading tori. For ``infinite disks, a global steady state can be reached, and the accretion variability has a dominant frequency ${sim}0.2Omega_{rm{b}}$ for $q_{rm{b}}>0.5$ and $Omega_{rm{b}}$ for $q_{rm{b}}<0.5$, ($Omega_{rm{b}}$ is the binary angular frequency). We find that the accretion ``eigenvalue $l_0$ -- the net angular momentum transfer from the disk to the binary per unit accreted mass -- is always positive and falls in the range ($0.65$-$0.85)a_{rm b}^2Omega_{rm{b}}$ (with $a_{rm{b}}$ the binary separation), depending weakly on the mass ratio and viscosity. This leads to binary expansion when $q_{rm{b}}gtrsim0.3$. Accretion from a finite torus can be separated into two phases: an initial transient phase, corresponding to the filling of the binary cavity, followed by a viscous pseudo-stationary phase, during which the torus viscously spreads and accretes onto the binary. In the viscous phase, the net torque on the binary per unit accreted mass is close to $l_0$, the value derived for ``infinite disks. We conclude that similar-mass binaries accreting from CBDs gain angular momentum and expand over long time scales. This result significantly impacts the coalescence of supermassive binary black holes and newly formed binary stars. We offer a word of caution against conclusions drawn from simulations of transient accretion onto empty circumbinary cavities.
Surface observations indicate that the speed of the solar meridional circulation in the photosphere varies in anti-phase with the solar cycle. The current explanation for the source of this variation is that inflows into active regions alter the global surface pattern of the meridional circulation. When these localized inflows are integrated over a full hemisphere, they contribute to the slow down of the axisymmetric poleward horizontal component. The behavior of this large scale flow deep inside the convection zone remains largely unknown. Present helioseismic techniques are not sensitive enough to capture the dynamics of this weak large scale flow. Moreover, the large time of integration needed to map the meridional circulation inside the convection zone, also masks some of the possible dynamics on shorter timescales. In this work we examine the dynamics of the meridional circulation that emerges from a 3D MHD global simulation of the solar convection zone. Our aim is to assess and quantify the behavior of meridional circulation deep inside the convection zone, where the cyclic large-scale magnetic field can reach considerable strength. Our analyses indicate that the meridional circulation morphology and amplitude are both highly influenced by the magnetic field, via the impact of magnetic torques on the global angular momentum distribution. A dynamic feature induced by these magnetic torques is the development of a prominent upward flow at mid latitudes in the lower convection zone that occurs near the equatorward edge of the toroidal bands and that peaks during cycle maximum. Globally, the dynamo-generated large-scale magnetic field drives variations in the meridional flow, in stark contrast to the conventional kinematic flux transport view of the magnetic field being advected passively by the flow.
A key component of solar interior dynamics is the meridional circulation (MC), whose poleward component in the surface layers has been well observed. Time-distance helioseismic studies of the deep structure of MC, however, have yielded conflicting inferences. Here, following a summary of existing results we show how a large center-to-limb systematics (CLS) in the measured travel times of acoustic waves affect the inferences through an analysis of frequency dependence of CLS, using data from the Helioseismic and Doppler Imager (HMI) onboard Solar Dynamics Observatory (SDO). Our results point to the residual systematics in travel times as a major cause of differing inferences on the deep structure of MC.
Recently, the vertical shear instability (VSI) has become an attractive purely hydrodynamic candidate for the anomalous angular momentum transport required for weakly ionized accretion disks. In direct three-dimensional numerical simulations of VSI turbulence in disks, a meridional circulation pattern was observed that is opposite to the usual viscous flow behavior. Here, we investigate whether this feature can possibly be explained by an anisotropy of the VSI turbulence. Using three-dimensional hydrodynamical simulations, we calculate the turbulent Reynolds stresses relevant for angular momentum transport for a representative section of a disk. We find that the vertical stress is significantly stronger than the radial stress. Using our results in viscous disk simulations with different viscosity coefficients for the radial and vertical direction, we find good agreement with the VSI turbulence for the stresses and meridional flow; this provides additional evidence for the anisotropy. The results are important with respect to the transport of small embedded particles in disks.