Development of a Free Surface Capability in a RANS Solver with Coupled Equations and Overset Grids
Licentiate thesis, 2013
This thesis is a result of investigations carried out during the implementation of a free surface capturing method in a RANS solver. The work consists of two distinct parts which complement each other, and the purpose is to improve the free surface interface sharpness and to reduce numerical diffusion of the gravity waves. The first part deals with numerical schemes and the second with local grid refinement and adaptivity near the interface. There is a vast amount of discretization schemes for convection equations. However, only a subset is found to be applicable to the water fraction transport equation. Here we investigate the performance of selected schemes applied to hydrodynamic problems. Moreover, an additional scheme is developed based on the experience gained during the work. This part of the work was carried out since there is very little information in the literature for problems beyond academic test cases such as bubble convection or two dimensional shape translation and rotation. The second part of the presented work concerns air-water interface refinements using fixed and adaptive overlapping grids. The adaptive refinement grid follows approximately the free surface interface and reduces the refinement thickness, thus reducing the required number of discretization cells. The grid alignment with the interface reduces further the diffusion of the water fraction in the transition region from air to water. The method is validated on 2D and 3D cases. The Duncan submerged hydrofoil test case as well as Wigley, Series 60 and KCS hulls are compared with experiments. A grid dependence study is presented for Wigley and KCS which shows that the code is robust and the deviations from the measurement data are within the expected accuracy of a CFD code for naval applications.