Computer Programs

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 ]

[ 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 |
---|---|---|---|

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.

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.

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 ]

[ top ]

[ top ]

[ top ]

[ 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.

- 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. GongUser'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 ]

[ top ]

[ top ]

[ top ]

[ top ]

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 |

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