Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe interaction of planets with a disc with MHD turbulence III: Flow morphology and conditions for gap formation in local and global simulations
104
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present the results of both global cylindrical disc simulations and local shearing box simulations of protoplanets interacting with a disc undergoing MHD turbulence with zero net flux magnetic fields. We investigate the nature of the disc response and conditions for gap formation. This issue is an important one for determining the type and nature of the migration of the protoplanet, with the presence of a deep gap being believed to enable slower migration. For both types of simulation we find a common pattern of behaviour for which the main parameter determining the nature of the response is $M_p R^3/(M_* H^3)$, with $M_p$, $M_*$, $R$, and $H$ being the protoplanet mass, the central mass, the orbital radius and the disc semi-thickness respectively. We find that as this parameter is increased towards 0.1, the presence of the protoplanet is first indicated by the appearance of the well known trailing wake which, although it may appear erratic on account of the turbulence, appears to be well defined. Once the above parameter exceeds a number around unity a gap starts to develop inside which the magnetic energy density tends to be concentrated in the high density wakes. This gap formation condition can be understood from simple dimensional considerations of the conditions for nonlinearity, and the balance of angular momentum transport due to Maxwell and Reynolds stresses with that due to tidal torques. An important result is that the basic flow morphology in the vicinity of the protoplanet is very similar in both the local and global simulations. This indicates that local shearing box simulations, which are computationally less demanding, capture much of the physics of disc-planet interaction. Thus they may provide a useful tool for studying the local interaction between forming protoplanets and turbulent, protostellar discs.
rate research
Read More
(Abridged) We present global disc and local shearing box simulations of planets interacting with a MHD turbulent disc. We examine the torque exerted by the disc on the embedded planets as a function of planet mass, and thus make a first study of orbital migration of planets due to interaction with turbulent discs. Global simulations were performed for a disc with H/R=0.07 and planet masses M_p=3,10,30 Earth masses, and 3 Jupiter masses. Shearing box runs were performed for values of (M_p/M_*)/(H/R)^3=0.1,0.3,1.0 and 2.0, M_* being the central mass. These allow embedded and gap forming planets to be examined. In all cases the instantaneous torque exerted on a planet showed strong fluctuations. In the embedded cases it oscillated between negative and positive values, and migration occurs as a random walk, unlike the usual type I migration. Running time averages for embedded planets over 20-25 orbital periods show that large fluctuations occur on longer time scales, preventing convergence of the average torque to well defined values, or even to a well defined sign. Fluctuations become relatively smaller for larger masses, giving better convergence, due to the planets perturbation of the disc becoming larger than the turbulence in its vicinity. Eventually gap formation occurs, with a transition to type II migration. The existence of significant fluctuations occurring in turbulent discs on long time scales is important for lower mass embedded protoplanets. If significant fluctuations occur on the longest disc evolutionary time scales, convergence of torque running averages for practical purposes will not occur, and the migration behaviour of low mass protoplanets considered as an ensemble would be very different from predictions of type I theory for laminar discs.
We present a global MHD simulation of a turbulent accretion disc interacting with a protoplanet of 5 Jupiter masses. The disc model had H/r=0.1,and a value of the Shakura & Sunyaev alpha ~ 0.005. The protoplanet opened a gap in the disc, with the interaction leading to inward migration on the expected time scale. Spiral waves were launched by the protoplanet and although they were diffused and dissipated through interaction with the turbulence, they produced an outward angular momentum flow which compensated for a reduced flux associated with the turbulence, so maintaining the gap. When compared with laminar disc models with the same estimated alpha, the gap was found to be deeper and wider indicating that the turbulent disc behaved as if it possessed a smaller alpha. This may arise for two reasons. First, the turbulence does not provide a source of constantly acting friction in the near vicinity of the planet that leads to steady mass flow into the gap region. Instead the turbulence is characterised by large fluctuations in the radial velocity, and time averaging over significant time scales is required to recover the underlying mass flow through the disc. Near the planet the disc material experiences high amplitude perturbations on time scales that are short relative to the time scale required for averaging. The disc response is thus likely to be altered relative to a Navier--Stokes model. Second, the simulation indicates that an ordered magnetic connection between the inner and outer disc can occur enabling angular momentum to flow out across the gap, helping to maintain it independently of the protoplanets tide. This type of effect may assist gap formation for smaller mass protoplanets which otherwise would not be able to maintain them.
Abridged: We study the properties of clumps formed in three-dimensional weakly magnetized magneto-hydrodynamic simulations of converging flows in the thermally bistable, warm neutral medium (WNM). We find that: (1) Similarly to the situation in the classical two-phase medium, cold, dense clumps form through dynamically-triggered thermal instability in the compressed layer between the convergent flows, and are often characterised by a sharp density jump at their boundaries though not always. (2) However, the clumps are bounded by phase-transition fronts rather than by contact discontinuities, and thus they grow in size and mass mainly by accretion of WNM material through their boundaries. (3) The clump boundaries generally consist of thin layers of thermally unstable gas, but these layers are often widened by the turbulence, and penetrate deep into the clumps. (4) The clumps are approximately in both ram and thermal pressure balance with their surroundings, a condition which causes their internal Mach numbers to be comparable to the bulk Mach number of the colliding WNM flows. (5) The clumps typically have mean temperatures 20 < T < 50 K, corresponding to the wide range of densities they contain (20 < n < 5000 pcc) under a nearly-isothermal equation of state. (6) The turbulent ram pressure fluctuations of the WNM induce density fluctuations that then serve as seeds for local gravitational collapse within the clumps. (7) The velocity and magnetic fields tend to be aligned with each other within the clumps, although both are significantly fluctuating, suggesting that the velocity tends to stretch and align the magnetic field with it. (8) The typical mean field strength in the clumps is a few times larger than that in the WNM. (9) The magnetic field strength has a mean value of B ~ 6 mu G ...
This chapter reviews the nature of turbulence in the Galactic interstellar medium (ISM) and its connections to the star formation (SF) process. The ISM is turbulent, magnetized, self-gravitating, and is subject to heating and cooling processes that control its thermodynamic behavior. The turbulence in the warm and hot ionized components of the ISM appears to be trans- or subsonic, and thus to behave nearly incompressibly. However, the neutral warm and cold components are highly compressible, as a consequence of both thermal instability in the atomic gas and of moderately-to-strongly supersonic motions in the roughly isothermal cold atomic and molecular components. Within this context, we discuss: i) the production and statistical distribution of turbulent density fluctuations in both isothermal and polytropic media; ii) the nature of the clumps produced by thermal instability, noting that, contrary to classical ideas, they in general accrete mass from their environment; iii) the density-magnetic field correlation (or lack thereof) in turbulent density fluctuations, as a consequence of the superposition of the different wave modes in the turbulent flow; iv) the evolution of the mass-to-magnetic flux ratio (MFR) in density fluctuations as they are built up by dynamic compressions; v) the formation of cold, dense clouds aided by thermal instability; vi) the expectation that star-forming molecular clouds are likely to be undergoing global gravitational contraction, rather than being near equilibrium, and vii) the regulation of the star formation rate (SFR) in such gravitationally contracting clouds by stellar feedback which, rather than keeping the clouds from collapsing, evaporates and diperses them while they collapse.
The stresses produced by magnetorotational turbulence can provide effective angular momentum transport in accretion disks. However, questions remain about the ability of simulated disks to reproduce observationally inferred stress-to-gas-pressure ratios. In this paper we present a set of high resolution global magnetohydrodynamic disk simulations which are initialised with different field configurations: purely toroidal, vertical field lines, and nested poloidal loops. A mass source term is included which allows the total disk mass to equilibrate in simulations with long run times, and also enables the impact of rapid mass injection to be explored. Notably different levels of angular momentum transport are observed during the early-time transient disk evolution. However, given sufficient time to relax, the different models evolve to a statistically similar quasi-steady state with a stress-to-gas-pressure ratio, $alpha sim 0.032-0.036$. The indication from our results is that {it steady, isolated} disks may be unable to maintain a large-scale magnetic field or produce values for the stress-to-gas-pressure ratio implied by some observations. Supplementary simulations exploring the influence of trapping magnetic field, injecting vertical field, and rapidly injecting additional mass into the disk show that large stresses ($alpha sim 0.1-0.25$) can be induced by these mechanisms. The simulations highlight the common late-time evolution and characteristics of turbulent disks for which the magnetic field is allowed to evolve freely. If the boundaries of the disk, the rate of injection of magnetic field, or the rate of mass replenishment are modified to mimic astrophysical disks, markedly different disk evolution occurs.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
