Abstract This paper describes RHALE, a multi-material arbitrary Lagrangian-Eulerian (MMALE) shock physics code. RHALE is the successor to CTH, Sandia's 3-D Eulerian shock physics code, and will be capable of solving problems that CTH cannot adequately address. We discuss the new Lagrangian capabilities of RHALE, which include arbitrary mesh connectivity, superior artificial viscosity, and improved equations of state. We also discuss some of the issues we have encountered in the choice of an axisymmetric element technology and our resolution of these issues. We discuss the MMALE algorithms that have been extended for arbitrary grids in both two and three dimensions and present the results of calculations that are of interest to the hypervelocity impact community. The MMALE addition to RHALE provides the accuracy of a Lagrangian code while allowing a calculation to proceed under very large distortions. Coupling an arbitrary quadrilateral or hexahedral grid to the MMALE solution facilitates modeling of complex shapes with a minimum number of computational cells. RHALE allows regions of a problem to be modeled with Lagrangian, Eulerian or ALE meshes. In addition, regions can switch from Lagrangian to ALE to Eulerian based on user input or mesh distortion. For ALE meshes, new node locations are determined with equipotential schemes. Element quantities are advected with donor, van Leer, or Super-B algorithms. Nodal quantities are advected with the second order SHALE or HIS algorithms. Currently, material interfaces are determined with the SLIC algorithm; however, both two and three dimensional versions of Youngs' interface tracker are being investigated. To facilitate the development of such a lengthy code, we choose to write in the C++ programming language. We feel that object-oriented programming techniques are superior to conventional programming techniques. However, we discuss a few of the efficiency problems we have encountered using these techniques and how we have addressed these problems.
[1]
David J. Benson,et al.
Momentum advection on a staggered mesh
,
1992
.
[2]
J. Z. Zhu,et al.
The finite element method
,
1977
.
[3]
Kent G. Budge,et al.
ALE Shock Calculations Using a Stabilized Serendipity Rezoning Scheme
,
1992
.
[4]
R. T. Barton,et al.
Rescaling of equipotential smoothing
,
1982
.
[5]
C. T. Vaughan,et al.
Massively Parallel Computing, C++ and Hydrocode Algorithms
,
1992
.
[6]
P. Woodward,et al.
SLIC (Simple Line Interface Calculation)
,
1976
.
[7]
D. Benson.
An efficient, accurate, simple ALE method for nonlinear finite element programs
,
1989
.
[8]
John K. Dukowicz,et al.
Vorticity errors in multidimensional Lagrangian codes
,
1992
.
[9]
Herbert Schildt.
Advanced C
,
1986
.
[10]
Bjarne Stroustrup,et al.
The Annotated C++ Reference Manual
,
1990
.
[11]
L. G. Margolin,et al.
A method for treating hourglass patterns
,
1987
.
[12]
L. G. Margolin,et al.
Remapping on the staggered mesh
,
1988
.