Computer Programs

NESC0865 PELE-IC.

NESC0865 PELE-IC.

PELE-IC, 2-D Eulerian Incompressible Hydrodynamic and Bubble Dynamic after LWR LOCA

NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROBLEM 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 AUTHOR, MATERIAL, CATEGORIES

[ top ]

1. NAME OR DESIGNATION OF PROGRAM: PELE-IC.

[ top ]

2. COMPUTERS

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
PELE-IC	NESC0865/01	Tested	01-DEC-1981

Machines used:

Package ID	Orig. computer	Test computer
NESC0865/01	CDC 7600	CDC 7600

[ top ]

3. DESCRIPTION OF PROBLEM OR FUNCTION

PELE-IC is a two-dimensional semi-implicit Eulerian hydrodynamics program for the solution of incompressible flow coupled to flexible structures. The code was developed to calculate fluid-structure interactions and bubble dynamics of a pressure-suppression system following a loss-of- coolant accident (LOCA). The fluid, structure, and coupling algorithms have been verified by calculation of benchmark problems and air and steam blowdown experiments. The code is written for both plane and cylindrical coordinates. The coupling algorithm is general enough to handle a wide variety of structural shapes. The concepts of void fractions and interface orientation are used to track the movement of free surfaces, allowing great versatility in following fluid-gas interfaces both for bubble definition and water surface motion without the use of marker particles.

[ top ]

4. METHOD OF SOLUTION

The solution strategy is to first solve the Navier-Stokes equations explicitly using values from the previous time-step. Since these values do not necessarily satisfy the continuity equation, the pressure field is iterated upon until the incompressibility condition for each computational cell is satisfied within prescribed limits.
The structural motion is computed by a finite element code from the applied pressure at the fluid-structure interface. The shell structure algorithm uses conventional thin-shell theory with transverse shear. The finite-element spatial discretization employs piecewise-linear interpolation functions and one-point quadrature applied to conical frustra. The Newmark implicit time integration method is used as a one-step module. The fluid code then uses the structure's position and velocity as boundary conditions. The fluid pressure field and the structure's response are corrected iteratively until the normal velocities of fluid and stucture are equal. The effects of steam condensation and oscillatory chugging on structures are calculated according to a simplified theory of condensation on a free surface due to C.S. Landram.

[ top ]

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM:

[ top ]

6. TYPICAL RUNNING TIME

Typical problem run times are from 5 to 40 minutes on the CDC7600. The running time for the six sample problems ranged from 32 to 172 CP seconds.

[ top ]

7. UNUSUAL FEATURES OF THE PROGRAM:

[ top ]

8. RELATED AND AUXILIARY PROGRAMS

The basic semi-implicit solution algorithm of SOLA (NESC Abstract 651) was used as a foundation for the development of the PELE-IC code.

[ top ]

9. STATUS

Package ID	Status date	Status
NESC0865/01	01-DEC-1981	Tested at NEADB

[ top ]

10. REFERENCES

- W.H. McMaster. D.M. Norris Jr., G.L. Goudreau, D.F. Quinones,
  E.Y. Gong, B. Moran, and N.A. Macken,
  Coupled Fluid-Structure method for Pressure Suppression Analysis,
  NUREG/CR-0607, May 1979.
- W.H. McMaster, E.Y. Gong, C.S. Landram, and D.F. Quinones,
  Fluid Structure Coupling Algorithm,
  Computers and Structures, Vol. 13, pp. 163-166, 1981.

NESC0865/01, included references:

- W. H. McMaster and E. Y. Gong
  User's Manual for PELE-IC: A Computer Code for Eulerian
  Hydrodynamics. UCRL-52609 (May 29, 1979).
- E. Y. Gong, E. E. Alexander, W. H. McMaster and D. F. Quinones
  PELE-IC Test Problems. UCRL-52835 (October 1, 1979).

[ top ]

11. MACHINE REQUIREMENTS

123000 (octal) words of SCM and 272000 (octal) words of ICM are used.

[ top ]

12. PROGRAMMING LANGUAGE(S) USED

Package ID	Computer language
NESC0865/01	FORTRAN-IV

[ top ]

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED: SCOPE 2.1, LTSS.

[ top ]

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

The
program relies on the contents of memory being preset to zero.

[ top ]

15. NAME AND ESTABLISHMENT OF AUTHOR

W.H. McMaster and E.Y. Gong
Lawrence Livermore National Laboratory
P.O. Box 808
Livermore, California 94550

[ top ]

16. MATERIAL AVAILABLE

NESC0865/01

File name	File description	Records
NESC0865_01.001	PELE-IC FORTRAN SOURCE	19532
NESC0865_01.002	PELE-IC UPDATE SOURCE	19582
NESC0865_01.003	PELE-IC SAMPLE PROBLEM 1,TAPE59	2
NESC0865_01.004	PELE-IC SAMPLE PROBLEM 1,TAPE9	28
NESC0865_01.005	PELE-IC SAMPLE PROBLEM 2,TAPE59	2
NESC0865_01.006	PELE-IC SAMPLE PROBLEM 2,TAPE9	67
NESC0865_01.007	PELE-IC SAMPLE PROBLEM 3,TAPE59	2
NESC0865_01.008	PELE-IC SAMPLE PROBLEM 3,TAPE9	52
NESC0865_01.009	PELE-IC SAMPLE PROBLEM 4,TAPE59	2
NESC0865_01.010	PELE-IC SAMPLE PROBLEM 4,TAPE9	50
NESC0865_01.011	PELE-IC SAMPLE PROBLEM 5,TAPE59	2
NESC0865_01.012	PELE-IC SAMPLE PROBLEM 5,TAPE9	67
NESC0865_01.013	PELE-IC SAMPLE PROBLEM 6,TAPE59	2
NESC0865_01.014	PELE-IC SAMPLE PROBLEM 6,TAPE9	54
NESC0865_01.015	JCL & INFORMATION	83
NESC0865_01.016	SAMPLE PROBLEM 1 OUTPUT	1797
NESC0865_01.017	SAMPLE PROBLEM 2 OUTPUT	1742
NESC0865_01.018	SAMPLE PROBLEM 3 OUTPUT	2370
NESC0865_01.019	SAMPLE PROBLEM 4 OUTPUT	1998
NESC0865_01.020	SAMPLE PROBLEM 5 OUTPUT	2015
NESC0865_01.021	SAMPLE PROBLEM 6 OUTPUT	2274

[ top ]

17. CATEGORIES

H. Heat Transfer and Fluid Flow

Keywords: blowdown, finite element method, fluid flow, hydrodynamics, incompressible flow, pressure.