A Lagrangian-Lagrangian model is developed using an in-house code to simulate bubble trajectory in two-phase bubbly flow around circular cylinder. Random Vortex Method (RVM) which is a Lagrangian approach is used for solving the liquid phase. The significance of RVM relative to other RANS/LES methods is its capability in directly modelling the turbulence. In RVM, turbulence is modeled by solving the vorticity transport equation and there is no need to use turbulence closure models. Another advantage of RVM relative to other CFD approaches is its independence from mesh generation. For the bubbles trajectory, equation of motion of bubbles which takes into account effect of different forces are coupled with the RVM. Comparison of the results obtained from current model with the experimental data confirms the validity of the model. Effect of different parameters including flow Reynolds number, bubble diameter and injection point on the bubbles' trajectory are investigated. Results show that increase in the Reynolds number reduces the rising velocity of the bubbles. Similar behavior is observed for the bubbles when their diameter was decreased. According to the analysis carried out, present Lagrangian-Lagrangian model solves the issues of mesh generation and turbulence modelling which exist in common two phase flow modelling schemes.
151*Corresponding author email address: shadema@uwindsor.ca along the wall at higher Stokes number, while at lower Stokes numbers a Gaussian-type profiles in the boundary layer was seen. In addition, three types of bubble trajectories were found in the boundary layer including sliding bubbles, bouncing bubbles and free-dispersion bubbles.In the field of numerical simulation, different methods have been developed including Eulerian-Eulerian, Eulerian-Lagrangian and Lagrangian-Lagrangian models to simulate various types of two phase bubbly flows. In Eulerian-Eulerian model one set of equations are solved for each phase which includes mass, momentum and energy equations. Basically, in this model there is no difference between the continuous and discrete phases and both phases are treated as continuous phases [5]. This model computes the interaction between two phases, however, it is computationally expensive. This method was widely used for unsteady gas-liquid flows in bubble columns. It was observed that this approach can predict the timeaveraged flow properties in good agreement with experiments [6,7,8]. In Eulerian-Lagrangian model, liquid phase is treated as continuum phase while bubbles are assumed as dispersed phase. Motion of continuous phase is investigated in Eulerian coordinate while bubbles motion is tracked in Lagrangian coordinate [9,10,11,12,13]. This model is also used for liquid-solid two phase flow problems. Although the models based on the Eulerian-Lagrangian approach are capable of predicting the time-averaged properties, the ability of these models to predict dynamic characteristics is not yet fully understood [14]. In Lagrangian-Lagrangian model which is the subject of c...