Nuclear Science and Techniques

《核技术》(英文版) ISSN 1001-8042 CN 31-1559/TL     2019 Impact factor 1.556

Nuclear Science and Techniques ›› 2018, Vol. 29 ›› Issue (8): 114 doi: 10.1007/s41365-018-0456-3

• NUCLEAR ENERGY SCIENCE AND ENGINEERING • Previous Articles     Next Articles

High-order fully implicit SIMPLE-based model for fully implicit simulation of upward two-phase flow

A. Hajizadeh 1 • H. Kazeminejad 2 • S. Talebi 3   

  1. 1 Research School of Reactor, Nuclear Science and Technology Research Institute, P.O. Box 11365-3486, Tehran, Islamic Republic of Iran
    2 Research School of Radiation Applications, Nuclear Science and Technology Research Institute, P.O. Box 11365-3486, Tehran, Islamic Republic of Iran
    3 Department of Energy Engineering and Physics, Amirkabir University of Technology (Tehran Polytechnic), 424 Hafez Avenue, P.O. Box 15875-4413, Tehran, Islamic Republic of Iran
  • Contact: S. Talebi E-mail:sa.talebi@aut.ac.ir
PDF ShareIt Export Citation
A. Hajizadeh, H. Kazeminejad, S. Talebi. High-order fully implicit SIMPLE-based model for fully implicit simulation of upward two-phase flow.Nuclear Science and Techniques, 2018, 29(8): 114     doi: 10.1007/s41365-018-0456-3

Abstract:

The drift-flux model has a practical importance in two-phase flow analysis. In this study, a finite volume solution is developed for a transient four-equation drift-flux model through the staggered mesh, leading to the development of a fully implicit discretization method. The main advantage of the fully implicit method is its unconditional stability. Newton’s scheme is a popular method of choice for the solution of a nonlinear system of equations arising from fully implicit discretization of field equations. However, the lack of convergence robustness and the construction of Jacobian matrix have created several difficulties for the researchers. In this paper, a fully implicit model is developed based on the SIMPLE algorithm for two-phase flow simulations. The drawbacks of Newton’s method are avoided in the developed model. Different limiter functions are considered, and the stabilized method is developed under steady and transient conditions. The results obtained by the numerical modeling are in good agreement with the experimental data. As expected, the results prove that the developed model is not restricted by any stability limit.

Key words: Fully implicit, Two-phase flow, Drift-flux model, Pressure-based algorithm