The computer programs for heat transfer analysis in a horizontal annular porous insulation layer

Two computor programs have been developed to analyze natural convective heat transfer in an annular porous insulation layer which has a horizontal axis. In the insulation layer, it is considered that the Darcy's law rules the fluid motion and the effect of density variation on natural convection is taken into account basis of Boussinesq approximation. The boundary condition at the outer surface of insulation layer is applicable to both constant temperature and given heat transfer conditions for the energy equation. In solving the equations the successive over relaxation (SOR) method was employed for both momentum and energy equations and up-wind difference-method was also applied to the convective term of energy equation to avoid the divergence of the solutions. The input data including some control parameters are briefly described with their formats and the output data are also explained. Elements of the programs including the main program and the fucntions of subroutines are described. The source program listing is appended. (author)