A WENO algorithm for the radiative transfer and ionized sphere at reionization

We show that the algorithm based on the weighted essentially nonoscillatory (WENO) scheme with anti-diffusive flux corrections can be used as a solver of the radiative transfer equations. This algorithm is highly stable and robust for solving problems with both discontinuities and smooth solution structures. We test this code with the ionized sphere around point sources. It shows that the WENO scheme can reveal the discontinuity of the radiative or ionizing fronts as well as the evolution of photon frequency spectrum with high accuracy on coarse meshes and for a very wide parameter space. This method would be useful to study the details of the ionized patch given by individual source in the epoch of reionization. We demonstrate this method by calculating the evolution of the ionized sphere around point sources in physical and frequency spaces. It shows that the profile of the fraction of neutral hydrogen and the ionized radius are sensitively dependent on the intensity of the source.