last modified: 01-DEC-1981 | catalog | categories | new | search |

NESC0865 PELE-IC.

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

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. Only liaison officers are authorised to submit online requests. Rules for requesters 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.