An improved SPH scheme for cosmological simulations
Authors:
A. M. Beck,
G. Murante,
A. Arth,
R. -S. Remus,
A. F. Teklu,
J. M. F. Donnert,
S. Planelles,
M. C. Beck,
P. Foerster,
M. Imgrund,
K. Dolag,
S. Borgani
Abstract:
We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-depend…
▽ More
We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.
△ Less
Submitted 20 October, 2015; v1 submitted 25 February, 2015;
originally announced February 2015.
A Bayesian method for pulsar template generation
Authors:
M. Imgrund,
D. J. Champion,
M. Kramer,
H. Lesch
Abstract:
Extracting Times of Arrival from pulsar radio signals depends on the knowledge of the pulsars pulse profile and how this template is generated. We examine pulsar template generation with Bayesian methods. We will contrast the classical generation mechanism of averaging intensity profiles with a new approach based on Bayesian inference. We introduce the Bayesian measurement model imposed and derive…
▽ More
Extracting Times of Arrival from pulsar radio signals depends on the knowledge of the pulsars pulse profile and how this template is generated. We examine pulsar template generation with Bayesian methods. We will contrast the classical generation mechanism of averaging intensity profiles with a new approach based on Bayesian inference. We introduce the Bayesian measurement model imposed and derive the algorithm to reconstruct a "statistical template" out of noisy data. The properties of these "statistical templates" are analysed with simulated and real measurement data from PSR B1133+16. We explain how to put this new form of template to use in analysing secondary parameters of interest and give various examples: We implement a nonlinear filter for determining ToAs of pulsars. Applying this method to data from PSR J1713+0747 we derive ToAs self consistently, meaning all epochs were timed and we used the same epochs for template generation. While the average template contains fluctuations and noise as unavoidable artifacts, we find that the "statistical template" derived by Bayesian inference quantifies fluctuations and remaining uncertainty. This is why the algorithm suggested turns out to reconstruct templates of statistical significance from ten to fifty single pulses. A moving data window of fifty pulses, taking out one single pulse at the beginning and adding one at the end of the window unravels the characteristics of the methods to be compared. It shows that the change induced in the classical reconstruction is dominated by random fluctuations for the average template, while statistically significant changes drive the dynamics of the proposed method's reconstruction. The analysis of phase shifts with simulated data reveals that the proposed nonlinear algorithm is able to reconstruct correct phase information along with an acceptable estimation of the remaining uncertainty.
△ Less
Submitted 14 January, 2015;
originally announced January 2015.