Numerical methods have become a powerful tool for research in astrophysics, but their utility
depends critically on the availability of suitable simulation codes. This calls for continuous …

Planets are expected to conclude their growth through a series of giant impacts: energetic,
global events that significantly alter planetary composition and evolution. Computer models …

Seismic images of Earth's interior have revealed two continent-sized anomalies with low
seismic velocities, known as the large low-velocity provinces (LLVPs), in the lowermost …

Numerical simulations have become one of the key tools used by theorists in all the fields of
astrophysics and cosmology. The development of modern tools that target the largest …

As the US Department of Energy (DOE) computing facilities began deploying petascale
systems in 2008, DOE was already setting its sights on exascale. In that year, DARPA …

We present a meshfree quadrature rule for compactly supported nonlocal integro-differential
equations (IDEs) with radial kernels. We apply this rule to develop a meshfree discretization …

This work details the coupling of a Smoothed Particle Hydrodynamics (SPH) fluid solver with
a general-purpose Differential Variational Inequality (DVI) based non-smooth multibody …

We present a new methodology for simulating self-gravitating general-relativistic fluids. In
our approach the fluid is modelled by means of Lagrangian particles in the framework of a …

We compare two state-of-the-art numerical codes to study the overall accuracy in modelling
the intergalactic medium and reproducing Lyman-α forest observables for DESI and high …

PySPH is an open-source, Python-based, framework for particle methods in general and
Smoothed Particle Hydrodynamics (SPH) in particular. PySPH allows a user to define a …