Shock-capturing particle hydrodynamics with reproducing kernels

Abstract

We present and explore a new shock-capturing particle hydrodynamics approach. Our starting point is a commonly used discretization of smoothed particle hydrodynamics. We enhance this discretization with Roe's approximate Riemann solver, we identify its dissipative terms, and in these terms, we use slope-limited linear reconstruction. All gradients needed for our method are calculated with linearly reproducing kernels that are constructed to enforce the two lowest-order consistency relations. We scrutinize our reproducing kernel implementation carefully on a "glass-like" particle distribution, and we find that constant and linear functions are recovered to machine precision. We probe our method in a series of challenging 3D benchmark problems ranging from shocks over instabilities to Schulz-Rinne-type vorticity-creating shocks. All of our simulations show excellent agreement with analytic/reference solutions.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…