Numerical modelling of 3D stratified free surface flows: a case study of sediment dumping

A three‐dimensional numerical model has been developed to simulate stratified flows with free surfaces. The model is based on the Reynolds‐averaged Navier–Stokes (RANS) equations with variable fluid density. The equations are solved in a transformed σ‐coordinate system with the use of operator‐splitting method (Int. J. Numer. Meth. Fluids 2002; 38 :1045–1068). The numerical model is validated against the one‐dimensional diffusion problem and the two‐dimensional density‐gradient flow. Excellent agreements are obtained between numerical results and analytical solutions. The model is then used to study transport phenomena of dumped sediments into a water body, which has been modelled as a strongly stratified flow. For the two‐dimensional problem, the numerical results compare well with experimental data in terms of mean particle falling velocity and spreading rate of the sediment cloud for both coarse and medium‐size sediments. The model is also employed to study the dumping of sediments in a three‐dimensional environment with the presence of free surface. It is found that during the descending process an annulus‐like cloud is formed for fine sediments whereas a plate‐like cloud for medium‐size sediments. The model is proven to be a good tool to simulate strongly stratified free surface flows. Copyright © 2005 John Wiley & Sons, Ltd.     

3D free‐surface flow computation using a RANSE/Fourier–Kochin coupling

A coupling method for numerical calculations of steady free‐surface flows around a body is presented. The fluid domain in the neighbourhood of the hull is divided into two overlapping zones. Viscous effects are taken in account near the hull using Reynolds‐averaged Navier–Stokes equations (RANSE), whereas potential flow provides the flow away from the hull. In the internal domain, RANSE are solved by a fully coupled velocity, pressure and free‐surface elevation method. In the external domain, potential‐flow theory with linearized free‐surface condition is used to provide boundary conditions to the RANSE solver. The Fourier–Kochin method based on the Fourier–Kochin formulation, which defines the velocity field in a potential‐flow region in terms of the velocity distribution at a boundary surface, is used for that purpose. Moreover, the free‐surface Green function satisfying this linearized free‐surface condition is used. Calculations have been successfully performed for steady ship‐waves past a serie 60 and then have demonstrated abilities of the present coupling algorithm. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulation of turbulent free surface flow with two‐equation k–ε eddy‐viscosity models

This paper presents a finite difference technique for solving incompressible turbulent free surface fluid flow problems. The closure of the time‐averaged Navier–Stokes equations is achieved by using the two‐equation eddy‐viscosity model: the high‐Reynolds k–ε (standard) model, with a time scale proposed by Durbin; and a low‐Reynolds number form of the standard k–ε model, similar to that proposed by Yang and Shih. In order to achieve an accurate discretization of the non‐linear terms, a second/third‐order upwinding technique is adopted. The computational method is validated by applying it to the flat plate boundary layer problem and to impinging jet flows. The method is then applied to a turbulent planar jet flow beneath and parallel to a free surface. Computations show that the high‐Reynolds k–ε model yields favourable predictions both of the zero‐pressure‐gradient turbulent boundary layer on a flat plate and jet impingement flows. However, the results using the low‐Reynolds number form of the k–ε model are somewhat unsatisfactory. Copyright © 2004 John Wiley & Sons, Ltd.     

A fully coupled Navier–Stokes solver for calculation of turbulent incompressible free surface flow past a ship hull

This paper deals with the calculation of free surface flow of viscous incompressible fluid around the hull of a boat moving with rectilinear motion. An original method used to avoid a large part of the theoretical problems connected with free surface boundary conditions in three‐dimensional Navier–Stokes–Reynolds equations is proposed here. The linearised system of convective equations for velocities, pressure and free surface elevation unknowns is discretised by finite differences and two methods to solve the fully coupled resulting matrix are presented here. The non‐linear convergence of fully coupled algorithm is compared with the velocity–pressure weakly coupled algorithm SIMPLER. Turbulence is taken into account through Reynolds decomposition and k–ε or k–ω model to close the equations. These two models are implemented without wall function and numerical calculations are performed up to the viscous sub‐layer. Numerical results and comparisons with experiments are presented on the Series 60 CB=0.60 ship model for a Reynolds number Rn=4.5×10⁶ and a Froude number Fn=0.316. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical study of breaking waves by a two‐phase flow model

A two‐phase flow model, which solves the flow in the air and water simultaneously, is presented for modelling breaking waves in deep and shallow water, including wave pre‐breaking, overturning and post‐breaking processes. The model is based on the Reynolds‐averaged Navier–Stokes equations with the k ?ε turbulence model. The governing equations are solved by the finite volume method in a Cartesian staggered grid and the partial cell treatment is implemented to deal with complex geometries. The SIMPLE algorithm is utilised for the pressure‐velocity coupling and the air‐water interface is modelled by the interface capturing method via a high resolution volume of fluid scheme. The numerical model is validated by simulating overturning waves on a sloping beach and over a reef, and deep‐water breaking waves in a periodic domain, in which good agreement between numerical results and available experimental measurements for the water surface profiles during wave overturning is obtained. The overturning jet, air entrainment and splash‐up during wave breaking have been captured by the two‐phase flow model, which demonstrates the capability of the model to simulate free surface flow and wave breaking problems.Copyright © 2011 John Wiley & Sons, Ltd.     

Interaction of an obliquely rising vortex ring with a free surface in a viscous fluid

A vortex ring, which approaches a free surface under various initial conditions and different values of flow parameters, is investigated by means of numerical solutions of the Navier-Stokes equations. The study focuses on the connection process of the vortex ring with the free surface at low Reynolds number. Discrepancies between a numerical solution and experimental observations found recently with regard to the shape of the reconnected vorticity tubes have been resolved. The crucial parameter, which determines the essential differences at the low Reynolds number selected, is the angle of inclination at which the vortex ring rises toward the free surface. At an angle of inclination of 20° circular vorticity tubes at the free surface evolve, while at an angle of 45° the reconnection takes place in the form of an almost circular sheet. At higher Reynolds number this form of reconnection may change to the other form through instability. The possibility of vortex-ring reflection at the free surface is discussed.

---

Sommario Usando soluzioni numeriche delle equazioni di Navier-Stokes, viene analizzato il comportamento di un vortice and anello che si avvicina ad una superficie libera, per diverse condizioni iniziali e valori dei parametri del flusso. Lo studio è incentrato sul processo di connessione del vortice adnanello con la superficie libera per bassi numeri di Reynolds. con riguardo alla forma dei tubi di vorticità riconnessi con la superficie, tutte le differenze tra osservazioni sperimentali e simulazioni numeriche sono state risolte. Per bassi numeri di Reynolds, il parametro cruciale dell' interazione è l'angolo di inclinazione tra vortice e superficie. Per un angolo di 20° si formano tubi circolari di vorticità alla superficie, mentre a 45° la riconnessione ha luogo in forma di strato sottile di vorticità. Per altri numeri di Reynolds la forma della riconnessione può cambiare in altre forme a causa dell'instabilità. Infine viene discussa la possibilità di riflessione del vortice ad anello alla superficie.

     

The motion of viscous liquid column with finite length in a vertical straight capillary tube

This paper deals with the motion of viscous liquid column with finite length and two free surfaces in a vertical straight capillary tube. It is assumed that fluid is Newtonian. Linearizing the boundary conditions, analytic expressions in the form of infinite series have been obtained for velocity, piessure and free surface at low Reynolds number. The numerical calculation is carried out for a set of cylinder’s length of water and blood. It has been revealed that there are considerable circulating currents at the upper and lower meniscuses. Its maximum velocity is about 57% of the average velocity of the mainstream. Iner-tial effect is also studied in this paper. Using the time-dependent method in finite difference techniques, numerical solution of the corresponding nonlinear equation at Re<24.5 is computed. Comparing it with analytic exact solution at low Reynolds number shows that inertial effect is negligible provided Re<24.5.     

Flow above the free end of a surface-mounted finite-height circular cylinder: A review

The wake of a surface-mounted finite-height circular cylinder and the associated vortex patterns are strongly dependent on the cylinder aspect ratio and the thickness of the boundary layer on the ground plane relative to the dimensions of the cylinder. Above a critical aspect ratio, the mean wake is characterized by streamwise tip vortex structures and Kármán vortex shedding from the sides of the cylinder. Below a critical aspect ratio, a unique mean wake structure is observed. Recent experimental studies in the literature that used phase-averaged techniques, as well as recent numerical simulations, have led to an improved physical understanding of the near-wake vortex flow patterns. However, the flow above the free end of the finite circular cylinder, and its relationship to the near wake, has not been systematically studied. The effects of aspect ratio and boundary layer thickness on the free-end flow field are also not completely understood, nor has the influence of Reynolds number on the free-end flow field been fully explored. Common features associated with the free end include separation from the leading edge, a mean recirculation zone containing a prominent cross-stream arch (or mushroom) vortex, and reattachment onto the free-surface. Other flow features that remain to be clarified include a separation bubble near the leading edge, one or two cross-stream vortices within this separation bubble, the origins of the streamwise tip or trailing vortices, and various critical points in the near-surface flow topology. This paper reviews the current understanding of the flow above the free end of a surface-mounted finite-height circular cylinder, with a focus on models of the flow field, surface oil flow visualization studies, pressure and heat flux distributions on the free-end surface, measurements of the local velocity field, and numerical simulations, found in the literature.     

A finite difference technique for solving the Newtonian jet swell problem

A finite difference technique that incorporates a numerical mapping has been successfully applied to analyse both planar and axisymmetric Newtonian jets. A pressure gradient equation and a free-surface slope equation have been derived for free-surface iteration. The computation of pressure inside the jet surface using the pressure gradient equation is stable and accurate at high Reynolds numbers. The free-surface slope equation is needed for updating the free surface and is applicable for jets with strong surface tension effects. The present development can simulate the Newtonian jets for Reynolds numbers as high as 2000 and capillary number as low as 10^?5. Numerical predictions by the present technique are close to the results of previous finite element simulations.     

Flow past a circular cylinder over a free surface: Interaction between the near wake and the free surface deformation

Air-flow around a circular cylinder placed above a free surface and liquid flow under the free surface were investigated experimentally in a wind/wave tunnel. The cylinder spanned the tunnel test-section and was oriented normal to the freestream direction. The main objective of this study was to investigate the interaction of the cylinder wake with the free surface. The flow structure was analyzed for various gap widths, H, between the cylinder and the free surface using a digital particle image velocimetry (PIV) system with a spatial resolution of 2048×2048 pixels. The Reynolds number based on the cylinder diameter was 3.3×10³. For each experimental condition, 400 instantaneous velocity fields were measured and ensemble-averaged to obtain spatial distributions of the mean velocity and turbulence statistics. The results showed that the cylinder near-wake inclined upward due to the influence of the free surface elevation. Vortices were shed, even at a small gap ratio of H/D=0.25, where D is the cylinder diameter. Strong jet-like flow appeared in the gap beneath the cylinder. At a gap ratio of H/D=0.50, the jet flow exhibited a quasi-periodic vibration with a period of 2–3 s. The free surface deformation was caused by the pressure difference in the air-flow immediately above it. As the gap ratio increased, the inclination angle of the wake and the height of the free surface elevation decreased gradually. The liquid flow under the free surface followed a convective flow motion, and the range of the convection depended on the gap width between the cylinder and the free surface.     

Finite element study of free surface turbulent flow in rotating channels

Under certain conditions of liquid flow through rotating channels, the Coriolis force can induce a free surface to be formed. This problem is of practical importance in a Coriolis wear tester, which is used for determining the sliding wear coefficient of wear materials in slurry handling equipment. A deforming Galerkin finite element method is presented for predicting two‐dimensional turbulent free surface mean flow in rotating channels. Reynolds‐averaged Navier–Stokes (RANS) equations are cast into weak(algebraic) form using primitive variables (velocity and pressure). Eddy viscosity is determined via a mixing length model. Velocity is interpolated biquadratically, while pressure is interpolated bilinearly. The kinematic condition is used to form the Galerkin residual for the free surface. The free surface is represented by Hermite polynomials of zeroeth order for continuity of position and slope. Combined Newton's iteration is used to simultaneously solve for the free surface and the field variables. Results of velocity and pressure fields, as well as the free surface are shown to converge with mesh‐size refinement. There is excellent respect for mass conservation. Results are presented for various values of Rossby number (Ro) and height‐based Reynolds number (Re_H). Parameter continuation in Ro and Re_H space is used to compute solutions at higher values of flow rate and angular velocity. Copyright © 2002 John Wiley & Sons, Ltd.     

Development of a PLIC-VOF method for the dynamic simulation of entry region flow in a laminar falling film

The present work describes a numerical procedure to simulate the development of hydrodynamic entry region in a gravity-driven laminar liquid film flow over an inclined plane. It provides a better insight into the physics of developing film in entry region. A novel numerical approach is proposed which has the potential to provide solutions for the complex physics of liquid film spreading on solid walls. The method employs an incompressible flow algorithm to solve the governing equations, a PLIC-VOF method to capture the free surface evolution and a continuum surface force (CSF) model to include the effect of surface tension. To account for the moving contact line on the solid substrate, a precursor film model based wall treatment is implemented. Liquid film flow has been simulated for the Reynolds number range of 5 ≤ Re ≤ 37.5, and the predicted results are found to agree well with the available analytical and experimental data.     

The solution of viscous incompressible jet flows using non‐staggered boundary fitted co‐ordinate methods

A new approach for the solution of the steady incompressible Navier–Stokes equations in a domain bounded in part by a free surface is presented. The procedure is based on the finite difference technique, with the non‐staggered grid fractional step method used to solve the flow equations written in terms of primitive variables. The physical domain is transformed to a rectangle by means of a numerical mapping technique. In order to design an effective free solution scheme, we distinguish between flows dominated by surface tension and those dominated by inertia and viscosity. When the surface tension effect is insignificant we used the kinematic condition to update the surface; whereas, in the opposite case, we used the normal stress condition to obtain the free surface boundary. Results obtained with the improved boundary conditions for a plane Newtonian jet are found to compare well with the available two‐dimensional numerical solutions for Reynolds numbers, up to Re=100, and Capillary numbers in the range of 0≤Ca<1000. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical simulations of laminar flow over a 3D backward-facing step

A numerical investigation of laminar flow over a three-dimensional backward-facing step is presented with comparisons with detailed experimental data, available in the literature, serving to validate the numerical results. The continuity constraint method, implemented via a finite element weak statement, was employed to solve the unsteady three-dimensional Navier–Stokes equations for incompressible laminar isothermal flow. Two-dimensional numerical simulations of this step geometry underestimate the experimentally determined extent of the primary separation region for Reynolds numbers Re greater than 400. It has been postulated that this disagreement between physical and computational experiments is due to the onset of three-dimensional flow near Re ≈ 400. This paper presents a full three-dimensional simulation of the step geometry for 100⩽ Re⩽ 800 and correctly predicts the primary reattachment lengths, thus confirming the influence of three-dimensionality. Previous numerical studies have discussed possible instability modes which could induce a sudden onset of three-dimensional flow at certain critical Reynolds numbers. The current study explores the influence of the sidewall on the development of three-dimensional flow for Re greater than 400. Of particular interest is the characterization of three-dimensional vortices in the primary separation region immediately downstream of the step. The complex interaction of a wall jet, located at the step plane near the sidewall, with the mainstream flow reveals a mechanism for the increasing penetration (with increasing Reynolds number) of three-dimensional flow structures into a region of essentially two-dimensional flow near the midplane of the channel. The character and extent of the sidewall-induced flow are investigated for 100⩽Re⩽ 800. © 1997 John Wiley & Sons, Ltd.     

Computation of turbulent free‐surface flows around modern ships

This paper presents the calculated results for three classes of typical modern ships in modelling of ship‐generated waves. Simulations of turbulent free‐surface flows around ships are performed in a numerical water tank, based on the FINFLO‐RANS SHIP solver developed at Helsinki University of Technology. The Reynolds‐averaged Navier–Stokes (RANS) equations with the artificial compressibility and the non‐linear free‐surface boundary conditions are discretized by means of a cell‐centred finite‐volume scheme. The convergence performance is improved with the multigrid method. A free surface is tracked using a moving mesh technology, in which the non‐linear free‐surface boundary conditions are given on the actual location of the free surface. Test cases recommended are a container ship, a US Navy combatant and a tanker. The calculated results are compared with the experimental data available in the literature in terms of the wave profiles, wave pattern, and turbulent flow fields for two turbulence models, Chien's low Reynolds number k–εmodel and Baldwin–Lomax's model. Furthermore, the convergence performance, the grid refinement study and the effect of turbulence models on the waves have been investigated. Additionally, comparison of two types of the dynamic free‐surface boundary conditions is made. Copyright © 2003 John Wiley& Sons, Ltd.     

Reynolds stress modelling of rectangular open‐channel flow

A Reynolds stress model for the numerical simulation of uniform 3D turbulent open‐channel flows is described. The finite volume method is used for the numerical solution of the flow equations and transport equations of the Reynolds stress components. The overall solution strategy is the SIMPLER algorithm, and the power‐law scheme is used to discretize the convection and diffusion terms in the governing equations. The developed model is applied to a flow at a Reynolds number of 77000 in a rectangular channel with a width to depth ratio of 2. The simulated mean flow and turbulence structures are compared with measured and computed data from the literature. The computed flow vectors in the plane normal to the streamwise direction show a small vortex, called inner secondary currents, located at the juncture of the sidewall and the free surface as well as the free surface and bottom vortices. This small vortex causes a significant increase in the wall shear stress in the vicinity of the free surface. A budget analysis of the streamwise vorticity is carried out. It is found that both production terms by anisotropy of Reynolds normal stress and by Reynolds shear stress contribute to the generation of secondary currents. Copyright © 2006 John Wiley & Sons, Ltd.     

Unsteady convective surface heat flux measurements on cylinder for CFD code validation

The surface convective heat transfer rates to a cylinder have been measured using platinum thin film gauges in a shock tube and the results have been used to validate two numerical codes. The investigations have been carried out at different incident shock Mach numbers and flow Reynolds numbers. The measured and simulated results give an insight into the transient flow fields around the model in the shock tube and are a valuable means of complementing the numerical and experimental techniques used in this study.Received: 27 August 2003, Accepted: 11 November 2003, Published online: 17 February 2004Correspondence to: T. Saito     

An implicit numerical algorithm for solving non‐hydrostatic free‐surface flow problems

Details are given of the development of a two‐dimensional vertical numerical model for simulating unsteady free‐surface flows, using a non‐hydrostatic pressure distribution. In this model, the Reynolds equations and the kinematic free‐surface boundary condition are solved simultaneously, so that the water surface elevation can be integrated into the solution and solved for, together with the velocity and pressure fields. An efficient numerical algorithm has been developed, deploying implicit parameters similar to those used in the Crank–Nicholson method, and generating a block tri‐diagonal algebraic system of equations. The model has been applied to simulate a range of unsteady flow problems involving relatively strong vertical accelerations. The results show that the numerical algorithm described is able to produce accurate predictions and is also easy to apply. Copyright © 2001 John Wiley & Sons, Ltd.