Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Athena Astrophysical MHD Code in Cylindrical Geometry
113
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A method for implementing cylindrical coordinates in the Athena magnetohydrodynamics (MHD) code is described. The extension follows the approach of Athenas original developers and has been designed to alter the existing Cartesian-coordinates code as minimally and transparently as possible. The numerical equations in cylindrical coordinates are formulated to maintain consistency with constrained transport, a central feature of the Athena algorithm, while making use of previously implemented code modules such as the Riemann solvers. Angular-momentum transport, which is critical in astrophysical disk systems dominated by rotation, is treated carefully. We describe modifications for cylindrical coordinates of the higher-order spatial reconstruction and characteristic evolution steps as well as the finite-volume and constrained transport updates. Finally, we present a test suite of standard and novel problems in one-, two-, and three-dimensions designed to validate our algorithms and implementation and to be of use to other code developers. The code is suitable for use in a wide variety of astrophysical applications and is freely available for download on the web.
rate research
Read More
The properties (temperature, density, chemical composition, velocity) of hot astrophysical plasma and the physical processes affecting them (heating/cooling, turbulence, shocks, acceleration) can be probed by high resolution X-ray spectroscopy, to be complemented by high spatial resolution imaging. The paper presents a status of the ESAs Advanced Telescope for High Energy Astrophysics (Athena) mission, particularly focusing on the science performance of its two focal plane instruments for the studies of extended X-ray sources: the Wide Field Imager (WFI) and the X-ray Integral Field Unit (X-IFU). This paper then provides a brief summary of the breakthroughs expected with Athena on the astrophysics of hot plasma, building on the vast heritage of the discoveries and revolutionary results obtained by Chandra and XMM-Newton in this field. As of November 12th, 2019, Athena successfully concluded its feasibility study, and has since then moved into the definition phase, with a launch date scheduled in the early 2030s.
We present a series of new, publicly available mock catalogs of X-ray selected active galactic nuclei (AGNs), non-active galaxies, and clusters of galaxies. They are based on up-to-date observational results on the demographic of extragalactic X-ray sources and their extrapolations.These mocks reach fluxes below 1E-20 erg s-1 cm-2 in the 0.5-2 keV band, i.e., more than an order of magnitude below the predicted limits of future deep fields, and therefore represent an important tool for simulating extragalactic X-ray surveys with both current and future telescopes. We use our mocks to perform a set of end-to-end simulations of X-ray surveys with the forthcoming Athena mission and with the AXIS probe, a sub-arcsecond resolution X-ray mission concept proposed to the Astro 2020 Decadal Survey. We find that these proposed, next generation surveys may transform our knowledge of the deep X-ray Universe. As an example, in a total observing time of 15,Ms, AXIS would detect ~225,000 AGNs and ~$50,000 non-active galaxies, reaching a flux limit f_0.5-2~5E-19 erg/s/cm2 in the 0.5-2 keV band, with an improvement of over an order of magnitude with respect to surveys with current X-ray facilities. Consequently, 90% of these sources would be detected for the first time in the X-rays. Furthermore, we show that deep and wide X-ray surveys with instruments like AXIS and Athena are expected to detect ~20,000 z>3 AGNs and ~250 sources at redshift z>6, thus opening a new window of knowledge on the evolution of AGNs over cosmic time and putting strong constraints on the predictions of theoretical models of black hole seed accretion in the early universe.
This paper describes the design and implementation of our new multi-group, multi-dimensional radiation hydrodynamics (RHD) code Fornax and provides a suite of code tests to validate its application in a wide range of physical regimes. Instead of focusing exclusively on tests of neutrino radiation hydrodynamics relevant to the core-collapse supernova problem for which Fornax is primarily intended, we present here classical and rigorous demonstrations of code performance relevant to a broad range of multi-dimensional hydrodynamic and multi-group radiation hydrodynamic problems. Our code solves the comoving-frame radiation moment equations using the M1 closure, utilizes conservative high-order reconstruction, employs semi-explicit matter and radiation transport via a high-order time stepping scheme, and is suitable for application to a wide range of astrophysical problems. To this end, we first describe the philosophy, algorithms, and methodologies of Fornax and then perform numerous stringent code tests, that collectively and vigorously exercise the code, demonstrate the excellent numerical fidelity with which it captures the many physical effects of radiation hydrodynamics, and show excellent strong scaling well above 100k MPI tasks.
We describe the magnetohydrodynamics (MHD) code CRONOS, which has been used in astrophysics and space physics studies in recent years. CRONOS has been designed to be easily adaptable to the problem at hand, where the user can expand or exchange core modules or add new functionality to the code. This modularity comes about through its implementation using a C++ class structure. The core components of the code include solvers for both hydrodynamical (HD) and MHD problems. These problems are solved on different rectangular grids, which currently support Cartesian, spherical, and cylindrical coordinates. CRONOS uses a finite-volume description with different approximate Riemann solvers that can be chosen at runtime. Here, we describe the implementation of the code with a view toward its ongoing development. We illustrate the codes potential through several (M)HD test problems and some astrophysical applications.
In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently make use of future exascale machines with their ever higher degree of parallel concurrency motivates the search for more efficient and more accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG) methods represent a promising class of methods in this regard, as they can be straightforwardly extended to arbitrarily high order while requiring only small stencils. Especially for applications involving comparatively smooth problems, higher-order approaches promise significant gains in computational speed for reaching a desired target accuracy. Here, we introduce our new astrophysical DG code TENET designed for applications in cosmology, and discuss our first results for 3D simulations of subsonic turbulence. We show that our new DG implementation provides accurate results for subsonic turbulence, at considerably reduced computational cost compared with traditional finite volume methods. In particular, we find that DG needs about 1.8 times fewer degrees of freedom to achieve the same accuracy and at the same time is more than 1.5 times faster, confirming its substantial promise for astrophysical applications.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
