We describe RAGE, the "Radiation Adaptive Grid Eulerian" radiation-hydrodynamics code, including its data structures, its parallelization strategy and performance, its hydrodynamic algorithm(s), its (gray) radiation diffusion algorithm, and some of the considerable amount of verification and validation efforts. The hydrodynamics is a basic Godunov solver, to which we have made significant improvements to increase the advection algorithm's robustness and to converge stiffnesses in the equation of state. Similarly, the radiation transport is a basic gray diffusion, but our treatment of the radiation-material coupling, wherein we converge nonlinearities in a novel manner to allow larger timesteps and more robust behavior, can be applied to any multi-group transport algorithm.
Abstract. We present calculations of the spatial and spectral distribution of the radio emission from a wide WR+OB collidingwind binary system based on high-resolution hydrodynamical simulations and solutions to the radiative transfer equation. We account for both thermal and synchrotron radio emission, free-free absorption in both the unshocked stellar wind envelopes and the shocked gas, synchrotron self-absorption, and the Razin effect. To calculate the synchrotron emission several simplifying assumptions are adopted: the relativistic particle energy density is a simple fraction of the thermal particle energy density, in equipartition with the magnetic energy density, and a power-law in energy. We also assume that the magnetic field is tangled such that the resulting emission is isotropic. The applicability of these calculations to modelling radio images and spectra of colliding-wind systems is demonstrated with models of the radio emission from the wide WR+OB binary WR 147. Its synchrotron spectrum follows a power-law between 5 and 15 GHz but turns down to below this at lower and higher frequencies. We find that while free-free opacity from the circum-binary stellar winds can potentially account for the low-frequency turnover, models that also include a combination of synchrotron self-absorption and Razin effect are favoured. We argue that the highfrequency turn down is a consequence of inverse-Compton cooling. We present our resulting spectra and intensity distributions, along with simulated MERLIN observations of these intensity distributions. From these we argue that the inclination of the WR 147 system to the plane of the sky is low. We summarise by considering extensions of the current model that are important for models of the emission from closer colliding wind binaries, in particular the dramatically varying radio emission of WR 140.
Abstract-Over the last few decades, rapid improvement of computer capabilities has allowed impact cratering to be modeled with increasing complexity and realism, and has paved the way for a new era of numerical modeling of the impact process, including full, three-dimensional (3D) simulations. When properly benchmarked and validated against observation, computer models offer a powerful tool for understanding the mechanics of impact crater formation. This work presents results from the first phase of a project to benchmark and validate shock codes. A variety of 2D and 3D codes were used in this study, from commercial products like AUTODYN, to codes developed within the scientific community like SOVA, SPH, ZEUS-MP, iSALE, and codes developed at U.S. National Laboratories like CTH, SAGE/RAGE, and ALE3D. Benchmark calculations of shock wave propagation in aluminum-on-aluminum impacts were performed to examine the agreement between codes for simple idealized problems. The benchmark simulations show that variability in code results is to be expected due to differences in the underlying solution algorithm of each code, artificial stability parameters, spatial and temporal resolution, and material models. Overall, the inter-code variability in peak shock pressure as a function of distance is around 10 to 20%. In general, if the impactor is resolved by at least 20 cells across its radius, the underestimation of peak shock pressure due to spatial resolution is less than 10%. In addition to the benchmark tests, three validation tests were performed to examine the ability of the codes to reproduce the time evolution of crater radius and depth observed in vertical laboratory impacts in water and two well-characterized aluminum alloys. Results from these calculations are in good agreement with experiments. There appears to be a general tendency of shock physics codes to underestimate the radius of the forming crater. Overall, the discrepancy between the model and experiment results is between 10 and 20%, similar to the inter-code variability.
Spectral and kinematic studies suggest that the nonthermal radio source Sgr A*, located at the center of the Milky Way, is a supermassive compact object with a mass ∼ 2 − 3 × 10 6 M ⊙ . Winds from nearby stars, located ≈ 0.06 pc to the east of Sgr A*, should, in the absence of any outflow from the putative black hole itself, be accreting onto this object. We report the results of the first 3D Bondi-Hoyle hydrodynamical numerical simulations of this process under the assumption that the Galactic center wind is generated by several different point sources (here assumed to be 10 pseudo-randomly placed stars). Our results show that the accretion rate onto the central object can be higher than in the case of a uniform flow since wind-wind shocks dissipate some of the bulk kinetic energy and lead to a higher capture rate for the gas. However, even for this highly non-uniform medium, most of the accreting gas carries with it a relatively low level of specific angular momentum, though large transient fluctuations can occur. Additionally, the post-bow-shock focusing of the gas can be substantially different than that for a uniform flow, but it depends strongly on the stellar spatial distribution. We discuss how this affects the morphology of the gas in the inner 0.15 pc of the Galaxy and the consequences for accretion disk models of Sgr A*.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.