NEA-0005 GAKER-KIRA.

1. NAME OR DESIGNATION OF PROGRAM:  GAKER-KIRA.

3. NATURE OF PHYSICAL PROBLEM SOLVED

Evaluation of the energy transfer matrix and its first, second and third Legendre components for protons bound in water or polyethylene, and for deuterons bound in heavy water. Further output, the isotropic scattering matrix with diagonal elements corrected for first order anisotropy, scattering and transport cross sections, nth moment of energy transfer mn(e) (n=1,2,3), and Maxwellian averages of these moments, of the transport mean free path and of the scattering cross sections.

4. METHOD OF SOLUTION

For down-scattering the elements of the matrix are calculated directly from the analytical formulae (BNL 5826, p 69 FF (1961)). For scattering in the same speed (energy) group, a linear average is taken of four terms, scattering from the typical speed point of that group to four equidistant points in that same interval. Up-scattering elements are computed from detailed balance. The Bessel functions that appear are calculated following an algorithm due to Stegun and Abramowitz (see references). The anisotropic version of GAKER, which is due to Koppel and Young (see  references) is also included in GAKER-KIRA, and the integration over the orientation angles is done by a 5-points Gauss quadrature.

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM:  A maximum number of 55 speed groups is permitted.

6. TYPICAL RUNNING TIME

Approximately 8 to 9 elements are calculated per second in the isotropic version, the anisotropic version is five times slower.

7. UNUSUAL FEATURES OF THE PROGRAM:  Inclusion of the anisotropic version.

8. RELATED AND AUXILIARY PROGRAMS:  GAKER.

10. REFERENCES

- H.C. Honeck:
  BNL 5826, p 69 FF (1961).
- H.C. Honeck:
  Trans.Am.Nucl.Soc.5,47 (1962).
- D.T. Goldman and F.D. Federighi:
  BNL 719, 1,100 (1962).
- Irene A. Stegun and M. Abramowitz:
  Math. Tables Aids Comput. 11, 255 (1957).
- J.U. Koppel and J.A. Young:
  Nucl. Sci. Eng 19, 412 (1964).

11. MACHINE REQUIREMENTS:  Standard CDC 3600.

13. OPERATING SYSTEM OR MONITOR UNDER WHICH PROGRAM IS EXECUTED:  FORTRAN 3600 under SCOPE 5.1.

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

ANY OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

15. NAME AND ESTABLISHMENT OF AUTHOR

     R.J.J. Stammler
     Institutt for Atomenergi
     P.O. Box 40,
     Kjeller
     Norway.

File name	File description	Records
NEA0005_01.001	INFORMATION	4
NEA0005_01.002	SOURCE PROGRAM (F4) EBCDIC	409
NEA0005_01.003	SAMPLE PROBLEM INPUT DATA	19
NEA0005_01.004	SAMPLE PROBLEM PRINTED OUTPUT	555

File name	File description	Records
NEA0005_02.001	SOURCE & DATA 360/65	429
NEA0005_02.002	OUTPUT	360