Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Computational experience with an explicit difference scheme for a one phase Stefan problem

Technical Report ·
OSTI ID:5217095

The results of a computational scheme of Rose (a method for calculating solutions of parabolic equations with free boundary, Math. Comp., 14: 249-256(1960)) are compared with the known analytic solution of a single-phase Stefan problem with constant boundary temperature. Rose's scheme was implemented in FORTRAN on a DEC PDP-10. Computations were performed for Stefan numbers St = 0.01, 0.1, 1, 10, and 100. Three relative mesh sizes (..gamma.. = DT/(DX)/sup 2/ = 1/2, 1/6, and 1/25) were used. For each St, the relative mesh size had very little effect on accuracy. However, taking ..gamma.. = 1/25 resulted in dramatically increased computing time. The relative errors of temperature were large near the phase change front. 11 figures, 5 tables.

Research Organization:
Oak Ridge National Lab., TN (USA)
DOE Contract Number:
W-7405-ENG-26
OSTI ID:
5217095
Report Number(s):
ORNL/CSD-62
Country of Publication:
United States
Language:
English