Analytical SN solutions in heterogeneous slabs using symbolic algebra computer programs

A modern symbolic algebra computer program, MAPLE, is used to compute solutions to the well-known analytical discrete ordinates, or SN, solutions in one-dimensional, slab geometry. Symbolic algebra programs compute the solutions with arbitrary precision and are free of spatial discretization error so they can be used to investigate new discretizations for one-dimensional slab, geometry SN methods. Pointwise scalar flux solutions are computed for several sample calculations of interest. Sample MAPLE command scripts are provided to illustrate how easily the theory can be translated into a working solution and serve as a complete tool capable of computing analytical SN solutions for mono-energetic, one-dimensional transport problems.

[1]  M. Yavuz.  Methods for computing SN eigenvalues and eigenvectors of slab geometry transport problems , 1998 .

[2]  E. Inönü,et al.  Orthogonality of a Set of Polynomials Encountered in Neutron‐Transport and Radiative‐Transfer Theories , 1970 .

[3]  J. P Hennart,et al.  A generalized nodal finite element formalism for discrete ordinates equations in slab geometry Part I: Theory in the continuous moment case , 1995 .

[4]  Edward W. Larsen.  On Numerical Solutions of Transport Problems in the Diffusion Limit , 1983 .

[5]  J. P. Hennart,et al.  A generalized nodal finite element formalism for discrete ordinates equations in slab geometry Part III: Numerical results , 1993 .

[6]  M. Yavuz.  A one-d simplified discrete-ordinates method with no spatial truncation error , 1995 .

[7]  G. C. Pomraning,et al.  Linear Transport Theory , 1967 .

[8]  Norman J. McCormick,et al.  A Review of Neutron Transport Approximations , 1982 .

[9]  E. Richard Cohen,et al.  Neutron Transport Theory , 1959 .

[10]  William H. Reed,et al.  New Difference Schemes for the Neutron Transport Equation , 1971 .

[11]  C. E. Siewert,et al.  A high-order spherical harmonics solution to the standard problem in radiative transfer , 1984 .

[12]  I. Chen,et al.  Discrete Eigenvalues of the Neutron Transport Equation with Anisotropic Scattering , 1979 .

[13]  Kenneth M. Case,et al.  Elementary solutions of the transport equation and their applications , 1960 .

[14]  C. Siewert,et al.  On computing the Chandrasekhar polynomials in high order and high degree , 1990 .

[15]  A. M. Voloschenko,et al.  Adaptive positive nodal method for transport equation , 1994 .

[16]  J. Mika,et al.  Neutron Transport with Anisotropic Scattering , 1961 .

[17]  H. Keller,et al.  Approximate Solutions of Transport Problems. II. Convergence and Applications of the Discrete Ordinate Method , 1960 .

[18]  James S. Warsa,et al.  p-adaptive numerical methods for particle transport , 1998 .

[19]  Marcos Pimenta de Abreu,et al.  On the spectrum of the one-speed slab-geometry discrete ordinates operator in neutron transport theory , 1998 .

[20]  James J. Duderstadt,et al.  Phase-space finite element methods applied to the first-order form of the transport equation , 1981 .

[21]  Bryan Minor,et al.  Exponential Characteristic Spatial Quadrature for Discrete Ordinates Radiation Transport in Slab Geometry , 1994 .