Monte Carlo Radiative Transfer Simulations
The simulations were performed in a 1.0^3 box, speed of light = 1.0.
The opacity 32^3 data was taken from Greg Bryan. Linear interpolation
technique was used to get the opacity at any point in the box.
The radiation field of
a point source in the box was simulated using Monte Carlo Method.
The time step was dt=0.005. The Point source was turned on at the bottom of the box
at the time T0=0.
The plates were exposed at the time tb=5.00, when the field
approached the equilibrium, and were accumulating photons
up to the time T1.
The source was emitting 2,000,000 photons per unit time.
High Resolution Pictures on 400 X 400 pixel plates:
- T1=128.505 Grey Scale images. At the time T1 the number of advancing photons in the simulation Box was 1,509,230;
- Ceiling plate[LIN][LOG]
27,137,107 total number of photons accumulated.
Statistics per grain: the darkest grain accumulated 8 photons,
the brightest accumulated 288 photons;
- Side plate 69,823,750 scattering events recorded, (LOG);
- T1=559.5 (altair time=11.9d) Grey Scale images. At the time T1 the number of advancing photons in the simulation Box was 1,509,916;
- Ceiling plate[LIN]
121,857,036 total number of photons accumulated.
Statistics per grain: the darkest grain accumulated 68 photons,
the brightest accumulated 1170 photons;
- Side plate 311,147,516 scattering events recorded, (LOG);
[back]