EPJ Web Conf.
Volume 247, 2021PHYSOR2020 – International Conference on Physics of Reactors: Transition to a Scalable Nuclear Future
|Number of page(s)||8|
|Section||Advanced Modelling and Simulation|
|Published online||22 February 2021|
NUMERICAL OPTIMIZATION OF A MULTIPHYSICS CALCULATION SCHEME
1 DEN, Service d’études des réacteurs et de mathématiques appliquées (SERMA), CEA, Université Paris-Saclay F-91191 Gif-sur-Yvette, France
2 LPSC - IN2P3-CNRS/UGA/Grenoble INP 53 rue des Martyrs, F-38026 Grenoble Cedex, France
3 DEN, Service de thermo-hydraulique et de mécanique des fluides (STMF), CEA, Université Paris-Saclay F-91191 Gif-sur-Yvette, France.
Published online: 22 February 2021
This work concerns the numerical optimization of a multiphysics calculation scheme. The considered application is a 5x5 Pressurized Water Reactor (PWR) assemblies mini-core surrounded by radial and axial reflectors. The scenario adopted for the analysis is steady-state nominal conditions and fission products set to the equilibrium concentration. The neutronics is modelled at the pin-cell scale and the thermal-hydraulics at the subchannel level. Depending on the scenario, the damped fixed-point algorithm might not be sufficiently robust or efficient enough. For this reason, a technique based on the partial convergence of the solvers is tested. In every multiphysic iteration, a maximum number of iterations is imposed for both the neutronics and the thermal-hydraulics solvers. In combination with that, the solver restarts from the results of the last calculation. In this way, if the method is convergent, the initialization progresses towards the fixed-point solution. The results show that such a technique improves both the robustness and the speed of the algorithm. Within this approach, the range of relaxation factors that makes the algorithm converge is significantly broadened and the importance of this parameter on the global performance is reduced. The computing time also decreases by a factor between 10 and 20. Furthermore, this gain does not strongly depend on the exact imposed maximum number of iterations. Some preliminary observations are also reported in respect with the application of such a technique to the Anderson acceleration method.
Key words: MULTIPHYSICS / PARTIAL-CONVERGENCE / APOLLO3® / FLICA4
© The Authors, published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.