Flow field analysis of submerged horizontal plate type breakwater

Submerged horizontal plate can be considered as a new concept breakwater. In order to reveal the wave elimination mechanism of this type breakwater, boundary element method is utilized to investigate the velocity field around plate carefully. The flow field analysis shows that the interaction between incident wave and reverse flow caused by submerged plate will lead to the formation of wave elimination area around both sides of the plate. The velocity magnitude of flow field has been reduced and this is the main reason of wave elimination.     

Study on scattering wave force of horizontal and vertical plate type breakwaters

The interaction between wave and horizontal and vertical plates is investigated by the boundary element method,and the relations of wave exciting force with plate thickness,submergence and length are obtained.It is found that:1) The efficient wave exciting force exists while plate submergence is less than 0.5 m,and the plate is very thin with order O(0.005 m).2) The maximum heave wave exciting force exists,and it is the main factor for surface and submerged horizontal plate while the roll force can be ignored.3) The maximum sway wave exciting force exists,it is the main factor for surface or submerged vertical plate,and the roll force is about 20 times of horizontal plate.     

Scattering of water wave by a submerged horizontal plate and a submerged permeable breakwater

Based on a two-dimensional linear water wave theory, this study develops the boundary element method (BEM) to examine normally incident wave scattering by a fixed, submerged, horizontal, impermeable plate and a submerged permeable breakwater in water of finite depth. Numerical results for the transmission coefficients are also presented. In addition, the numerical technique's accuracy is demonstrated by comparing the numerical results with previously published numerical and experimental ones. According to that comparison, the transmission coefficient relies not only on the submergence of the horizontal impermeable plate and the height of the permeable breakwater, but also on the distance between horizontal plate and permeable breakwater. Results presented herein confirm that the transmission coefficient is minimum for the distance approximately equal to four times the water depth.     

Wave motion over a submerged breakwater with an upper horizontal porous plate and a lower horizontal solid plate

This study examines the hydrodynamic performance of a modified two-layer horizontal-plate breakwater. The breakwater consists of an upper submerged horizontal porous plate and a lower submerged horizontal solid plate. By means of the matched eigenfunction expansion method, a linear analytical solution is developed for the interaction of water waves with the structure. Then the reflection coefficient, the transmission coefficient, the energy-loss coefficient and the wave forces acting on the plates are calculated. The numerical results obtained for limiting cases are exactly the same as previous predictions for a single submerged horizontal solid plate and a single submerged horizontal porous plate. Numerical results show that with a suitable geometrical porosity of the upper plate, the uplift wave forces on both plates can be controlled at a low level. Numerical results also show that the transmission coefficient will be always small if the dimensionless plate length (plate length versus incident wavelength) exceeds a certain moderate value. This is rather significant for practical engineering, as the incident wavelength varies over a wide range in practice. Moreover, it is found that the hydrodynamic performance of the present structure may be further enhanced if the lower plate is also perforated.     

Wave interaction with a perforated wall breakwater with a submerged horizontal porous plate

This study examines the hydrodynamic performance of a new perforated-wall breakwater. The breakwater consists of a perforated front wall, a solid back wall and a submerged horizontal porous plate installed between them. The horizontal porous plate enhances the stability and wave-absorbing capacity of the structure. An analytical solution based on linear potential theory is developed for the interaction of water waves with the new proposed breakwater. According to the division of the structure, the whole fluid domain is divided into three sub-domains, and the velocity potential in each domain is obtained using the matched eigenfunction method. Then the reflection coefficient and the wave forces and moments on the perforated front wall and the submerged horizontal porous plate are calculated. The numerical results obtained for limiting cases are exactly the same as previous predictions for a perforated-wall breakwater with a submerged horizontal solid plate [Yip, T.L., Chwang, A.T., 2000. Perforated wall breakwater with internal horiontal plate. Journal of Engineering Mechanics ASCE 126 (5), 533–538] and a vertical wall with a submerged horizontal porous plate [Wu, J.H., Wan, Z.P., Fang, Y., 1998. Wave reflection by a vertical wall with a horizontal submerged porous plate. Ocean Engineering 25 (9), 767–779]. Numerical results show that with suitable geometric porosity of the front wall and horizontal plate, the reflection coefficient will be always rather small if the relative wave absorbing chamber width (distance between the front and back walls versus incident wavelength) exceeds a certain small value. In addition, the wave force and moment on the horizontal plate decrease significantly with the increase of the plate porosity.     

没水倾斜板式防波堤消波性能分析

提出了一种将波浪中倾斜板问题等效化简为波浪中水平板单元组的方法,该方法建立在使用分离变量法求解水工结构边值问题的基础上,并使用伽辽金法精确求解连续边界条件,确定考虑衰减波态的速度势函数,从而求解没水倾斜板结构的消波性能。等效化简法计算精度于边界元法相当,且计算单元数量少、开销低。基于二维线性势波理论,对没水倾斜板式防波堤消波性能分析显示,没水板的倾斜角度、没水深度与板长是结构消波性能的控制因素:没水倾斜板防波堤的消波性能优于没水水平板防波堤的消波性能,随着没水板结构的倾斜角度增大,没水倾斜板结构的波浪透射系数显著减小,且长板优于短板,浅板优于深板;与前人的水槽实验对比显示,相对没水深度与波陡影响结构的消波性能,且波陡造成的波浪破碎贡献了显著的波能消耗。该结论对板式防波堤的结构配置、优化设计有重要意义。     

Wave motion over a breakwater system of a horizontal plate and a vertical porous wall

A two-dimensional analytical solution is presented to study the reflection and transmission of linear water waves propagating past a submerged horizontal plate and through a vertical porous wall. The velocity potential in each fluid domain is formulated using three sets of orthogonal eigenfunctions and the unknown coefficients are determined from the matching conditions. Wave elevations and hydrodynamic forces acting on the porous wall are computed. Reflection and transmission coefficients are presented to examine the performance of the breakwater system. The present analytical solutions are found in fairly good agreement with the available laboratory data. The results indicate that the plate length, the porous-effect, the gap between plate and porous wall, and the submerged depth of the plate all show a significant influence on the reflected and transmitted wave fields. It is also interesting to note that the submerged plate plays an important role in reducing the transmitted wave height, especially for long incident waves.     

一种起伏板防波堤最优支撑刚度计算方法

带有弹簧支撑的水平板结构,可在波浪作用下起伏运动进行消浪,作为防波堤具有广泛的应用前景,但如何设计支撑刚度仍是当前研究的重点课题。基于黏性流理论建立数值波浪水槽,开展规则波与起伏板防波堤相互作用的模拟,得到起伏板透反射系数随支撑刚度的变化曲线,并分析最小透射系数对应刚度下起伏板运动相位对消浪的影响。为快速得到最小透射系数对应的支撑刚度,根据弹簧—阻尼—质量模型,建立一种以运动相位为输入反推起伏板最优支撑刚度的计算方法,并将推导结果与数值结果进行比较。结果表明,该方法推导的最优支撑刚度值与数值最优值接近,且可省去遍历取最优值的步骤,可为起伏板防波堤的支撑刚度设计提供参考。     

波浪与起伏水平板防波堤相互作用数值模拟

利用自主研发的基于紧致插值曲线CIP(constrained interpolation profile)方法的数学模型,开展规则波与起伏水平板防波堤相互作用的数值模拟研究。模型在笛卡尔直角坐标下建立,以CIP方法为流场基本求解器,分步求解Navier-Stokes方程,利用高精度的流体体积类型的THINC/SW (tangent of hyperbola for interface capturing with slope weighting)方法重构自由液面,采用浸入边界IBM(immersed boundary method)方法处理波浪与起伏板防波堤的耦合作用问题,通过动量源项造波方法模拟波浪的产生。重点关注波浪的浅水变形和板两端涡旋脱落的非线性现象,分析不同潜深、波要素下的板周围流场分布、板的运动响应和波浪的反透射系数。结果表明:起伏水平板主要通过能量反射、板上浅水变形和板两端的涡脱落消能,能有效减小板后波高,具有作为防波堤的可行性。     

梳式沉箱翼缘板设计理论的研究

结合模型试验 ,分析了梳式沉箱翼缘板的静力和自振特性。研究了三种不同支承形式的翼缘板在最不利荷载下的内力分布 ,为梳式沉箱翼缘板的设计提供了依据     

Lateral and consolidation behaviors of seabed-type breakwater for very soft ground

A seabed-type of breakwater applicable to very soft ground without the need for soil improvement is newly developed. This type of soft-ground breakwater is expected to ensure sufficient lateral resistance and prevent excessive consolidation settlement due to self-weight of the breakwater. In this paper, lateral and consolidation behaviors of soft-ground breakwater were investigated by performing model tests and finite element simulations. The results revealed that the bottom wall and buoyant box, which are the main features of soft-ground breakwater, contribute to the increase in lateral resistance and to the control of the consolidation settlements, respectively, and that Terzaghi's consolidation theory could be conservatively adopted in deriving the consolidation settlements of soft-ground breakwater proposed herein.     

Minimization of drift force on a floating cylinder by optimizing the flexural rigidity of a concentric annular plate

The deployment of suitable configurations of mutually interacting floating bodies for efficiently controlling their hydrodynamic interactions towards the reduction of the wave drift forces and, thus, of the mooring lines’ loads, has, nowadays, gained a great scientific interest. In this paper, the hydrodynamic behaviour of a floating cylinder and a concentric annular flexible plate is analysed in the frequency domain aiming at the minimization of the drift forces acting on the cylinder by optimizing the flexural rigidity of the plate. The diffraction/radiation problem is solved using a higher-order boundary element method. The analysis is implemented assuming that both floating bodies oscillate freely in heave, while for the plate, flexible modes are, additionally, considered for describing its structural deformations. The required modes shapes are determined in vacuum (“dry” mode superposition approach) through analytical expressions. The flexural rigidity of the plate, D, is optimized at a specific wave number using a real-coded genetic algorithm. Initially, results are compared with numerical results of other investigators for the case of two rigid concentric floating cylinders. Next, extended results are presented, focusing on the effect of D, including its optimum value, on various physical quantities describing the behaviour of both the cylinder and the plate. Contrary to the isolated cylinder, the presence of the plate introduces sharp peaks in the variation pattern of the drift force of the cylinder, bounded at specific wave numbers, where resonance of the seiche mode of water motion in the annular cavity or of specific flexible modes of the plate occurs. However, by reducing D to its optimum value, the cylinder’s drift force obtains practically zero values at the target wave number, due to an efficient improvement of the wave field in the annular cavity around the cylinder. Moreover, a great reduction of the drift force compared to the isolated cylinder is achieved in the subsequent high frequency range.     

Spar平台垂荡板水动力系数计算与分析

Spar平台带有厚度很薄的垂荡板,使用边界元方法计算时存在着准奇异积分的数值计算问题。采用自适应高斯积分法对其准奇异积分问题进行处理,较少单元即可达到很高精度。在势流理论下对Spar平台垂荡板的水动力系数进行了计算,并与实验结果进行了对比,发现真实流体中垂荡板的附加质量约为势流理论下附加质量计算结果的1.1倍,而辐射阻尼在总阻尼中所占的百分比很小,此时粘性阻尼成为垂荡板的主要阻尼来源,在工程设计中不可忽略。     

多个淹没水平圆柱诱发高阶谐波特性的数值与试验研究

基于势流理论采用时域高阶边界元方法建立了模拟非线性波浪与淹没水平双、三圆柱作用的数值水槽模型,其中采用混合欧拉-拉格朗日方法更新自由水面,四阶Runga-Kutta方法进行时间步进。利用两点法分离得到潜体下游高阶自由波,进而研究潜体间距、布置形式等对高倍频自由波的影响。同时在物理水槽内开展多潜体高阶谐波物理模型试验,并将试验结果与数值结果进行对比,吻合较好。研究发现:淹没双圆柱下游高阶谐波幅值随潜体间距呈现周期性振荡变化,其振荡的重现距离约为波长一半。而淹没三圆柱下游的高阶谐波随第一个间距呈周期性振荡变化,随第二个间距的增大而减小。     

波浪与结构物作用分析的一种高阶边界元方法--自由项和柯西主值积分的直接数值计算

采用直接数值计算方法计算了势流问题高阶边界元方法中的自由项系数和柯西主值积分,建立了波浪与结构物作用的一种高阶边界元方法.通过算例研究了物体表面上固角系数的计算精度和不同网格剖分、不同阶高斯积分点对柯西主值积分的影响.对截断圆柱上的波浪作用力与解析解做了对比,发现本方法具有很高的计算精度,随网格的加密迅速收敛于解析解.     

部分反射直墙前二维物体的绕射和辐射问题

基于映像理论将部分反射直墙前物体的散射问题,等效于开敞水域中原物体的散射和关于直墙映像体散射的线性叠加进行求解。采用高阶边界元方法建立了部分反射直墙前二维任意形状物体波浪绕射和辐射问题的数值分析模型,通过与已发表的海底方箱和淹没圆柱结果的对比验证了数值模型的准确性。应用该模型研究了直墙反射系数幅值及相位、方箱与墙间距离等参数对水面方箱上波浪激振力、附加质量和辐射阻尼的影响。结果表明：直墙反射系数幅值越大,波浪激振力、附加质量和辐射阻尼的波动越大,附加质量在一些频率下出现负值;相位角的变化会改变波浪激振力、附加质量和辐射阻尼曲线的偏移,在低频区对升沉附加质量有显著影响;方箱距离直墙越远,方箱上的波浪激振力、附加质量和辐射阻尼随波数振荡的频率越快,峰值频率向低频侧移动。     

A Unique Solvable Higher Order BEM for Wave Diffraction and Radiation

- For the discretization of higher order elements, this paper presents a modified integral domain method to remove the irregular frequencies inherited in the integral equation of wave diffraction and radiation from a surface-piercing body. The set of over-determined linear equations obtained from the method is modified into a normal set of linear equations by superposing a set of linear equations with zero solutions. Numerical experiments have also been carried out to find the optimum choice of the size of the auxiliary domain and the discretization on it.     

An Efficient Model for Transient Surface Waves in Both Finite and Infinite Water Depth

A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropriate wave theory. At the downstream boundary, an artificial damping zone is used to prevent wave reflection back into the computational domain. Using the image Green function in the whole fluid domain, the integrations on the two lateral walls and bottom are excluded. The simulation results on extreme wave elevations in finite and infinite water-depths are compared with experimental results and second-order analytical solutions respectively. The wave kinematics is also discussed in the present study.     

An Efficient Model for Transient Surface Waves in Both Finite and Infinite Water Depths

A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropr...