Journal of Applied Mathematics
Volume 2011 (2011), Article ID 126309, 26 pages
doi:10.1155/2011/126309
Research Article

Consistent Algorithms Marching Along Characteristics for the Numerical Solution of the Boltzmann Equation

1Nuclear Engineering Program, Federal University of Rio de Janeiro, P.O. Box 68509, 21941-972 Rio de Janeiro, RJ, Brazil
2Nuclear Engineering Institute, National Commission for Nuclear Energy, P.O. Box 68550, 21941-972 Rio de Janeiro, RJ, Brazil

Received 30 June 2010; Revised 19 January 2011; Accepted 20 January 2011

Academic Editor: A. Gumel

Copyright © 2011 Nilson C. Roberty and Rogerio C. Nunes. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

We introduce algorithms marching over a polygonal mesh with elements consistent with the propagation directions of the particle (radiation) flux. The decision for adopting this kind of mesh to solve the one-speed Boltzmann transport equation is due to characteristics of the domain of the transport operator which controls derivatives only in the direction of propagation of the particles (radiation) flux in the absorbing and scattering media. This a priori adaptivity has the advantages that it formulates a consistent scheme which makes appropriate the application of the Lax equivalence theorem framework to the problem. In this work, we present the main functional spaces involved in the formalism and a description of the algorithms for the mesh generation and the transport equation solution. Some numerical examples related to the solution of a transmission problem in a high-contrast model with absorption and scattering are presented. Also, a comparison with benchmarks problems for source and reactor criticality simulations shows the compatibility between calculations with the algorithms proposed here and theoretical results.