Sensitivity analysis of numerical schemes in natural cooling flows for low power research reactors

Abstract

The advantages of using natural circulation (NC) as a cooling system, has prompted the worldwide development to investigate this phenomenon more than before. The interesting application of the NC in low power experimental facilities and research reactors, highlights the obligation of study in these laminar flows. The inherent oscillations of NC between hot source and cold sink in low Grashof numbers necessitates stability analysis of cooling flow with experimental or numerical schemes. For this type of analysis, numerical methods could be implemented to desired mass, momentum and energy equations as an efficient instrument for predicting the behavior of the flow field. In this work, using the explicit, implicit and Crank-Nicolson methods, the fluid flow parameters in a natural circulation experimental test loop are obtained and the sensitivity of solving approaches are discussed. In this way, at first, the steady state and transient results from explicit are obtained and compared with experimental data. The implicit and crank-Nicolson scheme is investigated in next steps and in subsequent this research is focused on the numerical aspects of instability prediction for these schemes. In the following, the assessment of the flow behavior with coarse and fine mesh sizes and time-steps has been reported and the numerical schemes convergence are compared. For more detail research, the natural circulation of fluid was modeled by ANSYS-CFX software and results for the experimental loop are shown. Finally, the stability map for rectangular closed loop was obtained with employing the Nyquist criterion.