Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCOSMOS: A Hybrid N-Body/Hydrodynamics Code for Cosmological Problems
109
0
0.0
(
0
)
Authors
P. M. Ricker
Ask ChatGPT about the research
No Arabic abstract
We describe a new hybrid N-body/hydrodynamical code based on the particle-mesh (PM) method and the piecewise-parabolic method (PPM) for use in solving problems related to the evolution of large-scale structure, galaxy clusters, and individual galaxies. The code, named COSMOS, possesses several new features which distinguish it from other PM-PPM codes. In particular, to solve the Poisson equation we have written a new multigrid solver which can determine the gravitational potential of isolated matter distributions and which properly takes into account the finite-volume discretization required by PPM. All components of the code are constructed to work with a nonuniform mesh, preserving second-order spatial differences. The PPM code uses vacuum boundary conditions for isolated problems, preventing inflows when appropriate. The PM code uses a second-order variable-timestep time integration scheme. Radiative cooling and cosmological expansion terms are included. COSMOS has been implemented for parallel computers using the Parallel Virtual Machine (PVM) library, and it features a modular design which simplifies the addition of new physics and the configuration of the code for different types of problems. We discuss the equations solved by COSMOS and describe the algorithms used, with emphasis on these features. We also discuss the results of tests we have performed to establish that COSMOS works and to determine its range of validity.
rate research
Read More
We present a newly developed cosmological hydrodynamics code based on weighted essentially non-oscillatory (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 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 accuracy and convergence rate of the code, we subject it to a number of typical tests including the Sod shock tube in multidimensions, the Sedov blast wave and formation of the Zeldovich pancake. These tests validate the WENO hydrodynamics with 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.
We describe a newly developed cosmological hydrodynamics code based on the weighted essentially non-oscillatory (WENO) schemes for hyperbolic conservation laws. High order finite difference WENO schemes are designed for problems with piecewise smooth solutions containing discontinuities, and have been successful in applications for problems involving both shocks and complicated smooth solution structures. We couple hydrodynamics based on the WENO scheme with standard Poisson solver - particle-mesh (PM) algorithm for evolving the self-gravitating system. A third order total variation diminishing (TVD) Runge-Kutta scheme has been used for time-integration of the system. We brief the implementation of numerical technique. The cosmological applications in simulating intergalactic medium and Ly$alpha$ forest in the CDM scenario are also presented.
Star clusters form via clustering star formation inside molecular clouds. In order to understand the dynamical evolution of star clusters in their early phase, in which star clusters are still embedded in their surrounding gas, we need an accurate integration of individual stellar orbits without gravitational softening in the systems including both gas and stars, as well as modeling individual stars with a realistic mass function. We develop a new tree-direct hybrid smoothed particle hydrodynamics/N-body code, ASURA+BRIDGE, in which stars are integrated using a direct N-body scheme or PeTar, a particle-particle particle-tree scheme code, without gravitational softening. In ASURA+BRIDGE, stars are assumed to have masses randomly drawn from a given initial mass function. With this code, we perform star-cluster formation simulations starting from molecular clouds without gravitational softening. We find that artificial dense cores in star-cluster centers due to the softening disappear when we do not use softening. We further demonstrate that star clusters are built up via mergers of smaller clumps. Star clusters formed in our simulations include some dynamically formed binaries with the minimum semi-major axes of a few au, and the binary fraction is higher for more massive stars.
We describe a novel N-body code designed for simulations of the central regions of galaxies containing massive black holes. The code incorporates Mikkolas algorithmic chain regularization scheme including post-Newtonian terms up to PN2.5 order. Stars moving beyond the chain are advanced using a fourth-order integrator with forces computed on a GRAPE board. Performance tests confirm that the hybrid code achieves better energy conservation, in less elapsed time, than the standard scheme and that it reproduces the orbits of stars tightly bound to the black hole with high precision. The hybrid code is applied to two sample problems: the effect of finite-N gravitational fluctuations on the orbits of the S-stars; and inspiral of an intermediate-mass black hole into the galactic center.
We describe the numerical code N-MODY, a parallel particle-mesh code for collisionless N-body simulations in modified Newtonian dynamics (MOND). N-MODY is based on a numerical potential solver in spherical coordinates that solves the non-linear MOND field equation, and is ideally suited to simulate isolated stellar systems. N-MODY can be used also to compute the MOND potential of arbitrary static density distributions. A few applications of N-MODY indicate that some astrophysically relevant dynamical processes are profoundly different in MOND and in Newtonian gravity with dark matter.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
