NESC9473 IMPSOR.

IMPSOR, 3-D Boundary Problems Solution for Thermal Conductivity Calculation

NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROGRAM OR FUNCTION, METHOD OF SOLUTION, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES OF THE PROGRAM, RELATED AND AUXILIARY PROGRAMS, STATUS, REFERENCES, MACHINE REQUIREMENTS, LANGUAGE, OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED, OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHORS, MATERIAL, CATEGORIES

[ top ]

1. NAME OR DESIGNATION OF PROGRAM: IMPSOR.

[ top ]

2. COMPUTERS

To submit a request, click below on the link of the version you wish to order. Only liaison officers are authorised to submit online requests. Rules for requesters are available here.

Program name	Package id	Status	Status date
IMPSOR	NESC9473/01	Tested	24-FEB-1994

Machines used:

Package ID	Orig. computer	Test computer
NESC9473/01	CRAY X-MP	CRAY EL98

[ top ]

3. DESCRIPTION OF PROGRAM OR FUNCTION

IMPSOR implements finite difference methods for multidimensional moving boundary problems with Dirichlet or Neumann boundary conditions. The geometry of the spatial domain is a rectangular parallelepiped with dimensions specified by the user. Dirichlet or Neumann boundary conditions may be specified on each face of the box independently. The user defines the initial and boundary conditions as well as the thermal and physical properties of the problem and several parameters for the numerical method, e.g. degree of implicitness, time-step size.

[ top ]

4. METHOD OF SOLUTION

The spatial domain is partitioned and the governing equation discretized, which yields a nonlinear system of equations at each time-step. This nonlinear system is solved using a successive overrelaxation (SOR) algorithm. For a given node, the previous iteration's temperature and thermal conductivity values are used for advanced points with current values at previous points. This constitutes a Gauss-Seidel iteration. Most of the computing time used by the numerical method is spent in the iterative solution of the nonlinear system. The SOR scheme employed is designed to accommodate vectorization on a Cray X-MP.

[ top ]

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM: Maximum of -, 70,000 nodes

[ top ]

6. TYPICAL RUNNING TIME:

NESC9473/01

The program was tested at NEA-DB on a CRAY-EL98 computer. By adjusting input parameters, CPU time could be limited to 19.5 seconds.

[ top ]

7. UNUSUAL FEATURES OF THE PROGRAM:

[ top ]

8. RELATED AND AUXILIARY PROGRAMS:

[ top ]

9. STATUS

Package ID	Status date	Status
NESC9473/01	24-FEB-1994	tested vectorized

[ top ]

10. REFERENCES:

NESC9473/01, included references:

- M.A. Williams and D.G. Wilson:
  IMPSOR - A Fully Vectorized FORTRAN Code for Three Dimensional
  Moving Boundary Problems with Dirichlet or Neumann Boundary
  Conditions
  ORNL-6393 (August 1987).

[ top ]

11. MACHINE REQUIREMENTS:

[ top ]

12. PROGRAMMING LANGUAGE(S) USED

Package ID	Computer language
NESC9473/01	FORTRAN

[ top ]

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:

NESC9473/01

The program was tested by NEA-DB under the UNICOS operating system.

[ top ]

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS:

[ top ]

15. NAME AND ESTABLISHMENT OF AUTHORS

D.G. Wilson and M.A. Williams
Oak Ridge National Laboratory

[ top ]

16. MATERIAL AVAILABLE

NESC9473/01

File name	File description	Records
NESC9473_01.001	Information file	41
NESC9473_01.002	JCL and control information	10
NESC9473_01.003	IMPSOR FORTRAN source	820
NESC9473_01.004	Sample problem input	18
NESC9473_01.005	Sample problem output 6	4857
NESC9473_01.006	Sample problem output 7	14180

[ top ]

17. CATEGORIES

H. Heat Transfer and Fluid Flow
P. General Mathematical and Computing System Routines

Keywords: Dirichlet problem, algorithms, finite difference method, temperature, thermal conductivity, three-dimensional.