Abstract
A method to treat the internal boundary condition in an optical diffusion calculation is proposed and is compared with the conventional methods. One of the existing internal boundary conditions is Haskel's method, which uses the effective reflection coefficient for partial currents. However, Haskel's method ignores incoming partial currents from the adjacent mesh in its derivation. As a result, the accuracy at the internal boundary is lower. This paper proposes a method to improve the accuracy by iteratively updating the effective reflection coefficient for partial current. The proposed method was applied to the benchmark calculations on a one-dimensional slab geometry and its accuracy was confirmed by comparing it with the reference solution obtained by the Monte Carlo code MCML, along with the previously proposed Haskel's method and Aronson's method. As a result, it was confirmed that the proposed method is more accurate than Haskel's method at the internal boundary and gives the same result as Aronson's method. The convergence of the effective reflection coefficient using iterative calculations in the proposed method was good.