A solution of the multi-group multi-regional diffusion equation in two-dimensional cylindrical (ρ-z) geometry is obtained in the form of a regionwise double series com-posed of Bessel and trigonometrical functions.
The diffusion equation is multiplied by weighting functions, which satisfy the homogeneous part of the diffusion equation, and the products are integrated over the region for obtaining the equations to determine the fluxes and their normal derivatives at the region boundaries. Multiplying the diffusion equation by each function of the set used for the flux expansion, then integrating the products, the coefficients of the double series of the flux inside each region are calculated using the boundary values obtained above.
Since the convergence of the series thus obtained is slow especially near the region boundaries, a method for improving the convergence has been developed. The double series of the flux is separated into two parts. The normal derivative at the region boundary of the first part is zero, and that of the second part takes the value which is obtained in the first stage of this method. The second part is replaced by a continuous function, and the flux is represented by the sum of the continuous function and the double series.
A sample critical problem of a two-group two-region system is numerically studied. The results show that the present method yields very accurately the flux integrals in each region with only a small number of expansion terms.
抄録全体を表示