Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDIAPHANE: a Portable Radiation Transport Library for Astrophysical Applications
80
0
0.0
(
0
)
Authors
Darren S. Reed
Ask ChatGPT about the research
No Arabic abstract
One of the most computationally demanding aspects of the hydrodynamical modelling of Astrophysical phenomena is the transport of energy by radiation or relativistic particles. Physical processes involving energy transport are ubiquitous and of capital importance in many scenarios ranging from planet formation to cosmic structure evolution, including explosive events like core collapse supernova or gamma-ray bursts. Moreover, the ability to model and hence understand these processes has often been limited by the approximations and incompleteness in the treatment of radiation and relativistic particles. The DIAPHANE project has focused in developing a portable and scalable library that handles the transport of radiation and particles (in particular neutrinos) independently of the underlying hydrodynamic code. In this work, we present the computational framework and the functionalities of the first version of the DIAPHANE library, which has been successfully ported to three different smoothed-particle hydrodynamic codes, GADGET2, GASOLINE and SPHYNX. We also present validation of different modules solving the equations of radiation and neutrino transport using different numerical schemes.
rate research
Read More
Astrophysical fluid flow studies often encompass a wide range of physical processes to account for the complexity of the system under consideration. In addition to gravity, a proper treatment of thermodynamic processes via continuum radiation transport and/or photoionization is becoming the state of the art. We present a major update of our continuum radiation transport module, MAKEMAKE, and a newly developed module for photoionization, SEDNA, coupled to the magnetohydrodynamics code PLUTO. These extensions are currently not publicly available; access can be granted on a case-by-case basis. We explain the theoretical background of the equations solved, elaborate on the numerical layout, and present a comprehensive test suite for radiation-ionization hydrodynamics. The grid-based radiation and ionization modules support static one-dimensional, two-dimensional, and three-dimensional grids in Cartesian, cylindrical, and spherical coordinates. Each module splits the radiation field into two components, one originating directly from a point source - solved using a ray-tracing scheme - and a diffuse component - solved with a three-dimensional flux-limited diffusion (FLD) solver. The FLD solver for the continuum radiation transport makes use of either the equilibrium one-temperature approach or the linearization two-temperature approach. The FLD solver for the photoionization module enables accounting for the temporal evolution of the radiation field from direct recombination of free electrons into hydrogens ground state as an alternative to the on-the-spot approximation. A brief overview of completed and ongoing scientific studies is given to explicitly illustrate the multipurpose nature of the numerical framework presented.
A long-standing problem of astrophysical research is how to simultaneously obtain spectra of thousands of sources randomly positioned in the field of view of a telescope. Digital Micromirror Devices, used as optical switches, provide a most powerful solution allowing to design a new generation of instruments with unprecedented capabilities. We illustrate the key factors (opto-mechanical, cryo-thermal, cosmic radiation environment,...) that constrain the design of DMD-based multi-object spectrographs, with particular emphasis on the IR spectroscopic channel onboard the EUCLID mission, currently considered by the European Space Agency for a 2017 launch date.
Searching for geometric objects that are close in space is a fundamental component of many applications. The performance of search algorithms comes to the forefront as the size of a problem increases both in terms of total object count as well as in the total number of search queries performed. Scientific applications requiring modern leadership-class supercomputers also pose an additional requirement of performance portability, i.e. being able to efficiently utilize a variety of hardware architectures. In this paper, we introduce a new open-source C++ search library, ArborX, which we have designed for modern supercomputing architectures. We examine scalable search algorithms with a focus on performance, including a highly efficient parallel bounding volume hierarchy implementation, and propose a flexible interface making it easy to integrate with existing applications. We demonstrate the performance portability of ArborX on multi-core CPUs and GPUs, and compare it to the state-of-the-art libraries such as Boost.Geometry.Index and nanoflann.
We describe URCHIN, a reverse ray tracing radiative transfer scheme optimised to model self-shielding from the post-reionisation ultraviolet (UV) background in cosmological simulations. The reverse ray tracing strategy provides several benefits over forward ray tracing codes including: (1) the preservation of adaptive density field resolution (2) completely uniform sampling of gas elements by rays; (3) the preservation of galilean invariance; (4) the ability to sample the UV background spectrum with hundreds of frequency bins; and (5) exact preservation of the input UV background spectrum and amplitude in optically thin gas. The implementation described here focuses on Smoothed Particle Hydrodynamics (SPH). However, the method can be applied to any density field representation in which resolution elements admit ray intersection tests and can be associated with optical depths. We characterise the errors in our implementation in stages beginning with comparison to known analytic solutions and ending with a realistic model of the z = 3 cosmological UV background incident onto a suite of spherically symmetric models of gaseous galactic halos.
Due to increase in computing power, high-order Eulerian schemes will likely become instrumental for the simulations of turbulence and magnetic field amplification in astrophysical fluids in the next years. We present the implementation of a fifth order weighted essentially non-oscillatory scheme for constrained-transport magnetohydrodynamics into the code WOMBAT. We establish the correctness of our implementation with an extensive number tests. We find that the fifth order scheme performs as accurately as a common second order scheme at half the resolution. We argue that for a given solution quality the new scheme is more computationally efficient than lower order schemes in three dimensions. We also establish the performance characteristics of the solver in the WOMBAT framework. Our implementation fully vectorizes using flattened arrays in thread-local memory. It performs at about 0.6 Million zones per second per node on Intel Broadwell. We present scaling tests of the code up to 98 thousand cores on the Cray XC40 machine Hazel Hen, with a sustained performance of about 5 percent of peak at scale.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
