NESC0586 PLENUM

1. NAME OR DESIGNATION OF PROGRAM:  PLENUM

3. DESCRIPTION OF PROBLEM OR FUNCTION

Bulk flow distribution in a cylindrical reactor coolant inlet plenum is calculated using potential flow theory.

4. METHOD OF SOLUTION

The Laplace equation for the velocity potential is solved for general, symmetric boundary conditions and one inlet nozzle by Gauss-Seidel relaxation (part 1). Nonsymmetric, dimensional solutions are then formed by superposition of the general solution (part 2).

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

The code is dimensioned for 8000 mesh points in an (r,theta,z) array. Dimensions can be expanded at the expense of run time and core size.

6. TYPICAL RUNNING TIME

Run time for part 1 is proportional to n**2,  where n is the number of mesh points. With 2299 mesh points CPU time is 42 seconds on a 360/195. Run time for part 2 is 1 second per problem. Sample problem running time is 2 minutes for part 1 and 1 minute for part 2.

7. UNUSUAL FEATURES OF THE PROGRAM

The code approximates bulk fluid motion through large cylindrical headers in circumstances where details of the flow, such as boundary layer effects, may be neglected. Three-dimensional effects are calculated, so that asymmetric fluid input and output problems may be calculated. Use of a linear field equation permits rapid solution by superposition of many combinations of input and output boundary conditions. First order effects of friction may be included by addition of constant coefficients in the finite difference formulations.

8. RELATED AND AUXILIARY PROGRAMS:

10. REFERENCE

S. D. Harris and J. W. Reece, Application of Potential
Flow Theory  to Reactor  Coolant Circulation,  Proceedings of  the
Eighth  Southeastern  Seminar on  Thermal  Sciences,  DP-MS-71-13,
March 23-24, 1972.

11. MACHINE REQUIREMENTS:  240K bytes of storage and an on-line punch and printer.

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:   OS/360.

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

PLENUM  is written as two separate programs. Part 1 provides a general dimensionless solution to the field equation. For a given header this code need only be executed once. Part 2 uses punched output from part 1 to determine solutions for a variety of user-specified boundary conditions by dimensionalization and superposition.

15. NAME AND ESTABLISHMENT OF AUTHOR

                 S. D. Harris
                 Savannah River Laboratory
                 E. I. duPont de Nemours and Company, Inc.
                 Aiken, South Carolina  29801

File name	File description	Records
NESC0586_01.001	PART 1 - SOURCE PROGRAM (F4)	312
NESC0586_01.002	PART 1 - SAMPLE PROBLEM DATA	3
NESC0586_01.003	PART 1 - SAMPLE PROBLEM PRINTED OUTPUT	510
NESC0586_01.004	PART 2 - SOURCE PROGRAM (F4)	220
NESC0586_01.005	INFORMATION	2
NESC0586_01.006	PART 2 - SAMPLE PROBLEM DATA	435
NESC0586_01.007	PART 2 - SAMPLE PROBLEM PRINTED OUTPUT	538