A finite volume solution of wave forces on submarine pipelines

In this paper, a numerical model is established for simulating the wave forces on a submarine pipeline. A set of two-dimensional Navier–Stokes equations is discretized numerically with a finite volume method in a moving mesh system. After each time step, the mesh is modified according to the changed wave surface boundary. The deffered correction second-order upwind scheme (SUDC) is adopted here to discretize the convective fluxes. The effects of the clearance between the pipeline and the seabed, water depth and wave height on wave forces are studied, respectively. The results by the numerical simulation agree well with the experimental data and theory value.     

Blockage effect of oscillatory flow past a fixed cylinder

A numerical study of the effect of the width of the computational domain on viscous oscillatory flow past a circular cylinder has been conducted, for Keulegan–Carpenter numbers ranging between 0.1 and 6 at a fixed frequency parameter equal to 50. The finite element method was used for the solution of the Navier–Stokes equations, in the formulation where the stream function and the vorticity are the field variables. Simulations for blockage ratios in the range between 0.10 and 0.50 were performed assuming frictionless flow at the outer boundaries, the blockage ratio being defined as the cylinder diameter divided by the width of the solution domain. The first set of simulations was carried out for a constant stream function along the horizontal boundaries. Then the procedure was repeated, for stream function values at the outer boundaries derived from the irrotational solution around a circular cylinder. This boundary condition relieves considerably the blockage effect on the flow pattern and on the drag coefficient of the in-line force.     

Three-dimensional numerical simulation of solitary wave run-up using the IB method

Although the finite difference method is computationally efficient, it is acknowledged to be inferior when dealing with flow-over on structures with a complex geometry because of its rectilinear grid system. Therefore, we developed a numerical procedure that can cope with flow over structures with complex shapes while, at the same time, retaining the simplicity and efficiency of a rectilinear grid system. We used the immersed boundary method, which involves application of immersed boundary forces at solid boundaries rather than conventional boundary conditions, to investigate wave interactions with coastal structures in a three-dimensional numerical wave tank by solving the Navier–Stokes equations for two-phase flows. We simulated the run-up of a solitary wave around a circular island. Maximum run-up heights were computed around the island and compared with available laboratory measurements and previous numerical results. The three-dimensional features of the run-up process were analyzed in detail and compared with those of depth-integrated equations models.     

Interaction of linear waves with an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth

The linear water wave scattering and radiation by an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth is investigated by use of the multipole expansion method. The diffracted and radiated potentials are expressed as a linear combination of infinite multipoles placed at the centre of each cylinder with unknown coefficients to be determined by the cylinder boundary conditions. Analytical expressions for wave forces, hydrodynamic coefficients, reflection and transmission coefficients and energies are derived. Comparisons are made between the present analytical results and those obtained by the boundary element method, and some examples are presented to illustrate the hydrodynamic behavior of multiple horizontal circular cylinders in a two-layer fluid. It is found that for two submerged circular cylinders the influence of the fluid density ratio on internal-mode wave forces is more appreciable than surface-mode wave forces, and the periodic oscillations of hydrodynamic results occur with the increase of the distance between two cylinders; for four submerged circular cylinders the influence of adding two cylinders on the wave forces of the former cylinders is small in low and high wave frequencies, but the influence is appreciable in intermediate wave frequencies.     

Interactions between nonlinear water waves and non-wall-sided 3D structures

Interactions between water waves and non-wall-sided cylinders are analyzed based on velocity potential theory with fully nonlinear boundary conditions on the free surface and the body surface. The finite element method (FEM) is adopted together with a 3D mesh generated through an extension of a 2D Delaunay grid on a horizontal plane along the depth. The linear matrix equation for the velocity potential is constructed by imposing the governing equation and boundary conditions through the Galerkin method and is solved through an iterative method. By imposing the gradient of the potential equal to the velocity, the Galerkin method is used again to obtain the velocity field in the fluid domain. Simulations are made for bottom mounted and truncated cylinders with flare in a numerical tank. Periodic waves and wave groups are generated by a piston type wave maker mounted on one end of the tank. Results are obtained for forces, wave profiles and wave runups. Further simulations are made for a cylinder with flare subjected to forced motion in otherwise still open water. Results are provided for surge and heave motion in different amplitudes, and for a body moving in a circular path in the horizontal plane. Comparisons are made in several cases with the results obtained from the second order solution in the time domain.     

Numerical simulation of solitary wave scattering by a circular cylinder array

The interaction of a solitary wave with an array of surface-piercing vertical circular cylinders is investigated numerically. The wave motion is modeled by a set of generalized Boussinesq equations. The governing equations are discretized using a finite element method. The numerical model is validated against the experimental data of solitary wave reflection from a vertical wall and solitary wave scattering by a vertical circular cylinder respectively. The predicted wave surface elevation and the wave forces on the cylinder agree well with the experimental data. The numerical model is then employed to study solitary wave scattering by arrays of two circular cylinders and four circular cylinders respectively. The effect of wave direction on the wave forces and the wave runup on the cylinders is quantified.     

Least Square Finite Element Method for Viscous Splitting of Unsteady Incompressible Navier–Stokes Equations

In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton's method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C~0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.     

Viscous Effects on Wave Forces on A Submerged Horizontal Circular Cylinder

Numerical simulations are carried out for wave action on a submerged horizontal circular cylinder by means of a viscous fluid model, and it is focused on the examination of the discrepancies between the viscous fluid results and the potential flow solutions. It is found that the lift force resulted from rotational flow on the circular cylinder is always in anti-phase with the inertia force and induces the discrepancies between the results. The influence factors on the magnitude of the lift force, especially the correlation between the stagnation-point position and the wave amplitude, and the effect of the vortex shedding are investigated by further examination on the flow fields around the cylinder. The viscous numerical calculations at different wave frequencies showed that the wave frequency has also significant influence on the wave forces. Under higher frequency and larger amplitude wave action, vortex shedding from the circular cylinder will appear and influence the wave forces on the cylinder substantially.     

三维物体上二阶波浪绕射的实时模拟

采用二阶时域理论对非线性波浪在任意三维物体周围的绕射问题进行了研究，对自由表面边界条件进行Taylor级数展开，应用摄动展开可以建立相应的边值问题，而且此边值问题的计算域不随时间变化，运用基于B－样条的边界元方法求解每一时刻的波浪场，二阶自由表面边界条件在时间上进行数值积分，在自由表面加了一个人工阻尼层以避免波浪的反射，速度势分解为已知的入射势和未知的散射势，初始条件采用二阶Stokes波浪场，通过加入物体表面边界条件，得到散射势在时间和空间上的发展，本文对圆柱所受规则波的二阶波浪力和波浪爬高进行了计算，数值结果表明此理论计算准确，效率高，数值稳定。     

Wave forces acting on a submerged horizontal circular cylinder in oblique waves and on a vertical cylinder in deep waves

This paper presents a method of estimating wave forces acting on a submerged horizontal circular cylinder fixed in oblique waves.The experiments show that drag and inertia coefficients in beam sea are available for calculating the wave forces in oblique waves.Wave forces exerted on a vertical circular cylinder in deep waves are also investigated.The experimental results show that wave forces acting on the vertical cylinder coincide approximately with hydrodynamic forces acting on a submerged circular cylinder in an oscillating fluid.     

A finite element model of interaction between viscous free surface waves and submerged cylinders

This paper describes the simulation of the flow of a viscous incompressible Newtonian liquid with a free surface. The Navier–Stokes equations are formulated using a streamline upwind Petrov–Galerkin scheme, and solved on a Q-tree-based finite element mesh that adapts to the moving free surface of the liquid. Special attention is given to fitting the mesh correctly to the free surface and solid wall boundaries. Fully non-linear free surface boundary conditions are implemented. Test cases include sloshing free surface motions in a rectangular tank and progressive waves over submerged cylinders.     

Wave deformation and vortex generation in water waves propagating over a submerged dike

The Navier–Stokes equations and the exact free surface boundary conditions are solved to simulate wave deformation and vortex generation in water waves propagating over a submerged dike. Incident waves are generated by a piston-type wavemaker set up in the computational domain. Numerical results are compared with experimental data in order to confirm the validity of the numerical model. The fast Fourier transform and a wave resolution technique are applied to decompose the transformed waves and the higher harmonics. Effects of different parameters on wave transformation and vortex generation are studied systematically. These parameters include the Ursell number, the Keulegan–Carpenter number, the water depth ratio, the Reynolds number, the length aspect ratio of the dike, and the type of dike.     

Wave forces on a large, horizontal submerged cylinder

A numerical boundary integral equation method combined with a non-linear time stepping procedure is used for the calculation of wave forces on a large, submerged, horizontal circular cylinder. As the method is based on potential theory, all computations are performed in the inertia dominated domain, that is, for small Keulegan-Carpenter numbers. Computations are carried out for the Eulerian mean current under wave trough level equal to zero. When the cylinder is moved towards the sea bed the computations show that the inertia coefficients increase significantly, which is associated with a blockage effect. Furthermore, the effect of the wave steepness is reduced when the submergence of the cylinder is increased. In the vicinity of the free water surface the vertical inertia coefficient is highly dependent upon the wave steepness, which tends to reduce it, whereas the horizontal inertia coefficient is only slightly dependent on the wave steepness. Computations are also carried out for cylinder diameters comparable with the wave length. Finally, inertia coefficients computed by the present method are compared with some analytical results by Ogilvie [(1963), First and second order forces on a cylinder submerged under a free surface. J. Fluid Mech. 16, 451–472]. As long as the assumptions leading to Ogilvie's theory are fulfilled (cylinder radius small compared to the wave length), the results are quite similar.     

Finite element analysis of flow control using porous media

A finite element model is established for simulating flow in and out a porous media. The extended Darcy equation inside the porous media and the Navier–Stokes equations in the fluid are coupled via the continuity condition at the interface between the two media. The model is firstly validated against the analytical and the numerical results available in literature. Then it is applied to simulate flow past a circular cylinder covered by a porous layer. The effect of the porous layer on the reduction of lift coefficient is investigated numerically. It is found that the lift reduction can be achieved by properly choosing the porous material. However, the amount of reduction greatly depends on the Reynolds number, the permeability and the Forchheimer coefficient.     

Nonlinear response of membranes to ocean waves using boundary and finite elements

The behavior of a highly deformable membrane to ocean waves was studied by coupling a nonlinear boundary element model of the fluid domain to a nonlinear finite element model of the membrane. The hydrodynamic loadings induced by water waves are computed assuming large body hydrodynamics and ideal fluid flow and then solving the transient diffraction/radiation problem. Either linear waves or finite amplitude waves can be assumed in the model and thus the nonlinear kinematic and dynamic free surface boundary conditions are solved iteratively. The nonlinear nature of the boundary condition requires a time domain solution. To implicitly include time in the governing field equation, Volterra's method was used. The approach is the same as the typical boundary element method for a fluid domain where the governing field equation is the starting point. The difference is that in Volterra's method the time derivative of the governing field equation becomes the starting point.The boundary element model was then coupled through an iterative process to a finite element model of membrane structures. The coupled model predicts the nonlinear interaction of nonlinear water waves with highly deformable bodies. To verify the coupled model a large scale test was conducted in the OH Hinsdale wave Research Laboratory at Oregon State University on a 3-ft-diameter fabric cylinder submerged in the wave tank. The model data verified the numerical prediction of the structure displacements and of the changes in the wave field.The boundary element model is an ideal modeling technique for modeling the fluid domain when the governing field equations is the Laplace equation. In this case the nonlinear boundary element model was coupled with a finite element model of membrane structures, but the model could have been coupled with other finite element models of more rigid structures, such as a pontoon floating breakwater.     

Time stepping solutions of the two-dimensional nonlinear wave radiation problem

The two-dimensional nonlinear time domain free surface flow problem is analysed using potential flow theory. The problem is solved by a time marching method. At each time step two numerical approaches are used. One is based on the boundary element method in the complex plane. The complex potential is assumed to vary linearly within each element and the solution is obtained by imposing the boundary conditions at the nodes of the elements. The other approach is based on the finite element formulation. Triangular elements and linear shape functions are used. The solution is obtained by the Galerkin method. Numerical results are obtained for the wave elevation generated by a vertical wave maker. Results are also provided for a circular cylinder oscillating below the free surface. For these cases the finite element method is found to provide substantially more efficient computations than the boundary element method using equivalent discretizations.     

A finite difference method for effective treatment of mild-slope wave equation subject to non-reflecting boundary conditions

A numerical model for coastal water wave motion that includes an effective method for treatment of non-reflecting boundaries is presented. The second-order one-way wave equation to approximate the non-reflecting boundary condition is found to be excellent and it ensures a very low level of reflection for waves approaching the boundary at a fairly wide range of the incidence angle. If the Newman approximation is adopted, the resulting boundary condition has a unique property to allow the free propagation of wave components along the boundary. The study is also based on a newly derived mild-slope wave equation system that can be easily made compatible to the one-way wave equation. The equation system is theoretically more accurate than the previous equations in terms of the mild-slope assumption. The finite difference method defined on a staggered grid is employed to solve the basic equations and to implement the non-reflecting boundary condition. For verification, the numerical model is then applied to three coastal water wave problems including the classical problem of plane wave diffraction by a vertical circular cylinder, the problem of combined wave diffraction and refraction over a submerged hump in the open sea, and the wave deformation around a detached breakwater. In all cases, the numerical results are demonstrated to agree very well with the relevant analytical solutions or with experimental data. It is thus concluded that the numerical model proposed in this study is effective and advantageous.     

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.     

Time-domain simulation of berthing problem between FPSO and shuttle tanker in waves

A time-domain simulation method based on potential flow model has been developed to investigate the berthing problem between two floating bodies in wave. The boundary value problem is formulated with respect to an earth-fixed coordinate system because the relative positions of the two vessels continuously change during the berthing operation. The classical finite element method is used to solve the Laplace equation in the fluid domain with moving boundary. The linearized free-surface boundary conditions are integrated in time by applying 4th-order Adams–Bashforth–Moulton method. A simple re-mesh algorithm with local and global mesh systems is introduced to update mesh by considering large horizontal movement of the berthing vessel. The developed numerical method is used to investigate the berthing problem between a FPSO and shuttle tanker in waves. The focus is on the wave-induced motion response during the berthing process. The characteristics of the motion responses in berthing operation are examined with various wave frequencies, berthing speeds and wave headings.