A High-Resolution Adaptive Moving Mesh Hydrodynamic Algorithm

An algorithm for simulating self-gravitating cosmological astrophysical fluids is presented. The advantages include a large dynamic range, parallelizability, high resolution per grid element, and fast execution speed. The code is based on a finite volume flux-conservative total variation diminishing (TVD) scheme for the shock-capturing hydro and an iterative multigrid solver for the gravity. The grid is a time-dependent field, whose motion is described by a generalized potential flow. Approximately constant mass per cell can be obtained, which provides all the advantages of a Lagrangian scheme. The grid deformation combined with appropriate limiting and smoothing schemes guarantees a regular and well-behaved grid geometry, in which nearest neighbor relationships remain constant. The full hydrodynamic fluid equations are implemented in the curvilinear moving grid, which allows for arbitrary fluid flow relative to the grid geometry. This combination retains all the advantages of the grid-based schemes including high speed per fluid element and a rapid gravity solver. The current implementation is described, and empirical simulation results are presented. Accurate execution speed calculations are given in terms of floating point operations per time step per grid cell. This code is freely available to the community.

[1]  J. Michael Owen,et al.  Adaptive smoothed particle hydrodynamics, with application to cosmology: Methodology , 1996 .

[2]  Joe F. Thompson,et al.  Numerical grid generation: Foundations and applications , 1985 .

[3]  O. C. Zienkiewicz,et al.  An adaptive finite element procedure for compressible high speed flows , 1985 .

[4]  Z. Xin,et al.  The relaxation schemes for systems of conservation laws in arbitrary space dimensions , 1995 .

[5]  William H. Press,et al.  The Cosmological constant , 1992 .

[6]  T. Mouschovias,et al.  Ambipolar diffusion and star formation: Formation and contraction of axisymmetric cloud cores. I. Formulation of the problem and method of solution , 1992 .

[7]  H. C. Yee,et al.  A class of high resolution explicit and implicit shock-capturing methods , 1989 .

[8]  Nickolay Y. Gnedin,et al.  Softened Lagrangian Hydrodynamics for Cosmology , 1995 .

[9]  J. W. Eastwood,et al.  On the clustering of particles in an expanding Universe , 1981 .

[10]  R. Cen.  A hydrodynamic approach to cosmology - Methodology , 1992 .

[11]  A. Evrard,et al.  A comparison of cosmological hydrodynamic codes , 1994, astro-ph/9404014.

[12]  Jeremiah P. Ostriker,et al.  A Cosmological Hydrodynamic Code Based on the Total Variation Diminishing Scheme , 1993 .

[13]  J. Hansen,et al.  Review of some adaptive node-movement techniques in finite-element and finite-difference solutions of partial differential equations , 1991 .

[14]  R. Trapp,et al.  Tornado-like vortexgenesis in a simplified numerical model , 1995 .

[15]  J. Pain,et al.  Gas Dynamics , 2018, Nature.

[16]  H. C. Yee,et al.  Implicit TVD schemes for hyperbolic conservation laws in curvilinear coordinates , 1987 .

[17]  A. Szalay,et al.  The statistics of peaks of Gaussian random fields , 1986 .

[18]  L. Hernquist,et al.  TREESPH: A Unification of SPH with the Hierarchical Tree Method , 1989 .

[19]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[20]  Jeremiah P. Ostriker,et al.  The universe in a box : thermal effects in the standard cold dark matter scenario , 1990 .

[21]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[22]  U. Pen.  A Linear Moving Adaptive Particle-Mesh N-Body Algorithm , 1995 .

[23]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[24]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[25]  R.A. Fiedler.  Optimization and Scaling of Shared-Memory and Message-Passing Implementations of the Zeus Hydrodynamics Algorithm , 1997, ACM/IEEE SC 1997 Conference (SC'97).

[26]  Guohong Xu.  Hydrodynamic and N-body schemes on an unstructured, adaptive mesh with applications to cosmological simulations , 1997 .

[27]  Richard Courant,et al.  Supersonic Flow And Shock Waves , 1948 .

[28]  A. Kashlinsky,et al.  Large-scale structure in the Universe , 1991, Nature.

[29]  S. White,et al.  Simulations of dissipative galaxy formation in hierarchically clustering universes – I: Tests of the code , 1993 .

[30]  T. Mouschovias,et al.  Ambipolar diffusion and star formation : formation and contraction of axisymmetric cloud cores. II: Results , 1993 .

[31]  Jeremiah P. Ostriker,et al.  A piecewise parabolic method for cosmological hydrodynamics , 1995 .

[32]  T. N. Stevenson,et al.  Fluid Mechanics , 2021, Nature.

[33]  C. Hirsch.  Numerical Computation of Internal and External Flows, Volume 2: Computational Methods for Inviscid and Viscous Flows , 1990 .

[34]  Brian H. Fiedler,et al.  A Fast Dynamic Grid Adaption Scheme for Meteorological Flows , 1993 .