last modified: 01-MAR-1965 | catalog | categories | new | search |

NEA-0005 GAKER-KIRA.

GAKER-KIRA, Energy Transfer of Protons in H2O or Polyethylene and Deuterons in D2O

top ]
1. NAME OR DESIGNATION OF PROGRAM:  GAKER-KIRA.
top ]
2. COMPUTERS
To submit a request, click below on the link of the version you wish to order. Only liaison officers are authorised to submit online requests. Rules for requesters are available here.
Program name Package id Status Status date
GAKER-KIRA NEA-0005/01 Tested 01-APR-1976
GAKER-KIRA NEA-0005/02 Tested 01-MAR-1965

Machines used:

Package ID Orig. computer Test computer
NEA-0005/01 IBM 370 series IBM 370 series
NEA-0005/02 CDC 3600 CDC 3600
top ]
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.
top ]
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.
top ]
5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM:  A maximum number of 55 speed groups is permitted.
top ]
6. TYPICAL RUNNING TIME

Approximately 8 to 9 elements are calculated per second in the isotropic version, the anisotropic version is five times slower.
top ]
7. UNUSUAL FEATURES OF THE PROGRAM:  Inclusion of the anisotropic version.
top ]
8. RELATED AND AUXILIARY PROGRAMS:  GAKER.
top ]
9. STATUS
Package ID Status date Status
NEA-0005/01 01-APR-1976 Tested at NEADB
NEA-0005/02 01-MAR-1965 Tested at NEADB
top ]
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).
top ]
11. MACHINE REQUIREMENTS:  Standard CDC 3600.
top ]
12. PROGRAMMING LANGUAGE(S) USED
Package ID Computer language
NEA-0005/01 FORTRAN-IV
NEA-0005/02 FORTRAN-IV
top ]
13. OPERATING SYSTEM OR MONITOR UNDER WHICH PROGRAM IS EXECUTED:  FORTRAN 3600 under SCOPE 5.1.
top ]
14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

ANY OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS
top ]
15. NAME AND ESTABLISHMENT OF AUTHOR

     R.J.J. Stammler
     Institutt for Atomenergi
     P.O. Box 40,
     Kjeller
     Norway.
top ]
16. MATERIAL AVAILABLE
NEA-0005/01
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
NEA-0005/02
File name File description Records
NEA0005_02.001 SOURCE & DATA 360/65 429
NEA0005_02.002 OUTPUT 360
top ]
17. CATEGORIES
  • A. Cross Section and Resonance Integral Calculations

Keywords: Bessel functions, cross sections, deuterons, energy transfer, heavy water, legendre polynomials, protons, scattering, transport theory.