A practical numerical method is described for the solution of transient heat conduction where the thermal properties are temperature dependent and changes of phase occur. The procedure involves the use of the explicit finite difference technique, for which a stability criterion is given, in conjunction with two integral transformations.
Get full access to this article
View all access options for this article.
References
1.
CarslawH. S.JaegerJ. C.Conduction of heat in solids1959, 2nd edition (Clarendon Press, Oxford).
2.
EyresN. R.HartreeD. R.InghamJ.JacksonR.SargantR. J.WagstaffJ. B.‘The calculation of variable heat flow in solids’, Phil. Trans. roy. Soc.1946–8, 240 (Series A), 1.
3.
PriceP. H.SlackM. R.‘The effect of latent heat on the numerical solution of the heat flow equation’, Brit. J. appl. Phys.1954, 5 (No. 8), 285.
4.
RichtmeyerR. D.Difference methods for initial value problems1957, chapters 1 and 6 (Interscience Publishers Inc.).
5.
CrankJ.NicolsonP.‘A practical method for numerical integration of solutions of partial differential equations of heat conduction type’, Proc. Camb. Phil. Soc.194743, 50.
6.
DusinberreG. M.‘A note on the implicit method for finite difference heat transfer calculations’, J. Heat Trans.196183C (No. 1), 92.
KellyK. K.‘Contributions to the data on theoretical metallurgy, XIII. High temperature heat content, heat capacity and entropy data for the elements and inorganic compounds’, Bur. Min. Bull.1960, 584.
11.
RyderJ. T.Building Research Station, Ministry of Technology, private communication.
