A novel Two-way coupled Euler-Lagrange mode, including bubble-bubble collision, coalescence, with variable bubble radius, and bubbles breakup was applied to simulate air-water bubbly flow in vertical pipes. This approach uses the Continuous Random Walks CRW models for creating the velocity fluctuations according to the given state of turbulent kinetic energy and Dissipation rate at the location of the bubbles which is solved by the turbulence model. This dissertation i) describes the development of the Euler-Lagrange Approach under study, ii) presents a study for the two-way coupling effect on both the continuous and dispersed phases properties, iii) studies the effect of both the lift force coefficient and bubble induced turbulence BIT relations on the gas void fraction distributions, iv) studies the effect of the bubbles coalescence and breakup on bubble sizes and gas void fraction distributions. And presents the results of the simulations performed under each of these considerations.The two-way coupling process takes the effect of the dispersed phase on the continuous one through inserting source terms in the conservation equations of momentum and Turbulence. Also it modifies the volume of the computational cells in the Euler solver available for the continuous phase according to the void fraction of each cell. During the two-way coupling process, some studies needed to be performed like the adjustment of the lift force coefficient and the relation due to the change of the liquid velocity profiles as a result of the two-way coupling.The bubble-bubble collision was applied in the two-way coupling process. It was found that considering the collision appears on the void fraction distribution only in the high velocity and high gas holdup cases with very small effect. The bubble-bubble coalescence was applied as a complementary part of the collision process using the film drainage model of Chesters (1991). This model compares the film drainage time with the contact time to calculate the coalescence efficiency. The coalescence model was tested first before applying the breakup, so the bubbles size increased only and this affected on the void fraction distribution badly. Then the breakup model of Martínez-Bazán (1999a, b) was applied to perform the equilibrium in the bubble sizes. These two processes of the coalescence and breakup found to consume long computational time, the reason that did not give us a chance for testing many cases with considering both coalescence and breakup.The main investigation point through the development of this work was the BIT kinetic energy term and its effect on the used CRW model. This was considered in nearly every phase of the model development to study the effect of the BIT under the different considerations of the bubbles interaction mechanisms. It could be concluded a final vi expression for the BIT relation that is used successfully with the CRW models under consideration.vii
ResumenUna nueva aproximación euleriana-lagarangiana, en su forma de acople en dos vías, para ...