RSICC CODE PACKAGE PSR-443

1. NAME AND TITLE

SALE3D: ICEd-ALE Treatment of 3-D Fluid Flow.

2. CONTRIBUTORS

Los Alamos National Laboratory, Los Alamos, New Mexico through the Energy Science and Technology Software Center, Oak Ridge, Tennessee.

3. CODING LANGUAGE AND COMPUTER

FORTRAN IV; CRAY1 (P00443CY00000).

4. NATURE OF PROBLEM SOLVED

SALE3D calculates three-dimensional fluid flow at all speeds, from the incompressible limit to highly supersonic. An implicit treatment of the pressure calculation similar to that in the Implicit Continuous-fluid Eulerian (ICE) technique provides this flow speed flexibility. In addition, the computing mesh may move with the fluid in a typical Lagrangian fashion, be held in an Eulerian manner, or move in some arbitrarily specified way to provide a continuous rezoning capability. This latitude results from use of an Arbitrary Lagrangian-Eulerian (ALE) treatment of the mesh. The partial differential equations solved are the Navier-Stokes equations and the mass and internal energy equations. The fluid pressure is determined from an equation of state and supplemented with an artificial viscous pressure for the computation of shock waves. The computing mesh consists of a three-dimensional network of arbitrarily shaped, six-sided deformable cells, and a variety of user-selectable boundary conditions are provided in the program.

5. METHOD OF SOLUTION

SALE-3D uses an ICEd-ALE technique, which combines the ICE method of treating flow speeds and the ALE mesh treatment to calculate three-dimensional fluid flow. The finite-difference approximations to the conservation of mass, momentum, and specific internal energy differential equations are solved in a sequence of time steps on a network of deformable computational cells. The basic hydrodynamic part of each cycle is divided into three phases: (1) an explicit solution of the Lagrangian equations of motion updating the velocity field by the effects of all forces, (2) an implicit calculation using a Newton-Raphson iterative scheme that provides time-advanced pressures and velocities, and (3) the addition of advective contributions for runs that are Eulerian or contain some relative motion of grid and fluid. A powerful feature of this three-phase approach is the ease with which different phases can be combined to suit the requirements of individual problems.

6. RESTRICTIONS OR LIMITATIONS

None noted.

7. TYPICAL RUNNING TIME

NESC executed the sample problem in 4 CP seconds on a Cray Y-MP/832.

8. COMPUTER HARDWARE REQUIREMENTS

430K (octal) words are required to execute the sample problem on a Cray Y-MP/832.

9. COMPUTER SOFTWARE REQUIREMENTS

CTSS (Cray 1), UNICOS5.0 (Cray Y-MP).

10. REFERENCES

a) Included with package:

C. Yuelys-Miksis, "SALE3D, NESC No. 1069, SALE3D Tape Description and Implementation Information," National Energy Software Center Note 90-35 (January 4, 1990).

A.A. Amsden and H.M. Ruppel, "SALE3D: A Simplified ALE Computer Program for Calculating Three-Dimensional Fluid Flow," NUREG/CR-2185 (LA-8905) (November 1981).

b) Background information:

A.A. Amsden, H.M. Ruppel, and C.W. Hirt, "SALE: A Simplified ALE Computer Program for Fluid Flow at All Speeds," LA-8095 (June 1980).

A.A. Amsden and C.W. Hirt, "YAQUI: An Arbitrary Lagrangian-Eulerian Computer Program for Fluid Flow at All Speeds," LA-5100 (March 1973).

11. CONTENTS OF CODE PACKAGE

Included are the referenced documents in (10.a) and a 3.5" DOS formatted diskette containing the source code and sample problem.

12. DATE OF ABSTRACT

September 1999.

KEYWORDS: FLUID DYNAMICS; COMPLEX GEOMETRY