This paper will present a hybrid solution algorithm for the two-phase flow equations coupled to wall heat conduction. The partial differential equations in the physical model are the same as in RELAP5. The hybrid solution algorithm couples two solution methods, the solution method currently employed by RELAP5 and an implicitly balanced solution method. The resulting hybrid solution method is both fast and accurate. Results will be presented that show when accuracy and CPU time are considered simultaneously that there are ranges when the hybrid solution algorithm is preferred over the RELAP5 solution method.
Issue Section:
Research Paper
1.
Mousseau, V., 2004, “A Hybrid Solution Method for the Two-Phase Fluid Flow Equations,” ICONE12-49536, Proceedings of ICONE 12 2004 12th International Conference on Nuclear Engineering, ASME, New York.
2.
Mousseau
, V.
, 2004
, “Implicitly Balanced Solution of the Two-phase Flow Equations Coupled to Nonlinear heat Conduction
,” J. Comput. Phys.
, 200
, pp. 104
–132
.3.
The RELAP5 code development team, 2001, “RELAP5/MOD3.3 Code Manual Volume I: Code Structure, System Models, and Solution Methods,” NUREG/CR-5535 ed., U.S. Nuclear Regulatory Commission, Washington, D.C., (http://www.edasolutions.com/RELAP5/manuals/index.htm)
4.
The RELAP5 code development team, 2002, “RELAP5-3D Code Manual Volume I: Code Structure, System Models, and Solution Methods,” INEEL-EXT-98-00834 Rev. 2.0 ed., Idaho National Engineering and Environmental Laboratory, Idaho Falls, (http://www.inel.gov/relap5/r5manuals.htm)
5.
Frepoli
, C.
, Mahaffy
, J.
, and Ohkawa
, K.
, 2003
, “Notes on the Implementation of a Fully-Implicit Numerical Scheme for a Two-Phase Three-Field Flow Model
,” Nucl. Eng. Des.
, 225
, pp. 191
–217
.6.
Thurgood, M., and George, T., 1983, “COBRA/TRAC—A Thermal-Hydraulic Code for Transient Analysis of Nuclear Reactor Vessels and Primary Coolant System,” NUREG/CR-3046, Vols. 1–4, U.S. Nuclear Regulatory Commission, Washington D.C.
7.
Knoll
, D.
, and Keyes
, D.
, 2003
, “Jacobian-free Newton-Krylov Methods: A Survey of Approaches and Applications
,” J. Comput. Phys.
, 193
, pp. 357
–397
.8.
Knoll
, D.
, Chacon
, L.
, Margolin
, L.
, and Mousseau
, V.
, 2003
, “On Balanced Approximations for Time Integration of Multiple Time Scale Systems
,” J. Comput. Phys.
, 185
, pp. 583
–611
.9.
Knoll
, D.
, Mousseau
, V.
, Chacon
, L.
, and Reisner
, J.
, 2004
, “Jacobian-free Newton-Krylov Methods for Accurate Time Integration of Stiff Wave Systems
,” J. Sci. Comput.
, in press.10.
Chan
, T.
, and Jackson
, K.
, 1984
, “Nonlinearly Preconditioned Krylov Subspace Methods for Discrete Newton Algorithms
,” SIAM J. Sci. Comput. (USA)
, 5
, pp. 533
–542
.11.
Brown
, P.
, and Saad
, Y.
, 1990
, “Hybrid Krylov Methods for Nonlinear Systems of Equations
,” SIAM J. Sci. Stat. Comput.
, 11
, pp. 450
–481
.12.
Mousseau, V., 2004, “Transitioning from Interpretive to Predictive in Thermal Hydraulic Codes,” 125514, Proceedings of the International Meeting on Updates in Best Estimate Methods in Nuclear Installation Safety Analysis, American Nuclear Society, La Grange Park, pp. 44–51.
13.
Dembo
, R.
, Eisenstat
, S.
, and Steihaug
, T.
, 1982
, “Inexact Newton Methods
,” SIAM J. Num. Anal.
, 19
, pp. 400
–408
.14.
Saad
, Y.
, and Schultz
, M.
, 1986
, “GMRES: A Generalized Minimal Residual Algorithm for Solving Non-Symmetric Linear Systems
,” SIAM J. Sci. Stat. Comput.
, 7
, pp. 856
–869
.15.
Brown
, P.
, and Hindmarsh
, A.
, 1986
, “Matrix-free Methods for Stiff Systems of ODE’s
,” SIAM J. Num. Anal.
, 23
, pp. 610
–638
.16.
Mortensen, G., 2004, personal communication, member of the NRC’s RELAP5 code development team at ISL.
17.
Riemke, R., 2004, personal communication, member of the DOE’s RELAP5-3D code development team at INEEL.
You do not currently have access to this content.
