Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model

Mohamed-Kamel Riahi (1) , Anne-Marie A.-M. Baudron (1) , Jean-Jacques Lautard (1) , Yvon Maday (1, 2, 3) , Julien Salomon (1, 4)

1 LLPR - Laboratoire de Logiciels pour la Physique des Réacteurs
2 LJLL - Laboratoire Jacques-Louis Lions
3 DAM - Division of Applied Mathematics
4 CEREMADE - CEntre de REcherches en MAthématiques de la DEcision

We propose in this paper a parareal in time algorithm for the resolution of the Langenbuch-Maurer-Werner (LMW) benchmark~\citep{LMWref} neutron diffusion transient model. The parallelization in time is prompt thanks to a the coarse solver, which is based on a stable coarse time step with a steady control rods model. The finite element implementation provides a good scalability of the algorithm, where the time discretization is based on a Euler implicit scheme. Numerical results show the efficiency of the parareal method on large light water reactor transient model.

Mots clés: kinetic of neutrons; parallel solver; parareal in time algorithm

Domaines: Equations aux dérivées partielles [math.AP]; Physique mathématique [math-ph]; Analyse numérique [math.NA]

Mohamed-Kamel Riahi, Anne-Marie A.-M. Baudron, Jean-Jacques Lautard, Yvon Maday, Julien Salomon. Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model. 2011. ⟨hal-00722802v2⟩