ESTS0486 FESH.

FESH, 2-D Multigroup Neutron Transport with Isotropic Scattering

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 ]

1. NAME OR DESIGNATION OF PROGRAM: FESH.

[ 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
FESH	ESTS0486/01	Arrived	17-APR-2001

Machines used:

Package ID	Orig. computer	Test computer
ESTS0486/01	CDC 6600

[ top ]

3. DESCRIPTION OF PROGRAM OR FUNCTION

FESH solves the spherical harmonics approximation to the second-order multigroup neutron transport equation in two dimensions. A three-dimensional spherical harmonics approximation of user-specified order is applied to the angular variable and two-dimensional spatial finite elements to the spatial variable. Scattering is assumed to be isotropic. Either multiplication eigenvalue problems or subcritical inhomogeneous source problems can be solved with several combinations of vacuum and reflective boundary conditions possible.

[ top ]

4. METHOD OF SOLUTION

FESH solves the multigroup transport equations (outer iterations) using a power iteration technique on the fission source with coarse mesh rebalance acceleration. The within-group (inner iteration) equation is solved by a block successive overrelaxation method with coarse mesh rebalance acceleration available as an option. Each block represents coupling and source information for each spatial mesh point.

[ top ]

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

Sources and scattering are assumed to be isotropic; up-scattering is not allowed. Flexible dimensioning permits combinations of problem parameters to occupy an array of 50,000 words. Most group-dependent data is stored in extended core storage (ECS).

[ top ]

6. TYPICAL RUNNING TIME

The time required for execution is problem-dependent. NESC executed the sample problem in 80 CP seconds on a CDC6600 and 20 CP seconds on a CDC7600.

[ top ]

7. UNUSUAL FEATURES OF THE PROGRAM

The unique feature is a rigorous application of transverse buckling which incorporates transverse flux anisotropy into the buckled leakage.

[ top ]

8. RELATED AND AUXILIARY PROGRAMS:

[ top ]

9. STATUS

Package ID	Status date	Status
ESTS0486/01	17-APR-2001	Masterfiled Arrived

[ top ]

10. REFERENCES:

ESTS0486/01, included references:

- R.N. Blomquist and E.E. Lewis:
  Instruction Manual for the Computer FESH - Finite Element
  Spherical Harmonics An x-y multigroup neutron transport computer
  code utilizing spatial finite element and angular spherical
  harmonic approximations to the even-parity form of the neutron
  transport equation
  (March 1978).

[ top ]

11. MACHINE REQUIREMENTS

200,000 (octal) words of central memory and 142,000 (octal) words of ECS are needed to execute the program. Central memory requirements for single precision arithmetic may be estimated as (45M+10) * (L(L+1)/2)**2 + NSP*(l(l+1)/2)+2) words, where M is the number of material regions, L is the order of the spherical harmonics approximation, and NSP is the number of spatial mesh points.

[ top ]

12. PROGRAMMING LANGUAGE(S) USED

Package ID	Computer language
ESTS0486/01	FORTRAN-IV

[ top ]

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED: SCOPE.

[ top ]

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

Because of the central memory required to load the program, the dimension of the Q array must be decreased before FESH can be loaded onto a CDC7600. NESC reduced the dimension of the Q array to 20,000 to execute the sample problem on a CDC7600. The dimensions of the Q and XLARG arrays should be altered to fit the user's specific problem.

[ top ]

15. NAME AND ESTABLISHMENT OF AUTHORS

          E.E. Lewis
      Northwestern University
      Evanston, IL (United States)

[ top ]

16. MATERIAL AVAILABLE

ESTS0486/01

source program   mag tapeFESH Source Program                        SRCTP
test-case data   mag tapeFESH Sample Problem Input                  DATTP
user's guide             Instruction Manual (March 1978)            WRKPT

[ top ]

17. CATEGORIES

C. Static Design Studies

Keywords: buckling, multigroup theory, spherical harmonics, three-dimensional, two-dimensional.