Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRAMSES-CH: A New Chemodynamical Code for Cosmological Simulations
133
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a new chemodynamical code - Ramses-CH - for use in simulating the self-consistent evolution of chemical and hydrodynamical properties of galaxies within a fully cosmological framework. We build upon the adaptive mesh refinement code Ramses, which includes a treatment of self-gravity, hydrodynamics, star formation, radiative cooling, and supernovae feedback, to trace the dominant isotopes of C, N, O, Ne, Mg, Si, and Fe. We include the contribution of Type Ia and II supernovae, in addition to low- and intermediate-mass asymptotic giant branch stars, relaxing the instantaneous recycling approximation. The new chemical evolution modules are highly flexible and portable, lending themselves to ready exploration of variations in the underpining stellar and nuclear physics. We apply Ramses-CH to the cosmological simulation of a typical L* galaxy, demonstrating the successful recovery of the basic empirical constraints regarding, [{alpha}/Fe]-[Fe/H] and Type Ia/II supernovae rates.
rate research
Read More
We present a new chemodynamical code based on the adaptive mesh refinement code RAMSES. The new code uses Eulerian hydrodynamics and N-body dynamics in a cosmological framework to trace the production and advection of several chemical species. It is the first such code to follow the self-consistent evolution of chemical elements in cosmological volumes while maintaining sub-kiloparsec resolution. The code will be used to simulate disk galaxies and explore the influence of chemical evolution models and star formation on galactic abundance ratios.
A new N-body and hydrodynamical code, called RAMSES, is presented. It has been designed to study structure formation in the universe with high spatial resolution. The code is based on Adaptive Mesh Refinement (AMR) technique, with a tree based data structure allowing recursive grid refinements on a cell-by-cell basis. The N-body solver is very similar to the one developed for the ART code (Kravtsov et al. 97), with minor differences in the exact implementation. The hydrodynamical solver is based on a second-order Godunov method, a modern shock-capturing scheme known to compute accurately the thermal history of the fluid component. The accuracy of the code is carefully estimated using various test cases, from pure gas dynamical tests to cosmological ones. The specific refinement strategy used in cosmological simulations is described, and potential spurious effects associated to shock waves propagation in the resulting AMR grid are discussed and found to be negligible. Results obtained in a large N-body and hydrodynamical simulation of structure formation in a low density LCDM universe are finally reported, with 256^3 particles and 4.1 10^7 cells in the AMR grid, reaching a formal resolution of 8192^3. A convergence analysis of different quantities, such as dark matter density power spectrum, gas pressure power spectrum and individual haloes temperature profiles, shows that numerical results are converging down to the actual resolution limit of the code, and are well reproduced by recent analytical predictions in the framework of the halo model.
We investigate how the stellar and gas-phase He abundances evolve as functions of time within simulated star-forming disc galaxies with different star formation histories. We make use of a cosmological chemodynamical simulation for galaxy formation and evolution, which includes star formation, as well as energy and chemical enrichment feedback from asymptotic giant branch stars, core-collapse supernovae, and Type Ia supernovae. The predicted relations between the He mass fraction, $Y$, and the metallicity, $Z$, in the interstellar medium of our simulated disc galaxies depend on the past galaxy star formation history. In particular, $dY/dZ$ is not constant and evolves as a function of time, depending on the specific chemical element that we choose to trace $Z$; in particular, $dY/dX_{text{O}}$ and $dY/dX_{text{C}}$ increase as functions of time, whereas $dY/dX_{text{N}}$ decreases. In the gas-phase, we find negative radial gradients of $Y$, due to the inside-out growth of our simulated galaxy discs as a function of time; this gives rise to longer chemical enrichment time scales in the outer galaxy regions, where we find lower average values for $Y$ and $Z$. Finally, by means of chemical evolution models, in the galactic bulge and inner disc, we predict steeper $Y$ versus age relations at high $Z$ than in the outer galaxy regions. We conclude that, for calibrating the assumed $Y$-$Z$ relation in stellar models, C, N, and C+N are better proxies for the metallicity than O, because they show steeper and less scattered relations.
We present a ray tracing code to compute integrated cosmological observables on the fly in AMR N-body simulations. Unlike conventional ray tracing techniques, our code takes full advantage of the time and spatial resolution attained by the N-body simulation by computing the integrals along the line of sight on a cell-by-cell basis through the AMR simulation grid. Moroever, since it runs on the fly in the N-body run, our code can produce maps of the desired observables without storing large (or any) amounts of data for post-processing. We implemented our routines in the RAMSES N-body code and tested the implementation using an example of weak lensing simulation. We analyse basic statistics of lensing convergence maps and find good agreement with semi-analytical methods. The ray tracing methodology presented here can be used in several cosmological analysis such as Sunyaev-Zeldovich and integrated Sachs-Wolfe effect studies as well as modified gravity. Our code can also be used in cross-checks of the more conventional methods, which can be important in tests of theory systematics in preparation for upcoming large scale structure surveys.
We present EvoL, the new release of the Padova N-body code for cosmological simulations of galaxy formation and evolution. In this paper, the basic Tree + SPH code is presented and analysed, together with an overview on the software architectures. EvoL is a flexible parallel Fortran95 code, specifically designed for simulations of cosmological structure formation on cluster, galactic and sub-galactic scales. EvoL is a fully Lagrangian self-adaptive code, based on the classical Oct-tree and on the Smoothed Particle Hydrodynamics algorithm. It includes special features such as adaptive softening lengths with correcting extra-terms, and modern formulations of SPH and artificial viscosity. It is designed to be run in parallel on multiple CPUs to optimize the performance and save computational time. We describe the code in detail, and present the results of a number of standard hydrodynamical tests.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
