skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Specular reflection treatment for the 3D radiative transfer equation solved with the discrete ordinates method

Abstract

The contribution of this paper relies in the development of numerical algorithms for the mathematical treatment of specular reflection on borders when dealing with the numerical solution of radiative transfer problems. The radiative transfer equation being integro-differential, the discrete ordinates method allows to write down a set of semi-discrete equations in which weights are to be calculated. The calculation of these weights is well known to be based on either a quadrature or on angular discretization, making the use of such method straightforward for the state equation. Also, the diffuse contribution of reflection on borders is usually well taken into account. However, the calculation of accurate partition ratio coefficients is much more tricky for the specular condition applied on arbitrary geometrical borders. This paper presents algorithms that calculate analytically partition ratio coefficients needed in numerical treatments. The developed algorithms, combined with a decentered finite element scheme, are validated with the help of comparisons with analytical solutions before being applied on complex geometries.

Authors:
 [1];  [1];  [1];  [2]
  1. Université de Nantes, LTN UMR CNRS 6607 (France)
  2. Sorbonne Universités, UPMC Université Paris 06, UMR 7598, inria de Paris, Laboratoire Jacques-Louis Lions, F-75005, Paris (France)
Publication Date:
OSTI Identifier:
22622276
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 334; Other Information: Copyright (c) 2017 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICAL METHODS AND COMPUTING; ALGORITHMS; ANALYTICAL SOLUTION; COMPARATIVE EVALUATIONS; COMPUTER CALCULATIONS; DISCRETE ORDINATE METHOD; EQUATIONS; FINITE ELEMENT METHOD; PARTITION; QUADRATURES; RADIANT HEAT TRANSFER; REFLECTION

Citation Formats

Le Hardy, D., Favennec, Y., E-mail: yann.favennec@univ-nantes.fr, Rousseau, B., and Hecht, F.. Specular reflection treatment for the 3D radiative transfer equation solved with the discrete ordinates method. United States: N. p., 2017. Web. doi:10.1016/J.JCP.2017.01.019.
Le Hardy, D., Favennec, Y., E-mail: yann.favennec@univ-nantes.fr, Rousseau, B., & Hecht, F.. Specular reflection treatment for the 3D radiative transfer equation solved with the discrete ordinates method. United States. doi:10.1016/J.JCP.2017.01.019.
Le Hardy, D., Favennec, Y., E-mail: yann.favennec@univ-nantes.fr, Rousseau, B., and Hecht, F.. Sat . "Specular reflection treatment for the 3D radiative transfer equation solved with the discrete ordinates method". United States. doi:10.1016/J.JCP.2017.01.019.
@article{osti_22622276,
title = {Specular reflection treatment for the 3D radiative transfer equation solved with the discrete ordinates method},
author = {Le Hardy, D. and Favennec, Y., E-mail: yann.favennec@univ-nantes.fr and Rousseau, B. and Hecht, F.},
abstractNote = {The contribution of this paper relies in the development of numerical algorithms for the mathematical treatment of specular reflection on borders when dealing with the numerical solution of radiative transfer problems. The radiative transfer equation being integro-differential, the discrete ordinates method allows to write down a set of semi-discrete equations in which weights are to be calculated. The calculation of these weights is well known to be based on either a quadrature or on angular discretization, making the use of such method straightforward for the state equation. Also, the diffuse contribution of reflection on borders is usually well taken into account. However, the calculation of accurate partition ratio coefficients is much more tricky for the specular condition applied on arbitrary geometrical borders. This paper presents algorithms that calculate analytically partition ratio coefficients needed in numerical treatments. The developed algorithms, combined with a decentered finite element scheme, are validated with the help of comparisons with analytical solutions before being applied on complex geometries.},
doi = {10.1016/J.JCP.2017.01.019},
journal = {Journal of Computational Physics},
number = ,
volume = 334,
place = {United States},
year = {Sat Apr 01 00:00:00 EDT 2017},
month = {Sat Apr 01 00:00:00 EDT 2017}
}
  • The Schur-decomposition for three-dimensional matrix equations is developed and used to directly solve the radiative discrete ordinates equations which are discretized by Chebyshev collocation spectral method. Three methods, say, the spectral methods based on 2D and 3D matrix equation solvers individually, and the standard discrete ordinates method, are presented. The numerical results show the good accuracy of spectral method based on direct solvers. The CPU time cost comparisons against the resolutions between these three methods are made using MATLAB and FORTRAN 95 computer languages separately. The results show that the CPU time cost of Chebyshev collocation spectral method with 3Dmore » Schur-decomposition solver is the least, and almost only one thirtieth to one fiftieth CPU time is needed when using the spectral method with 3D Schur-decomposition solver compared with the standard discrete ordinates method.« less
  • The S4 discrete-ordinates approximation is used to solve the radiative transfer equation in nonasymmetric (i.e., three-dimensional) cylindrical enclosures containing absorbing-emitting and scattering media, with and without the temperature profile known a priori. Because neither detailed experimental data nor predictions from a zone or Monte-Carlo model for three-dimensional cylindrical enclosures are available, cylindrical equivalents of three-dimensional rectangular enclosures, for which zone model predictions of radiative transfer are available, are used in model evaluation. Limited evaluation of the model shows that the discrete-ordinates method provides acceptable predictions of radiative transfer in nonaxisymmetric cylindrical enclosures. 12 refs.
  • A radiation code based on the method of lines (MOL) solution of the discrete ordinates method (DOM) for transient three-dimensional radiative heat transfer in rectangular enclosures for use in conjunction with a computational fluid dynamics (CFD) code based on the same approach was developed. Assessment of the predictive accuracy of the code by benchmarking its steady-state solutions against exact solutions on one- and three-dimensional test problems shows that the MOL solution of the DOM provides accurate and computationally efficient solutions for radiative heat fluxes and source terms and can be used with confidence in conjunction with CFD codes for transientmore » problems.« less
  • Advantages and disadvantages of modern discrete-ordinates finite-element methods for the solution of radiative transfer problems in meteorology, climatology, and remote sensing applications are evaluated. After the common basis of the formulation of radiative transfer problems in the fields of neutron transport and atmospheric optics is established, the essential features of the discrete-ordinates finite-element method are described including the limitations of the method and their remedies. Numerical results are presented for 1-D and 2-D atmospheric radiative transfer problems where integral as well as angular dependent quantities are compared with published results from other calculations and with measured data. These comparisons providemore » a verification of the discrete-ordinates results for a wide spectrum of cases with varying degrees of absorption, scattering, and anisotropic phase functions. Accuracy and computational speed are also discussed. Since practically all discrete-ordinates codes offer a builtin adjoint capability, the general concept of the adjoint method is described and illustrated by sample problems. Our general conclusion is that the strengths of the discrete-ordinates finite-element method outweight its weaknesses. We demonstrate that existing general-purpose discrete-ordinates codes can provide a powerful tool to analyze radiative transfer problems through the atmosphere, especially when 2-D geometries must be considered.« less
  • This paper introduces a powerful but simple methodology for solving the general equation of radiative transfer for scattering and/or absorbing one-dimensional systems. Existing methods, usually designed to handle specific boundary and energy equilibrium conditions, either provide crude estimates or involve intricate mathematical analysis coupled with numerical techniques. In contrast, the present scheme, which uses a discrete-ordinate technique to reduce the integro-differential equation to a system of ordinary differential equations, utilizes readily available software routines to solve the resulting set of coupled first-order ordinary differential equations as a two-point boundary value problem. The advantage of this approach is that the usermore » is freed from having to understand complicated mathematical analysis and perform extensive computer programming. Additionally, the software used is state of the art, which is less prone to numerical instabilities and inaccuracies. Any degree of scattering anisotropy and albedo can be incorporated along with different conditions of energy equilibrium or specified temperature distributions and boundary conditions. Examples are presented where the radiative transfer is computed by using different quadratures such as Gaussian, Lobatto, Fiveland, Chebyshev, and Newton-Cotes. Comparison with benchmark cases shows that in a highly forward scattering medium Gaussian quadrature provides the most accurate and stable solutions.« less