Linear stability analysis of nonaxisymmetric instabilities in self-gravitating polytropic disks

• Main Author: Hadley, Kathryn Z., 1955-
• Language: en_US
• Published: University of Oregon 2011
• Subjects: Hydrodynamics; Stability; Polytropic disks; Astrophysics; Astronomy
• Online Access: http://hdl.handle.net/1794/11253

xvii, 371 p. : col. ill. === An important problem in astrophysics involves understanding the formation of planetary systems. When a star-forming cloud collapses under gravity its rotation causes it to flatten into a disk. Only a small percentage of the matter near the rotation axis falls inward to create the central object, yet our Sun contains over 99% of the matter of our Solar System. We examine how global hydrodynamic instabilities transport angular momentum through the disk causing material to accrete onto the central star. We analyze the stability of polytropic disks in the linear regime. A power law angular velocity of power q is imposed, and the equilibrium disk structure is found through solution of the time-independent hydrodynamic equations via the Hachisu self-consistent field method. The disk is perturbed, and the time-dependent linearized hydrodynamic equations are used to evolve it. If the system is unstable, the characteristic growth rate and frequency of the perturbation are calculated. We consider modes with azimuthal e im[varphi] dependence, where m is an integer and [varphi] is the azimuthal angle. We map trends across a wide parameter space by varying m , q and the ratios of the star-to-disk mass M * /M d and inner-to-outer disk radius r - /r + . We find that low m modes dominate for small r - /r + , increasing to higher r - /r + as M * /M d increases, independent of q . Three main realms of behavior are identified, for M * << M d , M * [approximate] M d and M * >> M d , and analyzed with respect to the I, J and P mode types as discussed in the literature. Analysis shows that for M * << M d , small r - /r + disks are dominated by low m I modes, which give way to high m J modes at high r - /r + . Low m J modes dominate M * [approximate] M d disks for small r - /r + , while higher m I modes dominate for high r - /r + . Behavior diverges with q for M * >> M d systems with high q models approximating M * [approximate] M d characteristics, while low q models exhibit m = 2 I modes dominating where r - /r + < 0.60. === Committee in charge: Raymond Frey, Chairperson; James Imamura, Advisor; Robert Zimmerman, Member; Paul Csonka, Member; Alan Rempel, Outside Member