Summary
In Turek (Computing 54, 27–38, 1995), the concept of the generalized mean intensity has been proposed as a special numerical approach to the (linear) radiative transfer equations which can result in a significant reduction of the dimension of the discretized system, without eliminating any information for the specific intensities. Moreover, in combination with Krylov-space methods (CG, Bi-CGSTAB, etc.), robust and very efficient solvers as extensions of the classical approximate Λ-iteration have been developed. In this paper, the key tool is the combination of special renumbering techniques together with finite difference-like discretization strategies for the arising transport operators which are based on short-characteristic upwinding techniques of variable order and which can be applied to highly unstructured meshes with locally varying mesh widths. We demonstrate how such special upwinding schemes can be constructed of first order, and particularly of second order accuracy, always leading to lower triangular system matrices on general meshes. As a consequence, the global matrix assembling can be avoided (“on-the-fly”) so that the storage cost is almost optimal, and the solution of the corresponding convection-reaction subproblems for each direction can be obtained very efficiently. As a further consequence, this approach results in a direct solver in the case of no scattering, while in the case of non-vanishing absorption and scattering coefficients the resulting convergence rates for the full systems depend only on their ratio and the absolute size of these physical quantities, but only weakly on the grid size or mesh topology. We demonstrate these results via prototypical configurations and we examine the resulting accuracy and efficiency for different computational domains, meshes and problem parameters.
Similar content being viewed by others
References
Deo, N.: Graph theory with applications to engineering and computer science. Prentice Hall, Englewood Cliffs (1974), pp. 222–333
Führer, C.: A comparative study of finite element solvers for hyperbolic problems with applications to radiative transfer. Technical report SFB 359, 65, Universität Heidelberg (1993)
Hübner, T.: Spezielle Diskretisierungs- und Lösungsmethoden für Integro-Differentialgleichungen am Beispiel der Strahlungstransportgleichung. Diploma thesis, Dortmund (2005). www. mathematik.uni-dortmund.de/lsiii/static/showpdffile_Huebner2005.pdf
Johnson C. and Pitkäranta J. (1983). Convergence of a fully discrete scheme for two-dimensional neutron transport. SIAM J Numer Anal 20: 951–966
Kalkofen W. (1987). Numerical radiative transfer. Cambridge University Press, London
Klar, A., Seaid, M.: Efficient preconditioning of linear systems arising from the discretization of radiative transfer equation. In: Challenges in Scientific Computing, Berlin 2002. Lecture Notes in Computational Science and Engineering (2003)
Larsen E. W. (1982). Unconditionally stable diffusion-synthetic acceleration methods for the slab geometry discrete ordinates equation. Nuclear Sci Engng 82: 47
Larsen E.W. and Morel J.E. (1989). Asymptotic solutions of numerical transport problems in optically thick diffusive regimes II. J Comput Phys 83: 212
Meinköhn, E.: Modellierung dreidimensionaler Strahlungsfelder im frühen Universum. Ph.D thesis, Heidelberg (2002)
Mihalas D. and Weibel–Mihalas B. (1984). Foundations of radiation hydrodynamics. Oxford University Press, New York
Richling S., Meinköhn E., Kryzhevoi N. and Kanschat G. (2001). Radiative transfer with finite elements. Astronom Astrophs 380: 776–788
Schäfer, M., Turek, S. (with support by Durst, F., Krause, E., Rannacher, R.): Benchmark computations of laminar flow around cylinder. In: Hirschel, E. H. (ed.) Flow simulation with high-performance computers II. Notes on Numerical Fluid Mechanics 52, 547–566. Vieweg, Braunschweig (1996)
Seaïd M., Frank M., Klar A., Pinnau R. and Thömmes G. (2004). Efficient numerical methods for radiation in gas turbines. J Comp Appl Math 170(1): 217–239
Steiner, O.: A rapidly converging temperature correction procedure using operator perturbation. Astronom Astrophys 231, 278–288. Springer, Heidelberg (1990)
Turek S. (1993). An efficient solution technique for the radiative transfer equation. Impact Comput Sci Engng 5: 201–214
Turek S. (1995). A generalized mean intensity approach for the numerical solution of the radiative transfer equation. Computing 54: 27–38
Turek S., Runge A. and Becker Ch. (2001). The FEAST Indices – realistic evaluation of modern software components and processor technologies. Comput Math Appl 41: 1431–1464
Van der Vorst H. (1992). BI–CGSTAB: a fast and smoothly converging variant of BI–CG for the solution of non-symmetric linear systems. SIAM J Sci Stat Comput 13: 631–644
Väth, H. M.: Three-dimensional radiative transfer on a massively parallel computer. Astronom Astrophys 284, 319. Springer, Heidelberg (1994)
Author information
Authors and Affiliations
Corresponding author
Additional information
This paper is based on research which was supported in part by DFG grant No. TU102/16-1.
Rights and permissions
About this article
Cite this article
Hübner, T., Turek, S. An efficient and accurate short-characteristics solver for radiative transfer problems. Computing 81, 281–296 (2007). https://doi.org/10.1007/s00607-007-0250-2
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00607-007-0250-2