The computer programs for analysing two-dimensional permeation flow within horizontal thermal insulation

A computer program has been developed to analyze the behavior of the permeation flow occurring when partition plates are placed in a thermal insulation layer within a hot gas duct of a high-temperature gas-cooled reactor. This computer program simulates a single thermal insulation layer devided by partition plates and treats this layer as two-dimensional and homogeneous. In the thermal insulation layer, it is assumed that the D'arcy's law rules the fluid motion and also the temperatures of local gas and insulation are equal. In calculation the extrapolated Liebman method was employed for both momentum and energy equations and a centered difference method was applied to the equations to avoid the divergence of the solutions. Convergence of numerical solutions and stability of computation were successfully obtained when the state of convection within thermal insulation layer is in steady or quasi-steady states where forced convection is more dominant than natural convection. The input and output data are briefly described with their formats, and the elements of the program including the main program and the functions of the subroutines are appended. (author)