Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userArepo-MCRT: Monte Carlo Radiation Hydrodynamics on a Moving Mesh
112
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present Arepo-MCRT, a novel Monte Carlo radiative transfer (MCRT) radiation-hydrodynamics (RHD) solver for the unstructured moving-mesh code Arepo. Our method is designed for general multiple scattering problems in both optically thin and thick conditions. We incorporate numerous efficiency improvements and noise reduction schemes to help overcome efficiency barriers that typically inhibit convergence. These include continuous absorption and energy deposition, photon weighting and luminosity boosting, local packet merging and splitting, path-based statistical estimators, conservative (face-centered) momentum coupling, adaptive convergence between time steps, implicit Monte Carlo algorithms for thermal emission, and discrete-diffusion Monte Carlo techniques for unresolved scattering, including a novel advection scheme. We primarily focus on the unique aspects of our implementation and discussions of the advantages and drawbacks of our methods in various astrophysical contexts. Finally, we consider several test applications including the levitation of an optically thick layer of gas by trapped infrared radiation. We find that the initial acceleration phase and revitalized second wind are connected via self-regulation of the RHD coupling, such that the RHD method accuracy and simulation resolution each leave important imprints on the long-term behavior of the gas.
rate research
Read More
We present the public Monte Carlo photoionization and moving-mesh radiation hydrodynamics code CMacIonize, which can be used to simulate the self-consistent evolution of HII regions surrounding young O and B stars, or other sources of ionizing radiation. The code combines a Monte Carlo photoionization algorithm that uses a complex mix of hydrogen, helium and several coolants in order to self-consistently solve for the ionization and temperature balance at any given type, with a standard first order hydrodynamics scheme. The code can be run as a post-processing tool to get the line emission from an existing simulation snapshot, but can also be used to run full radiation hydrodynamical simulations. Both the radiation transfer and the hydrodynamics are implemented in a general way that is independent of the grid structure that is used to discretize the system, allowing it to be run both as a standard fixed grid code, but also as a moving-mesh code.
Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code textsc{Arepo} and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems textsc{Arepo} only reaches first-order convergence for its original formulation. This can be rectified by simple modifications we propose to the time integration scheme and the spatial gradient estimates of the code, both improving the accuracy of the code. We demonstrate that the new implementation is indeed second-order accurate under the $L^1$ norm, and in particular substantially improves conservation of angular momentum. Interestingly, whereas these improvements can significantly change the results of smooth test problems, we also find that cosmological simulations of galaxy formation are unaffected, demonstrating that the numerical errors eliminated by the new formulation do not impact these simulations. In contrast, simulations of binary stars followed over a large number of orbital times are strongly affected, as here it is particularly crucial to avoid a long-term build up of errors in angular momentum conservation.
Radiative transfer has a strong impact on the collapse and the fragmentation of prestellar dense cores. We present the radiation-hydrodynamics solver we designed for the RAMSES code. The method is designed for astrophysical purposes, and in particular for protostellar collapse. We present the solver, using the co-moving frame to evaluate the radiative quantities. We use the popular flux limited diffusion approximation, under the grey approximation (one group of photon). The solver is based on the second-order Godunov scheme of RAMSES for its hyperbolic part, and on an implicit scheme for the radiation diffusion and the coupling between radiation and matter. We report in details our methodology to integrate the RHD solver into RAMSES. We test successfully the method against several conventional tests. For validation in 3D, we perform calculations of the collapse of an isolated 1 M_sun prestellar dense core, without rotation. We compare successfully the results with previous studies using different models for radiation and hydrodynamics. We have developed a full radiation hydrodynamics solver in the RAMSES code, that handles adaptive mesh refinement grids. The method is a combination of an explicit scheme and an implicit scheme, accurate to the second-order in space. Our method is well suited for star formation purposes. Results of multidimensional dense core collapse calculations with rotation are presented in a companion paper.
We present a detailed comparison between the well-known SPH code GADGET and the new moving-mesh code AREPO on a number of hydrodynamical test problems. Through a variety of numerical experiments we establish a clear link between test problems and systematic numerical effects seen in cosmological simulations of galaxy formation. Our tests demonstrate deficiencies of the SPH method in several sectors. These accuracy problems not only manifest themselves in idealized hydrodynamical tests, but also propagate to more realistic simulation setups of galaxy formation, ultimately affecting gas properties in the full cosmological framework, as highlighted in papers by Vogelsberger et al. (2011) and Keres et al. (2011). We find that an inadequate treatment of fluid instabilities in GADGET suppresses entropy generation by mixing, underestimates vorticity generation in curved shocks and prevents efficient gas stripping from infalling substructures. In idealized tests of inside-out disk formation, the convergence rate of gas disk sizes is much slower in GADGET due to spurious angular momentum transport. In simulations where we follow the interaction between a forming central disk and orbiting substructures in a halo, the final disk morphology is strikingly different. In AREPO, gas from infalling substructures is readily depleted and incorporated into the host halo atmosphere, facilitating the formation of an extended central disk. Conversely, gaseous sub-clumps are more coherent in GADGET simulations, morphologically transforming the disk as they impact it. The numerical artefacts of the SPH solver are particularly severe for poorly resolved flows, and thus inevitably affect cosmological simulations due to their hierarchical nature. Our numerical experiments clearly demonstrate that AREPO delivers a physically more reliable solution.
In certain astrophysical systems the commonly employed ideal magnetohydrodynamics (MHD) approximation breaks down. Here, we introduce novel explicit and implicit numerical schemes of ohmic resistivity terms in the moving-mesh code AREPO. We include these non-ideal terms for two MHD techniques: the Powell 8-wave formalism and a constrained transport scheme, which evolves the cell-centred magnetic vector potential. We test our implementation against problems of increasing complexity, such as one- and two-dimensional diffusion problems, and the evolution of progressive and stationary Alfven waves. On these test problems, our implementation recovers the analytic solutions to second-order accuracy. As first applications, we investigate the tearing instability in magnetized plasmas and the gravitational collapse of a rotating magnetized gas cloud. In both systems, resistivity plays a key role. In the former case, it allows for the development of the tearing instability through reconnection of the magnetic field lines. In the latter, the adopted (constant) value of ohmic resistivity has an impact on both the gas distribution around the emerging protostar and the mass loading of magnetically driven outflows. Our new non-ideal MHD implementation opens up the possibility to study magneto-hydrodynamical systems on a moving mesh beyond the ideal MHD approximation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
