Computer Programs

NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROGRAM OR FUNCTION, METHOD OF SOLUTION, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES OF THE PROGRAM, RELATED AND AUXILIARY PROGRAMS, STATUS, REFERENCES, MACHINE REQUIREMENTS, LANGUAGE, OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED, OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHORS, MATERIAL, CATEGORIES

[ top ]

[ top ]

To submit a request, click below on the link of the version you wish to order. Rules for end-users are
available here.

Program name | Package id | Status | Status date |
---|---|---|---|

MUP-2 | IAEA0907/03 | Tested | 07-DEC-1989 |

Machines used:

Package ID | Orig. computer | Test computer |
---|---|---|

IAEA0907/03 | FACOM M-340-S | IBM 3090 |

[ top ]

3. DESCRIPTION OF PROGRAM OR FUNCTION

The complete program calculates the fast neutron data of medium-heavy nuclei. It includes: the total cross sections, elastic scattering cross section, nonelastic cross section, total inelastic cross section, (n,2n) and (n,3n), cross section, (n,n'x) cross section where x may be alpha, p, d, t and He3, (n,n') cross section for leaving the residual nucleus in from the first up to the 40th excited state and the continuum state, radiative capture cross section, (n,x) cross section where x may be p, d, t, He3 and alpha, (n,2alpha) and (n,2p) cross section, (n,px) cross section where x may be alpha, d, t and He3, (n,alphax) cross section where x may be d and t, the average cosine of the scattering angle (laboratory system) for elastic scattering, the average logarithmic energy decrement for elastic scattering, the average of the square of the logarithmic energy decrement for elastic scattering divided by twice the average logarithmic, the (n,x) cross section for leaving the residual nucleus in from the ground state up to be 17th excited state and in the continuum state where x may be p, d, t, He3 and alpha. The elastic scattering angular distribution and inelastic scattering angular distribution for leaving the residual nucleus in the first up to 40th excited state. The (n2n), (n,3n), (n,n'x) reactions, where x may be alpha, p, d, t

and He3 and the inelastic scattering secondary neutron spectra. All nuclear data are given for the natural elements as well as their isotopes and their output is according to ENDF/B-IV format. The incident neutron energy region, where 150 energy points at most may be included, is from 1 KeV to 20 MeV. The program includes the optical model, width fluctuation corrected Hauser-Feshbach formula and preequilibrium statistical theory based on the exciton model. The experimental direct reaction component of the inelastic scattering projectile exciting from the first up to the 4th excited state may be placed in the input. These are added to the calculated compound nucleus cross sections with a given weight. All nuclear level densities required are calculated according to the formulation of A. Cameron. The transmission coefficients or inverse cross sections of the particles used in statistical theory are calculated by the optical model. The gamma-ray transmission coefficients are calculated based on the giant dipole resonance model. The optical potentials used are Wood- Saxon for the real part, Wood Saxon and derivative Wood-Saxon for the imaginary part corresponding to the volume and surface absorptions respectively, and Thomas form for the spin-orbit part. For some special elements of which the experimental data are unknown, this program provides a microscopic optical potential calculation with skyrme forces for neutron or proton channels.

The complete program calculates the fast neutron data of medium-heavy nuclei. It includes: the total cross sections, elastic scattering cross section, nonelastic cross section, total inelastic cross section, (n,2n) and (n,3n), cross section, (n,n'x) cross section where x may be alpha, p, d, t and He3, (n,n') cross section for leaving the residual nucleus in from the first up to the 40th excited state and the continuum state, radiative capture cross section, (n,x) cross section where x may be p, d, t, He3 and alpha, (n,2alpha) and (n,2p) cross section, (n,px) cross section where x may be alpha, d, t and He3, (n,alphax) cross section where x may be d and t, the average cosine of the scattering angle (laboratory system) for elastic scattering, the average logarithmic energy decrement for elastic scattering, the average of the square of the logarithmic energy decrement for elastic scattering divided by twice the average logarithmic, the (n,x) cross section for leaving the residual nucleus in from the ground state up to be 17th excited state and in the continuum state where x may be p, d, t, He3 and alpha. The elastic scattering angular distribution and inelastic scattering angular distribution for leaving the residual nucleus in the first up to 40th excited state. The (n2n), (n,3n), (n,n'x) reactions, where x may be alpha, p, d, t

and He3 and the inelastic scattering secondary neutron spectra. All nuclear data are given for the natural elements as well as their isotopes and their output is according to ENDF/B-IV format. The incident neutron energy region, where 150 energy points at most may be included, is from 1 KeV to 20 MeV. The program includes the optical model, width fluctuation corrected Hauser-Feshbach formula and preequilibrium statistical theory based on the exciton model. The experimental direct reaction component of the inelastic scattering projectile exciting from the first up to the 4th excited state may be placed in the input. These are added to the calculated compound nucleus cross sections with a given weight. All nuclear level densities required are calculated according to the formulation of A. Cameron. The transmission coefficients or inverse cross sections of the particles used in statistical theory are calculated by the optical model. The gamma-ray transmission coefficients are calculated based on the giant dipole resonance model. The optical potentials used are Wood- Saxon for the real part, Wood Saxon and derivative Wood-Saxon for the imaginary part corresponding to the volume and surface absorptions respectively, and Thomas form for the spin-orbit part. For some special elements of which the experimental data are unknown, this program provides a microscopic optical potential calculation with skyrme forces for neutron or proton channels.

[ top ]

[ top ]

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

1) The calculations are restricted to medium-heavy nuclei, for which the fission reactions are absent.

2) The maximum number of the isotopes of the target nucleus is limited up to six.

3) The energy region of indicent neutrons is restricted to from 0.01 to 20MeV. In this region the maximum number of energy points is up to 150.

4) The maximum number of energy points of outgoing secondary neutrons is up to 100.

1) The calculations are restricted to medium-heavy nuclei, for which the fission reactions are absent.

2) The maximum number of the isotopes of the target nucleus is limited up to six.

3) The energy region of indicent neutrons is restricted to from 0.01 to 20MeV. In this region the maximum number of energy points is up to 150.

4) The maximum number of energy points of outgoing secondary neutrons is up to 100.

[ top ]

6. TYPICAL RUNNING TIME

The running time depends on the number of the open channels in the considered energy region and on the number of the isotopes of the target nucleus. For example, the running time for calculating a whole set of neutron data of the natural Sodium, for which the number of isotopes is only one and the third reaction process does not open, from 1MeV to 20MeV and for 79 energy points is 24.17 minutes on the M-340S.

The running time depends on the number of the open channels in the considered energy region and on the number of the isotopes of the target nucleus. For example, the running time for calculating a whole set of neutron data of the natural Sodium, for which the number of isotopes is only one and the third reaction process does not open, from 1MeV to 20MeV and for 79 energy points is 24.17 minutes on the M-340S.

IAEA0907/03

NEA-DB ran the test case included in this package on an IBM 3090 computer in 14 seconds fo CPU time.[ top ]

[ top ]

8. RELATED AND AUXILIARY PROGRAMS

MUP-2 is an extension and improvememt of MUP-1. AUJP is a related program, which automatically searches an optimum set of the neutron optical model parameters for medium-heavy nuclei by means of the CHI SQUARE method, this optimum set of the neutron optical model parameters is the input of MUP-2.

MUP-2 is an extension and improvememt of MUP-1. AUJP is a related program, which automatically searches an optimum set of the neutron optical model parameters for medium-heavy nuclei by means of the CHI SQUARE method, this optimum set of the neutron optical model parameters is the input of MUP-2.

[ top ]

10. REFERENCES

- C.M. Perey et al.,

Atomic Data and Nuclear Data Tables, 17, 3 (1976).

- Su Zong-Di et al.,

Physica Energiae Fortis et Physica Nuclears, 3,80 (1979).

- J.J. Griffrn,

Phys. Rev. Letters. 17,478 (1966).

- C.K. Cline,

Nuclear Physics, A 193,417 (1972).

- A. Gilbert et al.,

Can J. Phys. 43,1446 (1965).

- C.M. Perey et al.,

Atomic Data and Nuclear Data Tables, 17, 3 (1976).

- Su Zong-Di et al.,

Physica Energiae Fortis et Physica Nuclears, 3,80 (1979).

- J.J. Griffrn,

Phys. Rev. Letters. 17,478 (1966).

- C.K. Cline,

Nuclear Physics, A 193,417 (1972).

- A. Gilbert et al.,

Can J. Phys. 43,1446 (1965).

IAEA0907/03, included references:

- Cai Chonghai et al.:Directions for MUP2 Users (October 1986)

- Su Zong Di and P.E. Hodgson:

Comparison between Chinese Unified Program (MUP2) and

International Nuclear Model Programs.

[ top ]

11. MACHINE REQUIREMENTS

The tapes, discs and printer are necessary for output. Amount of memory: 4096K bytes.

The tapes, discs and printer are necessary for output. Amount of memory: 4096K bytes.

IAEA0907/03

Program execution on an IBM 3090 required 2764K bytes of main storage.[ top ]

IAEA0907/03

MVS/XA (IBM 3090).[ top ]

[ top ]

[ top ]

IAEA0907/03

File name | File description | Records |
---|---|---|

IAEA0907_03.001 | INFORMATION FILE | 59 |

IAEA0907_03.002 | MUP-2 FORTRAN SOURCE | 4092 |

IAEA0907_03.003 | JCL USED IN TESTING | 105 |

IAEA0907_03.004 | TEST CASE INPUT DATA | 142 |

IAEA0907_03.005 | TEST CASE OUTPUT DATA | 610 |

Keywords: cross sections, fast neutrons, nuclear reactions.