Aims. Our goal is to find a converged cooling limit for fragmentation in self-gravitating disks. This is especially interesting for the formation of planets, brown dwarfs or stars and the growth of black holes. While investigating the limit, we want to give a clear criterion for the state of convergence. Methods. We run two-dimensional shearingsheet simulations with the hydrodynamic package Fosite at high resolutions. Thereby resolution and limiters are altered. Subsequently, we investigate the spectra of important physical quantities at the length scales where fragmentation occurs. In order to avoid prompt fragmentation at high resolutions we start these simulations with a fully developed gravitoturbulent state obtained at a lower resolution. Results. We show nearly converged results for fragmentation with a critical cooling timescale $t_{\mathrm{crit}} \sim 10\,\Omega^{-1}$ . We can backtrace this claim by investigating the spectra of relevant physical variables at length scales around and below the pressure scale height. We argue that well behaved results cannot be expected if counteracting quantities are varying too much on these critical length scales, either by change of resolution or numerical method. A comparison of fragmentation behaviour with the related spectra reveals that simulations behave similar, if the spectra are converged to the length scales where self-gravity leads to instabilities. Observable deviations in the results obtained with different numerical setup are confined to scales below these critical length scales.

Due to the non-axisymmetric potential of the central bar, barred spiral galaxies form, in addition to their characteristic arms and bar, a variety of structures within the thin gas disk, like nuclear rings, inner spirals and dust-lanes. These structures in the inner kiloparsec are most important to explain and understand the rate of black hole feeding. The aim of this work is to investigate the influence of stellar bars in spiral galaxies on the thin self-gravitating gas disk. We focus on the accretion of gas onto the central supermassive black hole and its time-dependent evolution. We conduct multi-scale simulations simultaneously resolving the galactic disk and the accretion disk around the central black-hole. We vary in all simulations the initial gas disk mass. As additional parameter we choose either the gas temperature for isothermal simulations or the cooling timescale in case of non-isothermal simulations. Accretion is either driven by a gravitationally unstable or clumpy accretion disk or by energy dissipation in strong shocks. Most simulations show a strong dependence of the accretion rate at the outer boundary of the central accretion disk ($r< 300~\mathrm{pc}$) on the gas flow at kiloparsec scales. The final black hole masses reach up to $\sim 10^9 M_\odot$ after $1.6~\mathrm{Gyr}$. Our models show the expected influence of the Eddington limit and a decline in growth rate at the corresponding sub-Eddington limit.

Due to the non-axisymmetric potential of the central bar, barred spiral galaxies form, in addition to their characteristic arms and bar, a variety of structures within the thin gas disk, like nuclear rings, inner spirals and dust-lanes. In this first of two papers, we present a method to accurately simulate the gas flow within the galactic plane in the 2D finite volume software package FOSITE, which solves the transport equations for mass, momentum and energy, and apply it to this class of objects. To this extent, we introduced a new transport scheme for angular momentum and a very efficient pseudo-spectral Poisson solver. Moreover, we provide a simple and generally applicable method of how to take care of gravity in the energy equation.

We investigate exact solutions for isothermal shock problems in different one-dimensional geometries. These solutions are given as analytical expressions if possible, or are computed using standard numerical methods for solving ordinary differential equations. We test the numerical solutions against the analytical expressions to verify the correctness of all numerical algorithms. We use similarity methods to derive a system of ordinary differential equations (ODE) yielding exact solutions for power law density distributions as initial conditions. Further, the system of ODEs accounts for implosion problems (IP) as well as explosion problems (EP) by changing the initial or boundary conditions, respectively. Taking genuinely isothermal approximations into account leads to additional insights of EPs in contrast to earlier models. We neglect a constant initial energy contribution but introduce a parameter to adjust the initial mass distribution of the system. Moreover, we show that due to this parameter a constant initial density is not allowed for isothermal EPs. Reasonable restrictions for this parameter are given. Both, the (genuinely) isothermal implosion as well as the explosion problem are solved for the first time.

Context. It is still an open issue whether a self-gravitating accretion disk fragments. There are many different physical and numerical explanations for fragmentation, but simulations often show a non-convergent behavior for ever better resolution. Aims. We investigate the influence of different numerical limiters in Godunov type schemes on the fragmentation boundary in self- gravitating disks. Methods. We compare the linear and non-linear outcome in two-dimensional shearingsheet simulations using the VANLEER and the SUPERBEE limiter. Results. We show that choosing inappropriate limiting functions to handle shock-capturing in Godunov type schemes can lead to an overestimation of the surface density in regions with shallow density gradients. The effect amplifies itself on timescales comparable to the dynamical timescale even at high resolutions. This is exactly the environment, where clumps are expected to form. The effect is present without, but scaled up by, self-gravity and also does not depend on cooling. Moreover it can be backtracked to a well known effect called oversteepening. If the effect is also observed in the linear case, the fragmentation limit is shifted to larger values of the critical cooling timescale.

Context: A substantial fraction of protoplanetary disks forms around stellar binaries. The binary system generates a time-dependent non-axisymmetric gravitational potential, inducing strong tidal forces on the circumbinary disk. This leads to a change in basic physical properties of the circumbinary disk, which should in turn result in unique structures that are potentially observable with the current generation of instruments. Aims: The goal of this study is to identify these characteristic structures, to constrain the physical conditions that cause them, and to evaluate the feasibility to observe them in circumbinary disks. Methods: To achieve this, at first two-dimensional hydrodynamic simulations are performed. The resulting density distributions are post-processed with a 3D radiative transfer code to generate re-emission and scattered light maps. Based on these, we study the influence of various parameters, such as the mass of the stellar components, the mass of the disk and the binary separation on observable features in circumbinary disks. Results: We find that the Atacama Large (sub-)Millimetre Array (ALMA) as well as the European Extremely Large Telescope (E-ELT) are capable of tracing asymmetries in the inner region of circumbinary disks which are affected most by the binary-disk interaction. Observations at submillimetre/millimetre wavelengths will allow the detection of the density waves at the inner rim of the disk and the inner cavity. With the E-ELT one can partially resolve the innermost parts of the disk in the infrared wavelength range, including the disk's rim, accretion arms and potentially the expected circumstellar disks around each of the binary components.

Aims: We investigate the suitability of {\alpha}-viscosity in self-similar models for self-gravitating disks with a focus on active galactic nuclei (AGN) disks. Methods: We use a self-similar approach to simplify the partial differential equations arising from the evolution equation, which are then solved using numerical standard procedures. Results: We find a self-similar solution for the dynamical evolution of self-gravitating {\alpha}-disks and derive the significant quantities. In the Keplerian part of the disk our model is consistent with standard stationary {\alpha}-disk theory, and self-consistent throughout the self-gravitating regime. Positive accretion rates throughout the disk demand a high degree of self-gravitation. Combined with the temporal decline of the accretion rate and its low amount, the model prohibits the growth of large central masses. Conclusions: {\alpha}-viscosity cannot account for the evolution of the whole mass spectrum of super-massive black holes (SMBH) in AGN. However, considering the involved scales it seems suitable for modelling protoplanetary disks. Keywords: accretion, accretion disks, turbulence, hydrodynamics, methods: analytical

The importance of contact discontinuities in 2D isothermal flows has rarely been discussed, since most Riemann solvers are derived for 1D Euler equations. We present a new contact resolving approximate Riemann solver for the isothermal Euler equations and show its performance for several one- and two-dimensional test problems. The new solver extends the well-known HLL solver, while retaining its computational simplicity. The significant gain in resolution of vortices is displayed by a simulation of the K\'arm\'an vortex street. We discuss the loss of Galilean invariance and its implications for the resolution of contact discontinuities, which is experienced by all modern numerical schemes for hydrodynamics in non-moving grids.

A new one-dimensional, dynamical model is proposed for geometrically thin, self-gravitating viscous accretion discs. The vertically integrated equations are simplified using the slow accretion limit and the monopole approximation with a time-dependent central point mass to account for self-gravity and accretion. It is shown that the system of partial differential equations can be reduced to a single non-linear advection diffusion equation which describes the time evolution of angular velocity. In order to solve the equation three different turbulent viscosity prescriptions are considered. It is shown that for these parametrizations the differential equation allows for similarity transformations depending only on a single non-dimensional parameter. A detailed analysis of the similarity solutions reveals that this parameter is the initial power law exponent of the angular velocity distribution at large radii. The radial dependence of the self-similar solutions is in most cases given by broken power laws. At small radii the rotation law always becomes Keplerian with respect to the current central point mass. In the outer regions the power law exponent of the rotation law deviates from the Keplerian value and approaches asymptotically the value determined by the initial condition. It is shown that accretion discs with flatter rotation laws at large radii yield higher accretion rates. The methods are applied to self-gravitating accretion discs in active galactic nuclei. Fully self-gravitating discs are found to evolve faster than nearly Keplerian discs. The implications on supermassive black hole formation and Quasar evolution are discussed.