Steady free surface flows are of interest in the fields of marine and hydraulic engineering. Fitting methods are generally used to represent the free surface position with a deforming grid. Existing fitting methods tend to use time-stepping schemes, which is inefficient for steady flows. There also exists a steady iterative method, but that one needs to be implemented with a dedicated solver.Therefore a new method is proposed to efficiently simulate two-dimensional (2D) steady free surface flows, suitable for use in conjunction with black-box flow solvers. The free surface position is calculated with a quasi-Newton method, where the approximate Jacobian is constructed in a novel way by combining data from past iterations with an analytical model based on a perturbation analysis of a potential flow. The method is tested on two 2D cases: the flow over a bottom topography and the flow over a hydrofoil. For all simulations the new method converges exponentially and in few iterations. Furthermore, convergence is independent of the free surface mesh size for all tests.
K E Y W O R D Sfitting method, free surface flow, perturbation analysis, quasi-Newton
INTRODUCTIONThis paper considers the numerical simulation of steady free surface flow of incompressible, immiscible fluids with a large density difference, typically water and air. This type of flow is often encountered in the fields of marine and hydraulic engineering, for example to calculate ship hull resistance, analyze ship-wall interactions in narrow straight canals and study flow behavior in river confluences. These flows are governed by the incompressible Navier-Stokes equations, which are typically solved with computational fluid dynamics (CFD). The free surface causes additional difficulties for the flow computations: its position is unknown a priori, and needs to be determined as a component of the solution during the computation. Two classes of methods exist to represent the free surface, namely, surface capturing and surface fitting methods. * In surface capturing approaches, the computational mesh is not aligned with the interface, and the interface can intersect the mesh in an in principle arbitrary manner. Capturing methods are versatile in that they can handle complex phenomena such as wave breaking. Examples are the marker-and-cell method, 2 the volume-of-fluid method, 3 and the level-set method. 4 *Classification and terminology of these methods tends to be inconsistent in the literature; terminology is adopted from Wackers et al. 1Int J Numer Meth Fluids. 2020;92:785-801.wileyonlinelibrary.com/journal/fld