Fast non-overlapping Schwarz domain decomposition methods for solving the neutron diffusion equation

  • Authors:
  • Erell Jamelot;Patrick Ciarlet, Jr

  • Affiliations:
  • Commissariat í l'ínergie Atomique et aux ínergie Alternatives, CEA Saclay, 91191 Gif-sur-Yvette Cedex, France;POEMS Laboratory, CNRS-INRIA-ENSTA UMR 7231, ENSTA ParisTech 32, Boulevard Victor, 75739 Paris Cedex 15, France

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2013

Quantified Score

Hi-index 31.45

Visualization

Abstract

Studying numerically the steady state of a nuclear core reactor is expensive, in terms of memory storage and computational time. In order to address both requirements, one can use a domain decomposition method, implemented on a parallel computer. We present here such a method for the mixed neutron diffusion equations, discretized with Raviart-Thomas-Nedelec finite elements. This method is based on the Schwarz iterative algorithm with Robin interface conditions to handle communications. We analyse this method from the continuous point of view to the discrete point of view, and we give some numerical results in a realistic highly heterogeneous 3D configuration. Computations are carried out with the MINOS solver of the APOLLO3(R) neutronics code.