Grid Noise in Moving Mesh Codes: Fixing The Volume Inconsistency Problem

Abstract

Current Voronoi based moving mesh hydro codes suffer from "grid noise". We identify the cause of this noise as the volume inconsistency error, where the volume that is transferred between cells is inconsistent with the hydrodynamical calculations. As a result, the codes do not achieve second order convergence. In this paper we describe how a simple fix allows Voronoi based moving mesh codes to attain second order convergence. The fix is based on the understanding that the volume exchanged between cells should be consistent with the hydrodynamical calculations. We benchmark our fix with three test problems and show that it can significantly improve the computational accuracy. We also examine the effect of initial mesh initialization and present an improved model for the Green-Gauss based gradient estimator.

submitted to MNRAS. Comments welcome

References in corpus (8)

• E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• Improving the convergence properties of the moving-mesh code AREPO
• Subsonic turbulence in smoothed particle hydrodynamics and moving-mesh simulations
• TESS: A Relativistic Hydrodynamics Code on a Moving Voronoi Mesh
• A discontinuous Galerkin method for solving the fluid and MHD equations in astrophysical simulations
• Rich: Open Source Hydrodynamic Simulation on a Moving Voronoi Mesh
• High-Frequency Voronoi Noise Reduced by Smoothed Mesh Motion

Cited by in corpus (8)

• DISCO: A 3D Moving-Mesh Magnetohydrodynamics Code Designed for the Study of Astrophysical Disks
• Common Envelope Evolution on a Moving Mesh
• The Multi-Dimensional Structure of Radiative Shocks: Suppressed Thermal X-rays and Relativistic Ion Acceleration
• A Moving Mesh Hydrodynamic Solver for ChaNGa
• Moving and Reactive Boundary Conditions in Moving-Mesh Hydrodynamics
• Multidimensional upwind hydrodynamics on unstructured meshes using Graphics Processing Units I. Two-dimensional uniform meshes
• The Effect of Impact Parameter on Tidal Disruption Events
• Multigroup Radiation Diffusion on a Moving Mesh: Implementation in RICH and Application to Tidal Disruption Events