Yuanzineng kexue jishu (Feb 2022)
Simultaneous Solution for Whole High Temperature Gas-cooled Reactor Coupled System: Method Research and Code Development
Abstract
The numerical calculation for nuclear reactor coupled system is the fusion, integration and upgrade of the existing numerical technologies in various related fields. Generally, the whole nuclear power plant has several different coupling mechanisms inherently, such as the multiphysics coupling, multiscale coupling and multisystem coupling, which is an ultralargescale strong nonlinear coupling system. Moreover, each coupling mechanism has its own mathematic features. The combination of different coupling mechanisms is such complicated, which is still an open issue for nuclear engineering community. Simultaneous solution methods, such as Jacobianfree Newton Krylov (JFNK) method and NewtonKrylov (NK) method, are a promising choice for effectively solving such complex nonlinear largescale system. In the first part of this work, the recent progress of numerical calculation for nuclear coupling system worldwide is reviewed, especially, several existing coupling methods are summarized and compared. In the second part, it focuses on the development of simultaneous solution methods as well as the coupling platform for high temperature gas-cooled reactor (HTR) full coupled system by Institute of Nuclear and New Energy Technology (INET), Tsinghua University. Due to the multiphysics, multiscale, multicomponent, multicircuit and multimodule coupling features in HTR, the key technologies, such as the JFNK method with nonlinearelimination, scaling factor and steam generator model, were presented and discussed. In detail, the JFNK method with nonlinearelimination was firstly proposed to deal with the multiscale temperature coupling between the global pebblebed and the local fuel sphere. The scaling factor technology was used to transform the illposed Jacobian matrix into the wellposed one in the multiphysical coupling issue. For the steam generator model, the dynamical residual method is conducted for the twophase flow in the second circuit to avoid the multisolution issue when the phase generation or the phase disappearance occurs. Moreover, a unified framework was proposed to tackle the multi-level coupling structure in HTR system, after several attempts were made. In this framework, the whole coupling system is gradually decomposed into systemlevel issue, componentlevel issue and physical fieldlevel issue, and solved as a whole system by the simultaneous solution method. However, to completely solve the whole HTR nuclear power plant system, there are still many works to do, such as the fundamental algorithm research, program framework development and engineering applications.