Fully nonlinear simulations of wave resonance by an array of cylinders in vertical motions

The finite element method(FEM) is employed to analyze the resonant oscillations of the liquid confined within multiple or an array of floating bodies with fully nonlinear boundary conditions on the free surface and the body surface in two dimensions.The velocity potentials at each time step are obtained through the FEM with 8-node quadratic shape functions.The finite element linear system is solved by the conjugate gradient(CG) method with a symmetric successive overelaxlation(SSOR) preconditioner.The waves at the open boundary are absorbed by the combination of the damping zone method and the Sommerfeld-Orlanski equation.Numerical examples are given by an array of floating wedgeshaped cylinders and rectangular cylinders.Results are provided for heave motions including wave elevations,profiles and hydrodynamic forces.Comparisons are made in several cases with the results obtained from the second order solution in the time domain.It is found that the wave amplitude in the middle region of the array is larger than those in other places,and the hydrodynamic force on a cylinder increases with the cylinder closing to the middle of the array.     

Clear water scour around a finite array of cylinders

This experimental study presents clear-water scour and deposition patterns around hexagonal arrays of circular cylinders in steady flow conditions. Understanding the scour processes around such configurations could facilitate the design of several hydraulic and marine engineering structures, such as bridge piers and piles. The flow alteration caused by the examined porous obstacles depends on the solid volume fraction of the obstacles and on the angle of attack of the incoming flow, due to the limited number of cylinders constituting the array. Flume experiments with erodible bed were carried out for four array densities (solid volume fractions: 0.14, 0.20, 0.32 and 0.56) under three different orientations (regular, angled and staggered configurations). The scour/deposition characteristics were obtained by means of laser scanner and the results were compared to solid cylinders of equal circumambient diameter. Different angles of attack of the incoming flow lead to different blockage ratios, which have direct impact on the scour characteristics and deposition patterns. The arrays with the higher solid volume fraction generated scour/deposition patterns similar to solid cylinder, while in the arrays with the lower solid volume fractions, local scour around the individual small cylinders became evident. Finally, considering that the load bearing capacity of a pier basically depends on the area of its cross-section, a comparison of the maximum induced scour depth and volume by the cylinder arrays and the solid cylinder with equal solid cross-sectional area is presented, in order to introduce an alternative pier configuration that induces less scour. The results showed that the array of cylinders could generate 27% less scour volume and 22% less scour depth compared to its single solid cylinder counterpart.     

Solitary wave interactions with an array of two vertical cylinders

This paper addresses a numerical investigation of nonlinear waves interactions with an array of two surface-piercing vertical cylinders and the corresponding nonlinear hydrodynamic loads on each individual cylinder. The primary interest of this study is concentrated on the problem of three-dimensional scattering of solitary waves by cylinder arrays and the nonlinear interactions between scattered waves. The theoretical model adopted for simulation is the generalized Boussinesq two-equation model. The boundary-fitted coordinate transformation and multiple-grid technique are utilized here to simplify the computation domain and to facilitate the applications of the boundary conditions on the cylinder surfaces. The velocity potential, free-surface elevation and subsequent evolution of the scattered wave field are numerically evaluated. The hydrodynamic forces on each cylinder during wave impact are also determined. A study of the sheltering effect by the neighboring structures on wave loads is conducted. It is found that the presence of the neighboring cylinder has shown significant influence on the wave loads and the scattering of the primary incident waves. For two transversely arranged cylinders, the transverse force coefficient increases as the separation distance decreases.     

Interaction between waves and an array of floating porous circular cylinders

The present study theoretically as well as experimentally investigates the interaction between waves and an array of porous circular cylinders with or without an inner porous plate based on the linear wave theory.To design more effective floating breakwaters,the transmission rate of waves propagating through the array is evaluated.Each cylinder in the array is partly made of porous materials.Specifically,it possesses a porous sidewall and an impermeable bottom.In addition,an inner porous plate is horizontally fixed inside the cylinders.It dissipates the wave more effectively and eliminates the sloshing phenomenon.The approach suggested by Kagemoto and Yue(1986) is adopted to solve the multiple-scatter problem,while a hierarchical interaction theory is adopted to deal with hydrodynamic interactions among a great number of bodies,which efficiently saves computation time.Meanwhile,a series of model tests with an array of porous cylinders is performed in a wave basin to validate the theoretical work and the calculated results.The draft of the cylinders,the location of the inner porous plate,and the spacing between adjacent cylinders are also adjusted to investigate their effects on wave dissipation.     

Water wave interaction with an array of bottom-mounted surface-piercing porous cylinders

The interaction of water waves with arrays of bottom-mounted, surface-piercing circular cylinders is investigated theoretically. The sidewall of each cylinder is porous and thin. Under the assumptions of potential flow and linear wave theory, a semi-analytical solution is obtained by an eigenfunction expansion approach first proposed for impermeable cylinders by Spring and Monkmeyer (1974), and later simplified by Linton and Evans (1990). Analytical expressions are developed for the wave motion in the exterior and all interior fluid regions. Numerical results are presented which illustrate the effects of various wave and structural parameters on the hydrodynamic loads and the diffracted wave field. It is found that the porosity of the structures may result in a significant reduction in both the hydrodynamic loads experienced by the cylinders and the associated wave runup.     

Second-order resonance among an array of two rows of vertical circular cylinders

Experimental and theoretical works are conducted on the first-order and second-order wave interactions among an array of two rows of vertical circular cylinders.We show that both the first-order and the second-order resonance of dynamic free-surface displacement could in some cases be quite relevant among an array of two rows of vertical cylinders not only in theory but also in real phenomena. In some cases, on the other hand, we show that large free-surface displacement predicted theoretically is quite attenuated in reality probably due to viscous dissipation.     

Interaction of cnoidal waves with an array of vertical concentric porous cylinders

An array of large concentric porous cylinder arrays is mounted in shallow water exposed to cnoidal waves. The interactions between waves and cylinders are studied theoretically using an eigenfunction expansion approach. Semi-analytical solutions of hydrodynamic loads and wave run-up on each cylinder are obtained using first approximation to cnoidal waves. The square array configuration of four-legged identical concentric porous cylinder is investigated in present study. Numerical results reveal the variation of dimensionless wave force and wave run-up on individual cylinder with angle of incidence, porosity parameter, spacing between outer and inner cylinders, spacing between concentric porous cylinders and wave parameter. Different mechanism of wave force is found under different range of scattering parameter.     

多向不规则波浪作用下群墩结构上爬高的计算研究

波浪的方向分布对波浪的传播及其与工程结构物的作用都具有明显影响,目前现有的研究大多是基于单向波浪进行的。为了研究方向分布对群墩结构上的爬高影响,基于规则波浪与群墩作用的理论解,结合多向不规则波浪的造波方法,建立了多向不规则波浪与群墩作用的计算模型,同时进行了物理模型试验对模型的有效性进行了验证。系统地对群墩周围及表面上的波浪爬高进行了计算分析,结果表明,方向分布对波浪爬高具有较大的影响,且不同位置处的影响并不相同,在实际的工程设计中如果按照单向波浪计算,可能低估或者高估群墩周围的爬高。     

波浪对阻抗型圆柱群的绕射

应用柱面上阻抗型边界条件分析单柱和圆柱群上波面高度和总波浪力的分布规律。随着阻抗值减小波面高度明显减小,柱上总的波浪力也随之减小。由于柱间干扰和阻抗影响,波浪力在一定的波数尺寸范围内出现急降区,阻抗型柱的变化平缓,而刚性光滑壁面柱降低显著,说明阻抗越小柱间波浪散射强度越小。     

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.     

Simulation of flows around cylinders by a Lagrangian vortex scheme

A method is presented for calculating flows around cylinders of arbitrary shape, which may be in an array. Vorticity is discretised and is moved in a Lagrangian time-stepping procedure by the vortex-in-cell method. The zero-velocity surface condition is maintained at each time step by representing the surface by a vortex sheet through a boundary integral calculation. Vorticity diffusion is simulated through random walks. Results are presented for an isolated circular cylinder in steady and oscillatory incident flows.     

水下平板与波浪相互作用完全非线性数值模拟(英文)

利用完全非线性数值波浪水槽技术研究水下平板与波浪的相互作用。假定水下平板厚度极薄、刚性,位于有限水深并且非常接近自由水面。应用四阶龙格库塔方法追踪每一时刻的波面形状,采用阻尼层来吸收反射波以保证算法的稳定性,同时引入平滑和重组的方法抑制自由表面控制点的较高梯度。通过对波浪与浮动圆柱相互作用的数值模拟证实了数值波浪水槽方法的有效性,计算结果与线性理论吻合良好。在波浪数值水槽方法中引入造波板模拟波浪产生并与水下平板发生相互作用,应用傅立叶解析方法对波面变形、波浪力作了分析。结果表明在板非常接近自由水面的情况下会表现出现很强的非线性,揭示了线性理论的局限性。     

A quasi-cloaking phenomenon to reduce the wave drift force on an array of adjacent floating bodies

The mooring of offshore floating structures, such as offshore platforms, in large waves against drift forces and rotational moments is a challenging problem in offshore engineering. To accurately investigate such problems, called positioning problems, the time-averaged steady forces of the second order known as the wave drift forces must be taken into account. Fortunately, a cloaking phenomenon occurs under certain conditions and dramatically reduces the wave drift force acting on such a floating body, as previously reported by several researchers. In the diffraction problem of water waves, cloaking refers to the condition where there is no scattering in the form of radial outgoing waves. The reduction of wave drift force on a truncated cylinder with the occurrence of cloaking phenomenon has been numerically and experimentally confirmed. In this paper, the arrangement of several small circular cylinders at regular intervals in a circle concentric with a fixed floating body is considered as an effective means of reducing the wave drift force. Using a combination of a higher-order boundary element method (HOBEM) and wave interaction theory, the influences of the geometric parameters of the outer surrounding cylinders on the wave drift force and the total scattered-wave energy are systematically investigated and discussed. A quasi-cloaking phenomenon is first found and reported in the present study, which is beneficial and flexible for application in practical engineering. More than one quasi-cloaking trigger (where a trigger is an occurrence condition) can be found simply by varying the distance between the inner and outer floating bodies.     

Flow around circular and noncircular cylinders by a low-order panel method

Time-dependent cross-flow was studied around cylinders with circular and noncircular cross-sections. The numerical approach for the analysis was a low-order panel method based on constant source and dipole values along each panel. The method was previously used successfully for several applications, such as calculation of the added mass and damping coefficients. In simulating the viscous time-dependent flow around the cylinder, the time-dependent wake feature of the code was used. For the circular and D-cylinders, the results agreed well with the experiments. Suggestions for improving the results for T-cylinders with angle of attack are included.     

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.     

五圆柱体结构群绕流特性及互扰效应研究

基于Fluent流体计算平台,运用大涡模拟方法对亚临界雷诺数Re=3900下“X”形排列五圆柱体结构群三维绕流特性进行研究,主要分析来流攻角α与间距比L/D两个关键参数对五圆柱体结构群的尾流区三维涡结构演化与流体力系数的影响,并揭示其内在流动互扰机理。研究表明:来流攻角和间距比的变化对五圆柱体结构群流动控制及互扰效应的影响显著。在小间距比工况下,观察到柱体群间隙区域内流体高速流动的现象,导致五圆柱体之间的互扰作用十分强烈。间隙流对中间圆柱体和下游圆柱体有较强的冲击作用,对其表面的流体力分布特性有显著的影响。另外,大间距比工况下,当α=0°与L/D≥5.0工况时,柱体群尾流效应强于其间隙流效应。当α=22.5°与L/D=7.0时,位于下游与中间处的圆柱体流体绕流特性存在较大差异。而当α=45°与L/D≥6.0时,位于上游与中间处的圆柱体尾流区均会产生正负交替的漩涡结构。     

无限水深聚焦波完全非线性数值模拟

基于势流理论提出一种新的高阶边界元方法对无限水深的聚焦波浪进行完全非线性数值模拟.自由水面满足完全非线性边界条件,模拟波浪的非线性效果可以达到更高阶.利用镜像原理,建立一种全新的格林函数应用到无限水深的数值波浪水槽中,以致于两无限深水槽侧壁的积分可以被排除.为了产生相应的入射波和吸收出流波浪,一个由点源组成的造波装置被布置于计算域内,同时人工阻尼层被用来吸引出流波浪,由波浪聚焦的方法得到极限波浪.通过开展线性和完全非线性聚焦波浪的数值实验及与理论解对比,验证本数值模型可以用来模拟无限深水域的极限波浪,且在出流边界没有反射.     

不同倒角半径下方柱绕流的数值模拟及水动力特性研究

为了研究不同倒角半径对方柱绕流特性的影响,采用有限体积法,模拟了雷诺数Re为22 500、倒角半径为0.1D(D为方柱边长的长度)、0.2D和0.3D时方柱的绕流过程。方柱近壁面采用增强壁面函数,模型采用SST k–?湍流模型。根据模拟结果给出了不同倒角半径下方柱的流场涡量图以及阻力系数Cd和升力系数Cl;利用快速傅里叶变换法得到斯托罗哈数St。结果表明,倒角半径的增加改变了方柱的分离点,使得尾流区长度增加,旋涡尺度减小;Cd和Cl的振动幅值呈现先减小后增大的趋势,倒角半径为0.1D和0.2D时方柱受力较小,不存在倒角时方柱受力较大,倒角半径为0.3D时方柱受力最大;随着倒角半径的增加,柱体截面形式越接近圆形,斯托罗哈数逐渐增大,漩涡脱落频率更快。