Source author record

Patrick Michel

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

22 published item(s)

preprint2022arXiv

Dynamical Evolution of the Didymos-Dimorphos Binary Asteroid as Rubble Piles following the DART Impact

Previous efforts have modeled the Didymos system as two irregularly shaped rigid bodies, although it is likely that one or both components are in fact rubble piles. Here, we relax the rigid-body assumption to quantify how this affects the spin and orbital dynamics of the system following the DART impact. Given known fundamental differences between our simulation codes, we find that faster rigid-body simulations produce nearly the same result as rubble-pile models in scenarios with a moderate value for the momentum enhancement factor, $β$ ($β{\sim}3$) and an ellipsoidal secondary. This indicates that the rigid-body approach is likely adequate for propagating the post-impact dynamics necessary to meet DART Mission requirements. Although, if Dimorphos has a highly-irregular shape or structure, or if $β$ is unexpectedly large, then rubble-pile effects may become important. If Dimorphos's orbit and spin state are sufficiently excited, then surface particle motion is also possible. However, these simulations are limited in their resolution and range of material parameters, so they serve as a demonstration of principle, and Future work is required to fully understand the likelihood and magnitude of surface motion.

preprint2022arXiv

Effects of impact and target parameters on the results of a kinetic impactor: predictions for the Double Asteroid Redirection Test (DART) mission

The Double Asteroid Redirection Test (DART) spacecraft will impact into the asteroid Dimorphos on September 26, 2022 as a test of the kinetic impactor technique for planetary defense. The efficiency of the deflection following a kinetic impactor can be represented using the momentum enhancement factor, Beta, which is dependent on factors such as impact geometry and the specific target material properties. Currently, very little is known about Dimorphos and its material properties that introduces uncertainty in the results of the deflection efficiency observables, including crater formation, ejecta distribution, and Beta. The DART Impact Modeling Working Group (IWG) is responsible for using impact simulations to better understand the results of the DART impact. Pre-impact simulation studies also provide considerable insight into how different properties and impact scenarios affect momentum enhancement following a kinetic impact. This insight provides a basis for predicting the effects of the DART impact and the first understanding of how to interpret results following the encounter. Following the DART impact, the knowledge gained from these studies will inform the initial simulations that will recreate the impact conditions, including providing estimates for potential material properties of Dimorphos and Beta resulting from DARTs impact. This paper summarizes, at a high level, what has been learned from the IWG simulations and experiments in preparation for the DART impact. While unknown, estimates for reasonable potential material properties of Dimorphos provide predictions for Beta of 1-5, depending on end-member cases in the strength regime.

preprint2022arXiv

Predictions for the Dynamical States of the Didymos System before and after the Planned DART Impact

NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $β$ at a realistic scale, wherein ejecta from the impact provides an additional deflection to the target. Earth-based observations, the LICIACube spacecraft (to be detached from DART prior to impact), and ESA's follow-up Hera mission to launch in 2024, will provide additional characterization of the deflection test. Together Hera and DART comprise the Asteroid Impact and Deflection Assessment (AIDA) cooperation between NASA and ESA. Here the predicted dynamical states of the binary system upon arrival and after impact are presented. The assumed dynamically relaxed state of the system will be excited by the impact, leading to an increase in eccentricity and slight tilt of the orbit together with enhanced libration of Dimorphos with amplitude dependent on the currently poorly known target shape. Free rotation around the moon's long axis may also be triggered and the orbital period will experience variations from seconds to minutes over timescales of days to months. Shape change of either body due to cratering or mass wasting triggered by crater formation and ejecta may affect $β$ but can be constrained through additional measurements. Both BYORP and gravity tides may cause measurable orbital changes on the timescale of Hera's rendezvous.

preprint2022arXiv

Rotation-induced granular motion on the secondary component of binary asteroids: Application to the DART impact on Dimorphos

NASA's Double Asteroid Redirection Test (DART) mission will kinetically impact Dimorphos, the secondary component of the Didymos binary asteroid system, which will excite Dimorphos's dynamical state and lead to significant libration about the synchronous state and possibly chaotic non-principal axis rotation. Although this particular outcome is human caused, many other secondary components of binary systems are also prone to such exotic spin states. For a satellite in an excited spin state, the time-varying tidal and rotational environment can lead to significant surface accelerations. Depending on the circumstances, this mechanism may drive granular motion on the surface of the secondary. We modeled the dynamical evolution of a Didymos-like binary asteroid system using a fully coupled, three-dimensional simulation code. Then, we computed the time-varying gravitational and rotational accelerations felt over the entire surface resulting from the secondary's perturbed dynamical state. We find that an excited spin and orbit can induce large changes in the effective surface slope, potentially triggering granular motion and surface refreshment. However, for the case of the DART impact, this effect is highly dependent on many unknowns, such as Dimorphos's detailed shape, bulk density, surface geology, and the momentum transferred. Aside from the Didymos system and the DART mission, this effect also has important implications for binary systems in general.

preprint2021arXiv

Resurfacing processes constrained by crater distribution on Ryugu

Understanding the geological modification processes on asteroids can provide information concerning their surface history. Images of small asteroids from spacecraft show a depletion in terms of smaller craters. Seismic shaking was considered to be responsible for crater erasure and the main driver modifying the geology of asteroids via regolith convection or the Brazil nut effect. However, a recent artificial impact experiment on the asteroid Ryugu by the Japanese Hayabusa2 mission revealed minimal seismic activity. To investigate whether a seismic shaking model can reproduce the observed crater record, the crater distribution on Ryugu was analyzed using crater production functions under cohesionless conditions. Crater retention ages were estimated as a function of crater diameter for Ryugu, Itokawa, Eros, and Bennu using the crater size-frequency distribution and crater production function estimated for those asteroids. We found that the power-law indices "a" are inconsistent with diffusion processes (e.g., seismic shaking, a=2). This result suggests that seismic shaking models based on diffusion equations cannot explain the crater distribution on small asteroids. Alternative processes include surface flows, possibly at the origin of geomorphological and spectral features of Ryugu. We demonstrate that the vertical mixing of material at depths shallower than 1 m occurs over 10^3-10^5 yr by cratering and obliteration. The young surface age of Ryugu is consistent with the slow space weathering that results from cratering, as suggested in previous studies. The timescale (10^4-10^6 yr) required for resurfacing at depths of 2-4 m can be compared with the cosmic-ray exposure ages of returned samples to constrain the distribution of impactors that collide with Ryugu.

preprint2020arXiv

Collisional heating and compaction of small bodies: Constraints for their origin and evolution

The current properties of small bodies provide important clues to their origin and history. However, how much small bodies were processed by past collisions and to what extent they retain a record of processes that took place during the formation and early evolution of the Solar System is still poorly understood. Here we study the degree of collisional heating and compaction by analysing the large set of previous simulations of small body break-ups by Jutzi et al. (2019), which used porous targets of 50 - 400 km in diameter and investigated a large range of impact velocities, angles as well as energies. We find that the degree of impact processing is generally larger than found in previous studies which considered smaller objects (e.g. Jutzi et al., 2017; Schwartz et al., 2018). However, there is a clear dichotomy in terms of impact processing: the escaping material always experiences stronger heating than the material bound to the largest remnant. Assuming they originate from the same parent body, some of the observed differences between the recently visited asteroids Ryugu and Bennu may be explained by a different location of the material eventually forming these asteroids in the original parent body. Our results also provide constraints on the initial size of cometary nuclei.

preprint2020arXiv

Debiased albedo distribution for Near Earth Objects

We extend the most recent orbital and absolute magnitude Near Earth Object (NEO) model (Granvik et al., 2018) to provide a statistical description of NEO geometric albedos. Our model is calibrated on NEOWISE albedo data for the NEO population and reproduces these data very well once a simple model for the NEOWISE observational biases is applied. The results are consistent with previous estimates. There are about 1,000 NEOs with diameter D>1km and the mean albedo to convert absolute magnitude into diameter is 0.147. We do not find any statistically significant evidence that the albedo distribution of NEOs depends on NEO size. Instead, we find evidence that the disruption of NEOs at small perihelion distances found in Granvik et al. (2016) occurs preferentially for dark NEOs. The interval between km-sized bodies striking the Earth should occur on average once every 750,000 years. Low and high albedo NEOs are well mixed in orbital space, but a trend remains with higher albedo objects being at smaller semimajor axes and lower albedo objects more likely found at larger semimajor axes.

preprint2020arXiv

Spin-driven evolution of asteroids' top-shapes at fast and slow spins seen from (101955) Bennu and (162173) Ryugu

Proximity observations by OSIRIS-REx and Hayabusa2 provided clues on the shape evolution processes of the target asteroids, (101955) Bennu and (162173) Ryugu. Their oblate shapes with equatorial ridges, or the so-called top shapes, may have evolved due to their rotational conditions at present and in the past. Different shape evolution scenarios were previously proposed; Bennu's top shape may have been driven by surface processing, while Ryugu's may have been developed due to large deformation. These two scenarios seem to be inconsistent. Here, we revisit the structural analyses in earlier works and fill a gap to connect these explanations. We also apply a semi-analytical technique for computing the cohesive strength distribution in a uniformly rotating triaxial ellipsoid to characterize the global failure of top-shaped bodies. Assuming that the structure is uniform, our semi-analytical approach describes the spatial variations in failed regions at different spin periods; surface regions are the most sensitive at longer spin periods, while interiors fail structurally at shorter spin periods. This finding suggests that the shape evolution of a top shape may vary due to rotation and internal structure, which can explain the different evolution scenarios of Bennu's and Ryugu's top shapes. We interpret our results as the indications of top shapes' various evolution processes.

preprint2020arXiv

Tidal distortion and disruption of rubble-pile bodies revisited -- Soft-sphere discrete element analyses

During close approaches to planets or stars, the morphological and dynamical properties of rubble-pile small bodies can be modified, and some may catastrophically break up. This phenomenon is of particular interest for understanding the evolution and population of small bodies, and for making predictions for future encounters. Previous numerical explorations typically used methods that do not adequately represent the nature of rubble piles. The encounter outcomes and influence factors are still poorly constrained. Based on recent advances in modeling rubble piles, we aim to provide a better understanding of the tidal encounter processes through soft-sphere discrete element modeling (SSDEM) and to establish a database of encounter outcomes and the dependencies on encounter conditions and rubble-pile properties. We performed thousands of numerical simulations using the SSDEM implemented in the N-body code pkdgrav to study the dynamical evolution of rubble piles during close encounters with the Earth. The effects of encounter conditions, material strength, arrangement, and resolution of constituent particles are explored. Three typical tidal encounter outcomes are classified, i.e., deformation, mass shedding, and disruption, ranging from mild modifications to severe damages of the progenitor. The encounter speed and distance required for causing disruption events are much smaller than those predicted by previous studies, indicating a smaller creation rate of tidally disrupted small body populations. Extremely elongated fragments with axis ratios ~1:6 can be formed by moderate encounters. Our analyses of the spin-shape evolution of the largest remnants reveal reshaping mechanisms of rubble piles in response to tidal forces, consistent with stable rubble-pile configurations derived by continuum theory. A case study for SL9 suggests a low bulk density (0.2-0.3 g/cc) for its progenitor.

preprint2020arXiv

Validating N-body code Chrono for granular DEM simulations in reduced-gravity environments

The Discrete Element Method (DEM) is frequently used to model complex granular systems and to augment the knowledge that we obtain through theory, experimentation, and real-world observations. Numerical simulations are a particularly powerful tool for studying the regolith-covered surfaces of asteroids, comets, and small moons, where reduced-gravity environments produce ill-defined flow behaviors. In this work, we present a method for validating soft-sphere DEM codes for both terrestrial and small-body granular environments. The open-source code Chrono is modified and evaluated first with a series of simple two-body-collision tests, and then, with a set of piling and tumbler tests. In the piling tests, we vary the coefficient of rolling friction to calibrate the simulations against experiments with 1 mm glass beads. Then, we use the friction coefficient to model the flow of 1 mm glass beads in a rotating drum, using a drum configuration from a previous experimental study. We measure the dynamic angle of repose, the flowing layer thickness, and the flowing layer velocity for tests with different particle sizes, contact force models, coefficients of rolling friction, cohesion levels, drum rotation speeds and gravity levels. The tests show that the same flow patterns can be observed at Earth and reduced-gravity levels if the drum rotation speed and the gravity-level are set according to the dimensionless parameter known as the Froude number. Chrono is successfully validated against known flow behaviors at different gravity and cohesion levels, and will be used to study small-body regolith dynamics in future works.

preprint2016arXiv

Dealing with Uncertainties in Asteroid Deflection Demonstration Missions: NEOTwIST

Deflection missions to near-Earth asteroids will encounter non-negligible uncertainties in the physical and orbital parameters of the target object. In order to reliably assess future impact threat mitigation operations such uncertainties have to be quantified and incorporated into the mission design. The implementation of deflection demonstration missions offers the great opportunity to test our current understanding of deflection relevant uncertainties and their consequences, e.g., regarding kinetic impacts on asteroid surfaces. In this contribution, we discuss the role of uncertainties in the NEOTwIST asteroid deflection demonstration concept, a low-cost kinetic impactor design elaborated in the framework of the NEOShield project. The aim of NEOTwIST is to change the spin state of a known and well characterized near-Earth object, in this case the asteroid (25143) Itokawa. Fast events such as the production of the impact crater and ejecta are studied via cube-sat chasers and a flyby vehicle. Long term changes, for instance, in the asteroid's spin and orbit, can be assessed using ground based observations. We find that such a mission can indeed provide valuable constraints on mitigation relevant parameters. Furthermore, the here proposed kinetic impact scenarios can be implemented within the next two decades without threatening Earth's safety.

preprint2016arXiv

Ejecta Cloud from a Kinetic Impact on the Secondary of a Binary Asteroid: I. Mechanical Environment and Dynamic Model

An understanding of the post-impact dynamics of ejecta clouds are crucial to the planning of a kinetic impact mission to an asteroid, and also has great implications for the history of planetary formation. The purpose of this article to track the evolution of ejecta produced by AIDA mission, which targets for kinetic impact the secondary of near-Earth binary asteroid 65803 Didymos on 2022, and to feedback essential informations to AIDA's ongoing phase-A study. We present a detailed dynamic model for the simulation of an ejecta cloud from a binary asteroid that synthesizes all relevant forces based on a previous analysis of the mechanical environment. We apply our method to gain insight into the expected response of Didymos to the AIDA impact, including the subsequent evolution of debris and dust. The crater scaling relations from laboratory experiments are employed to approximate the distributions of ejecta mass and launching speed. The size composition of fragments is modeled with a power law fitted from observations of real asteroid surface. A full-scale demonstration is simulated using parameters specified by the mission. We report the results of the simulation, which include the computed spread of the ejecta cloud and the recorded history of ejecta accretion and escape. The violent period of the ejecta evolution is found to be short, and is followed by a stage where the remaining ejecta is gradually cleared. Solar radiation pressure proves to be efficient in cleaning dust-size ejecta, and the simulation results after two weeks shows that large debris on polar orbits (perpendicular to the binary orbital plane) has a survival advantage over smaller ejecta and ejecta that keep to low latitudes.

preprint2016arXiv

Small-body deflection techniques using spacecraft: techniques in simulating the fate of ejecta

We define a set of procedures to numerically study the fate of ejecta produced by the impact of an artificial projectile with the aim of deflecting an asteroid. Here we develop a simplified, idealized model of impact conditions that can be adapted to fit the details of specific deflection-test scenarios, such as what is being proposed for the AIDA project. Ongoing studies based upon the methodology described here can be used to inform observational strategies and safety conditions for an observing spacecraft. To account for ejecta evolution, the numerical strategies we are employing are varied and include a large N-Body component, a smoothed-particle hydrodynamics (SPH) component, and an application of impactor scaling laws. Simulations that use SPH-derived initial conditions show high-speed ejecta escaping at low angles of inclination, and very slowly moving ejecta lofting off the surface at higher inclination angles, some of which re-impacts the small-body surface. We are currently investigating the realism of this and other models' behaviors. Next steps will include the addition of solar perturbations to the model and applying the protocol developed here directly to specific potential mission concepts such as the proposed AIDA scenario.

preprint2015arXiv

Collisional Formation and Modeling of Asteroid Families

In the last decade, thanks to the development of sophisticated numerical codes, major breakthroughs have been achieved in our understanding of the formation of asteroid families by catastrophic disruption of large parent bodies. In this review, we describe numerical simulations of asteroid collisions that reproduced the main properties of families, accounting for both the fragmentation of an asteroid at the time of impact and the subsequent gravitational interactions of the generated fragments. The simulations demonstrate that the catastrophic disruption of bodies larger than a few hundred meters in diameter leads to the formation of large aggregates due to gravitational reaccumulation of smaller fragments, which helps explain the presence of large members within asteroid families. Thus, for the first time, numerical simulations successfully reproduced the sizes and ejection velocities of members of representative families. Moreover, the simulations provide constraints on the family dynamical histories and on the possible internal structure of family members and their parent bodies.

preprint2014arXiv

Low-speed Impact Simulations into Regolith in Support of Asteroid Sampling Mechanism Design I.: Comparison with 1-g Experiments

This study is carried out in the framework of sample-return missions to asteroids that use a low-speed projectile as the primary component of its sampling mechanism (e.g., JAXA's Hayabusa and Hayabusa2 missions). We perform numerical simulations of such impacts into granular materials using different projectile shapes under Earth's gravity. We then compare the amounts of ejected mass obtained in our simulations against what was found in experiments that used similar setups, which allows us to validate our numerical approach. For the targets, we consider 2 different monodisperse grain-diameter sizes: 5 mm and 3 mm. The impact speed of the projectile is 11 m s$^{-1}$ directed downward, perpendicular to the surface of the targets. Using an implementation of the soft-sphere discrete element method (SSDEM) in the $N$-Body gravity tree code PKDGRAV, previously validated in the context of low-speed impacts into sintered glass bead agglomerates, we find a noticeable dependence of the amount of ejected mass on the projectile shape. As found in experiments, in the case of the larger target grain size (5 mm), a conically shaped projectile ejects a greater amount of mass than do projectiles of other shapes, including disks and spheres. We then find that numerically the results are sensitive to the normal coefficient of restitution of the grains, especially for impacts into targets comprised of smaller grains (3 mm). We also find that static friction plays a more important role for impacts into targets comprised of the larger grains. As a preliminary demonstration, one of these considered setups is simulated in a microgravity environment. As expected, a reduction in gravity increases both the amount of ejected mass and the timescale of the impact process. Our methodology is also adaptable to the conditions of sampling mechanisms included in specific mission designs.

preprint2014arXiv

Numerical predictions of surface effects during the 2029 close approach of asteroid 99942 Apophis

Asteroid (99942) Apophis' close approach in 2029 will be one of the most significant small-body encounter events in the near future and offers a good opportunity for in situ exploration to determine the asteroid's surface properties and measure any tidal effects that might alter its regolith configuration. Resurfacing mechanics has become a new focus for asteroid researchers due to its important implications for interpreting surface observations, including space weathering effects. This paper provides a prediction for the tidal effects during the 2029 encounter, with an emphasis on whether surface refreshing due to regolith movement will occur. The potential shape modification of the object due to the tidal encounter is first confirmed to be negligibly small with systematic simulations, thus only the external perturbations are taken into account for this work (despite this, seismic shaking induced by shifting blocks might still play a weak role and we will look into this mechanism in future work). A two-stage approach is developed to model the responses of asteroid surface particles (the regolith) based on the soft-sphere implementation of the parallel N-body gravity tree code pkdgrav. A full-body model of Apophis is sent past the Earth on the predicted trajectory to generate the data of all forces acting at a target point on the surface. A sandpile constructed in the local frame is then used to approximate the regolith materials; all the forces the sandpile feels during the encounter are imposed as external perturbations to mimic the regolith's behavior in the full scenario. The local mechanical environment on the asteroid surface is represented in detail, leading to an estimation of the change in global surface environment due to the encounter. Typical patterns of perturbation are presented that depend on the asteroid orientation and sense of rotation at perigee.

preprint2014arXiv

Numerical Simulations of Collisional Disruption of Rotating Gravitational Aggregates: Dependence on Material Properties

Our knowledge of the strengths of small bodies in the Solar System is limited by our poor understanding of their internal structures, and this, in turn, clouds our understanding of the formation and evolution of these bodies. Observations of the rotational states of asteroids whose diameters are larger than a few hundreds of meters have revealed that they are dominated by gravity and that most are unlikely to be monoliths; however, there is a wide range of plausible internal structures. Numerical and analytical studies of shape and spin limits of gravitational aggregates and their collisional evolution show a strong dependence on shear strength. In order to study this effect, we carry out a systematic exploration of the dependence of collision outcomes on dissipation and friction parameters of the material components making up the bodies. We simulate the catastrophic disruption (leading to the largest remnant retaining 50% of the original mass) of km-size asteroids modeled as gravitational aggregates using pkdgrav, a cosmology N-body code adapted to collisional problems and recently enhanced with a new soft-sphere collision algorithm that includes more realistic contact forces. We find that for a range of three different materials, higher friction and dissipation values increase the catastrophic disruption threshold by about half a magnitude. Furthermore, we find that pre-impact rotation systematically increases mass loss on average, regardless of the target's internal configuration. Our results have important implications for the efficiency of planet formation via planetesimal growth, and also more generally to estimate the impact energy threshold for catastrophic disruption, as this generally has only been evaluated for non-spinning bodies without detailed consideration of material properties.

preprint2014arXiv

The Brazil-nut effect and its application to asteroids

Out of the handful of asteroids that have been imaged, some have distributions of blocks that are not easily explained. In this paper, we investigate the possibility that seismic shaking leads to the size sorting of particles in asteroids. In particular, we focus on the so-called Brazil Nut Effect (BNE) that separates large particles from small ones under vibrations. We study the BNE over a wide range of parameters by using the N-body code PKDGRAV, and find that the effect is largely insensitive to the coefficients of restitution, but sensitive to friction constants and oscillation speeds. Agreeing with the previous results, we find that convection drives the BNE, where the intruder rises to the top of the particle bed. For the wide-cylinder case, we also observe a "whale" effect, where the intruder follows the convective current and does not stay at the surface. We show that the non-dimensional critical conditions for the BNE agree well with previous studies. We also show that the BNE is scalable for low-gravity environments and that the rise speed of an intruder is proportional to the square root of the gravitational acceleration. Finally, we apply the critical conditions to observed asteroids, and find that the critical oscillation speeds are comparable to the seismic oscillation speeds that are expected from non-destructive impacts.

preprint2014arXiv

The Design Reference Asteroid for the OSIRIS-REx Mission Target (101955) Bennu

The Design Reference Asteroid (DRA) is a compilation of all that is known about the OSIRIS-REx mission target, asteroid (101955) Bennu. It contains our best knowledge of the properties of Bennu based on an extensive observational campaign that began shortly after its discovery, and has been used to inform mission plan development and flight system design. The DRA will also be compared with post-encounter science results to determine the accuracy of our Earth-based characterization efforts. The extensive observations of Bennu in 1999 has made it one of the best-characterized near-Earth asteroids. Many physical parameters are well determined, and span a number of categories: Orbital, Bulk, Rotational, Radar, Photometric, Spectroscopic, Thermal, Surface Analog, and Environment Properties. Some results described in the DRA have been published in peer-reviewed journals while others have been reviewed by OSIRIS-REx Science Team members and/or external reviewers. Some data, such as Surface Analog Properties, are based on our best knowledge of asteroid surfaces, in particular those of asteroids Eros and Itokawa. This public release of the OSIRIS-REx Design Reference Asteroid is a annotated version of the internal OSIRIS-REx document OREX-DOCS-04.00-00002, Rev 9 (accepted by the OSIRIS-REx project on 2014-April-14). The supplemental data products that accompany the official OSIRIS-REx version of the DRA are not included in this release. We are making this document available as a service to future mission planners in the hope that it will inform their efforts.

preprint2013arXiv

Numerical simulations of granular dynamics II. Particle dynamics in a shaken granular material

Surfaces of planets and small bodies of our Solar System are often covered by a layer of granular material that can range from a fine regolith to a gravel-like structure of varying depths. Therefore, the dynamics of granular materials are involved in many events occurring during planetary and small-body evolution thus contributing to their geological properties. We demonstrate that the new adaptation of the parallel N-body hard-sphere code pkdgrav has the capability to model accurately the key features of the collective motion of bidisperse granular materials in a dense regime as a result of shaking. As a stringent test of the numerical code we investigate the complex collective ordering and motion of granular material by direct comparison with laboratory experiments. We demonstrate that, as experimentally observed, the scale of the collective motion increases with increasing small-particle additive concentration. We then extend our investigations to assess how self-gravity and external gravity affect collective motion. In our reduced-gravity simulations both the gravitational conditions and the frequency of the vibrations roughly match the conditions on asteroids subjected to seismic shaking, though real regolith is likely to be much more heterogeneous and less ordered than in our idealised simulations. We also show that collective motion can occur in a granular material under a wide range of inter-particle gravity conditions and in the absence of an external gravitational field. These investigations demonstrate the great interest of being able to simulate conditions that are to relevant planetary science yet unreachable by Earth-based laboratory experiments.

preprint2013arXiv

Numerical simulations of granular dynamics. I. Hard-sphere discrete element method and tests

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular "atmosphere" that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate.

preprint2009arXiv

Fragment properties at the catastrophic disruption threshold: The effect of the parent body's internal structure

Numerical simulations of asteroid break-ups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by Q*D, which results in the escape of half of the target's mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize Q*D for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of Q*D, and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of Q*D on the impact speed of the projectile. In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. We then propose some power-law relationships between Q*D and both target's size and impact speed that can be used in collisional evolution models.