Fast and accurate frequency-dependent radiation transport for hydrodynamics simulations in massive star formation

Context: Radiative feedback plays a crucial role in the formation of massive stars. The implementation of a fast and accurate description of the proceeding thermodynamics in pre-stellar cores and evolving accretion disks is therefore a main effort in current hydrodynamics simulations. Aims: We introduce our newly implemented three-dimensional frequency dependent radiation transport algorithm for hydrodynamics simulations of spatial configurations with a dominant central source. Methods: The module combines the advantage of the speed of an approximate Flux Limited Diffusion (FLD) solver with the high accuracy of a frequency dependent first order ray-tracing routine. Results: We prove the viability of the scheme in a standard radiation benchmark test compared to a full frequency dependent Monte-Carlo based radiative transfer code. The setup includes a central star, a circumstellar flared disk, as well as an envelope. The test is performed for different optical depths. Considering the frequency dependence of the stellar irradiation, the temperature distributions can be described precisely in the optically thin, thick, and irradiated transition regions. Resulting radiative forces onto dust grains are reproduced with high accuracy. The achievable parallel speedup of the method imposes no restriction on further radiative (magneto-) hydrodynamics simulations. Conclusions: The proposed approximate radiation transport method enables frequency dependent radiation hydrodynamics studies of the evolution of pre-stellar cores and circumstellar accretion disks around an evolving massive star in a highly efficient and accurate manner.Comment: 16 pages, 11 figure

Accelerating Monte Carlo simulations with an NVIDIA® graphics processor

Modern graphics cards, commonly used in desktop computers, have evolved beyond a simple interface between processor and display to incorporate sophisticated calculation engines that can be applied to general purpose computing. The Monte Carlo algorithm for modelling photon transport in turbid media has been implemented on an NVIDIA® 8800gt graphics card using the CUDA toolkit. The Monte Carlo method relies on following the trajectory of millions of photons through the sample, often taking hours or days to complete. The graphics-processor implementation, processing roughly 110 million scattering events per second, was found to run more than 70 times faster than a similar, single-threaded implementation on a 2.67 GHz desktop computer

A public code for general relativistic, polarised radiative transfer around spinning black holes

Ray tracing radiative transfer is a powerful method for comparing theoretical models of black hole accretion flows and jets with observations. We present a public code, grtrans, for carrying out such calculations in the Kerr metric, including the full treatment of polarised radiative transfer and parallel transport along geodesics. The code is written in Fortran 90 and efficiently parallelises with OpenMP, and the full code and several components have Python interfaces. We describe several tests which are used for verifiying the code, and we compare the results for polarised thin accretion disc and semi-analytic jet problems with those from the literature as examples of its use. Along the way, we provide accurate fitting functions for polarised synchrotron emission and transfer coefficients from thermal and power law distribution functions, and compare results from numerical integration and quadrature solutions of the polarised radiative transfer equations. We also show that all transfer coefficients can play an important role in predicted images and polarisation maps of the Galactic center black hole, Sgr A*, at submillimetre wavelengths.Comment: 22 pages, 12 figures, submitted to MNRAS. code available at: github.com/jadexter/grtran

Embedded FIR filter design for real-time refocusing using a standard plenoptic video camera

Copyright 2014 Society of Photo-Optical Instrumentation Engineers and IS&T—The Society for Imaging Science and Technology. One print or electronic copy may be made for personal use only. Systematic reproduction and distribution, duplication of any material in this paper for a fee or for commercial purposes, or modification of the content of the paper are prohibited.A novel and low-cost embedded hardware architecture for real-time refocusing based on a standard plenoptic camera is presented in this study. The proposed layout design synthesizes refocusing slices directly from micro images by omitting the process for the commonly used sub-aperture extraction. Therefore, intellectual property cores, containing switch controlled Finite Impulse Response (FIR) filters, are developed and applied to the Field Programmable Gate Array (FPGA) XC6SLX45 from Xilinx. Enabling the hardware design to work economically, the FIR filters are composed of stored product as well as upsampling and interpolation techniques in order to achieve an ideal relation between image resolution, delay time, power consumption and the demand of logic gates. The video output is transmitted via High-Definition Multimedia Interface (HDMI) with a resolution of 720p at a frame rate of 60 fps conforming to the HD ready standard. Examples of the synthesized refocusing slices are presented