Researcher profile

Roland Speith

Roland Speith contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.EP astro-ph.SR physics.comp-ph physics.data-an physics.space-ph

Trust snapshot

Quick read

Trust 21 - Emerging

12works

0followers

5topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

preprint2016arXiv

A smooth particle hydrodynamics code to model collisions between solid, self-gravitating objects

Modern graphics processing units (GPUs) lead to a major increase in the performance of the computation of astrophysical simulations. Owing to the different nature of GPU architecture compared to traditional central processing units (CPUs) such as x86 architecture, existing numerical codes cannot be easily migrated to run on GPU. Here, we present a new implementation of the numerical method smooth particle hydrodynamics (SPH) using CUDA and the first astrophysical application of the new code: the collision between Ceres-sized objects. The new code allows for a tremendous increase in speed of astrophysical simulations with SPH and self-gravity at low costs for new hardware. We have implemented the SPH equations to model gas, liquids and elastic, and plastic solid bodies and added a fragmentation model for brittle materials. Self-gravity may be optionally included in the simulations and is treated by the use of a Barnes-Hut tree. We find an impressive performance gain using NVIDIA consumer devices compared to our existing OpenMP code. The new code is freely available to the community upon request.

preprint2016arXiv

Triggering Sublimation-Driven Activity of Main Belt Comets

It has been suggested that the comet-like activity of Main Belt Comets are due to the sublimation of sub-surface water-ice that has been exposed as a result of their surfaces being impacted by m-sized bodies. We have examined the viability of this scenario by simulating impacts between m-sized and km-sized objects using a smooth particle hydrodynamics approach. Simulations have been carried out for different values of the impact velocity and impact angle as well as different target material and water-mass fraction. Results indicate that for the range of impact velocities corresponding to those in the asteroid belt, the depth of an impact crater is slightly larger than 10 m suggesting that if the activation of MBCs is due to the sublimation of sub-surface water-ice, this ice has to exist no deeper than a few meters from the surface. Results also show that ice-exposure occurs in the bottom and on the interior surface of impact craters as well as the surface of the target where some of the ejected icy inclusions are re-accreted. While our results demonstrate that the impact scenario is indeed a viable mechanism to expose ice and trigger the activity of MBCs, they also indicate that the activity of the current MBCs is likely due to ice sublimation from multiple impact sites and/or the water contents of these objects (and other asteroids in the outer asteroid belt) is larger than the 5% that is traditionally considered in models of terrestrial planet formation providing more ice for sublimation. We present details of our simulations and discuss their results and implications.

preprint2015arXiv

Collisions of planetesimals and formation of planets

We present preliminary results of terrestrial planet formation using on the one hand classical numerical integration of hundreds of small bodies on CPUs and on the other hand -- for comparison reasons -- the results of our GPU code with thousands of small bodies which then merge to larger ones. To be able to determine the outcome of collision events we use our smooth particle hydrodynamics (SPH) code which tracks how water is lost during such events.

preprint2015arXiv

Planetary Systems and the Formation of Habitable Planets

As part of a national scientific network 'Pathways to Habitability' the formation of planets and the delivery of water onto these planets is a key question as water is essential for the development of life. In the first part of the paper we summarize the state of the art of planet formation - which is still under debate in the astronomical community - before we show our results on this topic. The outcome of our numerical simulations depends a lot on the choice of the initial distribution of planetesimals and planetary embryos after gas disappeared in the protoplanetary disk. We also take into account that some of these planetesimals of sizes in the order of the mass of the Moon already contained water; the quantity depends on the distance from the Sun - close-by bodies are dry, but starting from a distance of about 2 AU they can contain substantial amounts of water. We assume that the gas giants and terrestrial planets are already formed when we check the collisions of the small bodies containing water (in the order of a few percent) with the terrestrial planets. We thus are able to give an estimate of the respective contribution to the actual water content (of some Earth-oceans) in the mantle, in the crust and on the surface of Earth. In the second part we discuss in more detail how the formation of larger bodies after a collision may happen as the outcome depends on parameters like collision velocity, impact angle, and the materials involved. We present results obtained by SPH (Smooth Particle Hydrodynamics) simulations. We briefly describe this method and show different scenarios with respect to the formed bodies, possible fragmentation and the water content before and after the collision. In an appendix we discuss detection methods for extrasolar planets (close to 2000 such objects have been discovered so far).

preprint2014arXiv

Fragmentation of colliding planetesimals with water content

We investigate the outcome of collisions of Ceres-sized planetesimals composed of a rocky core and a shell of water ice. These collisions are not only relevant for explaining the formation of planetary embryos in early planetary systems, but also provide insight into the formation of asteroid families and possible water transport via colliding small bodies. Earlier studies show characteristic collision velocities exceeding the bodies' mutual escape velocity which - along with the distribution of the impact angles - cover the collision outcome regimes 'partial accretion', 'erosion', and 'hit-and-run' leading to different expected fragmentation scenarios. Existing collision simulations use bodies composed of strengthless material; we study the distribution of fragments and their water contents considering the full elasto-plastic continuum mechanics equations also including brittle failure and fragmentation.

preprint2014arXiv

The role of material strength in collisions -- Comparing solid body and hydrodynamic physics for simulating collisions of planetesimals with icy shells

Context. We investigate the effects of including material strength in multi-material planetesimal collisions. Aims. The differences between strengthless material models and including the full elasto-plastic model for solid bodies with brittle failure and fragmentation when treating collisions of asteroid-sized bodies as they occur frequently in early planetary systems are demonstrated. Methods. We study impacts of bodies of Ceres-mass with a solid rock impactor and a target with 30 weight-% water content as surface ice.} The initial impact velocities and impact parameters are varied between the escape velocity $v_\mathrm{esc}$ to about 6 $v_\mathrm{esc}$ and from head-on collisions to close fly-bys, respectively. We simulate the collisions using our own SPH code utilizing both strengthless material and the full elasto-plastic material model including brittle failure. Results. One of the most prominent differences is the higher degree of fragmentation and shattered debris clouds in the solid model. In most collision scenarios however, the final outcomes are very similar and differ primarily by the about one order of magnitude higher degree of fragmentation in the solid case. Also, the survivors tend to be of less mass in the solid case which also predicts a higher water loss than the strengthless hydro model. This may be an effect of the relatively low-energy impacts that cannot destroy the solid material instantly. As opposed to giant impacts we also observe an indication that some water ice gets transferred between the bodies.

preprint2013arXiv

Growth and fragmentation of centimetre-sized dust aggregates: the dependence on aggregate size and porosity

We carry out three-dimensional Smoothed Particle Hydrodynamics simulations of spherical homogeneous SiO2 dust aggregates to investigate how the mass and the porosity of the aggregates affects their ability to survive an impact at various different collision velocities (between 1 - 27.5m/s). We explore how the threshold velocities for fragmentation vary with these parameters. Crucially, we find that the porosity plays a part of utmost importance in determining the outcome of collisions. In particular, we find that aggregates with filling factors >37% are significantly weakened and that the velocity regime in which the aggregates grow is reduced or even non-existent (instead, the aggregates either rebound off each other or break apart). At filling factors less than ~37% we find that more porous objects are weaker but not as weak as highly compact objects with filling factors >37%. In addition we find that (for a given aggregate density) collisions between very different mass objects have higher threshold velocities than those between very similar mass objects. We find that fragmentation velocities are higher than the typical values of 1m/s and that growth can even occur for velocities as high as 27.5m/s. Therefore, while the growth of aggregates is more likely if collisions between different sized objects occurs or if the aggregates are porous with filling factor <37%, it may also be hindered if the aggregates become too compact.

preprint2013arXiv

SPH-based simulation of multi-material asteroid collisions

We give a brief introduction to smoothed particle hydrodynamics methods for continuum mechanics. Specifically, we present our 3D SPH code to simulate and analyze collisions of asteroids consisting of two types of material: basaltic rock and ice. We consider effects like brittle failure, fragmentation, and merging in different impact scenarios. After validating our code against previously published results we present first collision results based on measured values for the Weibull flaw distribution parameters of basalt.

preprint2012arXiv

Compression Behaviour of Porous Dust Agglomerates

The early planetesimal growth proceeds through a sequence of sticking collisions of dust agglomerates. Very uncertain is still the relative velocity regime in which growth rather than destruction can take place. The outcome of a collision depends on the bulk properties of the porous dust agglomerates. Continuum models of dust agglomerates require a set of material parameters that are often difficult to obtain from laboratory experiments. Here, we aim at determining those parameters from ab-initio molecular dynamics simulations. Our goal is to improveon the existing model that describe the interaction of individual monomers. We use a molecular dynamics approach featuring a detailed micro-physical model of the interaction of spherical grains. The model includes normal forces, rolling, twisting and sliding between the dust grains. We present a new treatment of wall-particle interaction that allows us to perform customized simulations that directly correspond to laboratory experiments. We find that the existing interaction model by Dominik & Tielens leads to a too soft compressive strength behavior for uni and omni-directional compression. Upon making the rolling and sliding coefficients stiffer we find excellent agreement in both cases. Additionally, we find that the compressive strength curve depends on the velocity with which the sample is compressed. The modified interaction strengths between two individual dust grains will lead to a different behaviour of the whole dust agglomerate. This will influences the sticking probabilities and hence the growth of planetesimals. The new parameter set might possibly lead to an enhanced sticking as more energy can be stored in the system before breakup.

preprint2011arXiv

The four-populations model: a new classification scheme for pre-planetesimal collisions

Within the collision growth scenario for planetesimal formation, the growth step from centimetre sized pre-planetesimals to kilometre sized planetesimals is still unclear. The formation of larger objects from the highly porous pre-planetesimals may be halted by a combination of fragmentation in disruptive collisions and mutual rebound with compaction. However, the right amount of fragmentation is necessary to explain the observed dust features in late T Tauri discs. Therefore, detailed data on the outcome of pre-planetesimal collisions is required and has to be presented in a suitable and precise format. We propose and apply a new classification scheme for pre-planetesimal collisions based on the quantitative aspects of four fragment populations: the largest and second largest fragment, a power-law population, and a sub-resolution population. For the simulations of pre-planetesimal collisions, we adopt the SPH numerical scheme with extensions for the simulation of porous solid bodies. By means of laboratory benchmark experiments, this model was previously calibrated and tested for the correct simulation of the compaction, bouncing, and fragmentation behaviour of macroscopic highly porous silica dust aggregates. It is shown that previous attempts to map collision data were much too oriented on qualitatively categorising into sticking, bouncing, and fragmentation events. We show that the four-populations model encompasses all previous categorisations and in addition allows for transitions. This is because it is based on quantitative characteristic attributes of each population such as the mass, kinetic energy, and filling factor. As a demonstration of the applicability and the power of the four-populations model, we utilise it to present the results of a study on the influence of collision velocity in head-on collisions of intermediate porosity aggregates.

preprint2010arXiv

Numerical Simulations of Highly Porous Dust Aggregates in the Low-Velocity Collision Regime

A highly favoured mechanism of planetesimal formation is collisional growth. Single dust grains, which follow gas flows in the protoplanetary disc, hit each other, stick due to van der Waals forces and form fluffy aggregates up to centimetre size. The mechanism of further growth is unclear since the outcome of aggregate collisions in the relevant velocity and size regime cannot be investigated in the laboratory under protoplanetary disc conditions. Realistic statistics of the result of dust aggregate collisions beyond decimetre size is missing for a deeper understanding of planetary growth. Joining experimental and numerical efforts we want to calibrate and validate a computer program that is capable of a correct simulation of the macroscopic behaviour of highly porous dust aggregates. After testing its numerical limitations thoroughly we will check the program especially for a realistic reproduction of various benchmark experiments. We adopt the smooth particle hydrodynamics (SPH) numerical scheme with extensions for the simulation of solid bodies and a modified version of the Sirono porosity model. Experimentally measured macroscopic material properties of silica dust are implemented. We calibrate and test for the compressive strength relation and the bulk modulus. SPH has already proven to be a suitable tool to simulate collisions at rather high velocities. In this work we demonstrate that its area of application can not only be extended to low-velocity experiments and collisions. It can also be used to simulate the behaviour of highly porous objects in this velocity regime to a very high accuracy.The result of the calibration process in this work is an SPH code that can be utilised to investigate the collisional outcome of porous dust in the low-velocity regime.

preprint2009arXiv

Towards a Dynamical Collision Model of Highly Porous Dust Aggregates

In the recent years we have performed various experiments on the collision dynamics of highly porous dust aggregates and although we now have a comprehensive picture of the micromechanics of those aggregates, the macroscopic understanding is still lacking. We are therefore developing a mechanical model to describe dust aggregate collisions with macroscopic parameters like tensile strength, compressive strength and shear strength. For one well defined dust sample material, the tensile and compressive strength were measured in a static experiment and implemented in a Smoothed Particle Hydrodynamics (SPH) code. A laboratory experiment was designed to compare the laboratory results with the results of the SPH simulation. In this experiment, a mm-sized glass bead is dropped into a cm-sized dust aggregate with the previously measured strength parameters. We determine the deceleration of the glass bead by high-speed imaging and the compression of the dust aggregate by x-ray micro-tomography. The measured penetration depth, stopping time and compaction under the glass bead are utilized to calibrate and test the SPH code. We find that the statically measured compressive strength curve is only applicable if we adjust it to the dynamic situation with a 'softness' parameter. After determining this parameter, the SPH code is capable of reproducing experimental results, which have not been used for the calibration before.