A Multiscale Finite Element Method for reaction-diffusion eigenproblems arising from neutronics

Abstract

We consider reaction-diffusion eigenproblems with oscillatory diffusion and reaction coefficients. The reaction coefficient magnitude is large: the corrector equation identified by periodic homogenization involves both the diffusion and the reaction operators. We study the numerical approximation of this problem using the Multiscale Finite Element Method (MsFEM). This now classical method is a finite element type method that performs a Galerkin approximation of the oscillatory problem on a specific, problem dependent, basis set. The basis functions are precomputed in an offline stage. Inspired by homogenization theory and using some filtering ideas, we show how to define these basis functions in order to obtain an efficient method. The comprehensive set of numerical experiments that we present, in periodic and non-periodic cases, for the scalar-valued version of the problem (which is then self-adjoint) and for the vector-valued version of the problem (which is then in general non self-adjoint), demonstrates the performance of the approach. Some theoretical arguments complement the numerical observations.