Multiple size group （MUSIG） model combined with a threedimensional twofluid model were em ployed to predict subcooled boiling flow of liquid nitrogen in a vertical upward tube. Based on the mechanism of boiling heat transfer, some important bubble model parameters were amended to be applicable to the modeling of liquid nitrogen. The distribution of different discrete bubble classes was demonstrated numerically and the distribu tion patterns of void fraction in the wallheated tube were analyzed. It was found that the average void fraction in creases nonlinearly along the axial direction with wall heat flux and it decreases with inlet mass flow rate and sub cooled temperature. The local void fraction exhibited a Ushape distribution in the radial direction. The partition of the wall heat flux along the tube was obtained. The results showed that heat flux consumed on evaporation is the leading part of surface heat transfer at the rear region of subcooled boiling. The turning point in the pressure drop curve reflects the instability of bubbly flow. Good agreement was achieved on the local heat transfer coefficient aalnst experimental measurements, which demonstrated the accuracy of the numerical model.
