TSU3

A Comparison of Radiative Transfer Codes


Numerous physical problems require a detailed understanding of the radiative transfer of photons into different environments, ranging from intergalactic and interstellar medium to stellar or planetary atmospheres. The solution of the complete radiative transfer equation is presently beyond our capabilities, but several groups are  attacking the problem by using different numerical schemes and approximations (i.e. ray-tracing,  Monte Carlo methods, local depth approximation...). At this stage, it appears worthwhile to make a comparison between these different approaches.

At the conference  The Physics of Galaxy Formation, Tsukuba, Japan, July 3-7 2000, Andrea Ferrara, in his concluding remarks, proposed to apply the available codes to a commonly defined radiative transfer problem . In honor of Tsukuba, guest city of the conference, this initiative has been called the TsuCube. During the workshop that will be hosted by CITA we plan to finalize the comparison between different codes.

Below you can find the necessary instruction for the proposed tests, while in the paper you can find more details.
Note that the codes that do not solve He chemistry should perform the Tests where He is included (Test 0 and Test 2) setting the initial H density to the total density (e.g. in Test 0 nH=1 cm-3).

The basic physics

- Test 0: rates and evolution of a single zone

Pure radiative transfer tests

- Test 1: pure-hydrogen isothermal HII region expansion

- Test 2: HII region expansion: the temperature state

- Test 3: I-front trapping in a dense clump and the formation of a shadow

- Test 4: Multiple sources in a cosmological density field

Radiative hydrodynamics tests

- Test 5: Classical HII region expansion

- Test 6: HII region expansion in 1/r2 density profile

- Test 7: Photoevaporation of a dense clump


Submission of tests results

To upload your files (we recommend to COMPRESS them) please go to:
www.cita.utoronto.ca/~iliev/workshop/form.html
which provides a form for file submission. Since all files you submit would go to the same directory, please follow the following naming conventions. The files should be called
Tnumber_grid_name.dat
and
Tnumber_front_name.dat
for the 3-D and the 1-D data respectively. Here, "name" is the name of your group or your code and "number" is the test number (e.g. T2_grid_iliev.dat).
Please note that the 1-D I-front positions and velocities for each test should be submitted as one file, in 3 columns (t,x_I,v_I), with units [Myr,kpc,km/s]. The x_I should be the linearly-interpolated I_front position inside the cell and v_I is finite-differenced from x_I.
For what concerns the format, all files should be provided as ASCII free-format single-precision (Real*4) files, i.e. the 3-D data should be readable using:
do i=1,128
  do j=1,128
   do k=1,128
    read(*,*,*) all requested quantities in the order listed
   end do
  end do
end do


Participants

The codes/groups that have so far joined the TsuCube are:
1.   CRASH: Cosmological RAdiative transfer Scheme for Hydrodynamics (A. Maselli, A. Ferrara & B. Ciardi)
2.   C2-Ray: Photon-conserving transport of ionizing radiation (G. Mellema, I. Iliev, P. Shapiro & M. Alvarez)
3.   OTVET: Optically Thin Variable Eddington Tensor (N. Gnedin & T. Abel)
4.   ART: Authentic Radiative Transfer with Discretized Long Beams (T. Nakamoto, H. Susa, K. Hiroi & M. Umemura)
5.   RSPH: SPH coupled with radiative transfer (H. Susa & M. Umemura)
6.   FLASH-HC: Hybrid Characteristics: 3D radiative transfer for parallelized AMR hydrodynamics codes (E.-J. Rijkhorst, T. Plewa, A. Dubey & G. Mellema)
7.   SimpleX: radiative transfer on unconstructed grids (J. Ritzerveld, V. Icke & E.-J. Rijkhorst)
8.   ZEUS-MP with radiative transfer (D. Whalen & M. Norman)
9.   FTTE: Fully Threaded Transport Engine (A. Razoumov)
10.   IFT: Ionization Front Tracking (M. Alvarez & P. Shapiro)
11.   CORAL (I. Iliev, A. Raga, G. Mellema, P. Shapiro)
 


This page is maintained by Benedetta Ciardi. Last update April 4 2006.