Pakal: A Three-dimensional Model to Solve the Radiative Transfer Equation

• Published in: The Astrophysical journal. Supplement series Vol. 188; no. 2; pp. 437 - 446
• Main Authors: De la Luz, Victor, Lara, Alejandro, Mendoza-Torres, J. E, Selhorst, Caius L
• Format: Journal Article
• Language: English
• Published: United States IOP Publishing 01.06.2010
• Subjects: ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; ATMOSPHERES; CHROMOSPHERE; CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ELEMENTS; FLUID MECHANICS; FLUIDS; GASES; HELIUM; HYDRODYNAMICS; HYDROGEN; MAGNETOHYDRODYNAMICS; MAIN SEQUENCE STARS; MECHANICS; NONMETALS; RADIATIONS; RARE GASES; SOLAR ATMOSPHERE; SOLAR RADIATION; STARS; STELLAR ATMOSPHERES; STELLAR RADIATION; SUN; THREE-DIMENSIONAL CALCULATIONS
• ISSN: 0067-0049; 1538-4365
• DOI: 10.1088/0067-0049/188/2/437

We present a new numerical model called 'Pakal' intended to solve the radiative transfer equation in a three-dimensional (3D) geometry, using the approximation for a locally plane-parallel atmosphere. Pakal uses pre-calculated radial profiles of density and temperature (based on hydrostatic, hydrodynamic, or MHD models) to compute the emission from 3D source structures with high spatial resolution. Then, Pakal solves the radiative transfer equation in a set of (3D) ray paths, going from the source to the observer. Pakal uses a new algorithm to compute the radiative transfer equation by using an intelligent system consisting of three structures: a cellular automaton; an expert system; and a program coordinator. The code outputs can be either two-dimensional maps or one-dimensional profiles, which reproduce the observations with high accuracy, giving detailed physical information about the environment where the radiation was generated and/or transmitted. We present the model applied to a 3D solar radial geometry, assuming a locally plane-parallel atmosphere, and thermal free-free radio emission from hydrogen-helium gas in thermodynamic equilibrium. We also present the convergence test of the code. We computed the synthetic spectrum of the centimetric-millimetric solar emission and found better agreement with observations (up to 10{sup 4} K at 20 GHz) than previous models reported in the literature. The stability and convergence test show the high accuracy of the code. Finally, Pakal can improve the integration time by up to an order of magnitude compared against linear integration codes.