We present an implicit numerical method to solve the time-dependent equations of radiation hydrodynamics (RHD) in axial symmetry assuming hydrostatic equilibrium perpendicular to the equatorial plane (1+1D) of a gaseous disk. The equations are formulated in conservative form on an adaptive grid and the corresponding fluxes are calculated by a spacial second order advection scheme. Self-gravity of the disk is included by solving the Possion equation. We test the resulting numerical method through comparison with a simplified analytical solution as well as through the long term viscous evolution of protoplanetary disk when due to viscosity matter is transported towards the central host star and the disk depletes. The importance of the inner boundary conditions on the structural behaviour of disks is demonstrated with several examples.

The collapse of interstellar gaseous clouds towards a protostar leads to the formation of accretion disks around the central star. Such disks can be dynamically stable if they settle in an axisymmetric state. In this letter, we investigate the long-term stability of astrophysical viscous disks around various protostars. We apply an implicit numerical code which solves the equations of radiation hydrodynamics and treats turbulence-induced viscosity according to the $\alpha$-viscosity model. We show how the viscosity is related to the disk mass. A stability criterion to determine the maximum disk mass can be formulated. We analyse such instabilities for a variety of radial points with different orbital distances from the host star and discuss the feedback on the disk in the event of an unstable protoplanetary disk. Additionally, we examine the critical disk-mass for disks with variable outer boundaries and compare them to observations of protostellar disks in the Upper Scorpius OB Association and near the Lupus complex. We derive an easily applicable method to obtain an estimate for maximum disk masses when the outer disk radius ins known.

Cosmic rays (CRs) are transported out of the galaxy by diffusion and advection due to streaming along magnetic field lines and resonant scattering off self-excited Magneto-Hydro-Dynamic (MHD) waves. Thus momentum is transferred to the plasma via the frozen-in waves as a mediator assisting the thermal pressure in driving a galactic wind. Galactic CRs (GCRs) are accelerated by shock waves generated in supernova remnants (SNRs), and they propagate from the disc into the halo. Therefore CR acceleration in the halo strongly depends on the inner disc boundary conditions. We performed hydrodynamical simulations of galactic winds in flux tube geometry appropriate for disc galaxies, describing the CR diffusive-advective transport in a hydrodynamical fashion (by taking appropriate moments of the Fokker-Planck equation) along with the energy exchange with self-generated MHD waves. Our time-dependent CR hydrodynamic simulations confirm that the evolution of galactic winds with feedback depends on the structure of the galactic halo. In case of a wind-structured halo, the wind breaks down after the last super nova (SN) has exploded. The mechanism described here offers a natural and elegant solution to explain the power-law distribution of CRs between the `knee' and the `ankle'. The transition will be naturally smooth, because the Galactic CRs accelerated at SN shocks will be `post-accelerated' by shocks generated at the inner boundary and travelling through the halo.

Planets that form early enough to be embedded in the circumstellar gas disk accumulate thick atmospheres of nebular gas. Models of these atmospheres need to specify the surface luminosity (i.e. energy loss rate) of the planet. This luminosity is usually associated with a continuous inflow of solid bodies, where the gravitational energy released from these bodies is the source of energy. However, if these bodies release energy in the atmosphere instead of at the surface, this assumption might not be justified. Our aim is to explore the interactions of infalling planetesimals with primordial atmospheres at an embedded phase of evolution. We investigate effects of atmospheric interaction on the planetesimals (mass loss) and the atmosphere (heating/cooling). We used atmospheric parameters from a snapshot of time-dependent evolution simulations for embedded atmospheres and simulated purely radial, infall events of siliceous planetesimals in a 1D, explicit code. We implemented energy transfer between friction, radiation transfer by the atmosphere and the body and thermal ablation; this gives us the possibility to examine the effects on the planetesimals and the atmosphere. We find that a significant amount of gravitational energy is indeed dissipated into the atmosphere, especially for larger planetary cores, which consequently cannot contribute to the atmospheric planetary luminosity. Furthermore, we examine that planetesimal infall events for cores, $M_\mathrm{C} > 2$M$_{\oplus}$, which actually result in a local cooling of the atmosphere; this is totally in contradiction with the classical model.