Related papers: The Material Point Method (MPM) for simulating hyp…
Planet formation simulations are capable of directly integrating the evolution of hundreds to thousands of planetary embryos and planetesimals, as they accrete pairwise to become planets. In principle such investigations allow us to better…
Long-term integrations of asteroid orbits with high-accuracy numerical integrators are essential for understanding dynamical evolution and ejection from the Solar System, but are computationally expensive. Here, we investigate the dynamical…
The constant potential molecular dynamics simulation method proposed by Siepmann and Sprik and reformulated later by Reed (SR-CPM) has been widely employed to investigate the metallic electrolyte/electrode interfaces, especially for…
Quantum computing has shown great potential in various quantum chemical applications such as drug discovery, material design, and catalyst optimization. Although significant progress has been made in quantum simulation of simple molecules,…
In this paper we present solutions to three short comings of Smoothed Particles Hydrodynamics (SPH) encountered in previous work when applying it to Giant Impacts. First we introduce a novel method to obtain accurate SPH representations of…
We use a smooth particle hydrodynamics method (SPH) to simulate colliding rocky and icy bodies from cm-scale to hundreds of km in diameter, in an effort to define self-consistently the threshold for catastrophic disruption. Unlike previous…
Low-thrust trajectories play a crucial role in optimizing scientific output and cost efficiency in asteroid belt missions. Unlike high-thrust transfers, low-thrust trajectories require solving complex optimal control problems. This…
Forming planetary systems are populated by large numbers of gravitationally interacting planetary bodies, spanning from massive giant planets to small planetesimals akin to present-day asteroids and comets. All these planetary bodies are…
We present a boundary-spheropolygon element method (BSEM), that combines the boundary integral method (BIM) and the spheropolygon-based discrete element method (SEM). The interaction between particles is simulated via the SEM, and the…
We develop a novel constitutive modeling approach for the analysis of fracture propagation in quasi-brittle materials using the Material Point Method. The kinematics of constitutive models is enriched with an additional mode of localized…
Characterizing microstructure-material-property relations calls for software tools which extract point-cloud- and continuum-scale-based representations of microstructural objects. Application examples include atom probe, electron, and…
Predicting the extreme hydrodynamic response of porous and architected lattice materials is a fundamental challenge in high energy density physics, where shock-induced pore collapse, baroclinic vorticity, and anomalous kinetic and…
Mechanical metamaterials utilize geometry to achieve exceptional mechanical properties, including those not typically possible for traditional materials. To achieve these properties, it is necessary to identify the proper structures and…
Path-integral molecular dynamics (PIMD) simulations are crucial for accurately capturing nuclear quantum effects in materials. However, their computational intensity and reliance on multiple software packages often limit their applicability…
Proper elements are quasi-integrals of motion, meaning that they can be considered constant over a certain timespan, and they permit to describe the long-term evolution with a few parameters. Near-Earth objects (NEOs) generally have a large…
We describe a newly developed hydrodynamic code for studying accretion disk processes. The numerical method uses a finite volume, nonlinear, Total Variation Diminishing (TVD) scheme to capture shocks and control spurious oscillations. It is…
We propose a novel projection method to treat near-incompressibility and volumetric locking in small- and large-deformation elasticity and plasticity within the context of higher order material point methods. The material point method is…
Estimates for asteroid masses are based on their gravitational perturbations on the orbits of other objects such as Mars, spacecraft, or other asteroids and/or their satellites. In the case of asteroid-asteroid perturbations, this leads to…
Simulation of fracturing processes in porous rocks can be divided into two main branches: (i) modeling the rock as a continuum which is enhanced with special features to account for fractures, or (ii) modeling the rock by a discrete (or…
Contact algorithm between different bodies plays an important role in solving collision problems. Usually it is not easy to be treated very well. Several ones for material point method were proposed by Bardenhangen, Brackbill, and…