Alternative scattering kernels for the first estimates of a reactor: diffusion length

Alternative scattering kernels for the first estimates of a reactor: diffusion length 1 corresp first-author Corresponding author E-mail: hakanozturk@osmaniye.edu.tr Corresponding author E-mail: hakanozturk@osmaniye.edu.tr hakanozturk@osmaniye.edu.tr 1 Osmaniye Korkut Ata University, Faculty of Arts and Sciences, Department of Physics, 80000 Osmaniye, Turkey. Osmaniye Korkut Ata University, Faculty of Arts and Sciences, Department of Physics, 80000 Osmaniye, Turkey. Diffusion length calculations of neutrons are performed using the Chebyshev polynomials of the second kind. The neutrons are assumed to move with constant energy in a uniform homogeneous slab. An alternative scattering kernel called an Anlı-Güngör (AG) phase function and a traditional Henyey-Greenstein (HG) phase function are used for the scattering function in the stationary neutron transport equation. First, analytic expressions and then numerical results are obtained for the diffusion length for various values of the scattering and cross-section parameters. Numerical results obtained from both scattering kernels for the diffusion length of the neutrons are given in tables side by side for comparison. The applicability of the method is easily demonstrated by these results. Scattering kernel Diffusion length Chebyshev polynomials Transport equation 1 Introduction 1 First estimates of a reactor are very important, to establish and operate it safely, as these studies provide details for the subsequent developments for the design of a reactor. In the solution algorithm of the approximate methods, higher order approximations are generally preferred to obtain a closer description of the real system. In these methods, the angular part of the neutron flux is expanded in terms of orthogonal polynomials. The arguments of these polynomials supply the definition interval of the cosine of the travelling direction of the neutron with respect to the x axis. However, it is still worth considering the first order approximation since it is valid for the first estimates of the system. The first order approximation of these polynomial expansion–based methods is called a diffusion approximation. It provides convenient results regarding the transport and energy spectrum of the neutrons if the number of secondary neutrons per collision, c , is close to unity [ 1 1 ]. As the Legendre polynomials (spherical harmonics) approximation is not the unique and valid method for all problems, the Chebyshev polynomials of both kinds can be used for the solutions in neutron transport theory because these polynomials have the same definition interval and are in the same polynomial family, i.e., Jacobi polynomials. In some recent studies, Chebyshev polynomials of the second kind, i.e., U N method, have been successfully applied to the one-dimensional transport equation in slab and spherical geometries for the solution of the diffusion length and criticality problems [ 2 2 - 5 5 2-5 ]. It is generally difficult to improve a solution algorithm for the solution of the transport equation because of its integro-differential form. Therefore, the integrand term, including the scattering function, is usually converted to a sum using an integral transformation. By doing so, this difficulty is overcome and the analytic, semi analytic, or numerical solution of the problem in transport theory is found. In those methods, the scattering function is expanded in terms of the orthogonal polynomials [ 6 6 , 7 7 ]. However, the structure of the equation allows only a few terms representing the degrees of anisotropy. Therefore, these approximate functions have serious restrictions on the exact solution of the transport equation. A scattering (phase) function, such as Henyey-Greenstein (HG), is preferred for the solution of problems in the field of particle and photon transport, radiative transfer, and light scattering [ 8 8 - 14 14 8-14 ]. An alternative scattering function derived by Anlı et al. is used for specifying the eigenvalue spectrum of the neutrons in slab geometry [ 15 15 ]. This scattering function has been successfully applied to the transport equation in a few studies, using Legendre and Chebyshev polynomials. Results consistent with the literature were obtained [ 16 16 , 17 17 ]. However, this present work is the first to consider the Chebyshev polynomials of the second kind and the HG and AG phase functions were used together for the solution of the transport equation. In this work, we apply the U N method in the calculation of the diffusion lengths for one-speed neutrons in a uniform homogeneous slab. The alternative scattering function previously derived and applied to transport equation [ 16 16 , 17 17 ] together with the well-known Henyey-Greenstein phase function is used. Therefore, by examining the derivations of the equations and analyzing the results given in the tables, the method can be extended to other problems in applied science and engineering. 2 U 1 approximation with Henyey-Greenstein (HG) Phase Function 1 The time-independent linear transport equation in one-dimensional geometry for one-speed neutrons without a source can be written as, <math display="block" id="M1"><mrow><mi>μ</mi><mfrac><mrow><mo>∂</mo><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac><mo>+</mo><msub><mi>σ</mi><mtext>T</mtext></msub><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><mrow><mn>2</mn><mi>π</mi></mrow></munderover><mrow><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mrow><mo>−</mo><mn>1</mn></mrow><mn>1</mn></munderover><mrow><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><msup><mi>μ</mi><mo>′</mo></msup><mo stretchy="false">)</mo><mtext> </mtext><msub><mi>σ</mi><mtext>S</mtext></msub><mo stretchy="false">(</mo><msub><mi>μ</mi><mn>0</mn></msub><mo stretchy="false">)</mo><mtext> </mtext></mrow></mrow></mstyle><mtext>d</mtext><msup><mi>μ</mi><mo>′</mo></msup><mtext> </mtext></mrow></mrow></mstyle><mtext>d</mtext><msup><mi>φ</mi><mo>′</mo></msup><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>−</mo><mi>a</mi><mo>≤</mo><mi>x</mi><mo>≤</mo><mi>a</mi><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>−</mo><mn>1</mn><mo>≤</mo><mi>μ</mi><mo>≤</mo><mn>1</mn><mo>,</mo></mrow></math> where <math id="M2"><mrow><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo></mrow></math> is the angular flux or flux density of neutrons at position x traveling in direction µ , μ μ′ is the direction cosine of the neutron after a collision, and σ T is the macroscopic total cross section [ 18 18 ]. µ 0 is the cosine of the scattering angle between the the directions of neutron velocity before <math id="M3"><mstyle mathvariant="bold" mathsize="normal"><mi>Ω</mi></mstyle></math> and after <math id="M4"><mstyle mathvariant="bold" mathsize="normal"><msup><mi>Ω</mi><mo>′</mo></msup></mstyle></math> a collision, and is obtained by the scalar product <math id="M5"><mrow><msub><mi>μ</mi><mn>0</mn></msub><mo>=</mo><mi>Ω</mi><mo>⋅</mo><msup><mi>Ω</mi><mo>′</mo></msup></mrow></math> ; <math display="block" id="M6"><mrow><msub><mi>μ</mi><mn>0</mn></msub><mo>=</mo><mi>μ</mi><msup><mi>μ</mi><mo>′</mo></msup><mo>+</mo><msqrt><mrow><mn>1</mn><mo>−</mo><msup><mi>μ</mi><mn>2</mn></msup></mrow></msqrt><msqrt><mrow><mn>1</mn><mo>−</mo><msup><msup><mi>μ</mi><mo>′</mo></msup><mn>2</mn></msup></mrow></msqrt><mi>cos</mi><mo stretchy="false">(</mo><mi>ϕ</mi><mo>−</mo><msup><mi>ϕ</mi><mo>′</mo></msup><mo stretchy="false">)</mo><mo>,</mo></mrow></math> and <math id="M7"><mrow><msub><mi>σ</mi><mtext>S</mtext></msub><mo stretchy="false">(</mo><msub><mi>μ</mi><mn>0</mn></msub><mo stretchy="false">)</mo></mrow></math> is the scattering kernel which describes all neutron interactions with the materials inside the system. One of the most popular scattering functions used in the solutions of the problems regarding radiative transfer is the Henyey-Greenstein (HG) phase function. In the first part of this study, the HG phase function is used for the scattering function in the transport equation; <math display="block" id="M8"><mrow><msub><mi>σ</mi><mi>S</mi></msub><msup><mrow/><mrow><mi>H</mi><mi>G</mi></mrow></msup><mo stretchy="false">(</mo><msub><mi>μ</mi><mn>0</mn></msub><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mi>σ</mi><mi>S</mi></msub><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>t</mi><mn>2</mn></msup><mo stretchy="false">)</mo></mrow><mrow><mn>4</mn><mi>π</mi><mtext> </mtext><msup><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><mn>2</mn><msub><mi>μ</mi><mn>0</mn></msub><mi>t</mi><mo>+</mo><msup><mi>t</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow></mfrac><mo>,</mo></mrow></math> where σ S is any non-negative coefficient and t represents all kinds of scattering, such as, isotropic, forward, backward, and anisotropic, and is in the range <math id="M9"><mrow><mo>−</mo><mn>1</mn><mo>≤</mo><mi>t</mi><mo>≤</mo><mn>1</mn></mrow></math> [ 19 19 ]. The integral of the HG phase function is needed on the right hand side of Eq. (1) and can be obtained with the addition theorem of the Legendre polynomials [ 15 15 ], <math display="block" id="M10"><mrow><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><mrow><mn>2</mn><mi>π</mi></mrow></munderover><mrow><msubsup><mi>σ</mi><mi>S</mi><mrow><mi>H</mi><mi>G</mi></mrow></msubsup><mo stretchy="false">(</mo><msub><mi>μ</mi><mn>0</mn></msub><mo stretchy="false">)</mo><mtext> </mtext><mtext>d</mtext><msup><mi>ϕ</mi><mo>′</mo></msup><mo>=</mo><mfrac><mrow><msub><mi>σ</mi><mi>S</mi></msub></mrow><mn>2</mn></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mi>∞</mi></munderover><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo><mtext> </mtext><msup><mi>t</mi><mi>n</mi></msup><msub><mi>P</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>μ</mi><mo stretchy="false">)</mo><msub><mi>P</mi><mi>n</mi></msub><mo stretchy="false">(</mo><msup><mi>μ</mi><mo>′</mo></msup><mo stretchy="false">)</mo></mrow></mstyle></mrow></mrow></mstyle><mo>.</mo></mrow></math> Using Eq. (4) in Eq. (1) with a definition of a dimensionless constant, σ T x / ν → x , we obtain; <math display="block" id="M11"><mrow><mi>μ</mi><mfrac><mrow><mo>∂</mo><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac><mo>+</mo><mi>v</mi><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mi>v</mi><mi>c</mi></mrow><mn>2</mn></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo><msup><mi>t</mi><mi>n</mi></msup><msub><mi>P</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>μ</mi><mo stretchy="false">)</mo><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mrow><mo>−</mo><mn>1</mn></mrow><mrow><mo>+</mo><mn>1</mn></mrow></munderover><mrow><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo>'</mo><mo stretchy="false">)</mo><msub><mi>P</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>μ</mi><mo>'</mo><mo stretchy="false">)</mo><mtext> </mtext><mtext>d</mtext><mi>μ</mi><mo>'</mo></mrow></mrow></mstyle></mrow></mstyle><mo>,</mo></mrow></math> where ν denotes the eigenvalues and <math id="M12"><mrow><mi>c</mi><mo>=</mo><msub><mi>σ</mi><mtext>S</mtext></msub><mo>/</mo><msub><mi>σ</mi><mtext>T</mtext></msub></mrow></math> . In the solution algorithm of Eq. (5), the neutron angular flux is expanded in terms of the Chebyshev polynomials of the second kind as [ 3 3 , 4 4 ], <math display="block" id="M13"><mrow><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>2</mn><mi>π</mi></mfrac><msqrt><mrow><mn>1</mn><mo>−</mo><msup><mi>μ</mi><mn>2</mn></msup></mrow></msqrt><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><mrow><msub><mi>Φ</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow></mstyle><mtext> </mtext><msub><mi>U</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>−</mo><mi>a</mi><mo>≤</mo><mi>x</mi><mo>≤</mo><mi>a</mi><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>−</mo><mn>1</mn><mo>≤</mo><mi>μ</mi><mo>≤</mo><mn>1.</mn></mrow></math> The orthogonality and recurrence relations of the Chebyshev polynomials of the second kind are needed for obtaining the U N moments of the angular flux [ 20 20 ], <math display="block" id="M14"><mrow><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mrow><mo>−</mo><mn>1</mn></mrow><mn>1</mn></munderover><mrow><msub><mi>U</mi><mi>n</mi></msub></mrow></mrow></mstyle><mo stretchy="false">(</mo><mi>μ</mi><mo stretchy="false">)</mo><msub><mi>U</mi><mi>m</mi></msub><mo stretchy="false">(</mo><mi>μ</mi><mo stretchy="false">)</mo><msqrt><mrow><mn>1</mn><mo>−</mo><msup><mi>μ</mi><mn>2</mn></msup></mrow></msqrt><mtext>d</mtext><mi>μ</mi><mo>=</mo><mrow><mo>{</mo><mtable columnalign="left"><mtr><mtd><mi>π</mi><mo>/</mo><mn>2</mn><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mi>n</mi><mo>=</mo><mi>m</mi></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mi>n</mi><mo>≠</mo><mi>m</mi></mtd></mtr></mtable></mrow></mrow></math> <math display="block" id="M15"><mrow><mn>2</mn><mi>μ</mi><msub><mi>U</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>U</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo stretchy="false">(</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>U</mi><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow></msub><mo stretchy="false">(</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>−</mo><mn>1</mn><mo>≤</mo><mi>μ</mi><mo>≤</mo><mn>1.</mn></mrow></math> In order to obtain the U N moments of the angular flux, Eq. (6) is substituted in Eq. (5). After multiplying by U m ( μ ) and integrating over μ in the interval <math id="M16"><mrow><mo>−</mo><mn>1</mn><mo>≤</mo><mi>μ</mi><mo>≤</mo><mn>1</mn></mrow></math> , <math display="block" id="M17"><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>Φ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><mi>x</mi></mrow></mfrac><mo>+</mo><mn>2</mn><mi>v</mi><mtext> </mtext><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mn>2</mn><mi>v</mi><mi>c</mi><mtext> </mtext><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>,</mo></mrow></math> <math display="block" id="M18"><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>Φ</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><mi>x</mi></mrow></mfrac><mo>+</mo><mfrac><mrow><mtext>d</mtext><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><mi>x</mi></mrow></mfrac><mo>+</mo><mn>2</mn><mi>v</mi><mtext> </mtext><msub><mi>Φ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mn>2</mn><mi>v</mi><mi>c</mi><mi>t</mi><mtext> </mtext><msub><mi>Φ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>.</mo></mrow></math> for n = 0 and n = 1, respectively. Equation (10), corresponding to of n = 1, is generally referred to as the diffusion approximation and the last contribution of the flux is negligible by setting <math id="M19"><mrow><mrow><mrow><mtext>d</mtext><msub><mi>Φ</mi><mrow><mi>N</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><mo>/</mo><mrow><mtext>d</mtext><mi>x</mi></mrow></mrow><mo>=</mo><mn>0</mn></mrow></math> , as for the spherical harmonics ( P N ) method [ 18 18 ]. Then for the U 1 approximation, the coupled Eqs. (9) and (10) can be solved for Φ 0 ( x ), <math display="block" id="M20"><mrow><mfrac><mrow><msup><mtext>d</mtext><mtext>2</mtext></msup><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><msup><mi>x</mi><mn>2</mn></msup></mrow></mfrac><mo>−</mo><mn>4</mn><msup><mi>v</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>c</mi><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>c</mi><mi>t</mi><mo stretchy="false">)</mo><mtext> </mtext><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mn>0.</mn></mrow></math> The diffusion length is defined as the inverse square root of the coefficient of Φ 0 ( x ) in Eq. (11), <math display="block" id="M21"><mrow><msup><mi>L</mi><mrow><mtext>HG</mtext></mrow></msup><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>v</mi><msqrt><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>c</mi><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>c</mi><mi>t</mi><mo stretchy="false">)</mo></mrow></msqrt></mrow></mfrac><mo>.</mo></mrow></math> 3 U 1 approximation with Anlı-Güngör (AG) Phase Function 1 In the second part of this study, the Anlı-Güngör (AG) phase function is used in the transport equation with the defined parameters in the previous section [ 15 15 ], <math display="block" id="M22"><mrow><msub><mi>σ</mi><mi>S</mi></msub><msup><mrow/><mrow><mi>A</mi><mi>G</mi></mrow></msup><mo stretchy="false">(</mo><msub><mi>μ</mi><mn>0</mn></msub><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mi>σ</mi><mi>S</mi></msub></mrow><mrow><mn>4</mn><mi>π</mi><mtext> </mtext><msup><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><mn>2</mn><msub><mi>μ</mi><mn>0</mn></msub><mi>t</mi><mo>+</mo><msup><mi>t</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow></mfrac></mrow></math> The integral of the AG phase function is needed and is replaced in the right hand side of the one-dimensional transport equation given in Eq. (1), <math display="block" id="M23"><mrow><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><mrow><mn>2</mn><mi>π</mi></mrow></munderover><mrow><msubsup><mi>σ</mi><mi>S</mi><mrow><mi>A</mi><mi>G</mi></mrow></msubsup><mo stretchy="false">(</mo><msub><mi>μ</mi><mn>0</mn></msub><mo stretchy="false">)</mo><mtext> </mtext><mtext>d</mtext><msup><mi>ϕ</mi><mo>′</mo></msup><mo>=</mo><mfrac><mrow><msub><mi>σ</mi><mi>S</mi></msub></mrow><mn>2</mn></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mi>∞</mi></munderover><mrow><msup><mi>t</mi><mi>n</mi></msup><msub><mi>P</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>μ</mi><mo stretchy="false">)</mo><msub><mi>P</mi><mi>n</mi></msub><mo stretchy="false">(</mo><msup><mi>μ</mi><mo>′</mo></msup><mo stretchy="false">)</mo></mrow></mstyle></mrow></mrow></mstyle><mo>.</mo></mrow></math> The transport equation with the integral transformation of the AG phase function given in Eq. (14) can then be written as, <math display="block" id="M24"><mrow><mi>μ</mi><mfrac><mrow><mo>∂</mo><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac><mo>+</mo><mi>ν</mi><mtext> </mtext><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mi>ν</mi><mtext> </mtext><mi>c</mi></mrow><mn>2</mn></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mi>∞</mi></munderover><mrow><msup><mi>t</mi><mi>n</mi></msup><msub><mi>P</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>μ</mi><mo stretchy="false">)</mo><msub><mi>α</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow></mstyle><mo>,</mo></mrow></math> where α n ( x ) is defined as <math display="block" id="M25"><mrow><msub><mi>α</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mrow><mo>−</mo><mn>1</mn></mrow><mn>1</mn></msubsup><mrow><msub><mi>P</mi><mi>n</mi></msub><mo stretchy="false">(</mo><msup><mi>μ</mi><mo>′</mo></msup><mo stretchy="false">)</mo><mrow><mo>[</mo><mrow><mfrac><mn>2</mn><mi>π</mi></mfrac><msqrt><mrow><mn>1</mn><mo>−</mo><msup><msup><mi>μ</mi><mo>′</mo></msup><mn>2</mn></msup></mrow></msqrt><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>∞</mi></munderover><mrow><msub><mi>Φ</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow></mstyle><msub><mi>U</mi><mi>n</mi></msub><mo stretchy="false">(</mo><msup><mi>μ</mi><mo>′</mo></msup><mo stretchy="false">)</mo></mrow><mo>]</mo></mrow><mtext> </mtext><mtext>d</mtext><msup><mi>μ</mi><mo>′</mo></msup></mrow></mrow></mstyle><mo>.</mo></mrow></math> The U N moments of the angular flux can easily be obtained by following the same procedure described for Eq. (9) and Eq. (10): <math display="block" id="M26"><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>Φ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><mi>x</mi></mrow></mfrac><mo>+</mo><mn>2</mn><mi>v</mi><mtext> </mtext><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mn>2</mn><mi>v</mi><mi>c</mi><mtext> </mtext><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>,</mo></mrow></math> <math display="block" id="M27"><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>Φ</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><mi>x</mi></mrow></mfrac><mo>+</mo><mfrac><mrow><mtext>d</mtext><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><mi>x</mi></mrow></mfrac><mo>+</mo><mn>2</mn><mi>v</mi><mtext> </mtext><msub><mi>Φ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>2</mn><mn>3</mn></mfrac><mi>v</mi><mi>c</mi><mi>t</mi><mi>Φ</mi><msub><mtext> </mtext><mn>1</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>.</mo></mrow></math> for n = 0 and n = 1, respectively. Then, Eqs. (17) and (18) can be solved simultaneously for Φ 0 ( x ) by following the same procedure for Eq. (11):: <math display="block" id="M28"><mrow><mfrac><mrow><msup><mtext>d</mtext><mtext>2</mtext></msup><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><msup><mi>x</mi><mn>2</mn></msup></mrow></mfrac><mo>−</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><msup><mi>v</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>c</mi><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>3</mn><mo>−</mo><mi>c</mi><mi>t</mi><mo stretchy="false">)</mo><mtext> </mtext><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mn>0.</mn></mrow></math> The diffusion length in the case of the AG phase function can then be obtained by <math display="block" id="M29"><mrow><msup><mi>L</mi><mrow><mtext>AG</mtext></mrow></msup><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>v</mi></mrow></mfrac><msqrt><mrow><mfrac><mn>3</mn><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>c</mi><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>3</mn><mo>−</mo><mi>c</mi><mi>t</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></msqrt><mo>.</mo></mrow></math> 4 Numerical Results 1 The diffusion approximation was performed for one-speed neutrons in a homogeneous finite slab using both conventional HG and AG phase functions. In the solution algorithm, the neutron angular flux was expanded in a series of the Chebyshev polynomials of the second kind. This problem has been investigated by Öztürk and Anlı using the AG phase function with Legendre and Chebyshev polynomials of the first kind [ 16 16 ]. Then, analytic expressions for the diffusion lengths were derived using the U 1 approximation with the HG and AG phase functions and the numerical results were calculated for various values of the cross-section parameters c and t . Finally, these were tabulated side by side for comparison. The parameter ν relating to the total macroscopic cross-section was assumed to be ν = 1 cm -1 . Maple software was used for all of the calculations. In Tables 1 Tables 1 and 2 2 , the diffusion lengths were calculated from a weakly absorbing medium ( c = 0) to a highly scattering medium ( c = 1) for t in the interval <math id="M30"><mrow><mo>−</mo><mn>1</mn><mo>≤</mo><mi>t</mi><mo>≤</mo><mn>1</mn></mrow></math> . The same results for the diffusion length were obtained for a weakly absorbing medium ( c = 0) in the case of both HG and AG phase functions (0.50000) for all values of the parameter t . As can be seen from these tables, the calculated diffusion lengths obtained from both HG and AG phase functions increased with increasing values of c and t . 1 1 t 2 1 c = 0.99 2 1 c = 0.98 2 1 c = 0.95 1 1   1 1 HG 1 1 AG 1 1 HG 1 1 AG 1 1 HG 1 1 AG 1 1 -1.00 1 1 3.5444 1 1 4.3356 1 1 2.5126 1 1 3.0696 1 1 1.6013 1 1 1.9487 1 1 -0.99 1 1 3.5533 1 1 4.3409 1 1 2.5188 1 1 3.0733 1 1 1.6052 1 1 1.9511 1 1 -0.75 1 1 3.7878 1 1 4.4767 1 1 2.6841 1 1 3.1686 1 1 1.7087 1 1 2.0101 1 1 -0.50 1 1 4.0893 1 1 4.6324 1 1 2.8964 1 1 3.2780 1 1 1.8412 1 1 2.0776 1 1 -0.25 1 1 4.4766 1 1 4.8057 1 1 3.1686 1 1 3.3995 1 1 2.0101 1 1 2.1525 1 1 0.00 1 1 5.0000 1 1 5.0000 1 1 3.5355 1 1 3.5355 1 1 2.2361 1 1 2.2361 1 1 0.25 1 1 5.7639 1 1 5.2200 1 1 4.0689 1 1 3.6894 1 1 2.5607 1 1 2.3302 1 1 0.50 1 1 7.0360 1 1 5.4718 1 1 4.9507 1 1 3.8653 1 1 3.0861 1 1 2.4373 1 1 0.75 1 1 9.8533 1 1 5.7639 1 1 6.8680 1 1 4.0689 1 1 4.1703 1 1 2.5607 1 1 0.99 1 1 35.4441 1 1 6.0935 1 1 20.4808 1 1 4.2982 1 1 9.1670 1 1 2.6988 1 1 1.00 1 1 50.0000 1 1 6.1085 1 1 25.0000 1 1 4.3086 1 1 10.0000 1 1 2.7050 1 1 t 2 1 c = 0.90 2 1 c = 0.80 2 1 c = 0.50 2 1 c = 0.00 1 1   1 1 HG 1 1 AG 1 1 HG 1 1 AG 1 1 HG 1 1 AG 1 1 HG 1 1 AG 1 1 -1.00 1 1 1.1471 1 1 1.3868 1 1 0.8333 1 1 0.9934 1 1 0.5774 1 1 0.6547 1 1 0.5000 1 1 0.5000 1 1 -0.99 1 1 1.1498 1 1 1.3884 1 1 0.8352 1 1 0.9945 1 1 0.5783 1 1 0.6551 1 1 0.5000 1 1 0.5000 1 1 -0.75 1 1 1.2217 1 1 1.4286 1 1 0.8839 1 1 1.0206 1 1 0.6030 1 1 0.6667 1 1 0.5000 1 1 0.5000 1 1 -0.50 1 1 1.3131 1 1 1.4744 1 1 0.9449 1 1 1.0502 1 1 0.6325 1 1 0.6794 1 1 0.5000 1 1 0.5000 1 1 -0.25 1 1 1.4286 1 1 1.5250 1 1 1.0206 1 1 1.0825 1 1 0.6667 1 1 0.6928 1 1 0.5000 1 1 0.5000 1 1 0.00 1 1 1.5811 1 1 1.5811 1 1 1.1180 1 1 1.1180 1 1 0.7071 1 1 0.7071 1 1 0.5000 1 1 0.5000 1 1 0.25 1 1 1.7961 1 1 1.6440 1 1 1.2500 1 1 1.1573 1 1 0.7559 1 1 0.7223 1 1 0.5000 1 1 0.5000 1 1 0.50 1 1 2.1320 1 1 1.7150 1 1 1.4434 1 1 1.2010 1 1 0.8165 1 1 0.7386 1 1 0.5000 1 1 0.5000 1 1 0.75 1 1 2.7735 1 1 1.7961 1 1 1.7678 1 1 1.2500 1 1 0.8944 1 1 0.7559 1 1 0.5000 1 1 0.5000 1 1 0.99 1 1 4.7891 1 1 1.8858 1 1 2.4515 1 1 1.3032 1 1 0.9950 1 1 0.7738 1 1 0.5000 1 1 0.5000 1 1 1.00 1 1 5.0000 1 1 1.8898 1 1 2.5000 1 1 1.3056 1 1 1.0000 1 1 0.7746 1 1 0.5000 1 1 0.5000 A comparison table for the calculated diffusion lengths obtained from the present U 1 and those obtained from the T 1 and P 1 approximations for selected values of the scattering parameters are given in Table 3 Table 3 . The diffusion lengths obtained from T 1 and P 1 approximations tabulated in Table 3 Table 3 are quoted from Öztürk and Anlı [ 16 16 , 17 17 ]. A general agreement, except for the case of highly scattering medium, was demonstrated between the results obtained from the HG and AG phase functions. However, in the case of a highly scattering medium and in particular for the forward peaked scattering of the neutrons, the diffusion lengths calculated from the HG phase function differed with the results calculated from the AG phase function. In addition, as expected the numerical results obtained from all methods were similar for all cases, since they all belong to the same polynomial family, i.e., Jacobi polynomials. 1 1 t 6 1 c = 0.99 6 1 c = 0.90 6 1 c = 0.50 1 1   3 1 HG 3 1 AG 3 1 HG 3 1 AG 3 1 HG 3 1 AG 1 1   1 1 U 1 1 1 T 1 [ 16 16 ] 1 1 P 1 [ 16 16 , 17 17 ] 1 1 U 1 1 1 T 1 [ 16 16 ] 1 1 P 1 [ 16 16 , 17 17 ] 1 1 U 1 1 1 T 1 [ 16 16 ] 1 1 P 1 [ 16 16 , 17 17 ] 1 1 U 1 1 1 T 1 [ 16 16 ] 1 1 P 1 [ 16 16 , 17 17 ] 1 1 U 1 1 1 T 1 [ 16 16 ] 1 1 P 1 [ 16 16 , 17 17 ] 1 1 U 1 1 1 T 1 [ 16 16 ] 1 1 P 1 [ 16 16 , 17 17 ] 1 1 -1.00 1 1 3.5444 1 1 5.0125 1 1 4.0927 1 1 4.3356 1 1 6.1314 1 1 5.0063 1 1 1.1471 1 1 1.6222 1 1 1.3245 1 1 1.3868 1 1 1.9612 1 1 1.6013 1 1 0.5774 1 1 0.8165 1 1 0.6667 1 1 0.6547 1 1 0.9258 1 1 0.7559 1 1 -0.99 1 1 3.5533 1 1 5.0251 1 1 4.1029 1 1 4.3409 1 1 6.1390 1 1 5.0125 1 1 1.1498 1 1 1.6261 1 1 1.3277 1 1 1.3884 1 1 1.9634 1 1 1.6031 1 1 0.5783 1 1 0.8179 1 1 0.6678 1 1 0.6551 1 1 0.9265 1 1 0.7565 1 1 -0.75 1 1 3.7878 1 1 5.3567 1 1 4.3737 1 1 4.4767 1 1 6.3309 1 1 5.1692 1 1 1.2217 1 1 1.7277 1 1 1.4107 1 1 1.4286 1 1 2.0203 1 1 1.6496 1 1 0.6030 1 1 0.8528 1 1 0.6963 1 1 0.6667 1 1 0.9428 1 1 0.7698 1 1 -0.50 1 1 4.0893 1 1 5.7831 1 1 4.7219 1 1 4.6324 1 1 6.5512 1 1 5.3490 1 1 1.3131 1 1 1.8570 1 1 1.5162 1 1 1.4744 1 1 2.0851 1 1 1.7025 1 1 0.6325 1 1 0.8944 1 1 0.7303 1 1 0.6794 1 1 0.9608 1 1 0.7845 1 1 -0.25 1 1 4.4766 1 1 6.3309 1 1 5.1691 1 1 4.8057 1 1 6.7963 1 1 5.5491 1 1 1.4286 1 1 2.0203 1 1 1.6496 1 1 1.5250 1 1 2.1567 1 1 1.7609 1 1 0.6667 1 1 0.9428 1 1 0.7698 1 1 0.6928 1 1 0.9798 1 1 0.8000 1 1 0.00 1 1 5.0000 1 1 7.0711 1 1 5.7735 1 1 5.0000 1 1 7.0711 1 1 5.7735 1 1 1.5811 1 1 2.2361 1 1 1.8257 1 1 1.5811 1 1 2.2361 1 1 1.8257 1 1 0.7071 1 1 1.0000 1 1 0.8165 1 1 0.7071 1 1 1.0000 1 1 0.8165 1 1 0.25 1 1 5.7639 1 1 8.1514 1 1 6.6556 1 1 5.2200 1 1 7.3821 1 1 6.0275 1 1 1.7961 1 1 2.5400 1 1 2.0739 1 1 1.6440 1 1 2.3250 1 1 1.8983 1 1 0.7559 1 1 1.0690 1 1 0.8729 1 1 0.7223 1 1 1.0215 1 1 0.8341 1 1 0.50 1 1 7.0360 1 1 9.9504 1 1 8.1244 1 1 5.4718 1 1 7.7382 1 1 6.3182 1 1 2.1320 1 1 3.0151 1 1 2.4618 1 1 1.7150 1 1 2.4254 1 1 1.9803 1 1 0.8165 1 1 1.1547 1 1 0.9428 1 1 0.7386 1 1 1.0445 1 1 0.8528 1 1 0.75 1 1 9.8533 1 1 13.9347 1 1 11.3776 1 1 5.7639 1 1 8.1514 1 1 6.6556 1 1 2.7735 1 1 3.9223 1 1 3.2026 1 1 1.7961 1 1 2.5400 1 1 2.0739 1 1 0.8944 1 1 1.2649 1 1 1.0328 1 1 0.7559 1 1 1.0690 1 1 0.8729 1 1 0.99 1 1 35.4441 1 1 50.1255 1 1 40.9273 1 1 6.0935 1 1 8.6175 1 1 7.0362 1 1 4.7891 1 1 6.7729 1 1 5.5300 1 1 1.8858 1 1 2.6669 1 1 2.1775 1 1 0.9950 1 1 1.4072 1 1 1.1490 1 1 0.7738 1 1 1.0944 1 1 0.8935 1 1 1.00 1 1 50.0000 1 1 70.7107 1 1 57.7350 1 1 6.1085 1 1 8.6387 1 1 7.0535 1 1 5.0000 1 1 7.0711 1 1 5.7735 1 1 1.8898 1 1 2.6726 1 1 2.1822 1 1 1.0000 1 1 1.4142 1 1 1.1547 1 1 0.7746 1 1 1.0954 1 1 0.8944 The diffusion lengths calculated from the HG phase function were very similar to the results obtained from the AG phase function if the neutrons moved in a weakly absorbing medium, i.e., c ~0. However, if the neutrons moved in a highly scattering medium, i.e., c ~1, the results for the diffusion lengths calculated by the HG phase function differed from those from the AG phase function. This situation was especially observed for t > 0.75, corresponding to the case of strongly forward scattering [ 15 15 ]. The discrepancies between the diffusion lengths were not observed when these calculations were obtained by using the AG phase function. This was also observed in the study by Öztürk and Anlı, which used the AG phase function with the Legendre and Chebyshev polynomials of the first kind [ 16 16 ]. In addition, the diffusion length corresponding to t = 0 is equivalent to the case of isotropic scattering. 5 Conclusion 1 In this work, the AG phase function was used in the scattering function of the transport equation to calculate the diffusion lengths for one-speed neutrons in a uniform homogeneous slab. The neutron angular flux was expanded in terms of the Chebyshev polynomials of the second kind and the diffusion lengths were calculated for various values of c and t . The HG and/or AG phase functions were examined to determine if they can be applied to other problems in science and engineering, and to further suggest new phase functions and methods to be developed for more accurate solutions of the transport equation. From the tabulated numerical results, we concluded that both the HG and AG phase functions can be applied to transport theory to calculate the diffusion lengths for one-speed neutrons in a slab, except in the case of a highly scattering medium with forward peaked scattering. In this case while the diffusion lengths obtained from the HG phase function were unexpectedly high, the diffusion lengths obtained from the AG phase function were realistic. Therefore, the AG phase function is more applicable than the HG phase function for this specific case. In addition, it was seen from the provided tables, that realistic results were always obtained from the AG phase function, independently of the scattering properties of the medium. References J.R. Lamarsh , A.J. Baratta , Introduction to nuclear engineering, 3rd edn . ( Prentice Hall, Inc ., 2001 ), pp. 52 - 110 Introduction to nuclear engineering, 3rd edn S. Yabushita , Tschebyscheff polynomials approximation method of the neutron transport equation . J. Math. Phys . 2 , 543 - 549 ( 1961 ). doi: 10.1063/1.1703739 http://doi.org/10.1063/1.1703739 Tschebyscheff polynomials approximation method of the neutron transport equation H. Öztürk , F. Anli , S. Güngör , Application of the UN method to the reflected critical slab problem for one-speed neutrons with forward and backward scattering . Kerntechnik 72 , 74 - 76 ( 2007 ). doi: 10.3139/124.100321 http://doi.org/10.3139/124.100321 Application of the UN method to the reflected critical slab problem for one-speed neutrons with forward and backward scattering H. Öztürk , The reflected critical slab problem for one-speed neutrons with strongly anisotropic scattering . Kerntechnik 73 , 66 - 74 ( 2008 ). doi: 10.3139/124.100532 http://doi.org/10.3139/124.100532 The reflected critical slab problem for one-speed neutrons with strongly anisotropic scattering M. Tıraş , H. Öztürk , A. Bülbül , T1 and U1 approximations to neutron transport equation in one-dimensional spherical geometry . Kerntechnik 79 , 163 - 166 ( 2014 ). doi: 10.3139/124.110415 http://doi.org/10.3139/124.110415 T1 and U1 approximations to neutron transport equation in one-dimensional spherical geometry D.C. Sahni , N.G. Sjöstrand , N.S. Garis , Criticality and time eigenvalues for one-speed neutrons in a slab with forward and backward scattering . J. Phys. D: Appl. Phys . 25 , 1381 - 1389 ( 1992 ). doi: 10.1088/0022-3727/25/10/001 http://doi.org/10.1088/0022-3727/25/10/001 Criticality and time eigenvalues for one-speed neutrons in a slab with forward and backward scattering C.E. Lee , M.P. Dias , Analytical solutions to the moment transport equations-I; one-group one-region slab and sphere criticality . Ann. Nucl. Energy 11 , 515 - 530 ( 1984 ). doi: 10.1016/0306-4549(84)90076-8 http://doi.org/10.1016/0306-4549(84)90076-8 Analytical solutions to the moment transport equations-I; one-group one-region slab and sphere criticality B. Hunter , Z. Guo , Phase-function normalization in the 3-D discrete-ordinates solution of radiative transfer—part I: Conservation of scattered energy and asymmetry factor . Numer. Heat Tr. Part B , 62 , 203 - 222 ( 2012 ). doi: 10.1080/10407790.2012.709163 http://doi.org/10.1080/10407790.2012.709163 Phase-function normalization in the 3-D discrete-ordinates solution of radiative transfer—part I: Conservation of scattered energy and asymmetry factor B. Hunter , Z. Guo , Applicability of phase-function normalization techniques for radiation transfer computation . Numer. Heat Tr. Part B , 67 , 1 - 24 ( 2015 ). doi: 10.1080/10407790.2014.949516 http://doi.org/10.1080/10407790.2014.949516 Applicability of phase-function normalization techniques for radiation transfer computation B. Hunter , Z. Guo , A new and simple technique to normalize the HG phase function for conserving scattered energy and asymmetry factor . Numer. Heat Tr. Part B , 65 , 195 - 217 ( 2014 ). doi: 10.1080/10407790.2013.849992 http://doi.org/10.1080/10407790.2013.849992 A new and simple technique to normalize the HG phase function for conserving scattered energy and asymmetry factor F. Zhang , J. Li , A note on double Henyey-Greenstein phase function . J. Quant. Specrosc. Rad. Transfer 184 , 40 - 43 ( 2016 ). doi: 10.1016/j.jqsrt.2016.06.016 http://doi.org/10.1016/j.jqsrt.2016.06.016 A note on double Henyey-Greenstein phase function A. Kaşkaş , C. Tezcan , An Application of Transport Theory in Optical Oceanography: The Estimation of the Apparent Optical Properties Using Henyey-Greenstein Phase Function . Transport Theor. Stat. Phys . 38 , 317 - 329 ( 2009 ). doi: 10.1080/00411450903372126 http://doi.org/10.1080/00411450903372126 An Application of Transport Theory in Optical Oceanography: The Estimation of the Apparent Optical Properties Using Henyey-Greenstein Phase Function A. Kaşkaş , C. Tezcan , M.Ç. Güleçyüz , An application of transport theory for optical oceanography: the isotropic point source . J. Opt. A-Pure Appl. Op . 8 , 683 - 688 ( 2006 ). doi: 10.1088/1464-4258/8/8/010 http://doi.org/10.1088/1464-4258/8/8/010 An application of transport theory for optical oceanography: the isotropic point source M. Roger , C. Caliot ,  N. Crouseilles , P.J. Coelho , A hybrid transport-diffusion model for radiative transfer in absorbing and scattering media . J. Comp. Phys . 275 , 346 - 362 ( 2014 ). doi: 10.1016/j.jcp.2014.06.063 http://doi.org/10.1016/j.jcp.2014.06.063 A hybrid transport-diffusion model for radiative transfer in absorbing and scattering media F. Anli , F. Yasa , S. Güngör , General eigenvalue spectrum in a one-dimensional slab geometry transport equation . Nucl. Sci. Eng . 150 , 72 - 77 ( 2005 ). General eigenvalue spectrum in a one-dimensional slab geometry transport equation H. Öztürk , F. Anli , Diffusion approximation for certain scattering parameters of the Anli-Güngör phase function . Kerntechnik 77 , 381 - 384 ( 2012 ). doi: 10.3139/124.110216 http://doi.org/10.3139/124.110216 Diffusion approximation for certain scattering parameters of the Anli-Güngör phase function H. Öztürk , Application of the Henyey-Greenstein and Anlı-Güngör phase functions for the solution of the neutron transport equation with Legendre polynomials: Reflected critical slab problem . Kerntechnik 78 , 447 - 453 ( 2013 ). doi: 10.3139/124.110328 http://doi.org/10.3139/124.110328 Application of the Henyey-Greenstein and Anlı-Güngör phase functions for the solution of the neutron transport equation with Legendre polynomials: Reflected critical slab problem B. Davison , Neutron transport theory . 1 st edn. ( London , Oxford University Press , 1958 ), pp. 116 - 144 Neutron transport theory L.G. Henyey , J.L. Greenstein , Diffuse radiation in the galaxy . Astrophys. J . 93 , 70 - 83 ( 1941 ). doi: 10.1086/144246 http://doi.org/10.1086/144246 Diffuse radiation in the galaxy G. Arfken , Mathematical methods for physicists . 3 rd edn. ( London , Academic Press, Inc ., 1985 ), pp. 817 - 879 Mathematical methods for physicists 10.1007/s41365-018-0380-6 This study was supported by the Academic Research Projects Unit of Osmaniye Korkut Ata University (No. OKÜBAP-2014-PT3-019). 2017-01-28 2017-07-03 2017-07-21 2018-02-19 2018-03-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2018 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-018-0380-6 http://dx.doi.org/10.1007/s41365-018-0380-6 . 2018 Nuclear Science and Techniques 03 29