---

In this article, two-phase slug flow is simulated numerically in a horizontal duct with rectangular cross-section using Volume Of Fluid (VOF) method. Conservation equations of mass, momentum and advection equation are solved in open source OpenFOAM code accompanying k-ω SST turbulence equations. Simulation is conducted based on the experimental results in the duct with rectangular cross-section. The results shows, due to Kelvin-Helmholtz (K-H) instability criteria slug initiation forms in the air-water interface during three dimensional turbulence modeling. Water level was increased slightly at interface in both numerical simulation and experiment. This level increase satisfies the K-H instability to generate a slug at interface. During slug initiation, the pressure behind slug is increased significantly. Big pressure gradient at the beginning of the slug in compare to the end of it causes the slug length to be increased as propagate along the duct. The numerical simulation of present research is capable of predicting the slug length accurately in accordance with experiment; however, the slug position with 22% inaccuracy was obtained. Comparison of the results with the numerical and experimental results of other researchers confirms higher accuracy of flow prediction in the present work.

---

In the present article, velocity and deformation of an air bubble have been considered in quiescent liquid at different consecutive slopes from 5 to 90 degrees in respect to horizontal condition. To establish these purposes, air-water two-phase flow has been simulated numerically by using volume of fluid method. The two-phase flow interface has been traced by using Piecewise Linear Interface Calculation (PLIC) method. Surface tension force was estimated by Continuum Surface Force (CSF) model. The simulation results show that maximum bubble velocity occurred at 45 degrees which is in agreement with the previous researchers result. Simulation of bubble movement was also continued to two consecutive slopes at different angles. At slope deviation location, a vortex was generated due to liquid movement governed by gravity forces. This vortex changes the bubble velocity as well as bubble shape. This vortex also reduces the bubble velocity and changes the bubble nose shape from sharp to flatten at deviation from low to high slope values. However, at deviations from high to low slope values, the bubble nose becomes more sharpened in addition to bubble velocity increase. The maximum average velocity of bubble movement at two consecutive slopes was obtained during the condition that the first and second slopes were set to 60 and 30 degrees, respectively.

---

In present study, volume of fluid method in OpenFOAM open source CFD package will be extended to consider phase change phenomena due to condensation process. Both phases (liquid – vapor) are incompressible and immiscible. Vapor phase is assumed in saturated temperature. Interface between two phases are tracked with color function volume of fluid (CF-VOF) method. ُSurface Tension is taken accounted by Continuous Surface Force (CSF) model and mass transfer occurs along interface is considered by Lee mass transfer model. Pressure-Velocity coupling will be solved with PISO algorithm in the collocated grid. This solver is validated with Stefan problem. In one dimensional Stefan problem, the desistance of interface motion from cold wall is compared by the analytical solution. Then condensate laminar liquid film flow over vertical plate is simulated in the presence of gravity. Numerical result shows calculated film thickness from numerical simulation is thinner than analytical solution. Also, it shows Nusselt number is a function of vapor specific heat which neglected in existing correlations, therefore analytical solution and experimental correlation should modified to consider this effect on the Nusselt Number.

---

One of obstacles in simulation of two phase flow is parasite currents. These currents cause unphysical distortion at interface which impairs interface capturing and numerical results. In present study, two methods (using Filter and s-CLSVOF) are implemented in OpenFOAM two phase flow solver called interFoam to reduce parasite current. 3 filters are added to color function volume of fluid (CF-VOF) method. These filters reduce parasite current in different ways, one smoothes color function, one smoothes curvature and the other one compresses the interface. The original and the modified solvers are tested with a quiescent bubble bench mark to investigate the effect of each filter on parasite currents. Then optimum arrangement of filters is compared with s-CLSVOF method and interFoam. Present study shows parasite current magnitude can be reduced at least up to 50% in the modified solvers. Also, the comparison of pressure jump from numerical results and analytical result with Young-Laplace equation shows modified solvers can predict pressure jump better than original solver. The pressure jump error is reduced up to 400% in the modified solvers. Also present study shows filters have better performance than s-CLVOF method and it can be considered as a suitable substitution of coupled methods.

---

In present study, impact of single bubble on an inclined wall and its movement are investigated by applying volume of fluid method (VOF) in OpenFOAM open source cfd package using a solver called interFoam. Both phases are incompressible and surface tension between two phases is estimated by CSF method. The effect of some parameters such as contact angle, wall slope and Bond and Morton dimensionless numbers on bubble shapes and velocity are studied. The numerical results show bubble velocity along wall increases with the increase of wall slope angle. The maximum bubble velocity happens at 50 degree. Three bubble regimes are recognized and introduced in this study named as: sliding, bouncing, and zigzagging based on wall slope. The bubble regime changes from sliding to bouncing when wall slope changes from 30 to 40 degrees. In constant Morton number, increment of Bond number increases both velocity and amplitude of fluctuations. In addition, an increment of Morton number in constant Bond number, decreases velocity and amplitude of fluctuations. Moreover, by increment of Morton number, the bubble motion will change from an accelerating motion to a constant velocity condition.