Solution of the finite slab criticality problem using an alternative phase function with the second kind of Chebyshev polynomials

Solution of the finite slab criticality problem using an alternative phase function with the second kind of Chebyshev polynomials 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. The critical size of a finite homogenous slab is investigated for one-speed neutrons using the alternative phase function (AG, Anlı-Güngör) in place of the scattering function of the transport equation. First, the neutron angular flux expanded in terms of the Chebyshev polynomials of second kind ( U N approximation) together with the AG phase function is applied to the transport equation to obtain a criticality condition for the system. Then, using various values of the scattering parameters, the numerical results for the critical half-thickness of the slab are calculated and they are tabulated in the tables together with the ones obtained from the conventional spherical harmonics ( P N ) method for comparison. They can be said to be in good accordance with each other. Criticality problem U N method neutron transport equation alternative phase function. 1 Introduction 1 The particle transport equation, which was first developed by Boltzmann for the kinetic theory of gases, is based on the conservation of the neutrons in a reactor system. The radiative transfer of stellar and planetary atmosphere and light scattering phenomenon are also related with the transport concept. Therefore, the description of the behavior of the neutrons in a reactor has a great importance for the first calculation and thus the construction and the operation of the reactor uneventfully. The number of fission neutrons is wanted to be constant in all types of reactors in order to obtain constant power and to control it safely. This situation of the reactor is defined as to be critical and the criticality of a fission system is one of the most important problems in neutron transport theory. Therefore, the critical size of a reactor can be said to be decided after the investigation of the criticality problem related with the system under consideration. As well known, the transport equation is an integro-differential equation and thus it is not easy to solve it analytically. The discrete ordinates ( S N ) and polynomial expansion based techniques are accepted to be the most common and powerful ones among the methods developed for the solution of the transport equation [ 1 1 - 4 4 1-4 ]. In about all methods, either the derivative of the neutron angular flux or the neutron scattering function presented in the integral part of the transport equation are treated by some approximations to simplify the solution of the equation. In some instances, using an approximated scattering function in the transport equation can be sufficient depending on the scope of the problem under consideration. However, these approximations can take the problems and thus the solutions more or less away from the real situations. Since the scattering cross-sections vary with the scattering angle incredibly, various difficulties occur when the scattering function is represented in terms of any polynomials. Therefore, instead of using approximate expressions, if an exact scattering model is used in place of the scattering function, one could obtain more realistic results representing the system better [ 1 1 , 2 2 ]. Henyey and Greenstein [ 5 5 ] had first developed an exact scattering function called the Henyey-Greenstein (HG) formula to verify the existence of the diffuse interstellar radiation. However, they did not explain physically in their paper why they used such a function in the radiative transport equation. Later in the following studies such as biomedical applications, significant discrepancies have been reported about using the Henyey-Greenstein formula [ 5 5 - 7 7 5-7 ]. In this paper, an alternative scattering function (Anlı-Güngör, AG) which was applied successfully to the criticality problem in neutron transport theory using Legendre and the Chebyshev polynomials of first kind is preferred [ 9 9 - 11 11 9-11 ]. However, instead of using the first kind of Chebyshev polynomials, the second kind of Chebyshev polynomials approximation ( U N method) is preferred to serve as an alternative solution method to the literature. This method has been successfully applied to the solutions of the problems of the transport equation in nearly a decade [ 13 13 , 15 15 ]. Therefore, this study can be thought of as the extension of the study previously carried out by Öztürk [ 11 11 ]. In the solution algorithm, the neutron angular flux is first expanded in terms of the Chebyshev polynomials of second kind then the AG phase function is inserted into the transport equation in place of the scattering function. At the end, applying the U N method to the resultant equation following the moment equations, a general expression for the criticality condition is obtained for one-speed neutrons. Then, the critical half-thicknesses of the slab are calculated using various scattering and collision parameters in the criticality condition. The numerical results are obtained with an increasing order of N = 1 to 9 and they are listed in Tables 1 Tables 1 , 2 2 , 3 3 and 4 4 . Finally, a comparison table including the results obtained from the present method and the results obtained from the conventional P N method is given. 1 1 t 1 1 U 1 1 1 U 3 1 1 U 5 1 1 U 7 1 1 U 9 1 1 -1 1 1 6.35409 1 1 7.29973 1 1 7.27243 1 1 7.27499 1 1 7.27369 1 1 -4/5 1 1 6.50873 1 1 7.48101 1 1 7.45336 1 1 7.45608 1 1 7.45448 1 1 -3/4 1 1 6.54923 1 1 7.52819 1 1 7.50059 1 1 7.50333 1 1 7.50195 1 1 -2/3 1 1 6.61848 1 1 7.60864 1 1 7.58120 1 1 7.58398 1 1 7.58252 1 1 -1/2 1 1 6.76399 1 1 7.77692 1 1 7.75001 1 1 7.75283 1 1 7.75102 1 1 -3/10 1 1 6.95230 1 1 7.99354 1 1 7.96732 1 1 7.97020 1 1 7.96833 1 1 -1/4 1 1 7.00194 1 1 8.05049 1 1 8.02441 1 1 8.02729 1 1 8.02529 1 1 -1/5 1 1 7.05270 1 1 8.10865 1 1 8.08269 1 1 8.08560 1 1 8.08387 1 1 0 1 1 7.26764 1 1 8.35447 1 1 8.32866 1 1 8.33159 1 1 8.32957 1 1 1/5 1 1 7.50408 1 1 8.62419 1 1 8.59787 1 1 8.60083 1 1 8.59894 1 1 1/4 1 1 7.56698 1 1 8.69585 1 1 8.66929 1 1 8.67226 1 1 8.66762 1 1 3/10 1 1 7.63154 1 1 8.76937 1 1 8.74251 1 1 8.74549 1 1 8.74048 1 1 1/2 1 1 7.90785 1 1 9.08370 1 1 9.05518 1 1 9.05821 1 1 9.05398 1 1 2/3 1 1 8.16329 1 1 9.37393 1 1 9.34341 1 1 9.34647 1 1 9.34066 1 1 3/4 1 1 8.30095 1 1 9.53021 1 1 9.49849 1 1 9.50158 1 1 9.49957 1 1 4/5 1 1 8.38706 1 1 9.62794 1 1 9.59543 1 1 9.59858 1 1 9.59680 1 1 1 1 1 8.76137 1 1 10.05244 1 1 10.01637 1 1 10.01973 1 1 10.01662 1 1 t 1 1 U 1 1 1 U 3 1 1 U 5 1 1 U 7 1 1 U 9 1 1 -1 1 1 4.36099 1 1 4.98773 1 1 4.95863 1 1 4.96141 1 1 4.95998 1 1 -4/5 1 1 4.46439 1 1 5.11058 1 1 5.08128 1 1 5.08421 1 1 5.08269 1 1 -3/4 1 1 4.49144 1 1 5.14236 1 1 5.11315 1 1 5.11610 1 1 5.11456 1 1 -2/3 1 1 4.53769 1 1 5.19641 1 1 5.16742 1 1 5.17040 1 1 5.16882 1 1 -1/2 1 1 4.63477 1 1 5.30895 1 1 5.28061 1 1 5.28362 1 1 5.28200 1 1 -3/10 1 1 4.76024 1 1 5.45300 1 1 5.42548 1 1 5.42852 1 1 5.42682 1 1 -1/4 1 1 4.79329 1 1 5.49075 1 1 5.46339 1 1 5.46644 1 1 5.46471 1 1 -1/5 1 1 4.82706 1 1 5.52926 1 1 5.50204 1 1 5.50509 1 1 5.50338 1 1 0 1 1 4.96992 1 1 5.69153 1 1 5.66454 1 1 5.66761 1 1 5.66581 1 1 1/5 1 1 5.12677 1 1 5.86883 1 1 5.84138 1 1 5.84447 1 1 5.84261 1 1 1/4 1 1 5.16846 1 1 5.91583 1 1 5.88814 1 1 5.89123 1 1 5.88928 1 1 3/10 1 1 5.21121 1 1 5.96399 1 1 5.93601 1 1 5.93911 1 1 5.93711 1 1 1/2 1 1 5.39393 1 1 6.16940 1 1 6.13975 1 1 6.14288 1 1 6.14085 1 1 2/3 1 1 5.56249 1 1 6.35836 1 1 6.32668 1 1 6.32984 1 1 6.32764 1 1 3/4 1 1 5.65317 1 1 6.45983 1 1 6.42692 1 1 6.43011 1 1 6.42782 1 1 4/5 1 1 5.70985 1 1 6.52319 1 1 6.48948 1 1 6.49270 1 1 6.49033 1 1 1 1 1 5.95572 1 1 6.79755 1 1 6.76016 1 1 6.76360 1 1 6.76113 1 1 t 1 1 U 1 1 1 U 3 1 1 U 5 1 1 U 7 1 1 U 9 1 1 -1 1 1 1.08842 1 1 1.16747 1 1 1.13348 1 1 1.13843 1 1 1.13597 1 1 -4/5 1 1 1.11017 1 1 1.19975 1 1 1.16634 1 1 1.17137 1 1 1.16891 1 1 -3/4 1 1 1.11584 1 1 1.20746 1 1 1.17423 1 1 1.17924 1 1 1.17679 1 1 -2/3 1 1 1.12549 1 1 1.22006 1 1 1.18717 1 1 1.19217 1 1 1.18972 1 1 -1/2 1 1 1.14565 1 1 1.24468 1 1 1.21255 1 1 1.21748 1 1 1.21505 1 1 -3/10 1 1 1.17147 1 1 1.27373 1 1 1.24257 1 1 1.24737 1 1 1.24496 1 1 -1/4 1 1 1.17823 1 1 1.28099 1 1 1.25005 1 1 1.25481 1 1 1.25241 1 1 -1/5 1 1 1.18512 1 1 1.28827 1 1 1.25754 1 1 1.26226 1 1 1.25986 1 1 0 1 1 1.21406 1 1 1.31772 1 1 1.28759 1 1 1.29212 1 1 1.28974 1 1 1/5 1 1 1.24547 1 1 1.34808 1 1 1.31805 1 1 1.32237 1 1 1.31998 1 1 1/4 1 1 1.25376 1 1 1.35587 1 1 1.32576 1 1 1.33002 1 1 1.32763 1 1 3/10 1 1 1.26222 1 1 1.36374 1 1 1.33352 1 1 1.33772 1 1 1.33532 1 1 1/2 1 1 1.29809 1 1 1.39623 1 1 1.36508 1 1 1.36905 1 1 1.36662 1 1 2/3 1 1 1.33071 1 1 1.42471 1 1 1.39214 1 1 1.39596 1 1 1.39347 1 1 3/4 1 1 1.34808 1 1 1.43949 1 1 1.40596 1 1 1.40972 1 1 1.40719 1 1 4/5 1 1 1.35887 1 1 1.44855 1 1 1.41434 1 1 1.41809 1 1 1.41553 1 1 1 1 1 1.40513 1 1 1.48641 1 1 1.44860 1 1 1.45245 1 1 1.44965 1 1 t 1 1 U 1 1 1 U 3 1 1 U 5 1 1 U 7 1 1 U 9 1 1 -1 1 1 0.32188 1 1 0.30549 1 1 0.27813 1 1 0.27642 1 1 0.27071 1 1 -4/5 1 1 0.32719 1 1 0.31680 1 1 0.29045 1 1 0.28985 1 1 0.28463 1 1 -3/4 1 1 0.32856 1 1 0.31926 1 1 0.29298 1 1 0.29252 1 1 0.28734 1 1 -2/3 1 1 0.33090 1 1 0.32308 1 1 0.29684 1 1 0.29654 1 1 0.29142 1 1 -1/2 1 1 0.33574 1 1 0.32980 1 1 0.30350 1 1 0.30335 1 1 0.29832 1 1 -3/10 1 1 0.34190 1 1 0.33655 1 1 0.31016 1 1 0.30997 1 1 0.30503 1 1 -1/4 1 1 0.34350 1 1 0.33805 1 1 0.31165 1 1 0.31143 1 1 0.30651 1 1 -1/5 1 1 0.34512 1 1 0.33948 1 1 0.31308 1 1 0.31281 1 1 0.30792 1 1 0 1 1 0.35191 1 1 0.34465 1 1 0.31826 1 1 0.31771 1 1 0.31294 1 1 1/5 1 1 0.35917 1 1 0.34906 1 1 0.32266 1 1 0.32170 1 1 0.31705 1 1 1/4 1 1 0.36107 1 1 0.35007 1 1 0.32363 1 1 0.32255 1 1 0.31793 1 1 3/10 1 1 0.36301 1 1 0.35105 1 1 0.32456 1 1 0.32335 1 1 0.31175 1 1 1/2 1 1 0.37113 1 1 0.35465 1 1 0.32772 1 1 0.32591 1 1 0.32135 1 1 2/3 1 1 0.37841 1 1 0.35736 1 1 0.32960 1 1 0.32714 1 1 0.32249 1 1 3/4 1 1 0.38224 1 1 0.35866 1 1 0.33024 1 1 0.32738 1 1 0.32260 1 1 4/5 1 1 0.38461 1 1 0.35943 1 1 0.33052 1 1 0.32739 1 1 0.32249 1 1 1 1 1 0.39462 1 1 0.36258 1 1 0.33076 1 1 0.32612 1 1 0.32021 2 U N method with AG phase function 1 The stationary transport equation for one energy group neutrons in a source free medium can be written as, <math display="block" id="M1"><mrow><mstyle mathvariant="bold" mathsize="normal"><mi>Ω</mi></mstyle><mo>⋅</mo><mo>∇</mo><mi>ψ</mi><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mstyle mathvariant="bold" mathsize="normal"><mi>Ω</mi></mstyle><mo stretchy="false">)</mo><mo>+</mo><msub><mi>σ</mi><mtext>T</mtext></msub><mi>ψ</mi><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mstyle mathvariant="bold" mathsize="normal"><mi>Ω</mi></mstyle><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><mrow><munder><mo>∫</mo><mstyle mathvariant="bold" mathsize="normal"><msup><mi>Ω</mi><mo>′</mo></msup></mstyle></munder><mrow><mi>ψ</mi><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mstyle mathvariant="bold" mathsize="normal"><msup><mi>Ω</mi><mo>′</mo></msup></mstyle><mo stretchy="false">)</mo><mtext> </mtext><msub><mi>σ</mi><mtext>S</mtext></msub><mo stretchy="false">(</mo><mstyle mathvariant="bold" mathsize="normal"><msup><mi>Ω</mi><mo>′</mo></msup></mstyle><mo>⋅</mo><mstyle mathvariant="bold" mathsize="normal"><mi>Ω</mi></mstyle><mo stretchy="false">)</mo><mtext> </mtext><mtext>d</mtext><mstyle mathvariant="bold" mathsize="normal"><msup><mi>Ω</mi><mo>′</mo></msup></mstyle></mrow></mrow></mstyle><mo>,</mo></mrow></math> where Ω′ and Ω are the unit vectors along the neutron velocity before and after a scattering collision, respectively. <math id="M2"><mrow><mi>c</mi><mo>=</mo><msub><mi>σ</mi><mtext>S</mtext></msub><mo>/</mo><msub><mi>σ</mi><mtext>T</mtext></msub></mrow></math> is the cross-section parameter know as the number of secondary neutrons per collision and σ T is the total macroscopic cross-section. ψ ( r , Ω ) is the neutron angular flux at position r traveling in direction Ω , and <math id="M3"><mrow><msub><mi>σ</mi><mtext>S</mtext></msub><mo stretchy="false">(</mo><mi>Ω</mi><mo>⋅</mo><msup><mi>Ω</mi><mo>′</mo></msup><mo stretchy="false">)</mo></mrow></math> is the scattering function [ 1 1 - 3 3 1-3 ]. Up to now, an appropriate scattering function representing the probabilities of the neutron interactions is found to be enough for the solutions of the problems in neutron transport theory because of its mathematical applicability. However, although using an approximate scattering function is usually accepted to be successful in many of the studies, the direct form of a scattering function like Henyey-Greenstein could be fascinating both in solution algorithm and in calculation of numerical results. On the other hand, the Henyey-Greenstein phase function is reported to be unsuccessful in some of the studies about radiative transfer [ 5 5 - 7 7 5-7 ]. Therefore, Anlı et al. constituted a new scattering kernel i.e., an AG phase function similar to Henyey-Greenstein formula and they first applied it to calculate the eigenvalue spectrum in the one-dimensional slab geometry transport equation [ 9 9 , 12 12 ]. In some recent studies, the AG phase function has been applied to diffusion equation and criticality problems in the transport theory using Legendre polynomials ( P N method) and Chebyshev polynomials of first kinds ( T N method) [ 10 10 , 11 11 ]. In this work, other than the previous studies about criticality calculations with defined scattering functions, the Chebyshev polynomials of second kind is preferred for use in the series expansion of the neutron angular flux. The AG phase function is used as the scattering function in the transport equation, and the critical half-thicknesses of the slab for one-speed neutrons are calculated for various values of the scattering parameters. Here, the detailed information about the P N method with the AG phase function does not needed to be given since it was mentioned in the previous work by the author [ 10 10 , 11 11 ]. The AG phase function, which is first used by Anlı et al. [ 9 9 ] in the determination of the eigenvalue spectrum, is then expressed as the scattering function, <math id="M4"><mrow><msub><mi>σ</mi><mtext>S</mtext></msub><mo stretchy="false">(</mo><mi>Ω</mi><mo>⋅</mo><msup><mi>Ω</mi><mo>′</mo></msup><mo stretchy="false">)</mo></mrow></math> , in neutron transport equation, i.e. Eq. (1), <math display="block" id="M5"><mrow><msub><mi>σ</mi><mtext>S</mtext></msub><msup><mrow/><mrow><mtext>AG</mtext></mrow></msup><mo stretchy="false">(</mo><msub><mi>μ</mi><mn>0</mn></msub><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mi>σ</mi><mtext>S</mtext></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><mo>,</mo></mrow></math> where σ S is any non-negative coefficient, t is the parameter representing all kinds of scattering (forward, backward and anisotropic, etc) and it is defined in the range of <math id="M6"><mrow><mo>−</mo><mn>1</mn><mo>≤</mo><mi>t</mi><mo>≤</mo><mn>1</mn></mrow></math> , and <math id="M7"><mrow><msub><mi>μ</mi><mn>0</mn></msub><mo>=</mo><mi>Ω</mi><mo>⋅</mo><msup><mi>Ω</mi><mo>′</mo></msup></mrow></math> is the cosine of the scattering angle [ 5 5 , 9 9 ], <math display="block" id="M8"><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> The neutron angular flux is used as in Ref. [ 13 13 ], <math display="block" id="M9"><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> When Eq. (2) is inserted on the right hand side of Eq. (1), the one-dimensional steady state transport equation can be written as, <math display="block" id="M10"><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><mtext> </mtext><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>=</mo><mtext> </mtext><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>d</mtext><msup><mi>μ</mi><mo>′</mo></msup><mtext> </mtext><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><mrow><mn>2</mn><mi>π</mi></mrow></munderover><mrow><mtext> </mtext><mfrac><mrow><msub><mi>σ</mi><mtext>S</mtext></msub></mrow><mrow><mn>4</mn><mi>π</mi><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><mtext>d</mtext><msup><mi>φ</mi><mo>′</mo></msup></mrow></mrow></mstyle></mrow></mrow></mstyle><mo>,</mo></mrow></math> with the free space boundary and symmetry conditions: <math display="block" id="M11"><mrow><mi>ψ</mi><mo stretchy="false">(</mo><mi>a</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math> <math display="block" id="M12"><mrow><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>=</mo><mi>ψ</mi><mo stretchy="false">(</mo><mo>−</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mi>μ</mi><mo>></mo><mn>0.</mn></mrow></math> The slab is assumed to be homogeneous expanding from x = - a to x = a . The integrand on the right hand side of Eq. (5) over <math id="M13"><mrow><mtext>d</mtext><msup><mi>ϕ</mi><mo>′</mo></msup></mrow></math> can be obtained using the addition theorem of the Legendre polynomials [ 9 9 ], <math display="block" id="M14"><mrow><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><mrow><mn>2</mn><mi>π</mi></mrow></munderover><mrow><mfrac><mrow><msub><mi>σ</mi><mtext>S</mtext></msub></mrow><mrow><mn>4</mn><mi>π</mi><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><mtext>d</mtext><msup><mi>φ</mi><mo>′</mo></msup><mo>=</mo><mfrac><mrow><msub><mi>σ</mi><mtext>S</mtext></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> Then, Eq. (5) can be rearranged using Eq. (7), <math display="block" id="M15"><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> A dimensionless space variable, such that σ T x / ν → x is defined in order to simplify the derivation of the equations and ν , is the eigenvalue. In the application procedure of the method, first the neutron angular flux, <math id="M16"><mrow><mi>ψ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo></mrow></math> , given in Eq. (4) is replaced in Eq. (8), and then the resultant equation is integrated over μ ∈ (-1,1) after multiplying it by U m ( μ ). The orthogonality and the recurrence relations of the Chebyshev polynomials of second kind are used during this procedure [ 13 13 - 15 15 13-15 ], <math display="block" id="M17"><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><mfrac><mi>π</mi><mn>2</mn></mfrac><mtext> </mtext><msub><mi>δ</mi><mrow><mi>n</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>,</mo></mrow></math> <math display="block" id="M18"><mrow><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><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><mn>0.</mn></mrow></math> By following the steps mentioned above, a general expression for the U N moments of the angular flux could not be reached in this study. However, individual expressions for n = 0,1,2,…,9 are obtained and they are; <math display="block" id="M19"><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>ν</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>ν</mi><mi>c</mi><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="M20"><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>ν</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>ν</mi><mi>c</mi><mi>t</mi><msub><mi>Φ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>,</mo></mrow></math> <math display="block" id="M21"><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>Φ</mi><mn>3</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>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>ν</mi><mtext> </mtext><msub><mi>Φ</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mo>−</mo><mn>2</mn><mi>ν</mi><mi>c</mi><mrow><mo>{</mo><mrow><mfrac><mn>1</mn><mrow><mn>15</mn></mrow></mfrac><mrow><mo>(</mo><mrow><msup><mi>t</mi><mn>2</mn></msup><mo>−</mo><mn>5</mn></mrow><mo>)</mo></mrow><msub><mi>Φ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>−</mo><mfrac><mrow><msup><mi>t</mi><mn>2</mn></msup></mrow><mn>5</mn></mfrac><msub><mi>Φ</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><mo>}</mo></mrow><mo>,</mo></mrow></math> <math display="block" id="M22"><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>Φ</mi><mrow><mn>10</mn></mrow></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>8</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>ν</mi><mtext> </mtext><msub><mi>Φ</mi><mn>9</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mo>−</mo><mn>2</mn><mi>ν</mi><mi>c</mi><mrow><mo>{</mo><mtable columnalign="left"><mtr><mtd><mrow><mo>(</mo><mrow><mfrac><mrow><mn>42</mn><msup><mi>t</mi><mn>9</mn></msup></mrow><mrow><mn>46189</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>112</mn><msup><mi>t</mi><mn>7</mn></msup></mrow><mrow><mn>109395</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mi>t</mi><mn>5</mn></msup></mrow><mrow><mn>429</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>40</mn><msup><mi>t</mi><mn>3</mn></msup></mrow><mrow><mn>3003</mn></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>5</mn><mi>t</mi></mrow><mrow><mn>99</mn></mrow></mfrac></mrow><mo>)</mo></mrow><msub><mi>Φ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mrow><mfrac><mrow><mn>48</mn><msup><mi>t</mi><mn>9</mn></msup></mrow><mrow><mn>20995</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>112</mn><msup><mi>t</mi><mn>7</mn></msup></mrow><mrow><mn>36465</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>28</mn><msup><mi>t</mi><mn>5</mn></msup></mrow><mrow><mn>2145</mn></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>100</mn><msup><mi>t</mi><mn>3</mn></msup></mrow><mrow><mn>3003</mn></mrow></mfrac></mrow><mo>)</mo></mrow><msub><mi>Φ</mi><mn>3</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mrow><mfrac><mrow><mn>9</mn><msup><mi>t</mi><mn>9</mn></msup></mrow><mrow><mn>1615</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>16</mn><msup><mi>t</mi><mn>7</mn></msup></mrow><mrow><mn>1105</mn></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>21</mn><msup><mi>t</mi><mn>5</mn></msup></mrow><mrow><mn>715</mn></mrow></mfrac></mrow><mo>)</mo></mrow><msub><mi>Φ</mi><mn>5</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mrow><mfrac><mrow><mn>8</mn><msup><mi>t</mi><mn>9</mn></msup></mrow><mrow><mn>323</mn></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>8</mn><msup><mi>t</mi><mn>7</mn></msup></mrow><mrow><mn>255</mn></mrow></mfrac></mrow><mo>)</mo></mrow><msub><mi>Φ</mi><mn>7</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>−</mo><mfrac><mrow><msup><mi>t</mi><mn>9</mn></msup></mrow><mrow><mn>19</mn></mrow></mfrac><msub><mi>Φ</mi><mn>9</mn></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mtd></mtr></mtable><mo>}</mo></mrow><mo>,</mo></mrow></math> where Φ -1 ( x ) = 0. A well-known solution [ 1 1 ], <math display="block" id="M23"><mrow><msub><mi>Φ</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>A</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mtext> </mtext><mi>exp</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>,</mo></mrow></math> is customarily used in Eqs. (11) in order to obtain analytic expressions of all A n ( ν )’s as follows, <math display="block" id="M24"><mrow><msub><mi>A</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><mn>2</mn><mi>ν</mi><msub><mi>A</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mn>2</mn><mi>ν</mi><mi>c</mi><msub><mi>A</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>,</mo></mrow></math> <math display="block" id="M25"><mrow><msub><mi>A</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>A</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><mn>2</mn><mi>ν</mi><msub><mi>A</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>2</mn><mn>3</mn></mfrac><mi>ν</mi><mi>c</mi><mi>t</mi><msub><mi>A</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>,</mo></mrow></math> <math display="block" id="M26"><mrow><mfrac><mrow><mn>2</mn><mi>ν</mi><mi>c</mi></mrow><mrow><mn>15</mn></mrow></mfrac><mrow><mo>(</mo><mrow><msup><mi>t</mi><mn>2</mn></msup><mo>−</mo><mn>5</mn></mrow><mo>)</mo></mrow><msub><mi>A</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>A</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><mn>2</mn><mi>ν</mi><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><mfrac><mrow><mi>c</mi><msup><mi>t</mi><mn>2</mn></msup></mrow><mn>5</mn></mfrac></mrow><mo>)</mo></mrow><msub><mi>A</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>A</mi><mn>3</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math> <math display="block" id="M27"><mtable columnalign="left"><mtr><mtd><mn>2</mn><mi>ν</mi><mi>c</mi><mrow><mo>(</mo><mrow><mfrac><mrow><mn>42</mn><msup><mi>t</mi><mn>9</mn></msup></mrow><mrow><mn>46189</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>112</mn><msup><mi>t</mi><mn>7</mn></msup></mrow><mrow><mn>109395</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mi>t</mi><mn>5</mn></msup></mrow><mrow><mn>429</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>40</mn><msup><mi>t</mi><mn>3</mn></msup></mrow><mrow><mn>3003</mn></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>5</mn><mi>t</mi></mrow><mrow><mn>99</mn></mrow></mfrac></mrow><mo>)</mo></mrow><msub><mi>A</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>+</mo><mn>2</mn><mi>ν</mi><mi>c</mi><mrow><mo>(</mo><mrow><mfrac><mrow><mn>48</mn><msup><mi>t</mi><mn>9</mn></msup></mrow><mrow><mn>20995</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>112</mn><msup><mi>t</mi><mn>7</mn></msup></mrow><mrow><mn>36465</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>28</mn><msup><mi>t</mi><mn>5</mn></msup></mrow><mrow><mn>2145</mn></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>100</mn><msup><mi>t</mi><mn>3</mn></msup></mrow><mrow><mn>3003</mn></mrow></mfrac></mrow><mo>)</mo></mrow><msub><mi>A</mi><mn>3</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>+</mo><mn>2</mn><mi>ν</mi><mi>c</mi><mrow><mo>(</mo><mrow><mfrac><mrow><mn>9</mn><msup><mi>t</mi><mn>9</mn></msup></mrow><mrow><mn>1615</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>16</mn><msup><mi>t</mi><mn>7</mn></msup></mrow><mrow><mn>1105</mn></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>21</mn><msup><mi>t</mi><mn>5</mn></msup></mrow><mrow><mn>715</mn></mrow></mfrac></mrow><mo>)</mo></mrow><msub><mi>A</mi><mn>5</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>+</mo><mn>2</mn><mi>ν</mi><mi>c</mi><mrow><mo>(</mo><mrow><mfrac><mrow><mn>8</mn><msup><mi>t</mi><mn>9</mn></msup></mrow><mrow><mn>323</mn></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>8</mn><msup><mi>t</mi><mn>7</mn></msup></mrow><mrow><mn>255</mn></mrow></mfrac></mrow><mo>)</mo></mrow><msub><mi>A</mi><mn>7</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>A</mi><mn>8</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><mn>2</mn><mi>ν</mi><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><mfrac><mrow><mi>c</mi><msup><mi>t</mi><mn>9</mn></msup></mrow><mrow><mn>19</mn></mrow></mfrac></mrow><mo>)</mo></mrow><msub><mi>A</mi><mn>9</mn></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>+</mo><msub><mi>A</mi><mrow><mn>10</mn></mrow></msub><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mn>0</mn><mo>,</mo></mtd></mtr></mtable></math> where A -1 ( ν , t ) = 0 and A 0 ( ν , t ) = 1. Eqs. (13) can also be written in a matrix form for an alternative solution algorithm, <math display="block" id="M28"><mrow><mo stretchy="false">[</mo><mi>M</mi><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mi>A</mi><mo stretchy="false">(</mo><mi>ν</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math> where M ( ν , t ) is ( N + 1) × ( N + 1) coefficient matrix and A ( ν , t ) = [ A 0 , A 1 ,…, A N ] T . It is possible to obtain non-trivial solutions for the discrete eigenvalues by equating the determinant of the coefficient matrix to zero, i.e. det[M( ν , t )] = 0. As well known in the P N approximation, the contribution of the ( N +1) th term to the neutron flux could be accepted as negligible. In addition, the Legendre and Chebyshev polynomials are the members of the Jacobi polynomials. Then by following the same procedure derived for the P N approximation, the discrete and continuum ν eigenvalues can be obtained by setting A N +1 ( ν , t ) = 0 for various values of c and t . A brief information about the profiles of the eigenvalues can be given as following: All roots of A N +1 ( ν ) = 0 are identical with the roots of U N +1 ( ν ) = 0 in the case of c = 0 and all values of t . When c = 1, one pair of the roots is <math id="M29"><mrow><mo>±</mo><mi>∞</mi><mi>i</mi></mrow></math> and the others are real lying in the interval [-1,1]. When 0 < c < 1, all roots are real and one pair of them is usually greater than 1. Finally when c > 1, one pair of the roots is purely imaginary and the others are real [ 9 9 , 19 19 ]. As an example for U 1 approximation, two linear algebraic equations, i.e. Eq. (13a) for n = 0 and Eq. (13b) for n = 1, are obtained and an analytic expression for the eigenvalues can easily be derived by setting A 2 ( ν , t ) = 0 in Eqs. (13a) and (13b), <math display="block" id="M30"><mrow><msub><mi>ν</mi><mi>k</mi></msub><mo>=</mo><mo>±</mo><mfrac><mn>1</mn><mn>2</mn></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> In other words, the same eigenvalues can be obtained using Eq. (14) by deriving a 2×2 matrix equation and then equating the determinant of the coefficient matrix to zero. Therefore, the general solution of the flux moments for odd numbers of N can be written with so computed discrete eigenvalues ν k for k = 1, …, N + 1, <math display="block" id="M31"><mrow><msub><mi>Φ</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mo stretchy="false">(</mo><mi>N</mi><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo><mo>/</mo><mn>2</mn></mrow></munderover><mrow><msub><mi>α</mi><mi>k</mi></msub><msub><mi>A</mi><mi>n</mi></msub></mrow></mstyle><mo stretchy="false">(</mo><msub><mi>ν</mi><mi>k</mi></msub><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mrow><mo>[</mo><mrow><mi>exp</mi><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>σ</mi><mtext>T</mtext></msub><mi>x</mi></mrow><mrow><msub><mi>v</mi><mi>k</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><msup><mrow><mo stretchy="false">(</mo><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow><mi>n</mi></msup><mi>exp</mi><mrow><mo>(</mo><mrow><mo>−</mo><mfrac><mrow><msub><mi>σ</mi><mtext>T</mtext></msub><mi>x</mi></mrow><mrow><msub><mi>ν</mi><mi>k</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow><mo>,</mo></mrow></math> where A n (- ν , t ) = (-1) n A n ( ν , t ) and α k ’s are the coefficients which can be determined from the physical boundary conditions of the system. The general solution Eq. (16) is constituted as the summation of all eigenvectors corresponding to each eigenvalues. 3 Boundary conditions and the critical system 1 The study of calculation of the eigenvalues of the problem representing the system under consideration can be said to be equivalent to find the critical size of that system. The values of the number of secondary neutrons per collision are very important to operate the reactor safely and decide whether it is critical or not. In a reactor, for each absorption collision the reactor cannot said to be critical when fewer neutrons are emitted than absorbed ( c < 1). However, a reactor may be subcritical or critical for a slab of finite thickness with c > 1 [ 3 3 ]. The angular neutron flux is continuous across material region boundaries except for the direction μ = 0 in slab geometries. Any finite sum of the Legendre polynomials is continuous over the range -1 ⪯ μ ⪯ 1 and, therefore, continuous at μ = 0. Then, the P N approximation in slab geometries is a rather poor representation of the angular flux near material boundaries. Although the Mark and Marshak boundary conditions are the most commonly used ones for the criticality problems, the Marshak boundary condition which is based on the condition of zero incoming current at the vacuum boundary is somewhat more accurate than the Mark condition, at least for small N [ 3 3 , 16 16 ]. Because of the reasons mentioned above, for the criticality of the slab, the Marshak boundary condition for U N approximation of odd order can be written as, <math display="block" id="M32"><mrow><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><mn>1</mn></munderover><mrow><mtext> </mtext><mi>ψ</mi><mo stretchy="false">(</mo><mi>a</mi><mo>,</mo><mo>−</mo><mi>μ</mi><mo stretchy="false">)</mo><mtext> </mtext><mtext> </mtext><msub><mi>U</mi><mi>k</mi></msub><mo stretchy="false">(</mo><mo>−</mo><mi>μ</mi><mo stretchy="false">)</mo><mtext> </mtext></mrow></mrow></mstyle><mtext>d</mtext><mi>μ</mi><mo>=</mo><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><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mtext> </mtext><mn>3</mn><mo>,</mo><mn>....</mn><mo>,</mo><mi>N</mi><mo>.</mo></mrow></math> In the application procedure of the boundary condition, first the neutron flux in Eq. (16) is replaced in Eq. (4) and then, the resulting equation is inserted into Eq. (17) with the parity relation of the Chebyshev polynomials of second kind U k (- μ ) = (-1) k U k ( μ ), <math display="block" id="M33"><mrow><mfrac><mn>2</mn><mi>π</mi></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><mrow><msup><mrow><mo stretchy="false">(</mo><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow><mrow><mi>n</mi><mo>+</mo><mi>k</mi></mrow></msup><mrow><mo>{</mo><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mo stretchy="false">(</mo><mi>N</mi><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo><mo>/</mo><mn>2</mn></mrow></munderover><mrow><msub><mi>α</mi><mi>k</mi></msub><msub><mi>A</mi><mi>n</mi></msub><mo stretchy="false">(</mo><msub><mi>ν</mi><mi>k</mi></msub><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mtext> </mtext><mrow><mo>[</mo><mrow><mi>exp</mi><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>σ</mi><mtext>T</mtext></msub><mi>a</mi></mrow><mrow><msub><mi>ν</mi><mi>k</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><msup><mrow><mo stretchy="false">(</mo><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow><mi>n</mi></msup><mi>exp</mi><mrow><mo>(</mo><mrow><mo>−</mo><mfrac><mrow><msub><mi>σ</mi><mtext>T</mtext></msub><mi>a</mi></mrow><mrow><msub><mi>ν</mi><mi>k</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></mstyle></mrow><mo>}</mo></mrow></mrow></mstyle><mtext> </mtext><msub><mi>I</mi><mrow><mi>n</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math> where I n,k is, <math display="block" id="M34"><mrow><msub><mi>I</mi><mrow><mi>n</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><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><mtext> </mtext><msub><mi>U</mi><mi>k</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> </mtext><mtext>d</mtext><mi>μ</mi><mo>=</mo><mrow><mo>{</mo><mtable columnalign="left"><mtr><mtd><mi>π</mi><mo>/</mo><mn>4</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><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><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><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><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><mtext> </mtext><mtext> </mtext><mtext> </mtext><mi>n</mi><mo>=</mo><mtext> </mtext><mi>k</mi><mo>,</mo></mtd></mtr><mtr><mtd><mfrac><mrow><mi>sin</mi><mrow><mo>[</mo><mrow><mo stretchy="false">(</mo><mi>n</mi><mo>−</mo><mi>k</mi><mo stretchy="false">)</mo><mi>π</mi><mo>/</mo><mn>2</mn></mrow><mo>]</mo></mrow></mrow><mrow><mn>2</mn><mo stretchy="false">(</mo><mi>n</mi><mo>−</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></mfrac><mtext> </mtext><mtext> </mtext><mo>+</mo><mfrac><mrow><mi>sin</mi><mrow><mo>[</mo><mrow><mo stretchy="false">(</mo><mi>n</mi><mo>+</mo><mi>k</mi><mo stretchy="false">)</mo><mi>π</mi><mo>/</mo><mn>2</mn></mrow><mo>]</mo></mrow></mrow><mrow><mn>2</mn><mo stretchy="false">(</mo><mi>n</mi><mo>+</mo><mi>k</mi><mo>+</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mfrac><mtext> </mtext><mtext> </mtext><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mi>n</mi><mo>≠</mo><mi>k</mi><mo>,</mo><mtext> </mtext></mtd></mtr></mtable></mrow></mrow></math> and <math display="block" id="M35"><mrow><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><mn>1</mn></munderover><mrow><msub><mi>U</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mo>−</mo><mi>μ</mi><mo stretchy="false">)</mo><mtext> </mtext></mrow></mrow></mstyle><msub><mi>U</mi><mi>k</mi></msub><mo stretchy="false">(</mo><mo>−</mo><mi>μ</mi><mo stretchy="false">)</mo><msqrt><mrow><mn>1</mn><mo>−</mo><msup><mi>μ</mi><mn>2</mn></msup></mrow></msqrt><mtext> </mtext><mtext>d</mtext><mi>μ</mi><mo>=</mo><msup><mrow><mo stretchy="false">(</mo><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow><mrow><mi>n</mi><mo>+</mo><mi>k</mi></mrow></msup><msub><mi>I</mi><mrow><mi>n</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>.</mo></mrow></math> Eq. (18) is referred to as the criticality condition and it can also be written in a matrix form, <math display="block" id="M36"><mrow><mrow><mo>[</mo><mrow><msubsup><mi>M</mi><mi>m</mi><mi>k</mi></msubsup><mo stretchy="false">(</mo><mi>a</mi><mo stretchy="false">)</mo></mrow><mo>]</mo></mrow><mtext> </mtext><msub><mi>β</mi><mi>k</mi></msub><mo>=</mo><mn>0</mn><mo>,</mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mi>m</mi><mo>,</mo><mi>k</mi><mo>=</mo><mtext> </mtext><mn>1</mn><mo>,</mo><mtext> </mtext><mn>2</mn><mo>,</mo><mtext> </mtext><mo>…</mo><mo>,</mo><mtext> </mtext><mrow><mo>(</mo><mrow><mi>N</mi><mo>+</mo><mtext> </mtext><mn>1</mn></mrow><mo>)</mo></mrow><mo>/</mo><mn>2</mn><mo>,</mo></mrow></math> where β k is the column vector comprising elements of <math id="M37"><mrow><msup><mrow><mrow><mo>[</mo><mrow><msub><mi>β</mi><mn>1</mn></msub><mo>,</mo><msub><mi>β</mi><mn>2</mn></msub><mo>,</mo><mn>...</mn><mo>,</mo><msub><mi>β</mi><mrow><mo stretchy="false">(</mo><mi>N</mi><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo><mo>/</mo><mn>2</mn></mrow></msub></mrow><mo>]</mo></mrow></mrow><mi>T</mi></msup></mrow></math> and <math id="M38"><mrow><msubsup><mi>M</mi><mi>m</mi><mi>k</mi></msubsup><mo stretchy="false">(</mo><mi>a</mi><mo stretchy="false">)</mo></mrow></math> is the coefficient matrix with ( N + 1)/2 × ( N + 1)/2 elements. Eq. (18) or (21) can be solved for a non-trivial solution when the coefficients β k ’s are nonzero or the determinant of the coefficient matrix is zero, i.e. <math id="M39"><mrow><mi>det</mi><mrow><mo>[</mo><mrow><msubsup><mi>M</mi><mi>m</mi><mi>k</mi></msubsup><mo stretchy="false">(</mo><mi>a</mi><mo stretchy="false">)</mo></mrow><mo>]</mo></mrow><mo>=</mo><mn>0</mn></mrow></math> . Since the eigenvalues were already calculated from Eq. (15), as final application by letting N = 1 in Eq. (18) or Eq. (21), an analytic solution for the critical half-thickness of the slab can easily be obtained for U 1 approximation, <math display="block" id="M40"><mrow><mi>a</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msub><mi>σ</mi><mtext>T</mtext></msub></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><mtext> </mtext><msup><mrow><mi>tanh</mi></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mtext> </mtext><mrow><mo>(</mo><mrow><mo>−</mo><mfrac><mn>8</mn><mrow><mn>3</mn><mi>π</mi></mrow></mfrac><msqrt><mrow><mfrac><mrow><mo stretchy="false">(</mo><mn>3</mn><mo>−</mo><mi>c</mi><mi>t</mi><mo stretchy="false">)</mo><mtext> </mtext></mrow><mrow><mn>3</mn><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>c</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></msqrt></mrow><mo>)</mo></mrow><mo>.</mo></mrow></math> 4 Numerical results 1 An application of the Chebyshev polynomials expansion ( U N approximation) is done for the critical slab problem for one-speed neutrons in a uniform homogeneous medium. In the method, the Chebyshev polynomials of second kind are used in the angular part of the neutron flux as it has been successfully applied before [ 13 13 , 15 15 ]. Contrary to previous approximation scattering functions, in order to get closer to accurate solution of the transport equation, the AG phase function is used. Various values of c and t are used in Eq. (13) or (14) to compute the discrete eigenvalues by setting A N +1 ( ν ) = 0. An analytic expression of the eigenvalues for U 1 approximation is obtained and it is given in Eq. (15). Various orders of the U N approximation with the AG phase function are applied to Eq. (18) or (21) and the numerical results for the critical half-thickness of the slab are tabulated in Tables 1 Tables 1 through 4 4 . A final table i.e. Table 5 Table 5 , has been needed to compare the results obtained from the present method with the ones already obtained from the conventional P N method in a previous study [ 11 11 ]. The Marshak boundary condition is used during the application of the criticality condition to the problem since it is accepted as to be more valid than the Mark for small N [ 3 3 , 16 16 ]. All calculations are carried out by means of the Maple software and the total macroscopic cross section is taken as to be its normalized value, σ T = 1 cm -1 . 1 1 t 2 1 c = 1.01 2 1 c = 1.20 2 1 c = 2.00 1 1 1 1 P 9 [11] 1 1 U 9 (present work) 1 1 P 9 [11] 1 1 U 9 (present work) 1 1 P 9 [11] 1 1 U 9 (present work) 1 1 -1 1 1 7.27723 1 1 7.27369 1 1 1.13840 1 1 1.13597 1 1 0.27052 1 1 0.27071 1 1 -4/5 1 1 7.45782 1 1 7.45448 1 1 1.17127 1 1 1.16891 1 1 0.28491 1 1 0.28463 1 1 -3/4 1 1 7.50499 1 1 7.50195 1 1 1.17906 1 1 1.17679 1 1 0.28775 1 1 0.28734 1 1 -2/3 1 1 7.58523 1 1 7.58252 1 1 1.19179 1 1 1.18972 1 1 0.29202 1 1 0.29142 1 1 -1/2 1 1 7.75318 1 1 7.75102 1 1 1.21664 1 1 1.21505 1 1 0.29922 1 1 0.29832 1 1 -3/10 1 1 7.96959 1 1 7.96833 1 1 1.24599 1 1 1.24496 1 1 0.30618 1 1 0.30503 1 1 -1/4 1 1 8.02650 1 1 8.02529 1 1 1.25333 1 1 1.25241 1 1 0.30769 1 1 0.30651 1 1 -1/5 1 1 8.08440 1 1 8.08387 1 1 1.26068 1 1 1.25986 1 1 0.30913 1 1 0.30792 1 1 0 1 1 8.33040 1 1 8.32957 1 1 1.29038 1 1 1.28974 1 1 0.31418 1 1 0.31294 1 1 1/5 1 1 8.60028 1 1 8.59894 1 1 1.32081 1 1 1.31998 1 1 0.31814 1 1 0.31705 1 1 1/4 1 1 8.67076 1 1 8.66762 1 1 1.32856 1 1 1.32763 1 1 0.31896 1 1 0.31793 1 1 3/10 1 1 8.74482 1 1 8.74048 1 1 1.33638 1 1 1.33532 1 1 0.31971 1 1 0.31175 1 1 1/2 1 1 9.05877 1 1 9.05398 1 1 1.36833 1 1 1.36662 1 1 0.32188 1 1 0.32135 1 1 2/3 1 1 9.34822 1 1 9.34066 1 1 1.39580 1 1 1.39347 1 1 0.32251 1 1 0.32249 1 1 3/4 1 1 9.50423 1 1 9.49957 1 1 1.40981 1 1 1.40719 1 1 0.32232 1 1 0.32260 1 1 4/5 1 1 9.60097 1 1 9.59680 1 1 1.41830 1 1 1.41553 1 1 0.32201 1 1 0.32249 1 1 1 1 1 10.02130 1 1 10.01662 1 1 1.45277 1 1 1.44965 1 1 0.31913 1 1 0.32021 In Tables 1 Tables 1 ， 2 2 , 3 3 , 4 4 , the critical half-thicknesses of the slab are listed for c = 1.01, 1.02, 1.20, and 2.00 and t is selected with increasing order from -1 to 1. One can easily test that the t = 0 case corresponds to isotropic scattering [ 9 9 , 13 13 , 15 15 ]. In other words, by replacing the case of t = 0 in Eqs. (15-18) one would obtain the equations necessary for calculating the eigenvalues and the critical half-thicknesses using the U N method in slab geometry for isotropic scattering [ 13 13 , 15 15 ]. It can be seen from the tables that in many cases, the critical half-thickness of the slab increases while t is increasing and c is decreasing. It was reported that the critical thickness of the slab can behave non-monotonic when neutrons tend to propagate in the forward direction. This is observed as first an increasing trend and then a decreasing trend with increasing forward anisotropy parameter according to the choice of the c . This behavior of the critical thickness is referred to as non-monotonic and it is observed in this study for t approaching to 1 and c = 2, especially in the case of higher order approximations with N > 5. That means the same non-monotonic behavior of the critical thickness as reported before [ 17 17 , 18 18 ] appears when the neutrons scatter in the forward peaked directions. This circumstance can be followed in Tables 4 Tables 4 and 5 5 by examining the values with t > ¼. However, since the criticality calculations are important especially for c near to 1, this anomaly for the AG phase function can be thought as negligible. It can be seen from these tables that this non-monotonic behavior has occurred when c = 2.00 (away from 1) and higher order approximation is used ( N = 9) which is pointed to as the advantage of the Marshak boundary condition against the Mark [ 3 3 ]. More discrepancies about the behavior of the critical thickness are reported as to be seen when using Henyey-Greenstein phase function [ 11 11 ]. 5 Conclusion 1 In this paper, the critical thickness of one-speed neutrons in a finite homogeneous slab is studied using U N approximation which is applied successfully in preceding studies [ 13 13 , 15 15 ]. As a second important application of this study, the AG phase function is chosen as the scattering kernel of the transport equation. The critical half-thicknesses of the slab are calculated numerically using increasing orders of the U N approximation up to N = 9 for both positive and negative values of the parameter t . While the positive values of t represent the forward peaked scattering of the neutrons, the negative values of it represent the backward peaked scattering of the neutrons. These are physically possible situations presented in a reactor. When a neutron interacts with a particle having a mass approximately equal to the mass of the interacting neutron, such as a hydrogen nucleus in a moderator, this interaction has a probability to end with a forward scattering. In a similar way, when a neutron interacts with a particle having a mass of greater than the mass of the interacting neutron, such as an Oxygen nucleus in a moderator, a nucleus in reactor material, or a daughter nucleus emitted from a fission reaction, this interaction has a probability to end with a backward scattering [ 9 9 ]. Therefore, this study can be evaluated as the calculation of the critical half-thickness of the slab for forward and backward scattering since both the positive and negative values of t are taken in Tables 1 Tables 1 , , 3 3 , 4 4 , 5 5 . In summary, one can easily assert from the derivation of the equations and the results already obtained here that the AG phase function is seen to be convenient for the solution of the problems in transport theory. Furthermore, the AG phase function with its easily applicable derivation can also be sufficient for other problems containing a phase function in particle or photon transport and in science and engineering. References B. Davison , Neutron transport theory . 1 st edn. ( London , Oxford University Press , 1958 ), pp. 116 - 144 Neutron transport theory K.M. Case , P.F. Zweifel , Linear transport theory , 1 st edn. ( Addison-Wesley Publishing Company , 1967 ), pp. 10 - 45 . Linear transport theory G.I. Bell , S. Glasstone , Nuclear reactor theory , 1 st edn. ( New York , VNR Company , 1972 ), pp. 33 - 72 . Nuclear reactor theory 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 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 L.O. Reynolds , N.J. McCormic , Approximate two parameter phase function for light scattering , J. Optical Soc. America ., 70 , 1206 - 1212 ( 1980 ). doi: 10.1364/JOSA.70.001206 http://doi.org/10.1364/JOSA.70.001206 Approximate two parameter phase function for light scattering M.M.R. Williams , A synthetic scattering kernel for particle transport in absorbing media with anisotropic scattering , J. Phys. D: Appl. Phys . 11 , 2455 - 2463 ( 1978 ). doi: 10.1088/0022-3727/11/18/004 http://doi.org/10.1088/0022-3727/11/18/004 A synthetic scattering kernel for particle transport in absorbing media with anisotropic scattering P. Liu , A new phase function approximating to mie scattering for radiative transport equations , Phys. Med. and Bio . 39 , 1025 - 1036 ( 1994 ). doi: 10.1088/0031-9155/39/6/008 http://doi.org/10.1088/0031-9155/39/6/008 A new phase function approximating to mie scattering for radiative transport equations F. Anlı , F. Yaşa , 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. Anlı , 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 F. Anlı , S. Güngör , Some useful properties of Legendre polynomials and its applications to neutron transport equation in slab geometry , Appl. Math. Mod . 31 , 727 - 733 ( 2007 ). doi: 10.1016/j.apm.2005.12.005 http://doi.org/10.1016/j.apm.2005.12.005 Some useful properties of Legendre polynomials and its applications to neutron transport equation in slab geometry 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 G. Arfken , Mathematical methods for physicists . 3 rd edn. ( London , Academic Press, Inc ., 1985 ), pp. 817 - 879 Mathematical methods for physicists H. Öztürk , F. Anlı , 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 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 C. Yıldız , Variation of the critical slab thickness with the degree of strongly anisotropic scattering in one-speed neutron transport theory , Ann. Nucl. Energy , 25 , 529 - 540 ( 1998 ). doi: 10.1016/S0306-4549(97)00114-X http://doi.org/10.1016/S0306-4549(97)00114-X Variation of the critical slab thickness with the degree of strongly anisotropic scattering in one-speed neutron transport theory D.C. Sahni , N.G. Sjöstrand , Non-monotonic variation of the criticality factor with the degree of anisotropy in one-speed neutron transport , Trans. Theory Statis. Phys . 20 , 339 - 349 ( 1991 ). Non-monotonic variation of the criticality factor with the degree of anisotropy in one-speed neutron transport F. Yaşa , F. Anlı , S. Güngör , Eigenvalue spectrum with Chebyshev polynomial approximation of the transport equation in slab geometry , J. Quant. Spectros. Rad. Transf . 97 , 51 - 57 ( 2006 ). doi: 10.1016/j.jqsrt.2004.12.017 http://doi.org/10.1016/j.jqsrt.2004.12.017 Eigenvalue spectrum with Chebyshev polynomial approximation of the transport equation in slab geometry 10.1007/s41365-019-0552-z 2018-05-15 2018-08-14 2018-09-01 2019-01-24 2019-02-01 2021-04-10 © 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. 2019 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-019-0552-z http://dx.doi.org/10.1007/s41365-019-0552-z . 2019 Nuclear Science and Techniques 02 30