Block structured adaptive mesh and time refinement for hybrid, hyperbolic + N-body systems

We present a new numerical algorithm for the solution of coupled collisional and collisionless systems, based on the block structured adaptive mesh and time refinement strategy (AMR). We describe the issues associated with the discretization of the system equations and the synchronization of the numerical solution on the hierarchy of grid levels. We implement a code based on a higher order, conservative and directionally unsplit Godunov's method for hydrodynamics; a symmetric, time centered modified symplectic scheme for collisionless component; and a multilevel, multigrid relaxation algorithm for the elliptic equation coupling the two components. Numerical results that illustrate the accuracy of the code and the relative merit of various implemented schemes are also presented.

[1]  Michael L. Norman.  Cosmological Simulations of X-ray Clusters: The Quest for Higher Resolution and Essential Physics , 2003 .

[2]  R. Klein,et al.  On the hydrodynamic interaction of shock waves with interstellar clouds. 1: Nonradiative shocks in small clouds , 1994 .

[3]  M. Minion.  A Projection Method for Locally Refined Grids , 1996 .

[4]  Francesco Miniati.  Glimm-Godunov's method for cosmic-ray-hydrodynamics , 2007, J. Comput. Phys..

[5]  R. Teyssier.  Cosmological hydrodynamics with adaptive mesh refinement - A new high resolution code called RAMSES , 2001, astro-ph/0111367.

[6]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[7]  Phillip Colella,et al.  A cell-centered adaptive projection method for the incompressible Navier-Stokes equations in three dimensions , 2007, J. Comput. Phys..

[8]  James Binney,et al.  Multi‐level adaptive particle mesh (MLAPM): a c code for cosmological simulations , 2001, astro-ph/0103503.

[9]  F. Pearce,et al.  Hydra: an Adaptive{mesh Implementation of P 3 M{sph , 1994 .

[10]  Richard H. White,et al.  The Hydrodynamic Behavior of Supernovae Explosions , 1964 .

[11]  Michael L. Norman,et al.  The Formation of the First Star in the Universe , 2001, Science.

[12]  Francesco Miniati.  COSMOCR: A numerical code for cosmic ray studies in computational cosmology , 2001 .

[13]  Phillip Colella,et al.  A modified higher order Godunov's scheme for stiff source conservative hydrodynamics , 2007, J. Comput. Phys..

[14]  Vicent Quilis,et al.  A new multidimensional adaptive mesh refinement hydro + gravity cosmological code , 2004 .

[15]  P. Colella,et al.  A Conservative Adaptive Projection Method for the Variable Density Incompressible Navier-Stokes Equations , 1998 .

[16]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[17]  Daniel F. Martin,et al.  Solving Poisson's Equation using Adaptive Mesh Renemen t , 1996 .

[18]  Richard I. Klein,et al.  Single and multiple star formation in turbulent molecular cloud cores , 2002 .

[19]  Richard I. Klein,et al.  Self-gravitational Hydrodynamics with Three-dimensional Adaptive Mesh Refinement: Methodology and Applications to Molecular Cloud Collapse and Fragmentation , 1998 .

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

[21]  B. Fryxell,et al.  On the Cellular Structure of Carbon Detonations , 2000 .

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

[23]  J. Saltzman,et al.  An unsplit 3D upwind method for hyperbolic conservation laws , 1994 .

[24]  S. White,et al.  A Universal Density Profile from Hierarchical Clustering , 1996, astro-ph/9611107.

[25]  R. Cen,et al.  Properties of Cosmic Shock Waves in Large-Scale Structure Formation , 2000, astro-ph/0005444.

[26]  Technology of China,et al.  A Hybrid Cosmological Hydrodynamic/N-Body Code Based on a Weighted Essentially Nonoscillatory Scheme , 2004 .

[27]  Janka,et al.  Nucleosynthesis and Clump Formation in a Core-Collapse Supernova. , 2000, The Astrophysical journal.

[28]  F. Pearce,et al.  Hydra: An Adaptive--Mesh Implementation of PPPM--SPH , 1994 .

[29]  G. Lake,et al.  Time stepping N body simulations , 1997, astro-ph/9710043.

[30]  A. S. Almgren,et al.  Low mach number modeling of type Ia supernovae. I. Hydrodynamics , 2005 .

[31]  M. Steinmetz,et al.  The Santa Barbara Cluster Comparison Project: A Comparison of Cosmological Hydrodynamics Solutions , 1999, astro-ph/9906160.

[32]  P. Colella.  Multidimensional upwind methods for hyperbolic conservation laws , 1990 .

[33]  G. Bryan,et al.  A Hybrid AMR Application for Cosmology and Astrophysics , 1997, astro-ph/9710187.

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

[35]  Richard I. Klein,et al.  Recent Advances in the Collapse and Fragmentation of Turbulent Molecular Cloud Cores: The Formation of Low Mass Stars , 2003 .

[36]  S. Cole,et al.  Merger rates in hierarchical models of galaxy formation , 1993 .

[37]  William Kahan,et al.  Unconventional Schemes for a Class of Ordinary Differential Equations-With Applications to the Korteweg-de Vries Equation , 1997 .

[38]  R. Walder,et al.  KNOTS, FILAMENTS, AND TURBULENCE IN RADIATIVE SHOCKS , 1998 .

[39]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[40]  Melott.  Comment on "Nonlinear gravitational clustering in cosmology." , 1986, Physical review letters.

[41]  Eric F Darve,et al.  Author ' s personal copy A hybrid method for the parallel computation of Green ’ s functions , 2009 .

[42]  Stephen L. W. McMillan,et al.  An O(N log N) integration scheme for collisional stellar systems , 1993 .

[43]  Y. Zel’dovich.  Gravitational instability: An Approximate theory for large density perturbations , 1969 .

[44]  Richard I. Klein,et al.  The Jeans Condition: A New Constraint on Spatial Resolution in Simulations of Isothermal Self-Gravitational Hydrodynamics , 1997 .

[45]  R. Klein,et al.  On the Hydrodynamic Interaction of Shock Waves with Interstellar Clouds. II. The Effect of Smooth Cloud Boundaries on Cloud Destruction and Cloud Turbulence , 2005, astro-ph/0511016.

[46]  Daniel F. Martin,et al.  A Cell-Centered Adaptive Projection Method for the Incompressible Euler Equations , 2000 .

[47]  P. Peebles.  Principles of Physical Cosmology , 1993 .

[48]  Michael Zingale,et al.  Low Mach Number Modeling of Type Ia Supernovae , 2005 .

[49]  P. Colella,et al.  A conservative three-dimensional Eulerian method for coupled solid-fluid shock capturing , 2002 .

[50]  Shaun Cole,et al.  Merger rates in hierarchical models of galaxy formation – II. Comparison with N-body simulations , 1994 .

[51]  Jean-Luc Vay,et al.  Mesh Refinement for Particle-In-Cell Plasma Simulations: Applications to - and benefits for - Heavy-Ion-Fusion , 2002 .