Introduction of Parallel GPGPU Acceleration Algorithms for the Solution of Radiative Transfer

General-purpose computing on graphics processing units (GPGPU) is a recent technique that allows the parallel graphics processing unit (GPU) to accelerate calculations performed sequentially by the central processing unit (CPU). To introduce GPGPU to radiative transfer, the Gauss-Seidel solution of the well-known expressions for 1-D and 3-D homogeneous, isotropic media is selected as a test case. Different algorithms are introduced to balance memory and GPU-CPU communication, critical aspects of GPGPU. Results show that speed-ups of one to two orders of magnitude are obtained when compared to sequential solutions. The underlying value of GPGPU is its potential extension in radiative solvers (e.g., Monte Carlo, discrete ordinates) at a minimal learning curve.

[1]  Önder Yιldιz,et al.  A Parallel Solution to the Radiative Transport in Three-Dimensional Participating Media , 2006 .

[2]  R. D. Skocypec,et al.  The use of high-performance computing to solve participating media radiative heat transfer problems-results of an NSF workshop , 1995 .

[3]  William F. Godoy,et al.  On the use of flux limiters in the discrete ordinates method for 3D radiation calculations in absorbing and scattering media , 2010, J. Comput. Phys..

[4]  M. E. Larsen.  The Exchange Factor Method: AN Alternative Zonal Formulation for Analysis of Radiating Enclosures Containing Participating Media. , 1983 .

[5]  Pedro J. Coelho,et al.  PARALLELIZATION OF THE DISCRETE ORDINATES METHOD , 1997 .

[6]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[7]  Bjarne Stroustrup,et al.  C++ Programming Language , 1986, IEEE Softw..

[8]  James Demmel,et al.  IEEE Standard for Floating-Point Arithmetic , 2008 .

[9]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[10]  R. Viskanta.  Radiation Transfer and Interaction of Convection with Radiation Heat Transfer , 1966 .

[11]  Abhilash J. Chandy,et al.  Parallelizing the Discrete Ordinates Method (DOM) for Three-Dimensional Radiative Heat Transfer Calculations using a Priority Queuing Technique , 2007 .

[12]  W. A. Fiveland,et al.  Three-dimensional radiative heat-transfer solutions by the discrete-ordinates method , 1988 .

[13]  John D. Owens,et al.  General Purpose Computation on Graphics Hardware , 2005, IEEE Visualization.

[14]  K. Liou,et al.  A monte carlo method for 3D thermal infrared radiative transfer , 2006 .

[15]  Robert Pincus,et al.  Computational Cost and Accuracy in Calculating Three-Dimensional Radiative Transfer: Results for New Implementations of Monte Carlo and SHDOM , 2009 .

[16]  U. R. Hanebutte,et al.  A Massively Parallel Discrete Ordinates Response Matrix Method for Neutron Transport , 1992 .

[17]  M. Modest.  Radiative heat transfer , 1993 .

[18]  Alireza Haghighat.  Angular parallelization of a curvilinear S sub n transport theory method , 1991 .

[19]  Uresh K. Vahalia.  UNIX Internals: The New Frontiers , 1995 .

[20]  R. F. Warming,et al.  Radiative transport and wall temperature slip in an absorbing planar medium , 1965 .

[21]  John R. Howell,et al.  Comparison of Monte Carlo Strategies for Radiative Transfer in Participating Media , 1998 .

[22]  Jens H. Krüger,et al.  A Survey of General‐Purpose Computation on Graphics Hardware , 2007, Eurographics.

[23]  T. S. West.  New Frontiers , 1968, Nature.

[24]  M. Modest.  CHAPTER 16 – THE METHOD OF DISCRETE ORDINATES (SN-APPROXIMATION) , 2003 .

[25]  A. Crosbie,et al.  Exact expressions for radiative transfer in a three-dimensional rectangular geometry☆ , 1982 .

[26]  Gautham Krishnamoorthy,et al.  PARALLEL COMPUTATIONS OF RADIATIVE HEAT TRANSFER USING THE DISCRETE ORDINATES METHOD , 2004 .

[27]  Bjarne Stroustrup,et al.  The C++ programming language (2nd ed.) , 1991 .