The Homogeneous B₁ Model as Polynomial Eigenvalue Problem

The homogeneous version of the B ₁ leakage model is a non-linear eigenvalue problem which is generally solved iteratively by a root-finding algorithm, combined to the supplementary eigenvalue problem of the multiplication factor. This problem is widely used for ordinary cross section preparation in reactor analysis. Our work approximates this problem with a polynomial eigenvalue problem, which can be easily transformed into an ordinary linear generalized eigenproblem of size equal to the initial one multiplied by the polynomial degree used for the approximation of a transcendental function. This procedure avoids recurring to numerical root-finding methods supported by extra eigenvalue problems. The solution of the fundamental buckling with increasing approximation order is compared to the reference value obtained by inverse iterations.