A numerical study of boiling flow instability of a reactor thermosyphon system

A numerical study has been carried out to investigate the boiling flow instability of a reactor thermosyphon system. The numerical model solves the conservation equations of mass, momentum and energy applicable to a two-fluid and three-field steam–water system using a finite difference technique. The computer code MONA was used for this purpose. The code was applied to the thermosyphon system of an EO (ethylene oxide) chemical reactor in which the heat released by a catalytic reaction is carried by boiling water under natural circulation conditions. The steady-state characteristics of the reactor thermosyphon system were predicted using the MONA code and conventional two-phase flow models in order to understand the model applicability for this type of thermosyphon system. The two-fluid model was found to predict the flow closest to the measured value of the plant. The stability behaviour of the thermosyphon system was investigated for a wide range of operating conditions. The effects of power, subcooling, riser length and riser diameter on the boiling flow instability were determined. The system was found to be unstable at higher power conditions which is typical for a Type II instability. However, with an increase in riser diameter, oscillations at low power were observed as well. These are classified as Type I instabilities. Stability maps were predicted for both Type I and Type II instabilities. Methods of improving the stability of the system are discussed.