Abstract
A semi-iterative method based on a nested application of Flexible Generalized Minimum Residual)FGMRES) was developed to solve the linear systems resulting from the application of the discretized two-phase hydrodynamics equations to nuclear reactor transient problems. The complex three-dimensional reactor problem is decomposed into simpler, more manageable problems which are then recombined sequentially by GMRES algorithms. Mathematically, the method consists of using an “inner” level GMRES to solve the preconditioner equation for an “outer” level GMRES. Applications were performed on practical, three-dimensional models of operating Pressurized Water Reactors (PWR). Serial and parallel applications were performed for a reactor model with two different details in the core representation. When appropriately tight convergence was enforced at each GMRES level, the results of the semi-iterative solver were in agreement with existing direct solution methods. For the larger model tested, the serial performance of GMRES was about a factor of 3 better than the direct solver and the parallel speedups were about 4 using 13 processors of the INTEL Paragon. Thus, for the larger problem over an order of magnitude reduction in the execution time was achieved indicating that the use of semi-iterative solvers and parallel computing can considerably reduce the computational load for practical PWR transient calculations.
Similar content being viewed by others
REFERENCES
Bank, R., and Rose, D. (1977). Marching Algorithms for Elliptic Boundary Value Problems. I: The Constant Coefficient Case, SIAM J. Numer. Anal. 14, 792-829.
de Sturler, E., and van der Vorst, H. A. (1993). Reducing the Effect of Global Communication in GMRES(m) and CG on Parallel Distributed Memory Computers, Tech. Report Preprint 832, Department of Mathematics, University Utrecht, Utrecht.
Duff, T. S., Erisman, A. M., Reid, J. K. (1986). Direct Methods for Sparse Matrices, Clarendon Press, New York.
Dutto, L. C. (1993). The Effect of Ordering on Preconditioned GMRES Algorithm, for Solving the Compressible Navier-Stokes Equations, Int'l. J. Numer. Methods Engng. 36, 457-497.
Hoffman, J. D. (1992). Numerical Methods for Engineers and Scientists, Mc-Graw Hill, New York.
Kavenoky, A. (1992). Supercomputing applied to nuclear reactors, Ann. Nucl. Energy, 19, 679-683.
McFadden, J., Paulsen, M. P., Gose, G. C., Peterson, C. E., McClure, J. A., Westacott, P. J. L. (1992). RETRAN-03: A Program for Transient Thermal-Hydraulic Analysis of Complex Fluid Flow Systems.
Murray, P. E. (1994). An Assessment of Advanced Numerical Methods for Two-Phase Fluid Flow, EGG-NRE-11262, Idaho National Engineering Laboratory.
Saad, Y. (1993). A Flexible InnerOuter Preconditioned GMRES Algorithm, SIAM J. Sci. Statist. Comput. 14, xx-xx.
Saad, Y., and Schultz, M. (1986). GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems, SIAM J. Sci. Statist. Comput. 7, xx-xx.
Shadid, J. N., and Tuminaro, R. S. (1994). A Comparison of Preconditioned Nonsymmetric Krylov Methods on a Large-Scale MIMD Machine, SIAM J. Sci. Comp. 15, 440-459.
Sonneveld, P. (1989). CGS, A Fast Lanczos-Type Solver for Nonsymmetric Linear Systems, SIAM J. Sci. Statist. Comput. 10, 36-52.
Turner, J. A., and Doster, J. (1991). Performance of conjugate gradient-like algorithms in transient two-phase subchannel analysis, Nuclear Thermal Hydraulics, ANS Winter Meeting.
van der Vorst, H., and Vuik, C. (1994). GMRESR: A family of nested GMRES methods, Numer. Linear Algebra Appl. 1, 369-386.
van der Vorst, H. A. (1992). Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems, SIAM J. Sci. Statist. Comput. 13, 631-644.
Woodruff, S. B. (1993). Some computational challenges of developing efficient parallel algorithms for data-dependent computation in thermal-hydraulics supercomputer applications. In High Performance Fortran Forum, pp. 1615-1621.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Janardhan, R., Downar, T. A Nested FGMRES Method for Parallel Calculation of Nuclear Reactor Transients. Journal of Scientific Computing 13, 65–93 (1998). https://doi.org/10.1023/A:1023256711255
Issue Date:
DOI: https://doi.org/10.1023/A:1023256711255