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.     

Fully three-dimensional direct numerical simulation of a plunging breakerSimulation numérique directe tridimensionnelle du déferlement plongeant

The scope of this paper is to show the results obtained for simulating three-dimensional breaking waves by solving the Navier–Stokes equations in air and water. The interface tracking is achieved by a Lax–Wendroff TVD scheme (Total Variation Diminishing), which is able to handle interface reconnections. We first present the equations and the numerical methods used in this work. We then proceed to the study of a three-dimensional plunging breaking wave, using initial conditions corresponding to unstable periodic sinusoidal waves of large amplitudes. We compare the results obtained for two simulations, a longshore depth perturbation has been introduced in the solution of the flow equations in order to see the transition from a two-dimensional velocity field to a fully three-dimensional one after plunging. Breaking processes including overturning, splash-up and breaking induced vortex-like motion beneath the surface are presented and discussed. To cite this article: P. Lubin et al., C. R. Mecanique 331 (2003).     

Analysis of efficient preconditioned defect correction methods for nonlinear water waves

Robust computational procedures for the solution of non‐hydrostatic, free surface, irrotational and inviscid free‐surface water waves in three space dimensions can be based on iterative preconditioned defect correction (PDC) methods. Such methods can be made efficient and scalable to enable prediction of free‐surface wave transformation and accurate wave kinematics in both deep and shallow waters in large marine areas or for predicting the outcome of experiments in large numerical wave tanks. We revisit the classical governing equations are fully nonlinear and dispersive potential flow equations. We present new detailed fundamental analysis using finite‐amplitude wave solutions for iterative solvers. We demonstrate that the PDC method in combination with a high‐order discretization method enables efficient and scalable solution of the linear system of equations arising in potential flow models. Our study is particularly relevant for fast and efficient simulation of non‐breaking fully nonlinear water waves over varying bottom topography that may be limited by computational resources or requirements. To gain insight into algorithmic properties and proper choices of discretization parameters for different PDC strategies, we study systematically limits of accuracy, convergence rate, algorithmic and numerical efficiency and scalability of the most efficient known PDC methods. These strategies are of interest, because they enable generalization of geometric multigrid methods to high‐order accurate discretizations and enable significant improvement in numerical efficiency while incuring minimal storage requirements. We demonstrate robustness using such PDC methods for practical ranges of interest for coastal and maritime engineering, that is, from shallow to deep water, and report details of numerical experiments that can be used for benchmarking purposes. Copyright © 2014 John Wiley & Sons, Ltd.     

A meshless method for numerical simulation of depth‐averaged turbulence flows using a k‐ϵ model

A three‐dimensional numerical model is developed to analyze free surface flows and water impact problems. The flow of an incompressible viscous fluid is solved using the unsteady Navier–Stokes equations. Pseudo‐time derivatives are introduced into the equations to improve computational efficiency. The interface between the two phases is tracked using a volume‐of‐fluid interface tracking algorithm developed in a generalized curvilinear coordinate system. The accuracy of the volume‐of‐fluid method is first evaluated by the multiple numerical benchmark tests, including two‐dimensional and three‐dimensional deformation cases on curvilinear grids. The performance and capability of the numerical model for water impact problems are demonstrated by simulations of water entries of the free‐falling hemisphere and cone, based on comparisons of water impact loadings, velocities, and penetrations of the body with experimental data. For further validation, computations of the dam‐break flows are presented, based on an analysis of the wave front propagation, water level, and the dynamic pressure impact of the waves on the downstream walls, on a specific container, and on a tall structure. Extensive comparisons between the obtained solutions, the experimental data, and the results of other numerical simulations in the literature are presented and show a good agreement. Copyright © 2015 John Wiley & Sons, Ltd.     

Semi‐coupled air/water immersed boundary approach for curvilinear dynamic overset grids with application to ship hydrodynamics

For many problems in ship hydrodynamics, the effects of air flow on the water flow are negligible (the frequently called free surface conditions), but the air flow around the ship is still of interest. A method is presented where the water flow is decoupled from the air solution, but the air flow uses the unsteady water flow as a boundary condition. The authors call this a semi‐coupled air/water flow approach. The method can be divided into two steps. At each time step the free surface water flow is computed first with a single‐phase method assuming constant pressure and zero stress on the interface. The second step is to compute the air flow assuming the free surface as a moving immersed boundary (IB). The IB method developed for Cartesian grids (Annu. Rev. Fluid Mech. 2005; 37 :239–261) is extended to curvilinear grids, where no‐slip and continuity conditions are used to enforce velocity and pressure boundary conditions for the air flow. The forcing points close to the IB can be computed and corrected under a sharp interface condition, which makes the computation very stable. The overset implementation is similar to that of the single‐phase solver (Comput. Fluids 2007; 36 :1415–1433), with the difference that points in water are set as IB points even if they are fringe points. Pressure–velocity coupling through pressure implicit with splitting of operators or projection methods is used for water computations, and a projection method is used for the air. The method on each fluid is a single‐phase method, thus avoiding ill‐conditioned numerical systems caused by large differences of fluid properties between air and water. The computation is only slightly slower than the single‐phase version, with complete absence of spurious velocity oscillations near the free surface, frequently present in fully coupled approaches. Validations are performed for laminar Couette flow over a wavy boundary by comparing with the analytical solution, and for the surface combatant model David Taylor Model Basin (DTMB) 5512 by comparing with Experimental Fluid Dynamics (EFD) and the results of two‐phase level set computations. Complex flow computations are demonstrated for the ONR Tumblehome DTMB 5613 with superstructure subject to waves and wind, including 6DOF motions and broaching in SS7 irregular waves and wind. Copyright © 2008 John Wiley & Sons, Ltd.     

A level set technique applied to unsteady free surface flows

An unsteady Navier–Stokes solver for incompressible fluid is coupled with a level set approach to describe free surface motions. The two‐phase flow of air and water is approximated by the flow of a single fluid whose properties, such as density and viscosity, change across the interface. The free surface location is captured as the zero level of a distance function convected by the flow field. To validate the numerical procedure, two classical two‐dimensional free surface problems in hydrodynamics, namely the oscillating flow in a tank and the waves generated by the flow over a bottom bump, are studied in non‐breaking conditions, and the results are compared with those obtained with other numerical approaches. To check the capability of the method in dealing with complex free surface configurations, the breaking regime produced by the flow over a high bump is analyzed. The analysis covers the successive stages of the breaking phenomenon: the steep wave evolution, the falling jet, the splash‐up and the air entrainment. In all phases, numerical results qualitatively agree with the experimental observations. Finally, to investigate a flow in which viscous effects are relevant, the numerical scheme is applied to study the wavy flow past a submerged hydrofoil. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical simulation of pollutant transport acted by wave for a shallow water sea bay

A numerical model of pollutant transport acted by water waves on a shallow‐water mild‐slope beach is established in this study. The numerical model is combined with a wave propagation model, a multiple wave‐breaking model, a wave‐induced current model and a pollutant convection–dispersion model. The wave propagation model is based on the higher‐order approximation of parabolic mild‐slope equation which can be used to simulate the wave refraction, diffraction and breaking in a large area of near‐shore zone combined with the wave‐breaking model. The wave‐induced current model is established using the concept of the radiation stress and considering the effect of bottom resistance caused by waves. The numerical model is verified by laboratory experiment results of regular and irregular waves over two mild beaches with different slopes. The numerical results agree well with experimental results. The numerical model has been applied in the near‐shore zone of Bohai bay in China. It is concluded that pollutant transport parallel to the shoreline due to the action of waves, which will induce serious pollution on the beach. Copyright © 2005 John Wiley & Sons, Ltd.     

Boussinesq cut‐cell model for non‐linear wave interaction with coastal structures

Boussinesq models describe the phase‐resolved hydrodynamics of unbroken waves and wave‐induced currents in shallow coastal waters. Many enhanced versions of the Boussinesq equations are available in the literature, aiming to improve the representation of linear dispersion and non‐linearity. This paper describes the numerical solution of the extended Boussinesq equations derived by Madsen and Sørensen (Coastal Eng. 1992; 15 :371–388) on Cartesian cut‐cell grids, the aim being to model non‐linear wave interaction with coastal structures. An explicit second‐order MUSCL‐Hancock Godunov‐type finite volume scheme is used to solve the non‐linear and weakly dispersive Boussinesq‐type equations. Interface fluxes are evaluated using an HLLC approximate Riemann solver. A ghost‐cell immersed boundary method is used to update flow information in the smallest cut cells and overcome the time step restriction that would otherwise apply. The model is validated for solitary wave reflection from a vertical wall, diffraction of a solitary wave by a truncated barrier, and solitary wave scattering and diffraction from a vertical circular cylinder. In all cases, the model gives satisfactory predictions in comparison with the published analytical solutions and experimental measurements. Copyright © 2007 John Wiley & Sons, Ltd.     

A new σ‐transform based Fourier‐Legendre‐Galerkin model for nonlinear water waves

This paper presents a new spectral model for solving the fully nonlinear potential flow problem for water waves in a single horizontal dimension. At the heart of the numerical method is the solution to the Laplace equation which is solved using a variant of the $\sigma$‐transform. The method discretizes the spatial part of the governing equations using the Galerkin method and the temporal part using the classical fourth‐order Runge‐Kutta method. A careful investigation of the numerical method's stability properties is carried out, and it is shown that the method is stable up to a certain threshold steepness when applied to nonlinear monochromatic waves in deep water. Above this threshold artificial damping may be employed to obtain stable solutions. The accuracy of the model is tested for: (i) highly nonlinear progressive wave trains, (ii) solitary wave reflection, and (iii) deep water wave focusing events. In all cases it is demonstrated that the model is capable of obtaining excellent results, essentially up to very near breaking.     

Numerical simulation of 3D free surface flows,with multiple incompressible immiscible phases. Applications to impulse waves

A numerical method for the solution to the density‐dependent incompressible Navier–Stokes equations modeling the flow of N immiscible incompressible liquid phases with a free surface is proposed. It allows to model the flow of an arbitrary number of liquid phases together with an additional vacuum phase separated with a free surface. It is based on a volume‐of‐fluid approach involving N indicator functions (one per phase, identified by its density) that guarantees mass conservation within each phase. An additional indicator function for the whole liquid domain allows to treat boundary conditions at the interface between the liquid domain and a vacuum. The system of partial differential equations is solved by implicit operator splitting at each time step: first, transport equations are solved by a forward characteristics method on a fine Cartesian grid to predict the new location of each liquid phase; second, a generalized Stokes problem with a density‐dependent viscosity is solved with a FEM on a coarser mesh of the liquid domain. A novel algorithm ensuring the maximum principle and limiting the numerical diffusion for the transport of the N phases is validated on benchmark flows. Then, we focus on a novel application and compare the numerical and physical simulations of impulse waves, that is, waves generated at the free surface of a water basin initially at rest after the impact of a denser phase. A particularly useful application in hydraulic engineering is to predict the effects of a landslide‐generated impulse wave in a reservoir. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical modeling of wave interactions with coastal structures by a constrained interpolation profile/immersed boundary method

A high‐order difference method based multiphase model is proposed to simulate nonlinear interactions between water wave and submerged coastal structures. The model is based on the Navier–Stokes equations using a constrained interpolation profile (CIP) method for the flow solver, and employs an immersed boundary method (IBM) for the treatment of wave–structure interactions. A more accurate interface capturing scheme, the volume of fluid/weighed line interface calculation (VOF/WLIC) scheme, is adopted as the interface capturing method. A series of computations are performed to verify the application of the model for simulations of fluid interaction with various structures. These problems include flow over a fixed cylinder, water entry of a circular cylinder and solitary waves passing various submerged coastal structures. Computations are compared with the available analytical, experimental and other numerical results and good agreement is obtained. The results of this study demonstrate the accuracy and applications of the proposed model to simulate the nonlinear flow phenomena and capture the complex free surface flow. Copyright © 2015 John Wiley & Sons, Ltd.     

A σ‐coordinate non‐hydrostatic model with embedded Boussinesq‐type‐like equations for modeling deep‐water waves

A σ‐coordinate non‐hydrostatic model, combined with the embedded Boussinesq‐type‐like equations, a reference velocity, and an adapted top‐layer control, is developed to study the evolution of deep‐water waves. The advantage of using the Boussinesq‐type‐like equations with the reference velocity is to provide an analytical‐based non‐hydrostatic pressure distribution at the top‐layer and to optimize wave dispersion property. The σ‐based non‐hydrostatic model naturally tackles the so‐called overshooting issue in the case of non‐linear steep waves. Efficiency and accuracy of this non‐hydrostatic model in terms of wave dispersion and nonlinearity are critically examined. Overall results show that the newly developed model using a few layers is capable of resolving the evolution of non‐linear deep‐water wave groups. Copyright © 2009 John Wiley & Sons, Ltd.     

Coupled ghost fluid/two‐phase level set method for curvilinear body‐fitted grids

A coupled ghost fluid/two‐phase level set method to simulate air/water turbulent flow for complex geometries using curvilinear body‐fitted grids is presented. The proposed method is intended to treat ship hydrodynamics problems. The original level set method for moving interface flows was based on Heaviside functions to smooth all fluid properties across the interface. We call this the Heaviside function method (HFM). The HFM requires fine grids across the interface. The ghost fluid method (GFM) has been designed to explicitly enforce the interfacial jump conditions, but the implementation of the jump conditions in curvilinear grids is intricate. To overcome these difficulties a coupled GFM/HFM method was developed in which approximate jump conditions are derived for piezometric pressure and velocity and pressure gradients based on exact continuous velocity and stress and jump in momentum conditions with the jump in density maintained but continuity of the molecular and turbulent viscosities imposed. The implementation of the ghost points is such that no duplication of memory storage is necessary. The level set method is adopted to locate the air/water interface, and a fast marching method was implemented in curvilinear grids to reinitialize the level set function. Validations are performed for three tests: super‐ and sub‐critical flow without wave breaking and an impulsive plunging wave breaking over 2D submerged bumps, and the flow around surface combatant model DTMB 5512. Comparisons are made against experimental data, HFM and single‐phase level set computations. The proposed method performed very well and shows great potential to treat complicated turbulent flows related to ship flows. Copyright © 2007 John Wiley & Sons, Ltd.     

A depth‐integrated non‐hydrostatic finite element model for wave propagation

This paper presents a two‐dimensional finite element model for simulating dynamic propagation of weakly dispersive waves. Shallow water equations including extra non‐hydrostatic pressure terms and a depth‐integrated vertical momentum equation are solved with linear distributions assumed in the vertical direction for the non‐hydrostatic pressure and the vertical velocity. The model is developed based on the platform of a finite element model, CCHE2D. A physically bounded upwind scheme for the advection term discretization is developed, and the quasi second‐order differential operators of this scheme result in no oscillation and little numerical diffusion. The depth‐integrated non‐hydrostatic wave model is solved semi‐implicitly: the provisional flow velocity is first implicitly solved using the shallow water equations; the non‐hydrostatic pressure, which is implicitly obtained by ensuring a divergence‐free velocity field, is used to correct the provisional velocity, and finally the depth‐integrated continuity equation is explicitly solved to satisfy global mass conservation. The developed wave model is verified by an analytical solution and validated by laboratory experiments, and the computed results show that the wave model can properly handle linear and nonlinear dispersive waves, wave shoaling, diffraction, refraction and focusing. Copyright © 2013 John Wiley & Sons, Ltd.     

Global volume conservation in unsteady free surface flows with energy absorbing far‐end boundaries

A wave absorption filter for the far‐end boundary of semi‐infinite large reservoirs is developed for numerical simulation of unsteady free surface flows. Mathematical model is based on finite volume solution of the Navier–Stokes equations and depth‐integrated continuity equation to track the free surface. The Sommerfeld boundary condition is applied at the far‐end of the truncated computational domain. A dissipation zone is formed by applying artificial pressure on water surface to dissipate the kinetic energy of the outgoing waves. The computational scheme is tested to verify the conservation of total fluid volume in the domain for long simulation durations. Combination of the Sommerfeld boundary and dissipation zone can effectively minimize reflections and prevent cumulative changes in total fluid volume in the domain. Solitary wave, nonlinear periodic waves and irregular waves are simulated to illustrate the numerical developments. Earthquake excited surface waves and nonlinear hydrodynamic pressures in a dam–reservoir are computed. Copyright © 2009 John Wiley & Sons, Ltd.     

A robust high‐resolution finite volume scheme for the simulation of long waves over complex domains

The propagation, runup and rundown of long surface waves are numerically investigated, initially in one dimension, using a well‐balanced high‐resolution finite volume scheme. A conservative form of the nonlinear shallow water equations with source terms is solved numerically using a high‐resolution Godunov‐type explicit scheme coupled with Roe's approximate Riemann solver. The scheme is also extended to handle two‐dimensional complex domains. The numerical difficulties related to the presence of the topography source terms in the model equations along with the appearance of the wet/dry fronts are properly treated and extended. The resulting numerical model accurately describes breaking waves as bores or hydraulic jumps and conserves volume across flow discontinuities. Numerical results show very good agreement with previously presented analytical or asymptotic solutions as well as with experimental benchmark data. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical simulation of free‐surface flow using the level‐set method with global mass correction

A new numerical method that couples the incompressible Navier–Stokes equations with the global mass correction level‐set method for simulating fluid problems with free surfaces and interfaces is presented in this paper. The finite volume method is used to discretize Navier–Stokes equations with the two‐step projection method on a staggered Cartesian grid. The free‐surface flow problem is solved on a fixed grid in which the free surface is captured by the zero level set. Mass conservation is improved significantly by applying a global mass correction scheme, in a novel combination with third‐order essentially non‐oscillatory schemes and a five stage Runge–Kutta method, to accomplish advection and re‐distancing of the level‐set function. The coupled solver is applied to simulate interface change and flow field in four benchmark test cases: (1) shear flow; (2) dam break; (3) travelling and reflection of solitary wave and (4) solitary wave over a submerged object. The computational results are in excellent agreement with theoretical predictions, experimental data and previous numerical simulations using a RANS‐VOF method. The simulations reveal some interesting free‐surface phenomena such as the free‐surface vortices, air entrapment and wave deformation over a submerged object. Copyright © 2009 John Wiley & Sons, Ltd.     

The numerical study of roll‐waves in inclined open channels and solitary wave run‐up

In this paper, we introduce a finite‐volume kinetic BGK scheme and its applications to the study of roll and solitary waves. The current scheme is based on the numerical solution of the gas‐kinetic Bhatnagar–Gross–Krook model in the flux evaluation across each cell interface. An intrinsic connection between the BGK model and time‐dependent, non‐linear, non‐homogeneous shallow‐water equations enables us to solve shallow‐water equations automatically with our kinetic scheme. The analytical solution, experimental measurements, and numerical calculations for problems associated with roll‐waves down an inclined open channel and solitary waves incident on a sloped beach are also presented. Copyright © 2005 John Wiley & Sons, Ltd.     

波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟

波浪破碎卷入气体易对建筑物受力产生压力振荡, 了解波浪作用下建筑物附近掺气水流的运动特性是精确计算建筑物受力的前提. 基于OpenFOAM开源程序包和修正速度入口造波方法建立三维数值波浪水槽, 模型采用S-A IDDES湍流模型进行湍流封闭, 并采用修正的VOF 方法捕捉自由液面, 数值模拟了规则波在1:10的光滑斜坡上与直立结构物的相互作用过程, 重点分析了结构物附近的水动力和掺气水流运动特性. 结果表明, 建立的数值模型能精确地捕捉波浪作用下直立结构物附近的自由液面的变化以及气泡输运过程, 较好地描述气体卷入所形成的气腔形态以及多气腔之间的融合、分裂等过程; 波浪与直立结构物相互作用产生强湍动掺气水流, 其运动过程十分复杂; 掺气流体输运过程中水气界面周围一直伴随着涡的存在, 其中, 气泡的分裂与周围正负涡量剪切作用密切相关, 且其输运轨迹主要受周围流场的影响; 研究揭示了结构物附近湍动能与掺气特性的关系, 发现波浪作用下直立结构物附近湍动能的分布与掺气水流特征参数(气泡数量、空隙率)整体呈现一定的线性关系.      

An anti‐diffusive volume of fluid method for interfacial fluid flows

A volume of fluid (VOF) method is developed combining a first‐order limited downwind scheme with higher order accurate schemes. The method is characterized by retaining a sharp fluid interface and a reduction in numerical diffusion near the interface, but avoids complicated geometrical reconstruction as occurs in most volume tracing algorithms. To demonstrate the accuracy and robustness of the method, a selection of numerical experiments are presented involving a pure advection problem, a water wave impact caused by a dam breaking and liquid sloshing in a partially filled tank. Copyright © 2011 John Wiley & Sons, Ltd.