0
RESEARCH PAPER

A Fully Implicit Hybrid Solution Method for a Two-Phase Thermal-Hydraulic Model

[+] Author and Article Information
Vincent A. Mousseau

Fluid Dynamics Group, T-3, M.S. B216, Los Alamos National Laboratory, Los Alamos, New Mexico 87545e-mail: vmss@lanl.gov

J. Heat Transfer 127(5), 531-539 (May 25, 2005) (9 pages) doi:10.1115/1.1865223 History: Received April 27, 2004; Revised December 23, 2004; Online May 25, 2005
Your Session has timed out. Please sign back in to continue.

References

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.
Mousseau,  V., 2004, “Implicitly Balanced Solution of the Two-phase Flow Equations Coupled to Nonlinear heat Conduction,” J. Comput. Phys., 200, pp. 104–132.
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)
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)
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.
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.
Knoll,  D., and Keyes,  D., 2003, “Jacobian-free Newton-Krylov Methods: A Survey of Approaches and Applications,” J. Comput. Phys., 193, pp. 357–397.
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.
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.
Chan,  T., and Jackson,  K., 1984, “Nonlinearly Preconditioned Krylov Subspace Methods for Discrete Newton Algorithms,” SIAM J. Sci. Comput. (USA), 5, pp. 533–542.
Brown,  P., and Saad,  Y., 1990, “Hybrid Krylov Methods for Nonlinear Systems of Equations,” SIAM J. Sci. Stat. Comput., 11, pp. 450–481.
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.
Dembo,  R., Eisenstat,  S., and Steihaug,  T., 1982, “Inexact Newton Methods,” SIAM J. Num. Anal., 19, pp. 400–408.
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.
Brown,  P., and Hindmarsh,  A., 1986, “Matrix-free Methods for Stiff Systems of ODE’s,” SIAM J. Num. Anal., 23, pp. 610–638.
Mortensen, G., 2004, personal communication, member of the NRC’s RELAP5 code development team at ISL.
Riemke, R., 2004, personal communication, member of the DOE’s RELAP5-3D code development team at INEEL.

Figures

Grahic Jump Location
Schematic diagram of the computational domain
Grahic Jump Location
Interfacial area versus vapor volume fraction
Grahic Jump Location
Liquid (solid) and vapor (dashed) wall area versus vapor volume fraction
Grahic Jump Location
Interfacial friction versus vapor volume fraction
Grahic Jump Location
Contact interfacial heat transfer (not associated with phase change) versus vapor volume fraction
Grahic Jump Location
Peak-clad temperature versus time for the high (dashed) and low (solid) velocity test problems
Grahic Jump Location
Vapor volume fraction versus time for the constant closures (dashed) and the varaible closures (solid) at a location near the center of the domain
Grahic Jump Location
Vapor velocity versus time for the constant closures (dashed) and the variable closures (solid) at a location near the center of the domain
Grahic Jump Location
Vapor volume fraction versus distance for the constant closures (dashed) and the variable closures (solid), thin lines initial, thick lines final
Grahic Jump Location
Vapor velocity versus distance for the constant closures (dashed) and the variable closures (solid), thin lines initial, thick lines final
Grahic Jump Location
High velocity test problem, error versus CFL number, squares hybrid, circles OSSI
Grahic Jump Location
High velocity test problem, error versus CPU time, squares hybrid, circles OSSI
Grahic Jump Location
Low velocity test problem, error versus CFL number, squares hybrid, circles OSSI
Grahic Jump Location
Low velocity test problem, error versus CPU time, squares hybrid, circles OSSI

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In