Section outline

  • Applied Reactor Physics
  • General information
  • Matlab scripts for Chapter 2
    • Return the Ferrer definition of the associated Legendre function
    • Real spherical harmonics as polynomials in direction cosines
    • Spherical harmonics difference coefficients (Alm to Hlm)
    • Doppler chi and psi functions

  • Matlab scripts for Chapter 3
    • Scripts related to the spherical harmonics method
    • Coefficients used with the collision probability method in 1D slab geometry
    • Coefficients used with the collision probability method in 1D cylindrical geometry
    • Scripts related to the collision probability method in 2D square pincell geometry
    • Scripts related to the discrete ordinates method
    • Scripts related to the method of characteristics (MOC)
    • Monte Carlo method and random walk of a neutron
  • Matlab scripts for Chapter 4
    • Dancoff factor calculation in an array of cylindrical pins
    • Set the coefficients of the Chebyshev Rational Approximation Method (CRAM) for exponentials

  • Matlab scripts for Chapter 5
    • Compute the coupling coefficients for the 1D analytic nodal method
    • Compute the coupling coefficients for the 2D analytic nodal method
    • Find the roots of the inhour equation
    • Solution of a non-linear equation using the Newton method
    • Half-interval (aka, bissection) method
  • Matlab scripts for Appendix A
    • Produce a Gauss tracking in 1D slab geometry

    • Produce a Gauss-Jacobi tracking in 1D curvilinear geometry
    • Compute the p11 component in 1D cylindrical geometry
    • Compute the p11 component in 1D spherical geometry
    • Produce a Gauss-Jacobi tracking in 2D pincell geometry
  • Matlab scripts for Appendix B
    • Exponential function of order n
    • Bickley function of order n
  • Matlab scripts for Appendix C
    • Solution of a tri-diagonal linear system using Gauss elimination
    • Cholesky factorization and linear system solution
    • QR factorization
    • Iterative solution of a linear system
    • Eigenvalue/eigenvector calculation
    • Find the solution of the fixed source eigenvalue equation (a-eval*b)*gamma=s using the inverse power method