اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدDART-Ray: a 3D ray-tracing radiative transfer code for calculating the propagation of light in dusty galaxies
43
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We present DART-Ray, a new ray-tracing 3D dust radiative transfer (RT) code designed specifically to calculate radiation field energy density (RFED) distributions within dusty galaxy models with arbitrary geometries. In this paper we introduce the basic algorithm implemented in DART-Ray which is based on a pre-calculation of a lower limit for the RFED distribution. This pre-calculation allows us to estimate the extent of regions around the radiation sources within which these sources contribute significantly to the RFED. In this way, ray-tracing calculations can be restricted to take place only within these regions, thus substantially reducing the computational time compared to a complete ray-tracing RT calculation. Anisotropic scattering is included in the code and handled in a similar fashion. Furthermore, the code utilizes a Cartesian adaptive spatial grid and an iterative method has been implemented to optimize the angular densities of the rays originated from each emitting cell. In order to verify the accuracy of the RT calculations performed by DART-Ray, we present results of comparisons with solutions obtained using the DUSTY 1D RT code for a dust shell illuminated by a central point source and existing 2D RT calculations of disc galaxies with diffusely distributed stellar emission and dust opacity. Finally, we show the application of the code on a spiral galaxy model with logarithmic spiral arms in order to measure the effect of the spiral pattern on the attenuation and RFED.
قيم البحث
اقرأ أيضاً
We describe a novel adaptive ray tracing scheme to solve the equation of radiative transfer around point sources in hydrodynamical simulations. The angular resolution adapts to the local hydrodynamical resolution and hence is of use for adaptive mesh
es as well as adaptive smooth particle hydrodynamical simulations. Recursive creation of rays ensures ease of implementation. The multiple radial integrations needed to solve the time dependent radiative transfer are sped up significantly using a quad-tree once the rays are cast. Simplifications advantageous for methods with one radiation source are briefly discussed. The suggested method is easily generalized to speed up Monte Carlo radiative transfer techniques. In summary a nearly optimal use of long characteristics is presented and aspects of its implementation and comparison to other methods are given.
We present a ray tracing code to compute integrated cosmological observables on the fly in AMR N-body simulations. Unlike conventional ray tracing techniques, our code takes full advantage of the time and spatial resolution attained by the N-body sim
ulation by computing the integrals along the line of sight on a cell-by-cell basis through the AMR simulation grid. Moroever, since it runs on the fly in the N-body run, our code can produce maps of the desired observables without storing large (or any) amounts of data for post-processing. We implemented our routines in the RAMSES N-body code and tested the implementation using an example of weak lensing simulation. We analyse basic statistics of lensing convergence maps and find good agreement with semi-analytical methods. The ray tracing methodology presented here can be used in several cosmological analysis such as Sunyaev-Zeldovich and integrated Sachs-Wolfe effect studies as well as modified gravity. Our code can also be used in cross-checks of the more conventional methods, which can be important in tests of theory systematics in preparation for upcoming large scale structure surveys.
EMMA is a cosmological simulation code aimed at investigating the reionization epoch. It handles simultaneously collisionless and gas dynamics, as well as radiative transfer physics using a moment-based description with the M1 approximation. Field qu
antities are stored and computed on an adaptive 3D mesh and the spatial resolution can be dynamically modified based on physically-motivated criteria. Physical processes can be coupled at all spatial and temporal scales. We also introduce a new and optional approximation to handle radiation : the light is transported at the resolution of the non-refined grid and only once the dynamics have been fully updated, whereas thermo-chemical processes are still tracked on the refined elements. Such an approximation reduces the overheads induced by the treatment of radiation physics. A suite of standard tests are presented and passed by EMMA, providing a validation for its future use in studies of the reionization epoch. The code is parallel and is able to use graphics processing units (GPUs) to accelerate hydrodynamics and radiative transfer calculations. Depending on the optimizations and the compilers used to generate the CPU reference, global GPU acceleration factors between x3.9 and x16.9 can be obtained. Vectorization and transfer operations currently prevent better GPU performances and we expect that future optimizations and hardware evolution will lead to greater accelerations.
We present mock optical images, broad-band and H$alpha$ fluxes, and D4000 spectral indices for 30,145 galaxies from the EAGLE hydrodynamical simulation at redshift $z=0.1$, modelling dust with the SKIRT Monte Carlo radiative transfer code. The modell
ing includes a subgrid prescription for dusty star-forming regions, with both the subgrid obscuration of these regions and the fraction of metals in diffuse interstellar dust calibrated against far-infrared fluxes of local galaxies. The predicted optical colours as a function of stellar mass agree well with observation, with the SKIRT model showing marked improvement over a simple dust screen model. The orientation dependence of attenuation is weaker than observed because EAGLE galaxies are generally puffier than real galaxies, due to the pressure floor imposed on the interstellar medium. The mock H$alpha$ luminosity function agrees reasonably well with the data, and we quantify the extent to which dust obscuration affects observed H$alpha$ fluxes. The distribution of D4000 break values is bimodal, as observed. In the simulation, 20$%$ of galaxies deemed `passive for the SKIRT model, i.e. exhibiting D4000 $> 1.8$, are classified `active when ISM dust attenuation is not included. The fraction of galaxies with stellar mass greater than $10^{10}$ M$_odot$ that are deemed passive is slightly smaller than observed, which is due to low levels of residual star formation in these simulated galaxies. Colour images, fluxes and spectra of EAGLE galaxies are to be made available through the public EAGLE database.
Thermal dust emission carries information on physical conditions and dust properties in many astronomical sources. Because observations represent a sum of emission along the line of sight, their interpretation often requires radiative transfer modell
ing. We describe a new radiative transfer program SOC for computations of dust emission and examine its performance in simulations of interstellar clouds with external and internal heating. SOC implements the Monte Carlo radiative transfer method as a parallel program for shared memory computers. It can be used to study dust extinction, scattering, and emission. We tested SOC with realistic cloud models and examined the convergence and noise of the dust temperature estimates and of the resulting surface brightness maps. SOC has been demonstrated to produces accurate estimates for dust scattering and for thermal dust emission. It performs well with both with CPUs and with GPUs, the latter providing up to an order of magnitude speed-up. In the test cases, ALI improved the convergence rates but also was sensitive to Monte Carlo noise. Run-time refinement of the hierarchical-grid models did not help in reducing the run times required for a given accuracy of solution. The use of a reference field, without ALI, works more robustly. It also allows the run time to be optimised if the number of photon packages is increased only as the iterations progress. The use of GPUs in radiative transfer computations should be investigated further.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
