Alot of phenomena in our life include multiple phases turbulent flow. While simulating such problems needs powerful interface tracking method to be involved in the simulation scheme, it is still challenging to find the method that can catch the interface curvature and to maintain volume conservation. Meanwhile, CVOF is a promising method that shows good results with high precision, and it was established with adopting Weymouth and Yue (WY) Eulerian-Implicit method which allows volume conservation. Neverthe- less, at high contrast in the density ratio between the fluids in the domain, the method tends to fail in adverting volume fraction in a consistent manner with the momentum advection. Therefore, a serious need to couple the advection of both fields in order to retain the advection of momentum with the same speed of volume fraction advection is taking place. Meanwhile, MGLT is single phase optimally developed code for DNS and turbulence simulations. In order to upgrade it to become capable of tackling two phase problems at high Reynolds numbers, the CVOF has to be integrated fully while preserving the consistency between the momentum and mass advection. In this thesis, we achieved integration between CVOF and MGLET to the limits of performing com- plex two phase flow problems. The outcome suggested that coupling both software fully empowered MGLET for two phase flow problem and at the same time it did not affect its computational performance in comparison to its original state. Meanwhile, multiple investigations to different benchmarks such as open-channel case, droplet of water in gaseous phase, Taylor-Green vortex in two phase flow, and Kelvin-Helmholtz instability that revealed that the coupling was successful, and yet it showed the limits of our in- tegration approach. Finally, the benchmark tests were informative regarding the areas of the coupling that need to be improved and it allowed us to figure the limits of our integration methodology.     «

Alot of phenomena in our life include multiple phases turbulent flow. While simulating such problems needs powerful interface tracking method to be involved in the simulation scheme, it is still challenging to find the method that can catch the interface curvature and to maintain volume conservation. Meanwhile, CVOF is a promising method that shows good results with high precision, and it was established with adopting Weymouth and Yue (WY) Eulerian-Implicit method which allows volume conservatio...     »