The subject of this thesis is the computation of bubbly flows. Bubbly flows occur for example in chemical reactors, boiling, fuel injectors and coating. The bubbles and the surrounding fluid are modeled directly without phase averaging. This is the most fundamental approach of describing multi-phase flows. Every bubble is modeled in great detail. Since the model is so demanding, an efficient numerical approach has to be taken. This work aims to develop an efficient, robust method for the direct numerical simulation of multi-phase flows. There exist two incompressible fluids (e.g. air and water) that are separated by an interface. The interface is a moving, internal boundary, where density and viscosity change discontinuously and surface tension forces act. When the governing equations are discretized, the question arises how to deal with this interface. Various methods have been put forward to treat moving boundary problems. The two methods that are of most interest are the so-called Volume-of-Fluid method and the Level-Set method. With both methods a coloring function is used to identify the individual phase. The interface is implicitly defined, which allows arbitrarily shaped interface. The advection of the interface is the key part of this thesis, for which the Level-Set approach is chosen. The major disadvantage of the Level-Set method is that the mass of each individual phase is not conserved when the interface is advected. Additional effort is necessary and in this research the Volume-of-Fluid function is used to conserve mass. The method is applied to a falling drop and a rising bubble in two and three dimensions, respectively. Merging of rising bubbles is studied for two aligned and two misaligned bubbles. Comparison with other numerical work and experimental data is made.Electrical Engineering, Mathematics and Computer Scienc