High Order Finite Element Method for the Lambda modes problem on hexagonal geometry

A High Order Finite Element Method to approximate the Lambda modes problem for reactors with hexagonal geometry has been developed. This method is based on the expansion of the neutron flux in terms of the modified Dubiner’s polynomials on a triangular mesh. This mesh is fixed and the accuracy of the method is improved increasing the degree of the polynomial expansions without the necessity of remeshing. The performance of method has been tested obtaining the dominant Lambda modes of different 2D reactor benchmark problems.