The number of extrasolar planet discoveries has increased dramatically over the last 15 years. Nearly 700 exoplanets have currently been observed through a variety of observation techniques. Most of the currently documented exoplanets differ greatly from the planets in our own Solar System, with various combinations of eccentric orbits, short orbital periods, and masses many times that of Jupiter. More recently, planets belonging to a new class of `distant gas giants' have also been discovered with orbits of 30 to 100 times that of Jupiter. The wide variety of different planet formation outcomes stem from a complex interplay between gravitational interactions, hydrodynamic interactions and competitive accretion among the planets that is not yet fully understood. Simulations performed using a series of modifications to an existing, widely used hydrodynamic code (FARGO) are presented. The main goal is to develop a more rigorous and robust gas accretion scheme that is valid and consistent for the ranges of exolanetary gas giant masses, eccentricities and semimajor axes that have been observed to better understand the mechanisms involved in their formation. The resulting scheme is a more robust and accurate prescription for gas accretion onto planetary cores in a manner that is mostly resolution independent and valid over a large range of masses (less than an Earth mass to multiple Jupiter masses). The modified scheme accounts for multiple, competing, dynamic accretion mechanisms (including atmospheric effects) and their associated time scales between an arbitrary number of protoplanets. This updated accretion scheme provides a means for exploring the entire formation process of gas giants out of a variety of initial conditions in a self-consistent manner. The modifications made to the code as well as simulation results will be discussed and explored.

In the case of cold accretion disk, coupling between charge neutral gas and magnetic field is too weak such that the magneto-rotational instability will be less effective or even stop working. In such a situation it is of prime interest to investigate the pure hydrodynamic turbulence and transport phenomenon. As the Reynolds number increases, the relative importance of the non-linear term in the hydrodynamic equation increases and in the case of accretion disk where molecular viscosity is too small the Reynolds number is large enough for the non-linear term to bring new effects. We investigate a scenario, the ‘weakly non-linear’ evolution of amplitude of linear mode when the flow is bounded by two parallel walls. The unperturbed flow is similar to plane Couette flow but with Coriolis force included in the hydrodynamic equation. Although there is no exponentially growing eigenmode, due to self-interaction the least stable eigenmode will grow in an intermediate phase. Later on this will lead to higher order non-linearity and plausible turbulence. Although the non-linear term in the hydrodynamic equation is energy conserving, within the weakly non-linear analysis it is possible to define a lower bound of the energy needed for flow to transform to turbulent phase. Such an unstable phase is possible only if the Reynolds number ≥ 103−4. In Chapter-2 we set up equation of amplitude for the hydrodynamic perturbation and study the effect of weak non-linear evolution of linear mode for general angular momentum distribution, where Keplerian disk is obtained as a special case. As we know that to explain observed hard X-rays the choice of Keplerian angular momentum profile is not adequate, we consider the sub-Keplerian regime of the disk. In Chapter-3 we assume that the cooling mechanism is dominated by bremsstrahlung process (without any strict knowledge of the magnetic field structure).We show that in a range of Shakura-Sunyaev viscosity 0.2 ≥ α ≥ 0.0005, flow behavior varies widely, particularly by means of the size of disk, efficiency of cooling and corresponding temperatures of ions and electrons. We also show that the disk around a rotating black hole is hotter compared to that around a Schwarzschild black hole, rendering a larger difference between ion and electron temperatures in the former case. We finally reproduce the observed luminosities(L) of two extreme cases—the under-fed AGNs and quasars and ultra-luminous X-ray sources at different combinations of mass accretion rate, ratio of specific heats, Shakura-Sunyaev viscosity parameter and Kerr parameter. In Chapter-4 we investigate the viscous two temperature accretion disk flows around rotating blackholes. We describe the global solution of accretion flows, unlike that in Chapter-3, with a sub-Keplerian angular momentum profile, by solving the underlying conservation equations including explicit cooling processes self-consistently. Bremsstrahlung, synchrotron and inverse comptonization of soft photons are considered as possible cooling mechanisms. We focus on the set of solutions for sub-Eddington, Eddington and super-Eddington mass accretion rates around Schwarzschild and Kerr black holes with a Kerr parameter 0.998. We analyse various phases of advection–general advective paradigm to radiatively inefficient paradigm. The solution may potentially explain the hard X-rays and γ-rays emitted from AGNs and X-ray binaries. We also compare the solutions for two different regimes of viscosity. We finally reproduce the observed luminosities of the under-fed AGNs and quasars, ultra-luminous X-ray sources at different combinations of input parameters such as mass accretion rate and ratio of specific heats.

京都大学 / 0048 / 新制・課程博士 / 博士(理学) / 甲第21572号 / 理博第4479号 / 新制||理||1643(附属図書館) / 京都大学大学院理学研究科物理学・宇宙物理学専攻 / (主査)准教授 野上 大作, 教授 長田 哲也, 教授 嶺重 慎 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DGAM

京都大学 / 0048 / 新制・課程博士 / 博士(理学) / 甲第22251号 / 理博第4565号 / 新制||理||1655(附属図書館) / 京都大学大学院理学研究科物理学・宇宙物理学専攻 / (主査)教授 嶺重 慎, 准教授 前田 啓一, 教授 長田 哲也 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DFAM

No description available.

Smoothed Particle Hydrodynamics (SPH) is analysed as the weighted residual method. In particular the analysis focuses on the collocation aspect of the method. Using Monte Carlo experiments we demonstrate that SPH is highly sensitive to node disorder, especially in its symmetrised energy and momentum conserving form. This aspect of the method is related to low [Beta] MHD instabilities observed by other authors. A remedy in the form of the Weighted Differences Method is suggested, which addresses this problem to some extent, but at a cost of losing automatic conservation of energy and momentum. ¶ The Weighted Differences Method is used to simulate propagation of Alfven and magnetosonic wave fronts in [Beta] = 0 plasma, and the results are compared with data obtained with the NCSA Zeus3D code with the Method of Characteristics (MOC) module. ¶ SPH is then applied to two interesting astrophysical situations: accretion on to a white dwarf in a compact binary system, which results in a formation of an accretion disk, and gravitational collapse of a magnetised vortex. Both models are 3 dimensional. ¶ The accretion disk which forms in the binary star model is characterised by turbulent flow: the Karman vortex street is observed behind the stream-disk interaction region. The shock that forms at the point of stream-disk interaction is controlled by the means of particle merges, whereas Monaghan-Lattanzio artificial viscosity is used to simulate Smagorinsky closure. ¶ The evolution of the collapsing magnetised vortex ends up in the formation of an expanding ring in the symmetry plane of the system. We observe the presence of spiralling inward motion towards the centre of attraction. That final state compares favourably with the observed qualitative and quantitative characteristics of the circumnuclear disk in the Galactic Centre. That simulation has also been verified with the NCSA Zeus3D run. ¶ In conclusions we contrast the result of our Monte Carlo experiments with the results delivered by our production runs. We also compare SPH and Weighted Differences against the new generation of conservative finite differences methods, such as the Godunov method and the Piecewise Parabolic Method. We conclude that although SPH cannot match the accuracy and performance of those methods, it appears to have some advantage in simulation of rotating flows, which are of special interest to astrophysics.

No description available.