CTF-R: A novel residual-based thermal hydraulic solver

Abstract The traditional scientific process has been revolutionized by the advent of computational modeling, but the nuclear industry generally uses “legacy codes,” which were developed early in the evolution of computers. One example of a legacy code, the thermal hydraulic subchannel code CTF, is modernized in this work through the development of a novel residual-based version, CTF-R. Unlike its predecessor, CTF-R is not limited by the strict computational limitations of the early 1980’s, and can therefore be designed such that it is inherently flexible and easy to use. A case study is examined to demonstrate how the flexibility of the code can be used to improve simulation results. In this example, the coupling between the solid and liquid fields is examined. Traditionally, this coupling is modeled explicitly, which imposes numerical stability limits on the time step size. These limits are derived and it is shown that they are removed when the coupling is made implicit. Further, the development of CTF-R will enable future improvements in next generation reactor modeling, numerical methods, and coupling to other codes. Through the further development of CTF-R and other residual-based codes, state-of-the-art simulation is possible.

[1]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[2]  J. Neumann,et al.  Numerical Integration of the Barotropic Vorticity Equation , 1950 .

[3]  S. T. Zalesak,et al.  High order “ZIP” differencing of convective terms , 1981 .

[4]  William A. Wieselquist,et al.  VERA Core Simulator Methodology for Pressurized Water Reactor Cycle Depletion , 2017 .

[5]  F. Dittus,et al.  Heat transfer in automobile radiators of the tubular type , 1930 .

[6]  Robert K. Salko,et al.  Implementation and assessment of wall friction models for LWR core analysis , 2018 .

[7]  Ephraim M Sparrow,et al.  Heat Transfer to Longitudinal Laminar Flow Between Cylinders , 1961 .

[8]  Vincent A. Mousseau,et al.  Implicitly balanced solution of the two-phase flow equations coupled to nonlinear heat conduction , 2004 .

[9]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[10]  Nathan Wayne Porter.  Development of a Novel Residual Formulation of CTF and Application of Parameter Estimation Techniques. , 2018 .

[11]  Roger P. Pawlowski,et al.  The Virtual Environment for Reactor Applications (VERA): Design and architecture☆ , 2016, J. Comput. Phys..

[12]  N. Sylvester,et al.  Explicit approximations to the solution of Colebrook's friction factor equation , 1982 .

[13]  H. H. Rosenbrock,et al.  Some general implicit processes for the numerical solution of differential equations , 1963, Comput. J..

[14]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[15]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[16]  Christopher Alan Dances.  Initial Residual Formulation of CTF , 2015 .

[17]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .