Numerical Convergence in Smoothed Particle Hydrodynamics

Abstract

We study the convergence properties of smoothed particle hydrodynamics (SPH) using numerical tests and simple analytic considerations. Our analysis shows that formal numerical convergence is possible in SPH only in the joint limit $N\to \infty$, $h\to 0$, and $N_{nb}\to \infty$, where $N$ is the total number of particles, $h$ is the smoothing length, and $N_{nb}$ is the number of neighbor particles within the smoothing volume used to compute smoothed estimates. Previous work has generally assumed that the conditions $N\to \infty$ and $h\to 0$ are sufficient to achieve convergence, while holding $N_{nb}$ fixed. We demonstrate that if $N_{nb}$ is held fixed as the resolution is increased, there will be a residual source of error that does not vanish as $N\to \infty$ and $h\to 0$. Formal numerical convergence in SPH is possible only if $N_{nb}$ is increased systematically as the resolution is improved. Using analytic arguments, we derive an optimal compromise scaling for $N_{nb}$ by requiring that this source of error balance that present in the smoothing procedure. For typical choices of the smoothing kernel, we find $N_{nb}\propto N^{1/2}$. This means that if SPH is to be used as a numerically convergent method, the required computational cost does not scale with particle number as $O(N)$, but rather as $O(N^{1+\delta })$, where $\delta \approx 1/2$, with a weak dependence on the form of the smoothing kernel.

arXiv commentThe revised version accepted in ApJ, with typos corrected and references added

References in corpus (30)

• The cosmological simulation code GADGET-2
• The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
• E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
• Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe
• Galaxies on FIRE (Feedback In Realistic Environments): Stellar Feedback Explains Cosmologically Inefficient Star Formation
• Properties of galaxies reproduced by a hydrodynamic simulation
• A model for cosmological simulations of galaxy formation physics
• Cosmological SPH simulations: The entropy equation
• Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
• Fundamental differences between SPH and grid methods
• The Aquila comparison Project: The Effects of Feedback and Numerical Methods on Simulations of Galaxy Formation
• A General Class of Lagrangian Smoothed Particle Hydrodynamics Methods and Implications for Fluid Mixing Problems
• Smoothed Particle Hydrodynamics in Astrophysics
• Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
• Inviscid SPH
• On Validating an Astrophysical Simulation Code
• Moving mesh cosmology: numerical techniques and global statistics
• A Necessary Condition for Individual Time-Steps in SPH Simulations
• Implementation of feedback in SPH: towards concordance of methods
• A Density Independent Formulation of Smoothed Particle Hydrodynamics
• Subsonic turbulence in smoothed particle hydrodynamics and moving-mesh simulations
• Resolving mixing in Smoothed Particle Hydrodynamics
• Following the flow: tracer particles in astrophysical fluid simulations
• SPHS: Smoothed Particle Hydrodynamics with a higher order dissipation switch
• SPHGal: Smoothed Particle Hydrodynamics with improved accuracy for Galaxy simulations
• Galaxy mergers on a moving mesh: a comparison with smoothed-particle hydrodynamics
• Calibrating an updated SPH scheme within GCD+
• Resolving high Reynolds numbers in SPH simulations of subsonic turbulence
• Direct Numerical Simulation of decaying two-dimensional turbulence in a no-slip square box using Smoothed Particle Hydrodynamics
• Equalizing resolution in smoothed-particle hydrodynamics calculations using self-adaptive sinc kernels

Cited by in corpus (37)

• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
• Simulating cosmic structure formation with the GADGET-4 code
• Improving the convergence properties of the moving-mesh code AREPO
• Accurate, Meshless Methods for Magneto-Hydrodynamics
• The impact of feedback on cosmological gas accretion
• Baryonic impact on the dark matter distribution in Milky Way-size galaxies and their satellites
• On the OVI Abundance in the Circumgalactic Medium of Low-Redshift Galaxies
• Entropically Damped Artificial Compressibility for SPH
• Smoothed particle hydrodynamics simulations of gas and dust mixtures
• SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code
• SPHYNX: an accurate density-based SPH method for astrophysical applications
• The Limited Accuracy of Linearized Gravity
• A new insight into the consistency of smoothed particle hydrodynamics
• Discriminating Topology in Galaxy Distributions using Network Analysis
• Primordial Earth mantle heterogeneity caused by the Moon-forming giant impact
• Simulations of core formation for frequent dark matter self-interactions
• The Effect of Protoplanetary Disk Cooling Times on the Formation of Gas Giant Planets by Gravitational Instability
• Observational signatures of outbursting protostars - I: From hydrodynamic simulations to observations
• The formation and evolution of star clusters in interacting galaxies
• Improved Performances in Subsonic Flows of an SPH Scheme with Gradients Estimated using an Integral Approach
• Techniques for second order convergent weakly-compressible smoothed particle hydrodynamics schemes without boundaries
• From large-scale to protostellar disk fragmentation into close binary stars
• Consistent SPH Simulations of Protostellar Collapse and Fragmentation
• Fast and accurate Voronoi density gridding from Lagrangian hydrodynamics data
• Conservative, density-based smoothed particle hydrodynamics with improved partition of the unity and better estimation of gradients
• Impetus: Consistent SPH calculations of 3D spherical Bondi accretion onto a black hole
• Frequency of the dark matter subhalo collisions and bifurcation sequence arising formation of dwarf galaxies
• Gas Cooling in Hydrodynamic Simulations with An Exact Time Integration Scheme
• The Formation of a Milky Way-sized Disk Galaxy 1. A Comparison of Numerical Methods
• Smoothed Particle Hydrodynamics Physically Reconsidered -- The Relation to Explicit Large Eddy Simulation and the Issue of Particle Duality
• Implicit Molecular Stresses in Weakly-Compressible Particle-Based Discretization Methods for Fluid Flow
• Novel Hydrodynamic Schemes Capturing Shocks and Contact Discontinuities and Comparison Study with Existing Methods
• Mixing Sinc kernels to improve interpolations in smoothed particle hydrodynamics without pairing instability
• Evidence on the incompatibility of smoothed particle hydrodynamics and eddy viscosity models for large eddy simulations
• Predictive Capability Maturity Quantification using Bayesian Network
• SPH modeling of natural convection around a heated horizontal cylinder: A comparison with experiments