Solution of the two dimensional diffusion and transport equations in a rectangular lattice with an elliptical fuel element using Fourier transform methods: One and two group cases

A method for solving the diffusion and transport equations in a rectangular lattice cell with an elliptical fuel element has been developed using a Fourier expansion of the neutron flux. The method is applied to a one group model with a source in the moderator. The cell flux is obtained and also the associated disadvantage factor. In addition to the one speed case, we also consider the two group equations in the cell which now become an eigenvalue problem for the lattice multiplication factor. The method of solution relies upon an efficient procedure to solve a large set of simultaneous linear equations and for this we use the IMSL library routines. Our method is compared with the results from a finite element code. The main drawback of the problem arises from the very large number of terms required in the Fourier series which taxes the storage and speed of the computer. Nevertheless, useful solutions are obtained in geometries that would normally require the use of finite element or analogous methods, for this reason the Fourier method is useful for comparison with that type of numerical approach. Extension of the method to more intricate fuel shapes, such as stars and cruciforms as well as superpositions of these, is possible.