Computational algorithms and their implementations are studied for the treatments of kernel function’s shadow effects, which occurs in the application of the BEM(boundary element method) for radiation heat transfer analysis. In this problem, surfaces are assumed to be diffuse and gray placed in two-dimensional enclosures with transparent medium. Self-blocking or third party blocking is possible when the radiosity from the BEM nodes cannot reach the destination points. Also the third party blocking can be the partial or total blocking case. Algorithms, which can accurately recognize the each blocking cases and reflect the shadow effects to the BEM kernel function, are studied in this paper. Effective implementation methods are presented, and their results are verified by the test problem