In this paper a new efficient algorithm for computation of radio wave ray trajectories is described. The algorithm is based on an original second-order difference scheme with a specific "length-conservation" property, which allows to resolve the ray shape even in the regions where its curvature is high. Besides the scheme, the algorithm includes a number of mechanisms ensuring its correctness and stability. The algorithm is intended for obtaining multi-pixel solar images in wide range of radio frequencies, and it is designed to be used in studies of the solar environment with modern high-resolution radio interferometers and radio telescopes such as the Murchison Widefield Array
@article{arxiv.1006.5635,
title = {Algorithm for Tracing Radio Rays in Solar Corona and Chromosphere},
author = {Leonid Benkevitch and Igor Sokolov and Divya Oberoi and Thomas Zurbuchen},
journal= {arXiv preprint arXiv:1006.5635},
year = {2010}
}
Comments
14 pages, 9 figures; author Thomas Zurbuchen added; equivalence of the ray equation to that of Walker-Haselgrove shown; derivation of the ray equation from Fermat's principle shown in detail