A Hybrid Cosmological Hydrodynamic/N-Body Code Based on a Weighted Essentially Nonoscillatory Scheme

• Published in: The Astrophysical journal Vol. 612; no. 1; pp. 1 - 13
• Main Authors: Feng, Long-Long, Shu, Chi-Wang, Zhang, Mengping
• Format: Journal Article
• Language: English
• Published: Chicago, IL IOP Publishing 01.09.2004; University of Chicago Press
• Subjects: Astronomy; Earth, ocean, space; Exact sciences and technology; Fundamental astronomy and astrophysics. Instrumentation, techniques, and astronomical observations; Mathematical procedures and computer techniques; Observation and data reduction techniques. Computer modeling and simulation; Algorithms; Digital simulation; Hydrodynamic model; N body system; Large-scale structure; Cosmological models; Finite difference method
• ISSN: 0004-637X; 1538-4357
• DOI: 10.1086/422513

We present a newly developed cosmological hydrodynamics code based on weighted essentially nonoscillatory (WENO) schemes for hyperbolic conservation laws. WENO is a higher order accurate finite difference scheme designed for problems with piecewise smooth solutions containing discontinuities and has been successfully applied for problems involving both shocks and complicated smooth solution structures. We couple hydrodynamics based on the WENO scheme with the standard Poisson solver-particle-mesh (PM) algorithm for evolving the self- gravitating system. The third-order low-storage total variation diminishing (TVD) Runge-Kutta scheme has been used for the time integration of the system. To test the accuracy and convergence rate of the code, we subject it to a number of typical tests, including the Sod shock tube in multiple dimensions, the Sedov blast wave, and formation of the Zel'dovich pancake. These tests validate the WENO hydrodynamics with a fast convergence rate and high accuracy. We also evolve a low-density flat cosmological model ( lambda CDM) to explore the validity of the code in practical simulations.