Resolution criteria to avoid artificial clumping in Lagrangian hydrodynamic simulations with a multi-phase interstellar medium
arXiv:2310.10721 · doi:10.1093/mnras/stad3935
Abstract
Large-scale cosmological galaxy formation simulations typically prevent gas in the interstellar medium (ISM) from cooling below K. This has been motivated by the inability to resolve the Jeans mass in molecular gas (>>) which would result in undesired artificial clumping. We show that the classical Jeans criteria derived for Newtonian gravity are not applicable in the simulated ISM if the spacing of resolution elements representing the dense ISM is below the gravitational force softening length and gravity is therefore softened and not Newtonian. We re-derive the Jeans criteria for softened gravity in Lagrangian codes and use them to analyse gravitational instabilities at and below the hydrodynamical resolution limit for simulations with adaptive and constant gravitational softening lengths. In addition, we define criteria for which a numerical runaway collapse of dense gas clumps can occur caused by over-smoothing of the hydrodynamical properties relative to the gravitational force resolution. This effect is illustrated using simulations of isolated disk galaxies with the smoothed particle hydrodynamics code Swift. We also demonstrate how to avoid the formation of artificial clumps in gas and stars by adjusting the gravitational and hydrodynamical force resolutions.
arXiv comment24 pages, 15 figures, accepted for publication in MNRAS, smaller updates to match published version
References in corpus (55)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Array Programming with NumPy
- A Universal Density Profile from Hierarchical Clustering
- The chemical composition of the Sun
- The cosmological simulation code GADGET-2
- Simulating the joint evolution of quasars, galaxies and their large-scale distribution
- The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
- The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
- Cosmological SPH simulations: A hybrid multi-phase model for star formation
- Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe
- Modeling feedback from stars and black holes in galaxy mergers
- Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
- Simulating Galaxy Formation with the IllustrisTNG Model
- First results from the IllustrisTNG simulations: the stellar mass content of groups and clusters of galaxies
- First results from the IllustrisTNG simulations: the galaxy color bimodality
- Simulating galaxy formation with black hole driven thermal and kinetic feedback
- GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
- Simba: Cosmological Simulations with Black Hole Growth and Feedback
- The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study
- First Results from the TNG50 Simulation: The evolution of stellar and gaseous disks across cosmic time
- Dancing in the dark: galactic properties trace spin swings along the cosmic web
- The Formation of a Star Cluster: Predicting the Properties of Stars and Brown Dwarfs
- Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
- Smoothed Particle Hydrodynamics and Magnetohydrodynamics
- On the relation between the Schmidt and Kennicutt-Schmidt star formation laws and its implications for numerical simulations
- The APOSTLE simulations: solutions to the Local Group's cosmic puzzles
- Simulating galactic outflows with thermal supernova feedback
- The Arepo public code release
- Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
- A General Class of Lagrangian Smoothed Particle Hydrodynamics Methods and Implications for Fluid Mixing Problems
- Simulating cosmic structure formation with the GADGET-4 code
- An energy-conserving formalism for adaptive gravitational force softening in SPH and N-body codes
- The APOSTLE project: Local Group kinematic mass constraints and simulation candidate selection
- Towards optimal softening in 3D N-body codes: I. Minimizing the force error
- Star formation and molecular hydrogen in dwarf galaxies: a non-equilibrium view
- SPHGal: Smoothed Particle Hydrodynamics with improved accuracy for Galaxy simulations
- Hydrodynamical simulations of the galaxy population: enduring successes and outstanding challenges
- Variable interstellar radiation fields in simulated dwarf galaxies: supernovae versus photoelectric heating
- FIREbox: Simulating galaxies at high dynamic range in a cosmological volume
- Radiative cooling rates, ion fractions, molecule abundances and line emissivities including self-shielding and both local and metagalactic radiation fields
- Energy equipartition between stellar and dark matter particles in cosmological simulations results in spurious growth of galaxy sizes
- Numerical convergence of simulations of galaxy formation: the abundance and internal structure of cold dark matter haloes
- Optimal softening for force calculations in collisionless N-body simulations
- Sphenix: Smoothed Particle Hydrodynamics for the next generation of galaxy formation simulations
- Optimal Smoothing for N-Body Codes
- Spurious heating of stellar motions in simulated galactic disks by dark matter halo particles
- Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
- Resolution requirements for simulating gravitational fragmentation using SPH
- Numerical convergence of hydrodynamical simulations of galaxy formation: the abundance and internal structure of galaxies and their cold dark matter haloes
- unyt: Handle, manipulate, and convert data with units in Python
- Optimal Choice of the Softening Length and Time-Step in N-body Simulations
- Chemical evolution of giant molecular clouds in simulations of galaxies
- Inconsistencies arising from the coupling of galaxy formation sub-grid models to Pressure-Smoothed Particle Hydrodynamics
- Testing the effect of resolution on gravitational fragmentation with Lagrangian hydrodynamic schemes
- Tests of subgrid models for star formation using simulations of isolated disk galaxies
Cited by in corpus (3)
- Hybrid-chimes: A model for radiative cooling and the abundances of ions and molecules in simulations of galaxy formation
- COLIBRE: calibrating subgrid feedback in cosmological simulations that include a cold gas phase
- WALLABY Pilot Survey & ASymba: Comparing HI Detection Asymmetries to the SIMBA Simulation