Publications

Computational physicists at Sandia National Laboratories have moved the Eulerian CTH code, and the arbitrary-Lagrangian-Eulerian ALEGRA code to distributed memory parallel computers. CTH is a three-dimensional solid mechanics code used for large-deformation, shock wave analysis. ALEGRA is a three-dimensional arbitrary Lagrangian-Eulerian solid-mechanics code used for coupled large-deformation, shock and structural mechanics problems. This paper discusses our experiences moving the codes to parallel computers, the algorithms we used and our experiences running the codes.

Computational physicists at Sandia National Laboratories have moved their production codes to distributed memory parallel computers. Such an effort required the development of parallel algorithms, parallel data bases and parallel support tools. The Eulerian CTH code was rewritten. Moving both ALEGRA and PRONTO to parallel computers required only a modest number of modifications. It involved restructuring the restart and graphics data bases to make them parallel and minimize the I/O to the parallel computer. It also involved developing mesh decomposition tools to divide a rectangular or arbitrary connectivity into sub-meshes. It also involved developing new visualization tools to process the very large, parallel data bases. This paper also discusses Sandia's experiences running these codes on its 1840 compute node Intel Paragon, 1024 processor nCUBE and networked stations.

Shock physics codes for armor analysis and the computers used to run them have improved dramatically over the last five years. Improved algorithms and material models allow accurate and efficient modeling of conventional armors. Desktop workstations routinely perform production two-dimensional calculations and massively parallel computers perform three-dimensional calculations. New developments in codes and computers promise improved accuracy, increased capabilities and faster simulation. This paper describes the current state of armor analysis codes and computers. Armor analysis code features fall into one of five categories: mesh, solution algorithm, material model, heuristic or computer. This paper discusses the current state of each of these categories and discusses the additional work needed.

This paper is an introductory discussion of stress pulse phenomena in simple solids and fluids. Stress pulse phenomena is a very rich and complex field that has been studied by many scientists and engineers. This paper describes the behavior of stress pulses in idealized materials. Inviscid fluids and simple solids are realistic enough to illustrate the basic behavior of stress pulses. Sections 2 through 8 deal with the behavior of pressure pulses. Pressure is best thought of as the average stress at a point. Section 9 deals with shear stresses which are most important in studying solids.

Continuum dynamics codes are categorized as Lagrangian or Eulerian according to the motion of the mesh. A Lagrangian code`s mesh moves with the material, so no mass flows between cells. An Eulerian code`s mesh is stationary, so mass flows between the cells. Eulerian codes have improved to the point where they are routinely used to solve a broad variety of large deformation solid and fluid dynamics problems ranging from air flow over an airplane wing to meteor impact on space structures. This presentation will concentrate on multi-fluid Eulerian codes capable of modeling transient were propagation in solids. These codes use a two-step process to integrate the physics across a time step. The first step, referred to as the Lagrangian step, integrates the physics on a Lagrangian mesh across the time step. The field values are then at the new time, but they are on the distorted Lagrangian mesh. The second step, referred to as the remap step, remaps the data on the distorted Lagrangian mesh back to the original Eulerian mesh thus completing one time step. The algorithms used in the first step are similar to those used in modern Lagrangian codes but they must be extended to handle multi-material cells. The algorithms used in the second step are complex and must be very carefully chosen to minimize errors. These algorithms include second-order, monotone advection equations to calculate the quantities flowing between cells. They also require algorithms that construct material interfaces inside multi-material cells. The strength and limitations of currently used numerical techniques will be discussed. New code development activities that combine the best features on both Lagrangian and Elueian codes will also be discussed. These new codes will employ the strengths of both technologies to address problems that cannot be adequately solved at this time.

Sandia National Laboratories is very active in developing multi-dimensional, multi-material shock wave physics codes. One example is the state-of-the-art, three-dimensional Eulerian code CTH which is used at numerous government and university sites. CTH is being ported to both Single Instruction Multiple Data (SIMD) and Multiple Instruction Multiple Data (MIMD) massively parallel computers. The next-generation arbitrary-Lagrangian-Eulerian code RHALE is under development. This paper will discuss these codes. CTH is an Eulerian code for modelling multi-dimensional, multi-material, large deformation, strong shock physics. Finite-volume numerical schemes are used with one-dimensional, two-dimensional and three-dimensional meshes. CTH has models for elastic-plastic materials, porous materials, high explosive detonation, fracture, and energy deposition. Several analytic equations of state are available including ideal gasses, Jones-Wilkins-Lee high explosive reaction products, Mie-Griineisen solids, and sophisticated multi-phase models that are valid for a very broad range of densities and temperatures. Tabular equations of state are also available. Second-order accurate advection schemes are used to minimize the dispersion found in Eulerian codes. Very large three-dimensional calculations may be run efficiently on a CRAY supercomputer because the code is highly vectorized and the data bases reside on the Solid State Disk (SSD). Data bases larger than one hundred million words are commonly used. Sophisticated color post-processing software was developed to aid in interpreting the results. Much of CTH has been ported to both SIMD and MIMD massively parallel computers. The two-dimensional version is running three times faster than a single CPU CRAY/YMP on the 16k node SIMD Connection Machine and five times faster on the 1024 node nCUBE2 MIMD computer. Both of the massively parallel computers can be expanded by a factor of four to eight yielding a system an order of magnitude faster than a CRAY. The next-generation, three-dimensional arbitrary-Lagrangian-Eulerian code RHALE is under development. Finite element techniques are used to integrate the physics through time. The mesh will move with the material (Lagrangian mesh) until the distortion becomes excessive and then the nodes are automatically repositioned to smooth the mesh and improve the accuracy (Eulerian mesh). Node motion occurs only where the distortion is excessive. This results in a code with the best features of both Lagrangian and Eulerian codes. Arbitrary-connectivity meshes are used to generate very complicated and sophisticated meshes. However, this dramatically increases the complexity of the Eulerian algorithms.