A continuous pseudo spectral method to compute the dominant Lambda modes of a nuclear power reactor core using a triangular mesh is studied. This method is based on the expansion of the neutronic flux on each node of the mesh in terms of the modified Dubiner's polynomials. Some coefficients of the expansion for the different nodes are selected in such a way that the solution for the neutronic flux is continuous by construction. Continuity conditions for the neutronic current are imposed on a set of collocation points on the faces of the triangles and on the internal vertices of the mesh, and the neutrons balance in each triangular node is taken into account by using an integral form of the diffusion equation. The performance of the proposed method has been tested obtaining the dominant modes of the benchmark problem IAEA reactor with reflector.
展开▼
