We present a solution to the long outstanding meter barrier problem in planet formation theory. As solids spiral inward due to aerodynamic drag, they will enter disk regions that are characterized by high temperatures, densities, and pressures. High
partial pressures of rock vapor can suppress solid evaporation, and promote collisions between partially molten solids, allowing rapid growth. This process should be ubiquitous in planet-forming disks, which may be evidenced by the abundant class of Systems with Tightly-packed Inner Planets (STIPs) discovered by the NASA Kepler mission.
Using analytic arguments and numerical simulations, we examine whether chondrule formation and the FU Orionis phenomenon can be caused by the burst-like onset of gravitational instabilities (GIs) in dead zones. At least two scenarios for bursting dea
d zones can work, in principle. If the disk is on the verge of fragmention, GI activation near $rsim4$ to 5 AU can produce chondrule-forming shocks, at least under extreme conditions. Mass fluxes are also high enough during the onset of GIs to suggest that the outburst is related to an FU Orionis phenomenon. This situation is demonstrated by numerical simulations. In contrast, as supported by analytic arguments, if the burst takes place close to $rsim1$ AU, then even low pitch angle spiral waves can create chondrule-producing shocks and outbursts. We also study the stability of the massive disks in our simulations against fragmentation and find that although disk evolution is sensitive to changes in opacity, the disks we study do not fragment, even at high resolution and even for extreme assumptions.
It is generally thought that protoplanetary disks embedded in envelopes are more massive and thus more susceptible to gravitational instabilities (GIs) than exposed disks. We present three-dimensional radiative hydrodynamics simulations of protoplane
tary disks with the presence of envelope irradiation. For a disk with a radius of 40 AU and a mass of 0.07 Msun around a young star of 0.5 Msun, envelope irradiation tends to weaken and even suppress GIs as the irradiating flux is increased. The global mass transport induced by GIs is dominated by lower-order modes, and irradiation preferentially suppresses higher-order modes. As a result, gravitational torques and mass inflow rates are actually increased by mild irradiation. None of the simulations produce dense clumps or rapid cooling by convection, arguing against direct formation of giant planets by disk instability, at least in irradiated disks. However, dense gas rings and radial mass concentrations are produced, and these might be conducive to accelerated planetary core formation. Preliminary results from a simulation of a massive embedded disk with physical characteristics similar to one of the disks in the embedded source L1551 IRS5 indicate a long radiative cooling time and no fragmentation. The GIs in this disk are dominated by global two and three-armed modes.
In this dissertation, I describe theoretical and numerical studies that address the three-dimensional behavior of spiral shocks in protoplanetary disks and the controversial topic of gas giant formation by disk instability. For this work, I discuss c
haracteristics of gravitational instabilities (GIs) in bursting and asymptotic phase disks; outline a theory for the three-dimensional structure of spiral shocks, called shock bores, for isothermal and adiabatic gases; consider convection as a source of cooling for protoplanetary disks; investigate the effects of opacity on disk cooling; use multiple analyses to test for disk stability against fragmentation; test the sensitivity of GI behavior to radiation boundary conditions; measure shock strengths and frequencies in GI-bursting disks; evaluate temperature fluctuations in unstable disks; and investigate whether spiral shocks can form chondrules when GIs activate. The numerical methods developed for these studies are discussed, including a radiation transport routine that explicitly couples the low and high optical depth regimes and a routine that models ortho and parahydrogen. Finally, I explore the hypothesis that chondrule formation and the FU Ori phenomenon are driven by GI activation in dead zones.
Recent three-dimensional radiative hydrodynamics simulations of protoplanetary disks report disparate disk behaviors, and these differences involve the importance of convection to disk cooling, the dependence of disk cooling on metallicity, and the s
tability of disks against fragmentation and clump formation. To guarantee trustworthy results, a radiative physics algorithm must demonstrate the capability to handle both the high and low optical depth regimes. We develop a test suite that can be used to demonstrate an algorithms ability to relax to known analytic flux and temperature distributions, to follow a contracting slab, and to inhibit or permit convection appropriately. We then show that the radiative algorithm employed by Mejia (2004) and Boley et al. (2006) and the algorithm employed by Cai et al. (2006) and Cai et al. (2007, in prep.) pass these tests with reasonable accuracy. In addition, we discuss a new algorithm that couples flux-limited diffusion with vertical rays, we apply the test suite, and we discuss the results of evolving the Boley et al. (2006) disk with this new routine. Although the outcome is significantly different in detail with the new algorithm, we obtain the same qualitative answers. Our disk does not cool fast due to convection, and it is stable to fragmentation. We find an effective $alphaapprox 10^{-2}$. In addition, transport is dominated by low-order modes.
