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.

In the framework of the coagulation scenario, kilometre-sized planetesimals form by subsequent collisions of pre-planetesimals of sizes from centimetre to hundreds of metres. Pre-planetesimals are fluffy, porous dust aggregates, which are inhomogeneous owing to their collisional history. Planetesimal growth can be prevented by catastrophic disruption in pre-planetesimal collisions above the destruction velocity threshold. We develop an inhomogeneity model based on the density distribution of dust aggregates, which is assumed to be a Gaussian distribution with a well-defined standard deviation. As a second input parameter, we consider the typical size of an inhomogeneous clump. These input parameters are easily accessible by laboratory experiments. For the simulation of the dust aggregates, we utilise a smoothed particle hydrodynamics (SPH) code with extensions for modelling porous solid bodies. The porosity model was previously calibrated for the simulation of silica dust, which commonly serves as an analogue for pre-planetesimal material. The inhomogeneity is imposed as an initial condition on the SPH particle distribution. We carry out collisions of centimetre-sized dust aggregates of intermediate porosity. We vary the standard deviation of the inhomogeneous distribution at fixed typical clump size. The collision outcome is categorised according to the four-population model. We show that inhomogeneous pre-planetesimals are more prone to destruction than homogeneous aggregates. Even slight inhomogeneities can lower the threshold for catastrophic disruption. For a fixed collision velocity, the sizes of the fragments decrease with increasing inhomogeneity. Pre-planetesimals with an active collisional history tend to be weaker. This is a possible obstacle to collisional growth and needs to be taken into account in future studies of the coagulation scenario.

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.

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.

Recent years have shown many advances in our knowledge of the collisional evolution of protoplanetary dust. Based on a variety of dust-collision experiments in the laboratory, our view of the growth of dust aggregates in protoplanetary disks is now supported by a deeper understanding of the physics involved in the interaction between dust agglomerates. However, the parameter space, which determines the collisional outcome, is huge and sometimes inaccessible to laboratory experiments. Very large or fluffy dust aggregates and extremely low collision velocities are beyond the boundary of today's laboratories. It is therefore desirable to augment our empirical knowledge of dust-collision physics with a numerical method to treat arbitrary aggregate sizes, porosities and collision velocities. In this article, we implement experimentally-determined material parameters of highly porous dust aggregates into a Smooth Particle Hydrodynamics (SPH) code, in particular an omnidirectional compressive-strength and a tensile-strength relation. We also give a prescription of calibrating the SPH code with compression and low-velocity impact experiments. In the process of calibration, we developed a dynamic compressive-strength relation and estimated a relation for the shear strength. Finally, we defined and performed a series of benchmark tests and found the agreement between experimental results and numerical simulations to be very satisfactory. SPH codes have been used in the past to study collisions at rather high velocities. At the end of this work, we show examples of future applications in the low-velocity regime of collisional evolution.

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.