In this study, we perform direct numerical simulations of horizontal turbulent channel flows, injecting multiple millimeter-sized bubbles to examine their dynamics Building upon our previous research on bubble dynamics in horizontal turbulent flows [1], the current study extends this analysis to the interaction with riblet surfaces and varying volumes of bubbles. We employ a Volume of Fluid (VOF) model to simulate the two immiscible fluids, water, and air, using the isoAdvector method[2] for interface capturing. the numerical results effectively show to suppress numerical diffusion at the interface, even during extended simulation periods. We established the simulation setup based on prior research[3], and the domain configurations were as follows: length (L) = 100 mm, width (W) = 70 mm, and a total average height (H) = 10 mm, accounting for the varying height due to the riblet shapes, with a flow velocity of 1 m/s. Riblet shapes were applied to the upper wall surface and varied along the span direction with configurations of 2 and 4 riblets, with incline angles of 4° and 8°, respectively. To maintain a consistent flow rate comparable to previous studies on flat wall surfaces, the riblet inclinations were calibrated based on their mid-height. This setup offers a valid basis for comparison with previous flat wall channel simulations. As illustrated in Figure 1, the bubbles tended to move to the topside of the riblet, hindering their movement outside of the riblet center. Additionally, bubbles were distributed in accordance with the angle and number of riblet.