Source author record

Graeme R. Sarson

Graeme R. Sarson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA Applications Biological Physics

Catalog footprint

What is connected

4works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2019arXiv

Modelling supernova driven turbulence

High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its high-order advection and time evolution algorithms. However, the high-order algorithms and lack of Riemann solver to follow shocks make it less well suited to handling high Mach number shocks, such as those produced by supernovae (SNe). Here, we outline methods required to enable the code to efficiently and accurately model SNe, using parameters that allow stable simulation of SN-driven turbulence, in order to construct a physically realistic galactic dynamo model. These include the resolution of shocks with artificial viscosity, thermal conductivity, and mass diffusion; the correction of the mass diffusion terms; and a novel generalization of the Courant condition to include all source terms in the momentum and energy equations. We test our methods with the numerical solution of the one-dimensional (1D) Riemann shock tube (Sod, J. Comput. Phys. 1978, 27), also extended to a 1D adiabatic shock with parameters and Mach number relevant to SN shock evolution, including shocks with radiative losses. We extend our test with the three-dimensional (3D) numerical simulation of individual SN remnant evolution for a range of ambient gas densities typical of the interstellar medium and compare these to the analytical solutions of Sedov-Taylor (adiabatic) and the snowplough and Cioffi, McKee and Bertschinger (Astrophys. J. 1988, 334) results incorporating cooling and heating processes. We show that our new timestep algorithm leads to linear rather than quadratic resolution dependence as the strength of the artificial viscosity varies, because of the corresponding change in the strength of interzone gradients.

preprint2013arXiv

Inference for population dynamics in the Neolithic period

We consider parameter estimation for the spread of the Neolithic incipient farming across Europe using radiocarbon dates. We model the arrival time of farming at radiocarbon-dated, early Neolithic sites by a numerical solution to an advancing wavefront. We allow for (technical) uncertainty in the radiocarbon data, lack-of-fit of the deterministic model and use a Gaussian process to smooth spatial deviations from the model. Inference for the parameters in the wavefront model is complicated by the computational cost required to produce a single numerical solution. We therefore employ Gaussian process emulators for the arrival time of the advancing wavefront at each radiocarbon-dated site. We validate our model using predictive simulations.

preprint2013arXiv

Supernova-regulated ISM: the effects of radiative cooling and thermal conductivity on the multi-phase structure

The hydrodynamic state of the interstellar medium (ISM) heated and randomly stirred by supernovae (SNe) is investigated. We use a three-dimensional non-ideal hydrodynamic ISM model in a domain extending 0.5 x 0.5 kpc horizontally and 2 kpc vertically to explore the relative importance of various physical and numerical effects on the multi-phase, turbulent ISM. We include both Type I and II SNe, the latter occurring only in dense regions. First we investigate the role of the thermal instability in the temperature range 300-6100 K, comparing results obtained for two different cooling functions, one susceptible to the instability, the other stable. The presence of thermal instability in the system is mainly visible as the tendency of the gas to avoid the relevant temperature range, as it quickly evolves towards either colder or warmer phases. Nevertheless, the formation of dense structures for both cooling functions appears to be dominated by expanding and colliding supernova remnants, rather than by the thermal instability. As we need to include a finite thermal conduction coefficient to resolve the thermal instability with our uniform grid, we also explore the effects of changing Prandtl number on the system. The purely divergent SN forcing is found to produce significant amounts of vorticity. The relative vorticity is around 0.6-0.7 for the highest Prandtl numbers explored, Pr=40, and is observed to diminish almost by a factor of two for the lowest Prandtl numbers studied, Pr=1. Rotation laws with angular velocity decreasing or increasing outwards are investigated, enabling us to separate the contributions to kinetic helicity due to rotation and shear. When angular velocity decreases outwards, these two contributions partly cancel each other, resulting in a smaller net helicity.

preprint2012arXiv

Bayesian inference for a wavefront model of the Neolithisation of Europe

We consider a wavefront model for the spread of Neolithic culture across Europe, and use Bayesian inference techniques to provide estimates for the parameters within this model, as constrained by radiocarbon data from Southern and Western Europe. Our wavefront model allows for both an isotropic background spread (incorporating the effects of local geography), and a localized anisotropic spread associated with major waterways. We introduce an innovative numerical scheme to track the wavefront, and use Gaussian process emulators to further increase the efficiency of our model, thereby making Markov chain Monte Carlo methods practical. We allow for uncertainty in the fit of our model, and discuss the inferred distribution of the parameter specifying this uncertainty, along with the distributions of the parameters of our wavefront model. We subsequently use predictive distributions, taking account of parameter uncertainty, to identify radiocarbon sites which do not agree well with our model. These sites may warrant further archaeological study, or motivate refinements to the model.