A fully hydrodynamic model for three‐dimensional,free‐surface flows

The hydrostatic pressure assumption has been widely used in studying water movements in rivers, lakes, estuaries, and oceans. While this assumption is valid in many cases and has been successfully used in numerous studies, there are many cases where this assumption is questionable. This paper presents a three‐dimensional, hydrodynamic model for free‐surface flows without using the hydrostatic pressure assumption. The model includes two predictor–corrector steps. In the first predictor–corrector step, the model uses hydrostatic pressure at the previous time step as an initial estimate of the total pressure field at the new time step. Based on the estimated pressure field, an intermediate velocity field is calculated, which is then corrected by adding the non‐hydrostatic component of the pressure to the estimated pressure field. A Poisson equation for non‐hydrostatic pressure is solved before the second intermediate velocity field is calculated. The final velocity field is found after the free surface at the new time step is computed by solving a free‐surface correction equation. The numerical method was validated with several analytical solutions and laboratory experiments. Model results agree reasonably well with analytical solutions and laboratory results. Model simulations suggest that the numerical method presented is suitable for fully hydrodynamic simulations of three‐dimensional, free‐surface flows. Copyright © 2003 John Wiley & Sons, Ltd.     

An implicit two‐dimensional non‐hydrostatic model for free‐surface flows

An implicit finite volume model in sigma coordinate system is developed to simulate two‐dimensional (2D) vertical free surface flows, deploying a non‐hydrostatic pressure distribution. The algorithm is based on a projection method which solves the complete 2D Navier–Stokes equations in two steps. First the pressure term in the momentum equations is excluded and the resultant advection–diffusion equations are solved. In the second step the continuity and the momentum equation with only the pressure terms are solved to give a block tri‐diagonal system of equation with pressure as the unknown. This system can be solved by a direct matrix solver without iteration. A new implicit treatment of non‐hydrostatic pressure, similar to the lower layers is applied to the top layer which makes the model free of any hydrostatic pressure assumption all through the water column. This treatment enables the model to evaluate both free surface elevation and wave celerity more accurately. A series of numerical tests including free‐surface flows with significant vertical accelerations and nonlinear behaviour in shoaling zone are performed. Comparison between numerical results, analytical solutions and experimental data demonstrates a satisfactory performance. Copyright © 2006 John Wiley & Sons, Ltd.     

An implicit three‐dimensional fully non‐hydrostatic model for free‐surface flows

An implicit method is developed for solving the complete three‐dimensional (3D) Navier–Stokes equations. The algorithm is based upon a staggered finite difference Crank‐Nicholson scheme on a Cartesian grid. A new top‐layer pressure treatment and a partial cell bottom treatment are introduced so that the 3D model is fully non‐hydrostatic and is free of any hydrostatic assumption. A domain decomposition method is used to segregate the resulting 3D matrix system into a series of two‐dimensional vertical plane problems, for each of which a block tri‐diagonal system can be directly solved for the unknown horizontal velocity. Numerical tests including linear standing waves, nonlinear sloshing motions, and progressive wave interactions with uneven bottoms are performed. It is found that the model is capable to simulate accurately a range of free‐surface flow problems using a very small number of vertical layers (e.g. two–four layers). The developed model is second‐order accuracy in time and space and is unconditionally stable; and it can be effectively used to model 3D surface wave motions. Copyright © 2004 John Wiley & Sons, Ltd.     

Space‐time NURBS‐enhanced finite elements for free‐surface flows in 2D

The accuracy of numerical simulations of free‐surface flows depends strongly on the computation of geometric quantities like normal vectors and curvatures. This geometrical information is additional to the actual degrees of freedom and usually requires a much finer discretization of the computational domain than the flow solution itself. Therefore, the utilization of a numerical method, which uses standard functions to discretize the unknown function in combination with an enhanced geometry representation is a natural step to improve the simulation efficiency. An example of such method is the NURBS‐enhanced finite element method (NEFEM), recently proposed by Sevilla et al. The current paper discusses the extension of the spatial NEFEM to space‐time methods and investigates the application of space‐time NURBS‐enhanced elements to free‐surface flows. Derived is also a kinematic rule for the NURBS motion in time, which is able to preserve mass conservation over time. Numerical examples show the ability of the space‐time NEFEM to account for both pressure discontinuities and surface tension effects and compute smooth free‐surface forms. For these examples, the advantages of the NEFEM compared with the classical FEM are shown. Copyright © 2015 John Wiley & Sons, Ltd.     

An adaptive numerical scheme for solving incompressible 2‐phase and free‐surface flows

In this paper, we present a numerical scheme for solving 2‐phase or free‐surface flows. Here, the interface/free surface is modeled using the level‐set formulation, and the underlying mesh is adapted at each iteration of the flow solver. This adaptation allows us to obtain a precise approximation for the interface/free‐surface location. In addition, it enables us to solve the time‐discretized fluid equation only in the fluid domain in the case of free‐surface problems. Fluids here are considered incompressible. Therefore, their motion is described by the incompressible Navier‐Stokes equation, which is temporally discretized using the method of characteristics and is solved at each time iteration by a first‐order Lagrange‐Galerkin method. The level‐set function representing the interface/free surface satisfies an advection equation that is also solved using the method of characteristics. The algorithm is completed by some intermediate steps like the construction of a convenient initial level‐set function (redistancing) as well as the construction of a convenient flow for the level‐set advection equation. Numerical results are presented for both bifluid and free‐surface problems.     

Three‐dimensional modeling of non‐hydrostatic free‐surface flows on unstructured grids

A three‐dimensional numerical model is presented for the simulation of unsteady non‐hydrostatic shallow water flows on unstructured grids using the finite volume method. The free surface variations are modeled by a characteristics‐based scheme, which simulates sub‐critical and super‐critical flows. Three‐dimensional velocity components are considered in a collocated arrangement with a σ‐coordinate system. A special treatment of the pressure term is developed to avoid the water surface oscillations. Convective and diffusive terms are approximated explicitly, and an implicit discretization is used for the pressure term to ensure exact mass conservation. The unstructured grid in the horizontal direction and the σ coordinate in the vertical direction facilitate the use of the model in complicated geometries. Solution of the non‐hydrostatic equations enables the model to simulate short‐period waves and vertically circulating flows. Copyright © 2015 John Wiley & Sons, Ltd.     

CFD simulation of two‐ and three‐dimensional free‐surface flow

The paper describes the implementation of moving‐mesh and free‐surface capabilities within a 3‐d finite‐volume Reynolds‐averaged‐Navier–Stokes solver, using surface‐conforming multi‐block structured meshes. The free‐surface kinematic condition can be applied in two ways: enforcing zero net mass flux or solving the kinematic equation by a finite‐difference method. The free surface is best defined by intermediate control points rather than the mesh vertices. Application of the dynamic boundary condition to the piezometric pressure at these points provides a hydrostatic restoring force which helps to eliminate any unnatural free‐surface undulations. The implementation of time‐marching methods on moving grids are described in some detail and it is shown that a second‐order scheme must be applied in both scalar‐transport and free‐surface equations if flows driven by free‐surface height variations are to be computed without significant wave attenuation using a modest number of time steps. Computations of five flows of theoretical and practical interest—forced motion in a pump, linear waves in a tank, quasi‐1d flow over a ramp, solitary wave interaction with a submerged obstacle and 3‐d flow about a surface‐penetrating cylinder—are described to illustrate the capabilities of our code and methods. Copyright © 2003 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.     

A marker‐and‐cell approach to free surface 2‐D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front‐tracking method. The velocity field is computed using a finite‐difference discretization of a modification of the Navier–Stokes equations. These equations together with the continuity equation are solved for the two‐dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright © 2012 John Wiley & Sons, Ltd.     

A two‐dimensional vertical non‐hydrostatic σ model with an implicit method for free‐surface flows

An implicit finite difference model in the σ co‐ordinate system is developed for non‐hydrostatic, two‐dimensional vertical plane free‐surface flows. To accurately simulate interaction of free‐surface flows with uneven bottoms, the unsteady Navier–Stokes equations and the free‐surface boundary condition are solved simultaneously in a regular transformed σ domain using a fully implicit method in two steps. First, the vertical velocity and pressure are expressed as functions of horizontal velocity. Second, substituting these relationship into the horizontal momentum equation provides a block tri‐diagonal matrix system with the unknown of horizontal velocity, which can be solved by a direct matrix solver without iteration. A new treatment of non‐hydrostatic pressure condition at the top‐layer cell is developed and found to be important for resolving the phase of wave propagation. Additional terms introduced by the σ co‐ordinate transformation are discretized appropriately in order to obtain accurate and stable numerical results. The developed model has been validated by several tests involving free‐surface flows with strong vertical accelerations and non‐linear waves interacting with uneven bottoms. Comparisons among numerical results, analytical solutions and experimental data show the capability of the model to simulate free‐surface flow problems. Copyright © 2004 John Wiley & Sons, Ltd.     

The solution of two‐dimensional free‐surface problems using automatic mesh generation

A new method is described for the iterative solution of two‐dimensional free‐surface problems, with arbitrary initial geometries, in which the interior of the domain is represented by an unstructured, triangular Eulerian mesh and the free surface is represented directly by the piecewise‐quadratic edges of the isoparametric quadratic‐velocity, linear‐pressure Taylor–Hood elements. At each time step, the motion of the free surface is computed explicitly using the current velocity field and, once the new free‐surface location has been found, the interior nodes of the mesh are repositioned using a continuous deformation model that preserves the original connectivity. In the event that the interior of the domain must be completely remeshed, a standard Delaunay triangulation algorithm is used, which leaves the initial boundary discretisation unchanged. The algorithm is validated via the benchmark viscous flow problem of the coalescence of two infinite cylinders of equal radius, in which the motion is due entirely to the action of capillary forces on the free surface. This problem has been selected for a variety of reasons: the initial and final (steady state) geometries differ considerably; in the passage from the former to the latter, large free‐surface curvatures—requiring accurate modelling—are encountered; an analytical solution is known for the location of the free surface; there exists a large body of literature on alternative numerical simulations. A novel feature of the present work is its geometric generality and robustness; it does not require a priori knowledge of either the evolving domain geometry or the solution contained therein. Copyright © 1999 John Wiley & Sons, Ltd.     

A new enrichment space for the treatment of discontinuous pressures in multi‐fluid flows

In this work, a new enrichment space to accommodate jumps in the pressure field at immersed interfaces in finite element formulations, is proposed. The new enrichment adds two degrees of freedom per element that can be eliminated by means of static condensation. The new space is tested and compared with the classical P₁ space and to the space proposed by Ausas et al (Comp. Meth. Appl. Mech. Eng., Vol. 199, 1019–1031, 2010) in several problems involving jumps in the viscosity and/or the presence of singular forces at interfaces not conforming with the element edges. The combination of this enrichment space with another enrichment that accommodates discontinuities in the pressure gradient has also been explored, exhibiting excellent results in problems involving jumps in the density or the volume forces. Copyright © 2011 John Wiley & Sons, Ltd.     

Conservation of circulation in SPH for 2D free‐surface flows

In this paper, we study how accurately the Smoothed Particle Hydrodynamics (SPH) scheme accounts for the conservation and the generation of vorticity and circulation, in a low viscosity, weakly compressible, barotropic fluid in the context of free‐surface flows. We consider a number of simple examples to clarify the processes involved and the accuracy of the simulations. The first example is a differentially rotating fluid where the integration path for the circulation becomes progressively more complicated, whereas the structure of the velocity field remains simple. The second example is the collision of two rectangular regions of fluid. We show that SPH accurately predicts the time variation of the circulation as well as the total vorticity for selected domains advected by the fluid. Finally, a breaking wave is considered. For such a problem we show how the dynamics of the vorticity generated by the breaking process is captured by the SPH model. Copyright © 2012 John Wiley & Sons, Ltd.     

Computation of free‐surface flows with local mesh adaptation

This paper describes a modern free‐surface capturing strategy implemented in an unstructured finite‐volume viscous flow solver that can handle moving grids composed of arbitrary‐shaped control volumes. An adaptive mesh strategy is fully integrated in the code making it a single tool for dynamically maintaining a prescribed density of grid points around the steady or unsteady interface between air and water. The whole adaptive procedure is described in detail. The efficiency of the overall approach is examined on two‐ and three‐dimensional hydrodynamic applications. The adaptive strategy achieves interesting gains in terms of computational and human efforts compared to single‐mesh computations. Copyright © 2005 John Wiley & Sons, Ltd.     

Semi‐implicit subgrid modelling of three‐dimensional free‐surface flows

In this paper, a semi‐implicit numerical model for two‐ and three‐dimensional free‐surface flows will be formulated in such a fashion as to intrinsically account for subgrid bathymetric details. It will be shown that with the proposed subgrid approach the model accuracy can be substantially improved without increasing the corresponding computational effort. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical solution of steady free‐surface flows by the adjoint optimal shape design method

Numerical solution of flows that are partially bounded by a freely moving boundary is of great importance in practical applications such as ship hydrodynamics. Free‐boundary problems can be reformulated into optimal shape design problems, which can in principle be solved efficiently by the adjoint method. In this work we investigate the suitability of the adjoint shape optimization method for solving steady free‐surface flows. The asymptotic convergence behaviour of the method is determined for free‐surface flows in 2D and 3D. It is shown that the convergence behaviour depends sensitively on the occurrence of critical modes. The convergence behaviour is moreover shown to be mesh‐width independent, provided that proper preconditioning is applied. Numerical results are presented for 2D flow over an obstacle in a channel. The observed convergence behaviour is indeed mesh‐width independent and conform the derived asymptotic estimates. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulation of the vertical structure of discontinuous flows

A numerical method to solve the Reynolds‐averaged Navier–Stokes equations with the presence of discontinuities is outlined and discussed. The pressure is decomposed into the sum of a hydrostatic component and a hydrodynamic component. The numerical technique is based upon the classical staggered grids and semi‐implicit finite difference methods applied for quasi‐ and non‐hydrostatic flows. The advection terms in the momentum equations are approximated in order to conserve mass and momentum following the principles recently developed for the numerical simulation of shallow water flows with large gradients. Conservation of these properties is the most important aspect to represent near local discontinuities in the solution, following from sharp bottom gradients or hydraulic jumps. The model is applied to reproduce the flow over a step where a hydraulic jump forms downstream. The hydrostatic pressure assumption fails to represent this type of flow mainly because of the pressure deviation from the hydrostatic values downstream the step. Fairly accurate results are obtained from the numerical model compared with experimental data. Deviation from the data is found to be inherent to the standard k–ε model implemented. Copyright © 2001 John Wiley & Sons, Ltd.     

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.     

Further experiences with computing non‐hydrostatic free‐surface flows involving water waves

A semi‐implicit, staggered finite volume technique for non‐hydrostatic, free‐surface flow governed by the incompressible Euler equations is presented that has a proper balance between accuracy, robustness and computing time. The procedure is intended to be used for predicting wave propagation in coastal areas. The splitting of the pressure into hydrostatic and non‐hydrostatic components is utilized. To ease the task of discretization and to enhance the accuracy of the scheme, a vertical boundary‐fitted co‐ordinate system is employed, permitting more resolution near the bottom as well as near the free surface. The issue of the implementation of boundary conditions is addressed. As recently proposed by the present authors, the Keller‐box scheme for accurate approximation of frequency wave dispersion requiring a limited vertical resolution is incorporated. The both locally and globally mass conserved solution is achieved with the aid of a projection method in the discrete sense. An efficient preconditioned Krylov subspace technique to solve the discretized Poisson equation for pressure correction with an unsymmetric matrix is treated. Some numerical experiments to show the accuracy, robustness and efficiency of the proposed method are presented. Copyright © 2004 John Wiley & Sons, Ltd.