This study mainly focused on the neutronics modeling of bubbles in bubbly flow in boiling water reactors. The bubble, ring and homogenous models were used for radial void fraction distribution. Effect of the bubble and ring models on the infinite multiplication factor and two-group flux distribution was investigated by comparing with the homogenous model. Square pitch unit cell geometry was used in the calculations. In the bubble model, spherical and non-spherical bubbles at random positions, sizes and shapes were produced by Monte Carlo method. The results show that there are significant differences among the proposed models from the viewpoint of physical interaction mechanism. For the fully-developed bubbly flow, k(inf) is overestimated in the ring model by about 720 +/- 6 pcm with respect to homogeneous model whereas underestimated in the bubble model by about -65 +/- 9 pcm with a standard deviation of 15 pcm. In addition, the ring model shows that the coolant must be separated into regions to properly represent the radial void distribution. Deviations in flux distributions principally occur in certain regions, such as corners. As a result, the bubble model in modeling the void fraction can be used in nuclear engineering calculations. (C) 2019 Korean Nuclear Society, Published by Elsevier Korea LLC.