Comparison of numerical schemes for computational simulation of liquid wavy film flow on vertical wall

Film flows are classified into no-wavy laminar, wavy laminar and turbulence along the Reynolds number or the flow stability. Since the wavy motions of the film flows are so intricate and nonlinear, the studies have largely been dependent upon the experimental way. The numerical approaches have been limited on the no-wavy flow regime. In this study, the film's full wavy motion is numerically investigated. The present numerical approach is based on the VOF interface capturing methods. The implicit and explicit schemes for the interface time marching are compared with each other. The HRIC (High Resolution Interface Capturing) and the geometry reconstruction method are tested as the treatment of convective terms at the liquid-air interface. The dependency of the grid density and the time step size are investigated in the aspect of the probability density function of the film thickness. All results are compared with Takamasa's experimental results.