Regular wave pressures and forces on submerged pipelines near a sloping boundary

The hydrodynamic pressures induced by regular waves around the circumference of a pipeline normal to the wave direction and near a rigid bed of slope 1:10 have been investigated in a wave flume. The pressures were integrated to obtain the force time history, from which the peak horizontal and vertical forces are evaluated. The maximum and root mean square horizontal and transverse force coefficients are correlated with the Keulegan–Carpenter (KC) number. The effect of the distance between the sloping bed and the pipeline on the force coefficients is discussed. The force coefficients are found to decrease with an increase in KC number and with the decrease in the relative clearance of the pipeline from the boundary. In addition, the reflection characteristics of the sloping bed in the presence of the pipeline as a function of surf similarity parameter and their comparison with the results from existing literature are also reported. The details of the model setup, experimental procedure, results and discussion are presented in this paper.     

Propagation of a solitary wave over rigid porous beds

The unsteady two-dimensional Navier–Stokes equations and Navier–Stokes type model equations for porous flows were solved numerically to simulate the propagation of a solitary wave over porous beds. The free surface boundary conditions and the interfacial boundary conditions between the water region and the porous bed are in complete form. The incoming waves were generated using a piston type wavemaker set up in the computational domain. Accuracy of the numerical model was verified by comparing the numerical results with the theoretical solutions. The main characteristics of the flow fields in both the water region and the porous bed were discussed by specifying the velocity fields. Behaviors of boundary layer flows in both fluid and porous bed regions were also revealed. Effects of different parameters on the wave height attenuation were studied and discussed. The results of this numerical model indicate that for the investigated incident wave as the ratio of the porous bed depth to the fluid depth exceeds 10, any further increase of the porous bed depth has no effect on wave height attenuation.     

Numerical simulation of wave overtopping of coastal structures using the non-linear shallow water equations

A one-dimensional high-resolution finite volume model capable of simulating storm waves propagating in the coastal surf zone and overtopping a sea wall is presented. The model (AMAZON) is based on solving the non-linear shallow water (NLSW) equations. A modern upwind scheme of the Godunov-type using an HLL approximate Riemann solver is described which captures bore waves in both transcritical and supercritical flows. By employing a finite volume formulation, the method can be implemented on an irregular, structured, boundary-fitted computational mesh. The use of the NLSW equations to model wave overtopping is computationally efficient and practically flexible, though the detailed structure of wave breaking is of course ignored. It is shown that wave overtopping at a vertical wall may also be approximately modelled by representing the wall as a steep bed slope. The AMAZON model solutions have been compared with analytical solutions and laboratory data for wave overtopping at sloping and vertical seawalls and good agreement has been found. The model requires more verification tests for irregular waves before its application as a generic design tool.     

A two-dimensional finite volume hydrodynamic model for coastal areas: Model development and validation

A two dimensional implicit finite volume scheme for solving the shallow-water equations is developed. The effects of the Coriolis force, surface wind stress, and waves are included. A non-uniform rectilinear forward staggered grid is used with Cartesian coordinates. The time integration is performed using the Euler implicit technique. The convective flux is treated using the deferred correction method. The viscous terms are discretized using a second order central difference approximation. The SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm is used for coupling the velocity components and the water elevation gradient for the water level correction. The system of equations is solved sequentially using the Strongly Implicit Procedure (SIP). To simulate wave driven current, a phase averaged wave model is used first to simulate wave transformation and calculate radiation stresses. The performance of the developed model is validated for different sources of external forces and different combinations of boundary conditions. The validation cases include tidal circulation in a harbor and wave induced currents behind a breakwater parallel to the coastline. The model is finally applied to simulate the flow pattern in a closed artificial lagoon and along the coastline near Damietta Port located along the Northern coast of Egypt. Results of the developed model agree well with the published results for the considered cases.     

Prediction of hydrodynamic forces on submarine pipelines using an improved Wake II Model

The hydrodynamic force model for prediction of forces on submarine pipelines as described includes flow history effect (wake effects) and time dependence in the force coefficients. The wake velocity correction is derived by using a closed-form solution to the linearized Navier–Stokes equations for oscillatory flow. This is achieved by assuming that the eddy viscosity in the wake is only time dependent and of a harmonic sinusoidal form. The forces predicted by the new Wake (Wake II) Model have been compared to Exxon Production Research Company Wake Model in terms of time histories (force shape) and magnitudes of peak forces. Overall, the model predictions by the Wake II Model are satisfactory and represent a substantial improvement over the predictions of the conventional models. The conventional force models representing adaptations of Morison's equation with ambient velocity and constant coefficients give predictions that are in poor agreement with the measurements especially for the lift force component. The Wake II Force Model can be used for submarine pipeline on-bottom stability design calculations for regular waves with various pipe diameters.     

The application of wave induced forces to a two-dimensional finite element long wave hydrodynamic model

This paper describes the modification of a two-dimensional finite element long wave hydrodynamic model in order to predict the net current and water levels attributable to the influences of waves. Tests examine the effects of the application of wave induced forces, including comparisons to a physical experiment. An example of a real river system is presented with comparisons to measured data, which demonstrate the importance of simulating the combined effects of tides and waves upon hydrodynamic behavior.     

Physical modeling of untrenched submarine pipeline instability

Wave-induced instability of untrenched pipeline on sandy seabed is a ‘wave–soil–pipeline’ coupling dynamic problem. To explore the mechanism of the pipeline instability, the hydrodynamic loading with U-shaped oscillatory flow tunnel is adopted, which is quite different from the previous experiment system. Based on dimensional analysis, the critical conditions for pipeline instability are investigated by altering pipeline submerged weight, diameter, soil parameters, etc. Based on the experimental results, different linear relationships between Froude number (Fr) and non-dimensional pipeline weight (G) are obtained for two constraint conditions. Moreover, the effects of loading history on the pipeline stability are also studied. Unlike previous experiments, sand scouring during the process of pipe’s losing stability is detected in the present experiments. In addition, the experiment results are compared with the previous experiments, based on Wake II model for the calculation of wave-induced forces upon pipeline. It shows that the results of two kinds of experiments are comparable, but the present experiments provide better physical insight of the wave–soil–pipeline coupling effects.     

A numerical model of wave overtopping on seadikes

This paper describes the development of a numerical model for wave overtopping on seadikes. The model is based on the flux-conservative form of the nonlinear shallow water equations (NLSW) solved with a high order total variation diminishing (TVD), Roe-type scheme. The goal is to reliably predict the hydrodynamics of wave overtopping on the dike crest and along the inner slope, necessary for the breach modelling of seadikes. Besides the mean overtopping rate, the capability of simulating individual overtopping events is also required. It is shown theoretically that the effect of wave breaking through the drastic motion of surface rollers in the surfzone is not sufficiently described by the conventional nonlinear shallow water equations, neglecting wave setup from the mean water level and thus markedly reducing the model predictive capacity for wave overtopping. This is significantly improved by including an additional source term associated with the roller energy dissipation in the depth-averaged momentum equation. The developed model has been validated against four existing laboratory datasets of wave overtopping on dikes. The first two sets are to validate the roller term performance in improving the model prediction of wave overtopping of breaking waves. The last two sets are to test the model performance under more complex but realistic hydraulic and slope geometric conditions. The results confirm the merit of the supplemented roller term and also demonstrate that the model is robust and reliable for the prediction of wave overtopping on seadikes.     

Nonlinear effect on inertia component of wave forces on a cylinder

Ｎｏｎｌｉｎｅａｒｅｆｆｅｃｔｏｎｉｎｅｒｔｉａｃｏｍｐｏｎｅｎｔｏｆｗａｖｅｆｏｒｃｅｓｏｎａｃｙｌｉｎｄｅｒ￥ＬｉＹａｎｂａｏａｎｄＳｏｎｇＲｅｎｇ（ＤｅｐａｒｔｍｅｎｔｏｆＨｙｄｒａｕｌｉｃＥｎｇｉｎｅｅｒｉｎｇ，ＴｉａｎｊｉｎＵｎｉｖｅｒｓｉ...     

A polar method for obtaining wave resonating quadruplets in finite depths

A polar method for obtaining wave resonating quadruplets {K₁, K₂, K₃, K₄} in the computation of nonlinear wave–wave interaction source term of the wave model is presented with results for both deep and finite water depths. The method first determines the end radial points of the locus equation for K₂, for each set of input wave vectors (K₁, K₃) on the symmetry. The locus of K₂ (and hence K₄) is then traced in the anti-clockwise direction starting with the maximum radial point on the line of symmetry. It is shown that when k₃>k₁, the number of points on the locus varies when the orientations of the input wave vectors are changed and reduces when the difference in the magnitude of the input wave vectors is increased. A significant advantage in this method is that the angular increment on the locus for K₂ can be kept constant.     

Finite volume scheme for the solution of 2D extended Boussinesq equations in the surf zone

In this paper, a hybrid finite volume-finite difference scheme is applied to study surf zone dynamics. The numerical model solves the 2DH extended Boussinesq equations proposed by Madsen and Sørensen (1992) where nonlinear and dispersive effects are both relevant whereas it solves NSWE equations where nonlinearity prevails. The shock-capturing features of the finite volume method allow an intrinsic representation of wave breaking and runup; therefore no empirical (calibration) parameters are necessary. Comparison with laboratory measurements demonstrates that the proposed model can accurately predict wave height decay and mean water level setup, for both regular and solitary wave breaking on a sloping beach. The model is also applied to reproduce two-dimensional wave transformation and breaking over a submerged circular shoal, showing good agreement with experimental data.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

A numerical study of three-dimensional wave interaction with a square cylinder

In this study, a three-dimensional numerical model is used to study the wave interaction with a vertical rectangular pile. The model employs the large eddy simulation (LES) method to model the effect of small-scale turbulence. The velocity and vorticity fields around the pile are presented and discussed. The drag and inertial coefficients are calculated based on the numerical computation. The calculated coefficients are found to be in a reasonable range compared with the experimental data. Additional analyses are performed to assess the relative importance of drag and initial effects, which could be quantified by the force-related Keulegan and Carpenter (KC) number: KC_f=UT/(4πL). Here U is the maximum fluid particle velocity, T the wave period and L the length of structure aligned with the wave propagation direction. For small KC_f, the effective drag coefficient is proportional to 1/KC_f, provided the wavelength is much longer than the structural length. When wavelength is comparable to the structure dimension, the effective drag coefficient would be reduced significantly due the cancellation of forces, which has been demonstrated by numerical results.     

圆弧面防波堤波浪力初步研究

圆弧面防波堤是在半圆形防波堤基础上开发的一种新型防波堤。首先通过与半圆形防波堤相同条件下的波浪试验,检验圆弧面防波堤的稳定性,并利用数值波浪水槽对圆弧面防波堤的水力特性进行初步研究,探求造成圆弧面防波堤与半圆形防波堤波浪力差别的主要原因。通过圆弧面防波堤与半圆形防波堤波浪力的对比试验,提出了圆弧面防波堤波浪力的简化计算方法,以半圆堤正向水平波浪力乘以一修正系数,在堤顶淹没情况下修正系数可取1.3,在堤顶出水情况下修正系数可取1.1。     

High-fidelity simulation of regular waves based on multi-moment finite volume formulation and THINC method

The performance of interFoam (a widely used solver within OpenFOAM package) in simulating the propagation of water waves has been reported to be sensitive to the temporal and spatial resolution. To facilitate more accurate simulations, a numerical wave tank is built based on a Navier–Stokes model, which employs the VPM (volume-average/point-value multi-moment) scheme as the fluid solver and the THINC/QQ method (THINC method with quadratic surface representation and Gaussian quadrature) for the free-surface capturing. Simulations of regular waves in an intermediate water depth are conducted and the results are assessed via comparing with the analytical solutions. The performance of the present model and interFoam solver in simulating the wave propagation is systematically compared in this work. The results clearly demonstrate that compared with interFoam solver, the present model significantly improves the dissipation properties of the propagating wave, where the waveforms as well as the velocity distribution can be substantially maintained while the waves propagating over long distances even with large time steps and coarse grids. It is also shown that the present model requires much less computation time to reach a given error level in comparison with interFoam solver.     

Boussinesq-type modelling using an unstructured finite element technique

A model for solving the two-dimensional enhanced Boussinesq equations is presented. The model equations are discretised in space using an unstructured finite element technique. The standard Galerkin method with mixed interpolation is applied. The time discretisation is performed using an explicit three-step Taylor–Galerkin method. The model is extended to the surf and swash zone by inclusion of wave breaking and a moving boundary at the shoreline. Breaking is treated by an existing surface roller model, but a new procedure for the detection of the roller thickness is devised. The model is verified using four test cases and the results are compared with experimental data and results from an existing finite difference Boussinesq model.     

系泊船非线性波浪力时域计算:二维模型

为找到具有工程实用价值的港口系泊船波浪力的时域计算方法,建立了在港口中存在系泊船时非线性波浪力时域计算的垂直二维耦合模型:用Boussinesq方程计算船的两侧的外域,用欧拉方程计算船底面下的内域,两域在交界面处的连接条件是流量连续和压力相等.将复平面内的边界元方法应用于所研究问题,对耦合模型进行了验证.进行了相关模型实验,实验结果与数值计算结果比较表明这两种数值计算模型都具有满意的精度,但耦合模型的计算效率要远远高于边界元方法的计算效率.本耦合模型的数学处理简单,可适用于工程计算.     

Numerical modeling of wave runup and forces on an idealized beachfront house

In this paper, a numerical wave model based on the incompressible Reynolds-averaged Navier–Stokes (RANS) and k–ε equations is used to estimate the impact of a solitary wave on an idealized beachfront house located at different elevations on a plane beach. The locations of the free surface are reconstructed by volume of fluid (VOF) method. The model is satisfactorily tested against the experimental data of wave runup, and the analytical solution of wave forces on vertical walls. The time histories of wave profiles, forces, and overturning moments on the idealized house are demonstrated and analyzed. The variations of wave forces and overturning moments with the elevation of the idealized beachfront house are also investigated.     

A coastal ocean model of semi-implicit finite volume unstructured grid

A two-dimensional coastal ocean model based on unstructured C-grid is built, in which the momentum equation is discretized on the faces of each cell, and the continuity equation is discretized on the cell. The model is discretized by semi-implicit finite volume method, in that the free surface is semi-implicit and the bottom friction is implicit, thereby removing stability limitations associated with the surface gravity wave and friction. The remaining terms in the momentum equations are discretized explicitly by integral finite volume method and second-order Adams-Bashforth method. Tidal flow in the polar quadrant with known analytic solution is employed to test the proposed model. Finally, the performance of the present model to simulate tidal flow in a geometrically complex domain is examined by simulation of tidal currents in the Pearl River Estuary.     

Nonlinear wave profiles of wave–wave interaction in a finite water depth by fixed point approach

In this paper, a superposition of two periodic wave profiles in a finite water depth was investigated. This paper is focused on the improvement of a wave profile on the linear superposition of two waves. This improvement was realized by introducing an iterative method, which was based on a fixed point approach. Application of the fixed point approach to the wave superposition made it possible to obtain a wave profile of wave–wave interaction. The improved result of the wave profile was in good agreement with that of the nonlinear perturbation solution of the second order. It was interesting that the improved result revealed the higher-order nonlinear frequencies for two interacting Stokes waves while Dalzell's solution by a perturbation method could not predict them.