NESC9599 TOPAZ3D.

1. NAME OR DESIGNATION OF PROGRAM:  TOPAZ3D.

3. DESCRIPTION OF PROGRAM OR FUNCTION

TOPAZ3D is a three-dimensional  implicit finite element computer code for heat transfer analysis. TOPAZ3D can be used to solve for the steady-state or transient temperature field on three-dimensional geometries. Material properties may be temperature-dependent and either, isotropic or orthotropic. A variety of time- and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation. By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functional representation of boundary conditions and internal heat generation. TOPAZ3D can solve problems of diffuse and specular band  radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.
  
This package is distributed by:
Energy Science and Technology Software Center
P.O. Box 62
1 Science.Gov Way
Oak Ridge, TN 37831
(865) 576-2606 TEL
(865) 576-6436 FAX
E-mail: ESTSC@osti.gov

4. METHOD OF SOLUTION

TOPAZ3D solves the differential equation of heat conduction in a three-dimensional solid by the finite element method. TOPAZ3D uses an eight-node trilinear hexahedral element for  spatial discretization of the geometry. The hexahedral element can degenerate to a six-node triangular prism and a four-node tetrahedron. These elements are integrated with a 2x2x2 Gauss quadrature rule, with temperature dependence of the properties accounted for at the Gauss point. Time integration is performed using a generalized trapezoidal method. Fixed point iteration with relaxation is used to satisfy equilibrium in nonlinear problems.

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

The phase change, slide surface, internal element, and bulk node features are not implemented.

6. TYPICAL RUNNING TIME:

7. UNUSUAL FEATURES OF THE PROGRAM:

8. RELATED AND AUXILIARY PROGRAMS

TOPAZ3D is an extension to three dimensions of the two-dimensional heat transfer code TOPAZ2D (NESC 9801). It is based to some extent on work done in developing the TACO3D code. TOPAZ3D has no general mesh generation capability. Rows of evenly spaced nodes and rows of sequential elements may be generated. For complex zoning, the mesh generation code and preprocessor INGRID (NESC 9649) should be used. The TAURUS (NESC 9908) interactive postprocessor can be used to produce temperature contour, flux contour, temperature-time history, and various geometry plots. Thermal stresses can be calculated using the solid mechanics code NIKE3D (NESC 9725) which reads the temperature state  data calculated by TOPAZ3D.

10. REFERENCES

- TOPAZ3D, NESC No. 9599.MPL2, Tape Description and Implementation
  Information for Multiple Machine Version of TOPAZ3D,
  National Energy Software Center 87-49, March 24, 1987.
- TOPAZ3D, NESC 9599.3090, TOPAZ3D IBM3090 Version Tape Description,   National Energy Software Center Note 88-76, June 8, 1988.

11. MACHINE REQUIREMENTS:

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:  VM/MVS (IBM3090).

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

TOPAZ3D is a vectorized code. The FORTRAN source as supplied for the Cray1,  DEC VAX11 version is set up for execution on the DEC VAX11. The Cray user must enable those source records with the identifier "CRAY" in  columns 73-76 and disable those with "VAX" identifiers. NESC packaged the contributed VAX executable image with the DEC VAX11 customized source as a separate machine version for distribution to  DEC VAX installations.

15. NAME AND ESTABLISHMENT OF AUTHORS

           A.B. Shapiro
           Mechanical Engineering Department
           Lawrence Livermore National Laboratory

   3090    G.A. Vitez
           IBM Palo Alto Science Center