RSICC CODE PACKAGE PSR-199
1. NAME AND TITLE
HEATING 7: Multidimensional, Finite-Difference Heat Conduction Analysis Code System, Versions 7.2i and 7.3.
AUXILIARY ROUTINES
H7MAP: Postprocessor code to print maps of node numbers and/or temperature distributions.
H7NODE: Returns node number at specified grid location.
H7TECPLOT: Reads HEATING plot file and produces input file for graphical postprocessing with TECPLOT, a proprietary software product of Tecplot, Inc., which is not included in this package.
DATA LIBRARY
H7MATLIB: Lawrence Livermore National Laboratory Material Properties Library
2. CONTRIBUTOR
Oak Ridge National Laboratory, Oak Ridge, Tennessee.
3. CODING LANGUAGE AND COMPUTER
HEATING 7.3: Fortran 90 for Pentium/Windows 2000/XP/Vista and IBM RS/6000.
HEATING 7.2i: Fortran 77 and C; UNIX Workstation or Mainframe
Both Windows and Unix versions are included in RSICC package P00199MNYCP06.
4. NATURE OF PROBLEM SOLVED
HEATING can solve steady-state and/or transient heat conduction problems in one-, two-, or three-dimensional Cartesian, cylindrical, or spherical coordinates. A model may include multiple materials, and the thermal conductivity, density, and specific heat of each material may be both time- and temperature-dependent. The thermal conductivity may also be anisotropic. Materials may undergo change of phase. Thermal properties of materials may be input or may be extracted from a material properties library. Heat- generation rates may be dependent on time, temperature, and position, and boundary temperatures may be time- and position-dependent. The boundary conditions, which may be surface-to-environment or surface-to-surface, may be specified temperatures or any combination of prescribed heat flux, forced convection, natural convection, and radiation. The boundary condition parameters may be time- and/or temperature-dependent. General graybody radiation problems may be modeled with user-defined factors for radiant exchange. The mesh spacing may be variable along each axis. HEATING uses a run-time memory allocation scheme to avoid having to recompile to match memory requirements for each specific problem. HEATING utilizes free-form input.
HEATING 7.2i and 7.3 are the most recent developments in a series of heat-transfer codes and obsolete all previous versions distributed by RSICC as SCA-1/HEATING5 and PSR-199/HEATING 6. No future modifications are planned for HEATING7.
In August 2007, the Windows package was updated with the addition of a verification report for HEATING7.3 and corresponding test case input and output files plus an executable that allows one to run in batch mode from the DOS prompt. This package replaced the previous Windows version. No changes were made to the Unix version.
5. METHOD OF SOLUTION
Three steady-state solution techniques are available: point-successive-overrelaxation iterative method with extrapolation, direct-solution (for one-dimensional or two-dimensional problems), and conjugate gradient. Transient problems may be solved using any one of several finite-difference schemes: Crank-Nicolson implicit, Classical Implicit Procedure (CIP), Classical Explicit Procedure (CEP), or Levy explicit method (which for some circumstances allows a time step greater than the CEP stability criterion.) The solution of the system of equations arising from the implicit techniques is accomplished by point-successive-overrelaxation iteration and includes procedures to estimate the optimum acceleration parameter.
6. RESTRICTIONS OR LIMITATIONS
All surfaces in a model must be parallel to one of the coordinate axes which makes modeling complex geometries difficult. Transient change of phase problems can only be solved with one of the explicit techniques - an implicit change-of-phase capability has not been implemented.
7. TYPICAL RUNNING TIME
Run time is dependent on many factors including among other things: number of nodes in the model, nonlinearities in the model, desired solution(s), and speed of the computer on which it is run. Therefore, it is not possible to specify a typical running time. The small sample problems supplied with this distribution required an average less than 2 cpu seconds each on a Pentium IV. The verification and reference test cases in the HEATING 7.3 distribution for Windows ran in ~10 seconds on a Pentium 4 3.0 GHz PC.
8. COMPUTER HARDWARE REQUIREMENTS
HEATING7.3 for PC runs under Windows operating systems. The Unix version of HEATING7.3 was tested only on IBM RS/6000. HEATING 7.2i ran on a variety of Unix workstations in the 1990s.
9. COMPUTER SOFTWARE REQUIREMENTS
HEATING7.3 for PC was tested under Windows XP, Windows 2000 and Windows Vista. Included HEATING 7.3 PC executables were created under Windows XP using Compaq Visual Fortran Version 6.6. HEATING7.3 for Unix was tested on IBM RS/6000 under AIX using C and XLF 3.2 compilers. The code was not tested on other platforms. Except for modifications to Subroutine MENSURAL in HEATING 7.3, the Unix versions have not been updated since the 1990's and will likely require modifications to run on current systems. HEATING 7.3 includes some Fortran 90 features.
Plotting with HEATING7 requires the proprietary software TECPLOT, which is not included in this package. If you wish to run H73TEC, you must purchase TECPLOT software from: Tecplot, Inc., phone: 800-763-7005, http://www.tecplot.com. Note that TECPLOT is not required to run the Heating7.3 code.
10. REFERENCES
10.a: included in RSICC document P199.PDF:
K. W. Childs, “Draft Input Description for HEATING 7.3" (February 26, 1997).
K. W. Childs, "HEATING 7.2 User's Manual," ORNL/TM-12262 (February 1993; update to Chapter 4, September 1998).
K. W. Childs, "Verification and Validation of HEATING7.3" ORNL Internal Report (March 2005).
K. W. Childs, “Note on Revisions: g, h, and i” (October 20, 1994).
10.b - Included in distribution files:
C. B. Bryan, K. W. Childs, and G. E. Giles, “HEATING6 Verification,” K/CSD/TM‑61, UC‑32 (December 1966).
RSICC, “README.1ST” (July 2004, revised August 2007).
RSICC, “READUNIX.TXT” (July 2002).
11. CONTENTS OF CODE PACKAGE
Included are the referenced documents and one CD-ROM with a WinZip file and a GNU compressed tar file for Unix. Both versions include source files, a material properties data library, and sample cases. In addition, the Windows version includes executables.
12. DATE OF ABSTRACT
February 1993, March 1993, November 1993, January 1995, June 1997, May 1998, September 1998, November 1999, July 2002, July 2004, August 2007.
KEYWORDS: COMPLEX GEOMETRY; CYLINDRICAL GEOMETRY; HEAT TRANSFER; MICROCOMPUTER; WORKSTATION