Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA Hybrid N-Body Code Incorporating Algorithmic Regularization and Post-Newtonian Forces
177
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
In this paper, we introduce PoMiN, a lightweight $N$-body code based on the post-Minkowskian $N$-body Hamiltonian of Ledvinka et. al., which includes general relativistic effects up to first order in Newtons constant $G$, and all orders in the speed of light $c$. PoMiN is written in C and uses a fourth-order Runge-Kutta integration scheme. PoMiN has also been written to handle an arbitrary number of particles (both massive and massless), with a computational complexity that scales as $O(N^2)$. We describe the methods we used to simplify and organize the Hamiltonian, and the tests we performed (convergence, conservation, and analytical comparison tests) to validate the code.
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 determine the gravitational interaction between two compact bodies up to the sixth power in Newtons constant GN, in the static limit. This result is achieved within the effective field theory approach to General Relativity, and exploits a manifest factorization property of static diagrams which allows to derive static post Newtonian (PN) contributions of (2n+1)-order in terms of lower order ones. We recompute in this fashion the 1PN and 3PN static potential, and present the novel 5PN contribution.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
