A simple numerical technique for turbulent flows with free surfaces

A simple technique is presented for the numerical solution of two‐dimensional time‐dependent flows, either laminar or turbulent, involving multiple free surfaces of arbitrary configuration. The governing equations are the Reynolds equations for incompressible fluids with Boussinesq closure, the k ‐ and ϵ‐equations and an additional equation describing the fluid configuration. This technique can potentially describe the propagation, deformation and overturning of pre‐breaking waves and the mean flow, surface configuration and turbulence field after breaking. The properties of the method are illustrated by several calculational examples. The main parts of the algorithm are optimized for vector processing in a form suitable for installation in supercomputer facilities.