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.     

Development of a generalized multi‐layer model for 3‐D simulation of free surface flows

A mathematical model was developed for three‐dimensional (3‐D) simulation of free surface flows. In this model, the flow depth is divided into a number of layers and shallow water equations are integrated in each layer to derive the hydrodynamic equations. To give a general form to this model, each layer is assumed to be non‐horizontal with varying thickness in the flow domain. A non‐orthogonal curvilinear coordinate system is employed in the model, to allow for flexibility in dealing with the irregular geometry of natural watercourses. Due to the similarity in governing equations, two‐dimensional (2‐D) depth averaged programs can be developed into a multi‐layer model. The development for a depth averaged program and its numerical scheme is described in this paper. Experimental data and semi‐analytical solutions are used to evaluate the performance of the model. Three different cases of open channel flow are tested: 1‐flow in a straight open channel, 2‐the flow development region in a channel, and 3‐flow in a meandering channel. It is shown that the model has the capability to predict velocity distribution and secondary flows in complex 3‐D flow conditions. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Three‐dimensional numerical modelling of free surface flows with non‐hydrostatic pressure

A three‐dimensional numerical model is developed for incompressible free surface flows. The model is based on the unsteady Reynolds‐averaged Navier–Stokes equations with a non‐hydrostatic pressure distribution being incorporated in the model. The governing equations are solved in the conventional sigma co‐ordinate system, with a semi‐implicit time discretization. A fractional step method is used to enable the pressure to be decomposed into its hydrostatic and hydrodynamic components. At every time step one five‐diagonal system of equations is solved to compute the water elevations and then the hydrodynamic pressure is determined from a pressure Poisson equation. The model is applied to three examples to simulate unsteady free surface flows where non‐hydrostatic pressures have a considerable effect on the velocity field. Emphasis is focused on applying the model to wave problems. Two of the examples are about modelling small amplitude waves where the hydrostatic approximation and long wave theory are not valid. The other example is the wind‐induced circulation in a closed basin. The numerical solutions are compared with the available analytical solutions for small amplitude wave theory and very good agreement is obtained. Copyright © 2002 John Wiley & Sons, Ltd.     

A weakly compressible MPS method for modeling of open‐boundary free‐surface flow

A mesh‐free particle method, based on the moving particle semi‐implicit (MPS) interaction model, has been developed for the simulation of two‐dimensional open‐boundary free‐surface flows. The incompressibility model in the original MPS has been replaced with a weakly incompressible model. The effect of this replacement on the efficiency and accuracy of the model has been investigated. The new inflow–outflow boundary conditions along with the particle recycling strategy proposed in this study extend the application of the model to open‐boundary problems. The final model is able to simulate open‐boundary free surface flow in cases of large deformation and fragmentation of free surface. The models and proposed algorithms have been validated and applied to sample problems. The results confirm the model's efficiency and accuracy. Copyright © 2009 John Wiley & Sons, Ltd.     

Comparing analytical and numerical solution of nonlinear two and three‐dimensional hydrostatic flows

New test cases for frictionless, three‐dimensional hydrostatic flows have been derived from some known analytical solutions of the two‐dimensional shallow water equations. The flow domain is a paraboloid of revolution and the flow is determined by the initial conditions, the nonlinear advective terms, the Coriolis acceleration and by the hydrostatic pressure. Wetting and drying is also included. Some specific properties of the exact solutions are discussed under different hypothesis and relative importance of the forcing terms. These solutions are proposed for testing the stability, the accuracy and the efficiency of numerical models to be used for simulating environmental hydrostatic flows. The computed solutions obtained with a semi‐implicit finite difference—finite volume algorithm on unstructured grid are compared with the corresponding analytical solutions in both two and three space dimension. Excellent agreement are obtained for the velocity and for the resulting water surface elevation. Comparison of the computed inundation area also shows a good agreement with the analytical solution with degrading accuracy observed when the inundation area becomes relatively large and for long simulation time. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

A high‐resolution wetting and drying algorithm for free‐surface hydrodynamics

A new wetting and drying algorithm for numerical modeling free‐surface flows is proposed and analyzed. A well structured, mildly nonlinear system for the discrete water surface elevation is derived from the governing differential equations by requiring a correct mass balance in wet areas as well as in the region of transition from wet to dry and from dry to wet. Existence and uniqueness of the numerical solution, along with a convergence analysis of an iterative scheme for the mildly nonlinear system, is provided. The present algorithm is devised to use high‐resolution bathymetric data at subgrid level. The resulting model is quite efficient, does not require a threshold value for minimal water depth, does not produce un‐physical negative water depths and generates accurate results with relatively coarse mesh and large time step size. These features are illustrated on a severe test‐case with known analytical solution. Copyright © 2008 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.     

Depth‐averaged two‐dimensional curvilinear explicit finite analytic model for open‐channel flows

A depth‐averaged two‐dimensional model has been developed in the curvilinear co‐ordinate system for free‐surface flow problems. The non‐linear convective terms of the momentum equations are discretized based on the explicit–finite–analytic method with second‐order accuracy in space and first‐order accuracy in time. The other terms of the momentum equations, as well as the mass conservation equation, are discretized by the finite difference method. The discretized governing equations are solved in turn, and iteration in each time step is adopted to guarantee the numerical convergence. The new model has been applied to various flow situations, even for the cases with the presence of sub‐critical and supercritical flows simultaneously or sequentially. Comparisons between the numerical results and the experimental data show that the proposed model is robust with satisfactory accuracy. Copyright © 2000 John Wiley & Sons, Ltd.     

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.     

An improved technique for solving two‐dimensional shallow water problems

A numerical model for solving the 2D shallow water equations is proposed herewith. This model is based on a finite volume technique in a generalized co‐ordinate system, coupled with a semi‐implicit splitting algorithm in which a Helmholtz equation is used for the surface elevation. Several benchmark problems have proven the good accuracy of this method in complex geometries. Nevertheless, several numerical perturbations were noted in the surface elevation. After finding the origin, a new numerical technique is suggested, to avoid these perturbations. Several severe tests are proposed to validate this technique. Copyright © 1999 John Wiley & Sons, Ltd.     

A fully non‐linear model for three‐dimensional overturning waves over an arbitrary bottom

An accurate three‐dimensional numerical model, applicable to strongly non‐linear waves, is proposed. The model solves fully non‐linear potential flow equations with a free surface using a higher‐order three‐dimensional boundary element method (BEM) and a mixed Eulerian–Lagrangian time updating, based on second‐order explicit Taylor series expansions with adaptive time steps. The model is applicable to non‐linear wave transformations from deep to shallow water over complex bottom topography up to overturning and breaking. Arbitrary waves can be generated in the model, and reflective or absorbing boundary conditions specified on lateral boundaries. In the BEM, boundary geometry and field variables are represented by 16‐node cubic ‘sliding’ quadrilateral elements, providing local inter‐element continuity of the first and second derivatives. Accurate and efficient numerical integrations are developed for these elements. Discretized boundary conditions at intersections (corner/edges) between the free surface or the bottom and lateral boundaries are well‐posed in all cases of mixed boundary conditions. Higher‐order tangential derivatives, required for the time updating, are calculated in a local curvilinear co‐ordinate system, using 25‐node ‘sliding’ fourth‐order quadrilateral elements. Very high accuracy is achieved in the model for mass and energy conservation. No smoothing of the solution is required, but regridding to a higher resolution can be specified at any time over selected areas of the free surface. Applications are presented for the propagation of numerically exact solitary waves. Model properties of accuracy and convergence with a refined spatio‐temporal discretization are assessed by propagating such a wave over constant depth. The shoaling of solitary waves up to overturning is then calculated over a 1:15 plane slope, and results show good agreement with a two‐dimensional solution proposed earlier. Finally, three‐dimensional overturning waves are generated over a 1:15 sloping bottom having a ridge in the middle, thus focusing wave energy. The node regridding method is used to refine the discretization around the overturning wave. Convergence of the solution with grid size is also verified for this case. Copyright © 2001 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.     

A new fully non‐hydrostatic 3D free surface flow model for water wave motions

A new fully non‐hydrostatic model is presented by simulating three‐dimensional free surface flow on a vertical boundary‐fitted coordinate system. A projection method, known as pressure correction technique, is employed to solve the incompressible Euler equations. A new grid arrangement is proposed under a horizontal Cartesian grid framework and vertical boundary‐fitted coordinate system. The resulting model is relatively simple. Moreover, the discretized Poisson equation for pressure correction is symmetric and positive definite, and thus it can be solved effectively by the preconditioned conjugate gradient method. Several test cases of surface wave motion are used to demonstrate the capabilities and numerical stability of the model. Comparisons between numerical results and analytical or experimental data are presented. It is shown that the proposed model could accurately and effectively resolve the motion of short waves with only two layers, where wave shoaling, nonlinearity, dispersion, refraction, and diffraction phenomena occur. Copyright © 2010 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.     

An efficient finite difference scheme for free‐surface flows in narrow rivers and estuaries

This paper presents a free‐surface correction (FSC) method for solving laterally averaged, 2‐D momentum and continuity equations. The FSC method is a predictor–corrector scheme, in which an intermediate free surface elevation is first calculated from the vertically integrated continuity equation after an intermediate, longitudinal velocity distribution is determined from the momentum equation. In the finite difference equation for the intermediate velocity, the vertical eddy viscosity term and the bottom‐ and sidewall friction terms are discretized implicitly, while the pressure gradient term, convection terms, and the horizontal eddy viscosity term are discretized explicitly. The intermediate free surface elevation is then adjusted by solving a FSC equation before the intermediate velocity field is corrected. The finite difference scheme is simple and can be easily implemented in existing laterally averaged 2‐D models. It is unconditionally stable with respect to gravitational waves, shear stresses on the bottom and side walls, and the vertical eddy viscosity term. It has been tested and validated with analytical solutions and field data measured in a narrow, riverine estuary in southwest Florida. Model simulations show that this numerical scheme is very efficient and normally can be run with a Courant number larger than 10. It can be used for rivers where the upstream bed elevation is higher than the downstream water surface elevation without any problem. Copyright © 2003 John Wiley & Sons, Ltd.     

A fully coupled numerical technique for 2D bed morphology calculations

A fully coupled two‐dimensional subcritical and/or supercritical, viscous, free‐surface flow numerical model is developed to calculate bed variations in alluvial channels. Vertically averaged free‐surface flow equations in conjunction with sediment transport equation are numerically solved using an explicit finite‐volume scheme using transformed grid in order to handle complex geometry fluvial problems. Convergence is accelerated with use of a multi‐grid technique. Firstly the capabilities of the proposed method are demonstrated by analyzing subcritical and supercritical hydrodynamic flows. Thereafter, an analysis of one‐ and two‐dimensional flows is performed referring to aggradation and scouring. For all reported test cases the computed results compare reasonably well with measurements as well as with other numerical solutions. The method is stable, reliable and accurate handling a variety of sediment transport equations with rapid changes of sediment transport at the boundaries. Copyright © 2002 John Wiley & Sons, Ltd.     

An implicit three‐dimensional numerical model to simulate transport processes in coastal water bodies

A three‐dimensional baroclinic numerical model has been developed to compute water levels and water particle velocity distributions in coastal waters. The numerical model consists of hydrodynamic, transport and turbulence model components. In the hydrodynamic model component, the Navier–Stokes equations are solved with the hydrostatic pressure distribution assumption and the Boussinesq approximation. The transport model component consists of the pollutant transport model and the water temperature and salinity transport models. In this component, the three‐dimensional convective diffusion equations are solved for each of the three quantities. In the turbulence model, a two‐equation k–ϵ formulation is solved to calculate the kinetic energy of the turbulence and its rate of dissipation, which provides the variable vertical turbulent eddy viscosity. Horizontal eddy viscosities can be simulated by the Smagorinsky algebraic sub grid scale turbulence model. The solution method is a composite finite difference–finite element method. In the horizontal plane, finite difference approximations, and in the vertical plane, finite element shape functions are used. The governing equations are solved implicitly in the Cartesian co‐ordinate system. The horizontal mesh sizes can be variable. To increase the vertical resolution, grid clustering can be applied. In the treatment of coastal land boundaries, the flooding and drying processes can be considered. The developed numerical model predictions are compared with the analytical solutions of the steady wind driven circulatory flow in a closed basin and of the uni‐nodal standing oscillation. Furthermore, model predictions are verified by the experiments performed on the wind driven turbulent flow of an homogeneous fluid and by the hydraulic model studies conducted on the forced flushing of marinas in enclosed seas. Copyright © 2000 John Wiley & Sons, Ltd.