VADER: A Flexible, Robust, Open-Source Code for Simulating Viscous Thin Accretion Disks

The evolution of thin axisymmetric viscous accretion disks is a classic problem in astrophysics. While models based on this simplified geometry provide only approximations to the true processes of instability-driven mass and angular momentum transport, their simplicity makes them invaluable tools for both semi-analytic modeling and simulations of long-term evolution where two- or three-dimensional calculations are too computationally costly. Despite the utility of these models, the only publicly-available frameworks for simulating them are rather specialized and non-general. Here we describe a highly flexible, general numerical method for simulating viscous thin disks with arbitrary rotation curves, viscosities, boundary conditions, grid spacings, equations of state, and rates of gain or loss of mass (e.g., through winds) and energy (e.g., through radiation). Our method is based on a conservative, finite-volume, second-order accurate discretization of the equations, which we solve using an unconditionally-stable implicit scheme. We implement Anderson acceleration to speed convergence of the scheme, and show that this leads to factor of $\sim 5$ speed gains over non-accelerated methods in realistic problems, though the amount of speedup is highly problem-dependent. We have implemented our method in the new code Viscous Accretion Disk Evolution Resource (VADER), which is freely available for download from https://bitbucket.org/krumholz/vader/ under the terms of the GNU General Public License.Comment: 58 pages, 13 figures, accepted to Astronomy & Computing; this version includes more discussion, but no other changes; code is available for download from https://bitbucket.org/krumholz/vader

Bulge Globular Clusters in Spiral Galaxies

There is now strong evidence that the metal-rich globular clusters (GC) near the center of our Galaxy are associated with the Galactic bulge rather than the disk as previously thought. Here we extend the concept of bulge GCs to the GC systems of nearby spiral galaxies. In particular, the kinematic and metallicity properties of the GC systems favor a bulge rather than a disk origin. The number of metal-rich GCs normalized by the bulge luminosity is roughly constant (i.e. bulge S_N ~ 1) in nearby spirals, and this value is similar to that for field ellipticals when only the red (metal--rich) GCs are considered. We argue that the metallicity distributions of GCs in spiral and elliptical galaxies are remarkably similar, and that they obey the same correlation of mean GC metallicity with host galaxy mass. We further suggest that the metal-rich GCs in spirals are the direct analogs of the red GCs seen in ellipticals. The formation of a bulge/spheroidal stellar system is accompanied by the formation of metal-rich GCs. The similarities between GC systems in spiral and elliptical galaxies appear to be greater than the differences.Comment: 5 pages, Latex, 2 figures, 1 table, Accepted for publication in ApJ Letter

Modelling end-pumped solid state lasers

The operation dynamics of end-pumped solid-state lasers are investigated by means of a spatially resolved numerical rate-equation model and a time-dependent analytical thermal model. The rate-equation model allows the optimization of parameters such as the output coupler transmission and gain medium length, with the aim of improving the laser output performance. The time-dependent analytical thermal model is able to predict the temperature and the corresponding induced thermal stresses on the pump face of quasi-continuous wave (qcw) end-pumped laser rods. Both models are found to be in very good agreement with experimental results