Skip to Main content Skip to Navigation
Conference papers

A Domain Decomposition Method Applied to Large Eigenvalue Problems in Neutron Physics

Maxime Barrault 1 Bruno Lathuilière 1 Pierre Ramet 2, 3 Jean Roman 2, 3
2 SCALAPPLIX - Algorithms and high performance computing for grand challenge applications
CNRS - Centre National de la Recherche Scientifique : UMR5800, École Nationale Supérieure d'Électronique, Informatique et Radiocommunications de Bordeaux (ENSEIRB), Inria Bordeaux - Sud-Ouest, Université Sciences et Technologies - Bordeaux 1, Université Bordeaux Segalen - Bordeaux 2
Abstract : The simulation of the neutron transport inside a nuclear reactor leads to the computation of the lowest eigen pair of a simplified transport operator. This computation is done by a power inverse algorithm accelerated by a Chebyshev polynomials based process. At each iteration, a large linear system is solved inexactly by a block Gauss-Seidel algorithm. For our applications, one Gauss-Seidel iteration is already sufficient to ensure the right convergence of the inverse power algorithm. For the approximate resolution of the linear system at each inverse power iteration, we propose a non overlapping domain decomposition based on the introduction of Lagrange multipliers in order to: - get a parallel algorithm, which allows to circumvent memory consumption problem and to reduce the computational time; - deal with different numerical approximations in each subdomain; - minimize the code modifications in our industrial solver. When the Chebyshev acceleration process is switched off, the method performs well on up to 100 processors for an industrial test case. It exhibits a good efficiency which allows us to realize some computations beyond the reach of standard workstations. Besides, we study the efficiency of the Chebyshev acceleration process in our domain decomposition method.
keyword : DDM
Document type :
Conference papers
Complete list of metadata

https://hal.inria.fr/inria-00346025
Contributor : Pierre Ramet <>
Submitted on : Wednesday, December 10, 2008 - 6:39:53 PM
Last modification on : Thursday, February 11, 2021 - 2:38:02 PM

Identifiers

  • HAL Id : inria-00346025, version 1

Collections

Citation

Maxime Barrault, Bruno Lathuilière, Pierre Ramet, Jean Roman. A Domain Decomposition Method Applied to Large Eigenvalue Problems in Neutron Physics. PMAA'08, 2008, Neuchatel, Switzerland. ⟨inria-00346025⟩

Share

Metrics

Record views

121