Jonathan’s Rust Adventures and the Quest for the Numerically Stable Soft-Sphere Integrator

In this thesis I will implement a numerically stable soft-sphere collision model using Hooke’s Law as the restoring force. This model allows choosing a desired coefficient of restitution and maximum penetration depth in order to generate appropriate values for the spring constant and damping. Then I will explore the applicability of various languages to the field of kD-tree based N-body simulations, concluding that Rust is competitive in both the time and memory usage to C/C++, which are the languages traditionally used for high-performance applications. Additionally, I found that the usage of higher-level languages like Java come at significant time or memory usage costs, making them of little value in the field of large astrophysics simulations. Lastly, as soft-sphere simulations require much smaller time steps for integration compared to hard-spheres, I implemented and benchmarked a priority-queue based adaptive time step system to only take the small time steps when needed. The benchmarks showed a significant speedup compared soft-sphere simulations without this adaptive time step method at roughly 8x faster. In direct particle-particle testing, the adaptive time step showed that with as little as 15 steps during a collision, which the adaptive time step method allows me to set, stable collision handling can be achieved across a spectrum of impact speeds and particle size ratios

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.Comment: 19 pages, 11 figures, accepted for publication in Advances in Space Research, Special Issue: Asteroids & Space Debri

The Development and Implementation of Numerical Tools for Investigation into the Granular Dynamics of Solid Solar System Bodies

The work advanced in this thesis joins together the disciplines of planetary science and granular physics. Grain dynamics have played a prominent role in the evolution of our Solar System from planetesimal formation billions of years ago to the surface processes that take place today on terrestrial planets, moons, and small bodies. Recent spacecraft images of small Solar System bodies provide strong evidence that the majority of these bodies are covered in regolith. This regolith ranges in size from the fine powder found on the Moon to large rocks and boulders, like the 27 m Yoshinodai boulder on the small asteroid, Itokawa. Accordingly, the processes that take place on the solid bodies of the Solar System vary widely based upon the material properties of the regolith and the gravitational environments on their surfaces. An understanding of granular dynamics is also critical for the design and operations of landers, sampling devices and rovers to be included in space missions. Part of my research is concerned with the development of numerical tools that have the ability to provide explanations for the types of processes that our spacecraft have observed. Granular processes on Earth are incredibly complex and varied, and constitute an enormous field of study on their own, with input taken from across the broad disciplines of engineering and the physical sciences. In micro-gravity, additional forces, which on Earth are relevant only to micron-size particles or smaller, are expected to become important for material up to the size of large rocks, adding further complexity. The numerical tools developed in this work allow for the simulation of grains using an adaptation of the Soft-Sphere Discrete Element Method (SSDEM) along with implementations of cohesive forces between particles into an existing parallel gravity tree code

Modelling of Bulk Material Flow Properties

Flow property testing is important in the design of handling equipment for bulk solids and the Jenike shear tester is a common and reliable method for ensuring flow from hoppers under the force of gravity alone. The Jenike shear test procedure has various stages and is known to have issues with operator dependency, questions have also been raised regarding the stress state within the shear cell. Currently the data for high pressure flow functions, which are applicable to large capacity storage facilities, are extrapolated from low pressure test data using a 3-parameter equation. Very little literature is provided to support the use of the 3-parameter equation and most shear testing devices are limited to major consolidation stresses of 100 kPa

Modelling the behaviour of granular material on the surface of asteroids and under different gravity conditions (e.g., Mars, the Moon)

This thesis, at the interface between the scientific disciplines of planetary science and granular physics, has two key components, both of which intend to increase our understanding of granular dynamics in varying gravitational conditions. The dynamics of granular materials are involved in the evolution of solid planets and small bodies in our Solar System, whose surfaces are generally covered with regolith. Understanding granular dynamics is also critical for the design and/or operations of landers, sampling devices and rovers to be included in space missions. The first component of this thesis is the validation of the hard-sphere discrete element method implementation in the N-body code pkdgrav to model the dynamics of granular material. By direct comparison with results from laboratory experiments, it is demonstrated that the hard-sphere discrete element method implementation in pkdgrav is valid for modelling granular material in dilute regimes and is capable of reproducing the complex dynamical behaviour of a specific dense system as well. The second component is focussed on the AstEx parabolic flight experiment. This experiment, with the aim of characterising the response of granular material to rotational shear forces in a microgravity environment, was designed, constructed, flown and the data were analysed as part of this thesis. It was found that the effect of constant shearing on a granular material in a direction perpendicular to the gravity field is not strongly influenced by gravity. The AstEx experiment has demonstrated, for the first time, that the efficiency of granular convection may decrease in the presence of a weak gravitational field, similar to that on the surface of small bodies. The first measurements of transient weakening of granular material after shear reversal in microgravity are also reported. Results suggest that the force contact network may be weaker in microgravity, although the influence of any change in the contact network is felt by the granular material over much larger distances. This may have important implications for our interpretation of asteroid surfaces. Continued advancement of our understanding of granular materials in varying gravitational conditions requires futher experiments and the development of the soft-sphere discrete element method implementation in pkdgrav in order to model the granular regimes that are inaccessbile to the hard-sphere implementation

Analysis and modeling of structure formation in granular and fluid-solid flows

Granular and multiphase flows are encountered in a number of industrial processes with particular emphasis in this manuscript given to the particular applications in cement pumping, pneumatic conveying, fluid catalytic cracking, CO2 capture, and fast pyrolysis of bio-materials. These processes are often modeled using averaged equations that may be simulated using computational fluid dynamics. Closure models are then required that describe the average forces that arise from both interparticle interactions, e.g. shear stress, and interphase interactions, such as mean drag. One of the biggest hurdles to this approach is the emergence of non-trivial spatio-temporal structures in the particulate phase, which can significantly modify the qualitative behavior of these forces and the resultant flow phenomenology. For example, the formation of large clusters in cohesive granular flows is responsible for a transition from solid-like to fluid-like rheology. Another example is found in gas-solid systems, where clustering at small scales is observed to significantly lower in the observed drag. Moreover, there remains the possibility that structure formation may occur at all scales, leading to a lack of scale separation required for traditional averaging approaches. In this context, several modeling problems are treated 1) first-principles based modeling of the rheology of cement slurries, 2) modeling the mean solid-solid drag experienced by polydisperse particles undergoing segregation, and 3) modeling clustering in homogeneous gas-solid flows. The first and third components are described in greater detail. In the study on the rheology of cements, several sub-problems are introduced, which systematically increase in the number and complexity of interparticle interactions. These interparticle interactions include inelasticity, friction, cohesion, and fluid interactions. In the first study, the interactions between cohesive inelastic particles was fully characterized for the first time. Next, kinetic theory was used to predict the cooling of a gas of such particles. DEM was then used to validate this approach. A study on the rheology of dry cohesive granules with and without friction was then carried out, where the physics of different flow phenomenology was exhaustively explored. Lastly, homogeneous cement slurry simulations were carried out, and compared with vane-rheometer experiments. Qualitative agreement between simulation and experiment were observed. Lastly, the physics of clustering in homogeneous gas-solid flows is explored in the hopes of gaining a mechanistic explanation of how particle-fluid interactions lead to clustering. Exact equations are derived, detailing the evolution of the two particle density, which may be closed using high-fidelity particle-resolved direct numerical simulation. Two canonical gas-solid flows are then addressed, the homogeneously cooling gas-solid flow (HCGSF) and sedimenting gas-solid flow (SGSF). A mechanism responsible for clustering in the HCGSF is identified. Clustering of plane-wave like structures is observed in the SGSF, and the exact terms are quantified. A method for modeling the dynamics of clustering in these systems is proposed, which may aid in the prediction of clustering and other correlation length-scales useful for less expensive computations

Publications of the Jet Propulsion Laboratory July 1965 through July 1966

Bibliography on Jet Propulsion Laboratory technical reports and memorandums, space programs summary, astronautics information, and literature searche

Marshall Space Flight Center Research and Technology Report 2019

Today, our calling to explore is greater than ever before, and here at Marshall Space Flight Centerwe make human deep space exploration possible. A key goal for Artemis is demonstrating and perfecting capabilities on the Moon for technologies needed for humans to get to Mars. This years report features 10 of the Agencys 16 Technology Areas, and I am proud of Marshalls role in creating solutions for so many of these daunting technical challenges. Many of these projects will lead to sustainable in-space architecture for human space exploration that will allow us to travel to the Moon, on to Mars, and beyond. Others are developing new scientific instruments capable of providing an unprecedented glimpse into our universe. NASA has led the charge in space exploration for more than six decades, and through the Artemis program we will help build on our work in low Earth orbit and pave the way to the Moon and Mars. At Marshall, we leverage the skills and interest of the international community to conduct scientific research, develop and demonstrate technology, and train international crews to operate further from Earth for longer periods of time than ever before first at the lunar surface, then on to our next giant leap, human exploration of Mars. While each project in this report seeks to advance new technology and challenge conventions, it is important to recognize the diversity of activities and people supporting our mission. This report not only showcases the Centers capabilities and our partnerships, it also highlights the progress our people have achieved in the past year. These scientists, researchers and innovators are why Marshall and NASA will continue to be a leader in innovation, exploration, and discovery for years to come

16th European Symposium on Comminution and Classification: book of extended abstracts

Extended abstracts from the 16th European Symposium on Comminution and Classification, ESCC 2019 held at the University of Leeds, 2-4 September 2019. Based on the abstracts received, the symposium was structured in the following themes: fundamentals of size reduction, innovations in milling and classification, nanomilling, mechano-chemistry and solid state transformations, pharmaceuticals and foods, attrition and wear, and related modelling. The notable number of abstracts received on modelling made it possible to divide them in sub-themes: mechanistic, population balance, discrete element and coupling with computational fluid dynamics

Robust event-driven particle tracking in complex geometries

The authors gratefully acknowledge the support of the German Research Foundation (DFG) through Grants PO 472/20 and SFB-814. FPS and ZISC are thanked for support. We would like to thank Prapanch Nair for helpful discussions.Peer reviewedPostprin