We studied the effects of rotation on standing accretion shock instability (SASI) by performing three-dimensional hydrodynamics simulations. Taking into account a realistic equation of state and neutrino heating/cooling, we prepared a spherically symmmetric and steady accretion flow through a standing shock wave onto a proto-neutron star (PNS). When the SASI entered the nonlinear phase, we imposed uniform rotation on the flow advecting from the outer boundary of the iron core, whose specific angular momentum was assumed to agree with recent stellar evolution models. Using spherical harmonics in space and Fourier decompositions in time, we performed mode analysis of the nonspherical deformed shock wave to observe rotational effects on the SASI in the nonlinear phase. We found that rotation imposed on the axisymmetric SASI did not make any spiral modes and hardly affected sloshing modes, except for steady l=2, m=0 modes. In contrast, rotation imposed on the non-axisymmetric flow increased the amplitude of spiral modes so that some spiral flows accreting on the PNS were more clearly formed inside the shock wave than without rotation. The amplitudes of spiral modes increased significantly with rotation in the progressive direction.

We consider the effect of a supernova (SN) explosion in a very massive binary that is expected to form in a portion of Population III stars with the mass higher than 100$M_\odot$. In a Population III binary system, a more massive star can result in the formation of a BH and a surrounding accretion disc. Such BH accretion could be a significant source of the cosmic reionization in the early universe. However, a less massive companion star evolves belatedly and eventually undergoes a SN explosion, so that the accretion disc around a BH might be blown off in a lifetime of companion star. In this paper, we explore the dynamical impact of a SN explosion on an accretion disc around a massive BH, and elucidate whether the BH accretion disc is totally demolished or not. For the purpose, we perform three-dimensional hydrodynamic simulations of a very massive binary system, where we assume a BH of $10^3 M_{\odot}$ that results from a direct collapse of a very massive star and a companion star of $100 M_{\odot}$ that undergoes a SN explosion. We calculate the remaining mass of a BH accretion disc as a function of time. As a result, it is found that a significant portion of gas disc can survive through three-dimensional geometrical effects even after the SN explosion of a companion star. Even if the SN explosion energy is higher by two orders of magnitude than the binding energy of gas disc, about a half of disc can be left over. The results imply that the Population III BH accretion disc can be a long-lived luminous source, and therefore could be an important ionizing source in the early universe.

We have studied non-axisymmetric standing accretion shock instability, or SASI, by 3D hydrodynamical simulations. This is an extention of our previous study on axisymmetric SASI. We have prepared a spherically symmetric and steady accretion flow through a standing shock wave onto a proto-neutron star, taking into account a realistic equation of state and neutrino heating and cooling. This unperturbed model is supposed to represent approximately the typical post-bounce phase of core-collapse supernovae. We then have added a small perturbation (~1%) to the radial velocity and computed the ensuing evolutions. Not only axisymmetric but non-axisymmetric perturbations have been also imposed. We have applied mode analysis to the non-spherical deformation of the shock surface, using the spherical harmonics. We have found that (1) the growth rates of SASI are degenerate with respect to the azimuthal index m of the spherical harmonics Y_l^m, just as expected for a spherically symmetric background, (2) nonlinear mode couplings produce only m=0 modes for the axisymmetric perturbations, whereas m=!0 modes are also generated in the non-axisymmetric cases according to the selection rule for the quadratic couplings, (3) the nonlinear saturation level of each mode is lower in general for 3D than for 2D because a larger number of modes are contributing to turbulence in 3D, (4) low l modes are dominant in the nonlinear phase, (5) the equi-partition is nearly established among different m modes in the nonlinear phase, (6) the spectra with respect to l obey power laws with a slope slightly steeper for 3D, and (7) although these features are common to the models with and without a shock revival at the end of simulation, the dominance of low l modes is more remarkable in the models with a shock revival.

We perform two-dimensional hydrodynamic simulation on the gas stripping by radiation drag from an interstellar cloud moving in uniform radiation fields. To properly include relativistic radiation drag, the radiation hydrodynamic equation is solved with taking into account the dilution of radiation fields by optical depth of the cloud. As a result, it is found that the optically-thin surface layers are effectively stripped by radiation drag from an optically-thick gas cloud, and simultaneously stripped gas loses momentum. The momentum loss time-scale is found to be on the order of $10^{8}$ years under intensive radiation fields which are expected in the early phase of galaxy evolution. The present results show that the radiation drag is an effective mechanism to extract angular momentum from interstellar medium and allows it to accrete onto the galactic center. The mass accretion driven by radiation drag may lead to the formation of a central supermassive black hole.

Two-dimensional numerical simulations of an accretion flow in a close binary system are performed by solving the Euler equations with radiative transfer. In the present study, the specific heat ratio is assumed to be constant while radiative cooling effect is included as a non-adiabatic process. The cooling effect of the disc is considered by discharging energy in the vertical directions from the top and bottom surfaces of the disc. We use the flux-limited diffusion approximation to calculate the radiative heat flux values. Our calculations show that a disc structure appears and the spiral shocks are formed on the disc. These features are similar to that observed in the case of an adiabatic gas with a lower specific heat ratio, $\gamma=1.01$. It is found that when radiative cooling effect is accounted for, the mass of the disc becomes larger than that assuming $\gamma=5/3$, and smaller than that assuming $\gamma=1.01$. It is concluded that employing an adiabatic gas with a lower specific heat ratio is almost a valid assumption for simulating accretion disc with radiative cooling effect.