Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRAMSES-CH: A New Chemodynamics Code
116
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
We trace the formation and advection of several elements within a cosmological adaptive mesh refinement simulation of an L* galaxy. We use nine realisations of the same initial conditions with different stellar Initial Mass Functions (IMFs), mass limits for type-II and type-Ia supernovae (SNII, SNIa) and stellar lifetimes to constrain these sub-grid phenomena. Our code includes self-gravity, hydrodynamics, star formation, radiative cooling and feedback from multiple sources within a cosmological framework. Under our assumptions of nucleosynthesis we find that SNII with progenitor masses of up to 100 Msun are required to match low metallicity gas oxygen abundances. Tardy SNIa are necessary to reproduce the classical chemical evolution knee in [O/Fe]-[Fe/H]: more prompt SNIa delayed time distributions do not reproduce this feature. Within our framework of hydrodynamical mixing of metals and galaxy mergers we find that chemical evolution is sensitive to the shape of the IMF and that there exists a degeneracy with the mass range of SNII. We look at the abundance plane and present the properties of different regions of the plot, noting the distinct chemical properties of satellites and a series of nested discs that have greater velocity dispersions, are more alpha-rich and metal poor with age.
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 have implemented non-ideal Magneto-Hydrodynamics (MHD) effects in the Adaptive Mesh Refinement (AMR) code RAMSES, namely ambipolar diffusion and Ohmic dissipation, as additional source terms in the ideal MHD equations. We describe in details how we have discretized these terms using the adaptive Cartesian mesh, and how the time step is diminished with respect to the ideal case, in order to perform a stable time integration. We have performed a large suite of test runs, featuring the Barenblatt diffusion test, the Ohmic diffusion test, the C-shock test and the Alfven wave test. For the latter, we have performed a careful truncation error analysis to estimate the magnitude of the numerical diffusion induced by our Godunov scheme, allowing us to estimate the spatial resolution that is required to address non-ideal MHD effects reliably. We show that our scheme is second-order accurate, and is therefore ideally suited to study non-ideal MHD effects in the context of star formation and molecular cloud dynamics.
High levels of deuterium fractionation of $rm N_2H^+$ (i.e., $rm D_{frac}^{N_2H^+} gtrsim 0.1$) are often observed in pre-stellar cores (PSCs) and detection of $rm N_2D^+$ is a promising method to identify elusive massive PSCs. However, the physical and chemical conditions required to reach such high levels of deuteration are still uncertain, as is the diagnostic utility of $rm N_2H^+$ and $rm N_2D^+$ observations of PSCs. We perform 3D magnetohydrodynamics simulations of a massive, turbulent, magnetised PSC, coupled with a sophisticated deuteration astrochemical network. Although the core has some magnetic/turbulent support, it collapses under gravity in about one freefall time, which marks the end of the simulations. Our fiducial model achieves relatively low $rm D_{frac}^{N_2H^+} sim 0.002$ during this time. We then investigate effects of initial ortho-para ratio of $rm H_2$ ($rm OPR^{H_2}$), temperature, cosmic ray (CR) ionization rate, CO and N-species depletion factors and prior PSC chemical evolution. We find that high CR ionization rates and high depletion factors allow the simulated $rm D_{frac}^{N_2H^+}$ and absolute abundances to match observational values within one freefall time. For $rm OPR^{H_2}$, while a lower initial value helps the growth of $rm D_{frac}^{N_2H^+}$, the spatial structure of deuteration is too widespread compared to observed systems. For an example model with elevated CR ionization rates and significant heavy element depletion, we then study the kinematic and dynamic properties of the core as traced by its $rm N_2D^+$ emission. The core, undergoing quite rapid collapse, exhibits disturbed kinematics in its average velocity map. Still, because of magnetic support, the core often appears kinematically sub-virial based on its $rm N_2D^+$ velocity dispersion.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
