Search

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., Comment: for details check https://authors.elsevier.com/tracking/article/details.do?aid=107437&jid=COMPHY&surname=Ragossnig

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., Comment: Not yet published

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.

Young terrestrial planets, when they are still embedded in a circumstellar disk, accumulate an atmosphere of nebula gas. The evolution and eventual evaporation of the protoplanetary disk affect the structure and dynamics of the planetary atmosphere. These processes, combined with other mass loss mechanisms, such as thermal escape driven by extreme ultraviolet and soft X-ray radiation (XUV) from the young host star, determine how much of the primary atmosphere, if anything at all, survives into later stages of planetary evolution. Our aim is to explore the structure and the dynamic outflow processes of nebula-accreted atmospheres in dependency on changes in the planetary environment. We integrate stationary hydrostatic models and perform time-dependent dynamical simulations to investigate the effect of a changing nebula environment on the atmospheric structure and the timescales on which the protoatmosphere reacts to these changes. We find that the behavior of the atmospheres strongly depends on the mass of the planetary core. For planets of about Mars-mass the atmospheric structure, and in particular the atmospheric mass, changes drastically and on very short timescales whereas atmospheres around higher mass planets are much more robust and inert.

Despite all advances in multi-dimensional hydrodynamics, investigations of stellar evolution and stellar pulsations still depend on one-dimensional computations. The present work devises an alternative to the mixing length theory or turbulence models usually adopted for the modelling of convective transport in such studies. Assuming that the largest convective patterns generate the majority of convective transport, the convective velocity field is described using two parallel radial columns to represent up- and downstream flows. Horizontal exchange in the form of fluid flow and radiation over their connecting interface couples the two columns and allows a simple circulating motion. The main parameters of this convective description have a straightforward geometrical meaning, namely the diameter of the columns (corresponding to the size of the convective cells) and the ratio of cross section between up- and downdrafts. For this geometrical setup, the time-dependent solution of the equations of radiation hydrodynamics is computed from an implicit scheme which has the advantage of being not affected by the Courant-Friedrichs-Lewy time step limit. In order to demonstrate the approach, results for the example of convection zones in Cepheids are presented.