Astrophysical Smooth Particle Hydrodynamics

The paper presents a detailed review of the smooth particle hydrodynamics (SPH) method with particular focus on its astrophysical applications. We start by introducing the basic ideas and concepts and thereby outline all ingredients that are necessary for a practical implementation of the method in a working SPH code. Much of SPH’s success relies on its excellent conservation properties and therefore the numerical conservation of physical invariants receives much attention throughout this review. The self-consistent derivation of the SPH equations from the Lagrangian of an ideal fluid is the common theme of the remainder of the text. We derive a modern, Newtonian SPH formulation from the Lagrangian of an ideal fluid. It accounts for changes of the local resolution lengths which result in corrective, so-called “grad-h-terms”. We extend this strategy to special relativity for which we derive the corresponding grad-h equation set. The variational approach is further applied to a general-relativistic fluid evolving in a fixed, curved background space-time. Particular care is taken to explicitely derive all relevant equations in a coherent way.

💡 Research Summary

The paper provides a comprehensive review of Smooth Particle Hydrodynamics (SPH) with a strong emphasis on its astrophysical applications. It begins by motivating the need for a mesh‑free, Lagrangian method in astrophysics, where complex, time‑varying geometries, strong shocks, large dynamic ranges, and strict conservation requirements are common. The authors then lay out the basic “vanilla” SPH formulation: they discretize the fluid Lagrangian, introduce the kernel interpolation, and derive the standard density estimate, pressure gradient, and artificial viscosity terms. Practical aspects such as kernel choice, smoothing length selection, adaptive resolution, time‑integration schemes (Leapfrog, predictor‑corrector, adaptive Runge‑Kutta), and best‑practice recommendations are discussed in detail.

The core of the review is the variational derivation of modern SPH equations. Starting from a particle Lagrangian L = Σ_i m_i (½ v_i² – u_i), the authors perform a systematic variation while allowing the smoothing length h_i to depend on particle coordinates. This yields the so‑called “grad‑h” correction terms, which naturally arise from ∂W/∂h contributions. These terms correct the pressure gradient and energy equations in regions where the resolution changes rapidly, thereby preserving energy and momentum to machine precision even in highly non‑uniform particle distributions. The variational approach also removes any ambiguity in symmetrizing the force law, guaranteeing that the resulting discretization inherits the exact conservation properties of the continuous fluid equations.

Building on this foundation, the paper extends SPH to relativistic regimes. In the special‑relativistic case, the fluid 4‑velocity and specific enthalpy are introduced, and the Lagrangian is written in terms of the stress‑energy tensor contracted with the 4‑velocity. The resulting equations contain Lorentz factors and pressure terms that are symmetrized with the same grad‑h corrections, ensuring accurate shock handling at relativistic speeds. For general relativity, the authors assume a fixed background metric g_{μν} and incorporate its dependence into the particle Lagrangian. Variation then produces additional terms involving the Christoffel symbols, which act as gravitational forces on the particles. This formulation allows SPH to model flows in strong gravitational fields (e.g., neutron star mergers, accretion onto black holes) while still conserving the relativistic energy‑momentum tensor.

Beyond the core methodology, the review surveys a wide range of physics modules that have been coupled to SPH in practice: self‑gravity (tree codes, Fast Multipole Methods, particle‑mesh hybrids), various equations of state (polytropic, tabulated nuclear EOS), solid‑mechanics extensions for planetary impacts, physical viscosity, thermal conduction, nuclear reaction networks, chemistry, radiation transport (flux‑limited diffusion, Monte‑Carlo, ray‑tracing), neutrino physics, magnetohydrodynamics (including divergence‑cleaning and vector‑potential formulations), and sub‑grid models for star formation, feedback, and cosmic rays. For each module the authors cite key references and discuss the numerical challenges involved.

The final sections summarize the strengths of SPH—its intrinsic adaptivity, exact conservation, and ease of adding new physics—while acknowledging limitations such as sensitivity to particle disorder, the need for careful tuning of artificial viscosity, and difficulties in maintaining ∇·B = 0 in MHD. The authors suggest future directions: higher‑order kernels, GPU‑accelerated gravity solvers, improved divergence‑cleaning schemes, and fully dynamical spacetime evolution coupled to SPH.

In essence, the paper serves as both a pedagogical guide and a technical reference, showing how a variationally derived SPH framework can be systematically extended from Newtonian to special‑ and general‑relativistic hydrodynamics, and how it can be enriched with the myriad physical processes required for realistic astrophysical simulations.