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

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

Interaction of ocean waves with floating and submerged horizontal flexible structures in three-dimensions

A three-dimensional general mathematical hydroelastic model dealing with the problem of wave interaction with a floating and a submerged flexible structure is developed based on small amplitude wave theory and linear structural response. The horizontal floating and submerged flexible structures are modelled with a thin plate theory. The linearized long wave equations based on shallow water approximations are derived and results are compared. Three-dimensional Green’s functions are derived using fundamental source potentials in water of finite and infinite depths. The expansion formulae associated with orthogonal mode-coupling relations are derived based on the application of Fourier transform in finite and infinite depths in case of finite width in three-dimensions. The usefulness of the expansion formula is demonstrated by analysing a physical problem of surface gravity wave interaction with a moored finite floating elastic plate in the presence of a finite submerged flexible membrane in three-dimensions. The numerical accuracy of the method is demonstrated by computing the complex values of reflected wave amplitudes for different modes of oscillation and mooring stiffness. Further, the effect of compressive force and modes of oscillations on a free oscillation hydroelastic waves in a closed channel of finite width and length for floating and submerged elastic plate system is analysed.     

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.     

内孤立波作用下竖直圆柱体的横向力研究

本文对内孤立波作用下圆柱体的水动力特性进行了数值研究。对不同直径、不同水深比、不同浸没深度的圆柱进行了水平方向受力分析,从多个因素出发,研究了内孤立波对水下竖直圆柱的水平作用情况。通过分析水平作用力和力矩的时程曲线,得到内孤立波行进过程中对于圆柱体的横向作用力和力矩与时间的关系。通过对比不同工况下内孤立波橫向作用力幅值,得到内孤立波作用力幅值与圆柱直径、水深比以及浸没深度的关系。结果表明,内孤立波横向作用力(力矩)幅值与圆柱直径呈线性相关关系,但与水深比相关关系不明显,随着圆柱浸没深度的增加,内孤立波横向作用力幅值呈现出先增大后减小的趋势。     

Hydroelastic Response of VLFS with An Attached Submerged Horizontal Solid/Porous Plate Under Wave Action

The hydroelastic responses of a submerged horizontal solid/porous plate attached at the front of a very large rectangular floating structure(VLFS) under wave action has been investigated in the context of linear water wave theory. Darcy's law is adopted to represent energy dissipation in pores. It is assumed that the porous plates are made of material with very fine pores so that the normal velocity across the perforated porous is linearly associated with the pressure drop. In the analytic method, the eigenfunction expansion-matching method(EEMM) for multiple domains is applied to solve the hydrodynamic problem and the elastic equation of motion is solved by the modal expansion method. The performance of the proposed submerged horizontal solid/porous plate can be significantly enhanced by selecting optimal design parameters, such as plate length, horizontal position, submerged depth and porosity. It is concluded that good damping effect can be achieved through installation of solid and porous plate.Porous plate has better damping effect at low frequencies, while solid plate has better damping effect at high frequencies. The optimal ratio of plate length to water depth is 0.25-0.375, and the optimal ratio of submerged depth to water depth is 0.09-0.181.     

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 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.     

Interaction mechanisms among waves,currents and a submerged plate

Wave-induced loads on a submerged plate, representative of submerged breakwater, coastal-bridge deck and a certain type of wave energy converter, in a uniform current are investigated in this study using fully nonlinear numerical wave tanks (NWTs) based on potential flow theory. The coupling effect of wave and current is explored, and the underlying interaction mechanisms of the hydrodynamic forces are described. The presence of a background current modifies the frequency dispersion. It produces changes of the water-surface elevation, and also has an effect on wave-induced loads. Depending on the nonlinearity, higher harmonic wave components are generated above the submerged plate. These contribute to the wave forces. It is found that the horizontal and the vertical force, hence the moment, are affected in the opposite way by the currents. The Doppler shifted effect dominates the vertical force and the moment on the plate. Whereas, the Doppler shifted effect and the generation of higher wave harmonics play opposite roles on the horizontal forces. The contribution of 2^nd order harmonics is found to be up to 30% of the linear component. The current-induced drag force, represented by the advection term ρU∂φ/∂x in the pressure equation, is found to lead to a decrease in the moment for the most range of wavelengths considered, and an increase in the moment for a small range of longer waves.     

Hydrodynamic Performance of the Arc-Shaped Bottom-Mounted Breakwater

The problem of the hydrodynamic interaction with the arc-shaped bottom-mounted breakwaters is investigated theoretically. The breakwater is assumed to be rigid, thin, impermeable and vertically located in a finite water depth. The fluid domain is divided into two sub-regions of inner and outer by an auxiliary circular interface. Linear theory is assumed and the eigenfunction expansion approach is used to determine the wave field. In order to examine the validity of the theoretical model, the analytical solutions are compared to agree well with published results with the same parameters. Numerical results including wave amplitude, surge pressure, and wave force are presented with different model parameters. The major factors including wave parameters, structure configuration, and water depth that affect the surge pressure, wave forces, and wave amplitudes are discussed and illustrated by some graphs and cloud maps.     

A new approximate analytic solution for water wave scattering by a submerged horizontal porous disk

This study gives a new approximate analytic solution for water wave scattering by a submerged horizontal porous disk in the context of the linear potential theory. The solution is based on the domain decomposition method. The velocity potentials are determined by two different approaches. One approach is to adopt decompositions for velocity potentials, and the other is to expand the vertical derivative of the velocity potential on the porous disk along the radial direction. Hence the velocity potentials are determined by the matched eigenfunction expansions. Differing from previous solutions with respect to the porous disk, the present solution needs no complex dispersion relations. Thus the new solution is easier for numerical implementation. According to numerical examples, the convergence of the present solution is satisfactory. In addition, the present predictions of the wave surface elevation and the vertical wave force on the disk agree very well with previous results by different approaches. The present solution can also be extended to other structures involving disks, such as a fish cage, a porous disk with finite thickness, and a submerged elastic disk.     

A Finite Element Solution of Wave Forces on Submerged Horizontal Circular Cylinders

In this paper, a numerical model is established for estimating the wave forces on a submerged horizontal circular cylinder. For predicting the wave motion, a set of two-dimensional Navier-Stokes equations is solved numerically with a finite element method. In order to track the moving non-linear wave surface boundary, the Navier-Stokes equations are discretized in a moving mesh system. After each computational time step, the mesh is modified according to the changed wave surface boundary. In order to stabilize the numerical procedure, a three-step finite element method is applied in the time integration. The water sloshing in a tank and wave propagation over a submerged bar are simulated for the first time to validate the present model. The computational results agree well with the analytical solution and the experimental data.Finally, the model is applied to the simulation of interaction between waves and a submerged horizontal circular cylinder.The effects of the KC number and the cylinder depth on the wave forces are studied.     

Forces and moment on a horizontal plate due to wave scattering

Wave reflection and transmission from a fixed horizontal plate have been widely studied but theoretical solutions are only available for certain limiting conditions. A general solution for this wave scattering problem is presented using the finite-element method, covering the whole range of relative depth ratio from shallow to deep water limits and submergence depth ratio from the water surface to the bed. Existing long-wave solutions for the surface plate and the submerged plate have been extended to obtain the hydrodynamic forces and overturning moment exerted on the plate. Results from the finite-element program compare well with these solutions. Variations of the reflection coefficient, wave forces and moment, with the plate width to wave length ratio, relative depth ratio and submergence depth ratio are discussed.     

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

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

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.     

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.     

Surface gravity wave interaction with elastic bottom

In the present paper, a hydroelastic model is developed to deal with surface gravity wave interaction with an elastic bed based on the small amplitude water wave theory and plate deflection in finite water depth. The elastic bottom bed is modelled as a thin elastic plate and is based on the Euler-Bernoulli beam equation. The wave characteristics in the presence of the elastic bed is analyzed in both the cases of deep and shallow water waves. Further, the linearized long wave equation is generalized to include bottom flexibility. A generalized expansion formula for the velocity potential is derived to deal with the boundary value problems associated with surface gravity waves having an elastic bed. The utility of the expansion formula is illustrated by demonstrating specific physical problems which will play significant role in the analysis of wave structure interaction problems. Behavior of the wave spectra are discussed in the case of closed basin having a free surface and an elastic bottom topography.     

MPS Method for Interaction Between Solitary Waves and Submerged Horizontal Plate

The interaction between structure and wave is a typical phenomenon in naval architecture and ocean engineering. In this paper, numerical simulation is carried out to study the interaction between a two-dimensional submerged, fixed,horizontal rigid plate and solitary wave with our in-house meshless particle CFD solver MLParticle-SJTU. First, the in-house CFD solver is verified by experimental results conducted at the State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology. During the verification, the plate is submerged under water and the solitary wave with a given amplitude is generated by a piston-type wave maker. Free surface elevation of the wave and the pressure impacting on the plate is recorded and compared with experimental data respectively. The predicted pressure and surface elevation agree well with the experimental results. Then in order to further investigate factors affecting wave-structure interaction, wave height, submerged depth and plate length are analyzed.     

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.     

浅浸没水平圆柱浮子波浪力计算研究

以海漂垃圾收集装置浮式围栏的圆柱浮子为研究对象,基于圆柱浮子的浅浸没特性,改进了Morison方程,并结合物理模型试验,对波浪与浅浸没水平圆柱浮子作用问题进行探讨。结果表明,改进的Morison方程可对水平圆柱浮子的波浪力做精确预测,并揭示了波幅、浸没深度以及周期对圆柱浮子波浪力的影响规律。对于浅浸没的圆柱浮子,所受水平波浪力随波高的增大而增大,随周期的增大而衰减到某一特定值。波浪力的正向大小分布要大于负向,而正负向大小的差异主要受浸没深度的影响。     

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.