数值波浪水槽消波器的改进与检验

数值波浪水槽技术是一种新兴的水波动力学数值模拟技术, 它能够实现水波现象的精确模拟。作为数值波浪水槽的一个重要功能, 消波技术被用于消除传入波浪在水槽末端的反射效应, 以防止反射波对有效实验区域的污染, 从而保证特定要求的水波实验的精确性。目前被广泛采用的消波技术可分为阻尼消波区和主动消波器两类。     

波浪与带窄缝方箱作用共振现象的数值模拟

By introducing a source term into the Laplace equation, a two-dimensional fully nonlinear time-domain numerical wave flume (NWF) is developed to investigate the resonance induced by the interaction bet...     

波浪与带窄缝多箱体作用共振现象的模拟研究

针对波浪与带有窄缝多箱体结构作用产生的流体共振问题,建立了基于域内源造波技术的二维非线性时域数值波浪水槽模型,其中自由水面满足完全非线性运动学和动力学边界条件,窄缝内流体引入人工阻尼来等效由于涡旋运动和流动分离引起的黏性耗散,计算域边界采用高阶边界元进行离散。通过模拟三箱体间两窄缝内相对波高变化,并与已发表的数值与实验结果对比,验证了本模型的准确性。同时通过大量的数值计算,分析了箱体数量对窄缝内水体共振频率、共振波高以及对结构反射波高和透射波高的影响。     

用ＯｐｅｎＦＯＡＭ 实现数值水池造波和消波

基于OpenFOAM求解器interDyMFoam,开发实现数值水池造波(包括推板和摇板造波)和阻尼消波。所编写的造波边界条件,可依据实验造波理论将各种造波模式植入其中,从而实现各种波的数值造波。首先进行了线性波的数值造波实验,通过结合阻尼消波段的应用,造波水池可提供稳定的线性波。还进行了瞬时极限波和有限振幅的数值造波实验,与实验数据或同类数值结果吻合很好,进一步验证了的数值造波和消波方式的可靠性。     

集成波能转换功能的方箱—挡浪板式浮式防波堤水动力特性研究

浮式防波堤与振荡浮子式波浪能转换装置集成是一种较为合理的波浪能开发利用方式,基于方箱式浮式防波堤—波浪能转换集成系统和幕帘式防波堤的研究成果,提出了一种新型方箱—垂直挡浪板式浮式防波堤—波浪能转换集成系统,建立数学模型对该集成系统的水动力特性和能量输出特性进行研究。模型基于N-S方程,采用紧致插值曲线(CIP)方法结合浸没边界法(IBM)求解。运用数值模型探究在一定波浪条件下,动力输出系统(PTO)阻尼力的大小以及挡浪板对集成系统的水动力特性和能量转换特性的影响,得到如下结论:集成系统的俘获宽度比随PTO阻尼力的增大呈现先增大后减小的趋势,在阻尼力F_(PTO)=150 N时达到最大;相对于方箱型集成系统,增设0.1 m挡浪板后可使其最大俘获宽度比η_e提高33%左右;此外,集成系统的俘获宽度比随挡浪板长度增加而增大,增长趋势逐渐变缓,在挡浪板长度S_p=0.5 m时达到最大,此时俘获宽度比η_e=0.563 1。     

完全非线性深水波的数值模拟

基于势流理论,并结合深水波质点运动从水面向下呈e指数衰减的特性,建立了完全非线性数值变深水槽模型,通过实时模拟活塞式造波机运动来产生波浪.采用时域高阶边界元法进行模拟,利用混合欧拉-拉格朗日方法和四阶Runge-Kutta方法追踪流体瞬时水面,应用镜像格林函数消除了水槽两个侧面的积分,在水槽末端布置人工阻尼层来消除反射...     

Experimental and numerical investigation of viscous effects on solitary wave propagation in a wave tank

Two-dimensional depth-averaged Boussinesq-type equations were presented with the consideration of slowly varying bathymetry and effects of bottom viscous boundary layer. These Boussinesq-type equations were written in terms of the horizontal velocity components evaluated at an arbitrary elevation in the water depth and the free surface displacement. The leading order effects of the bottom boundary layer were represented by a convolution integral in the depth-integrated continuity equation. To test the validity of the theory, a set of laboratory experiments was performed to measure the viscous damping and shoaling of a solitary wave propagating in a wave tank. The time histories of the free surface profiles were measured at several locations along the centerline of the flume. To compare these laboratory data with theoretical results, the two-dimensional Boussinesq-type equations were integrated across the wave tank, resulting in a set of one-dimensional equations, while the side-wall boundary layers were properly considered. The agreement between the experimental data and numerical results was very good. The bottom shear stress formula was also given and its impact on the sediment transport rate was discussed.     

Simulation of Fully Nonlinear 3-D Numerical Wave Tank

A fully nonlinear numerical wave tank (NWT) has been simulated by use of a three-dimensional higher order bouodary element method (HOBEM) in the time domain. Within the frame of potential flow and the adoption of simply Rankine source, the resulting boundary integral equation is repeatedly solved at each time step and the fully nonlinear free surface boundary conditions are integrated with time to update its position and boundary values. A smooth technique is also adopted in order to eliminate the possible saw-tooth numerical instabilities. The incident wave at the uptank is given as theoretical wave in this paper. The outgoing waves are absorbed inside a damping zone by spatially varying artificial damping on the free surface at the wave tank end. The numerical results show that the NWT developed by these approaches has a high accuracy and good numerical stability.     

不规则波浪数值水槽的造波和阻尼消波

利用MAC(marker and cell)法直接数值求解连续方程和N-S方程.为模拟不规则波长时段连续造波及消除波浪遇结构物后形成的二次反射,采用了源造波法.对开敞边界,采用了海绵阻尼消波和Sommerfeld条件相结合的处理方式.     

Time-domain simulation of wave–structure interaction based on multi-transmitting formula coupled with damping zone method for radiation boundary condition

Based on the Rankine source, this paper proposed a time-domain method for analyzing the three-dimensional wave–structure interaction problem in irregular wave. A stable integral form of the free-surface boundary condition (IFBC) is employed to update the velocity potential on the free surface. A multi-transmitting formula, with an artificial wave speed, is used to eliminate the wave reflection for radiation condition on the artificial boundary. An effective multi-transmitting formula, coupled with damping zone method, is further used to analyze the irregular wave diffraction at the artificial boundary. We investigate hydrodynamic forces on floating structure and compare our solution to the frequency-domain solution. It is shown that long time simulation can be done with high stability and the numerical results agree well with the solution obtained under the frequency domain. The efficiency of the proposed multi-transmitting formula and the coupled methods for radiation boundary make them promising candidates in studying the irregular water wave problem in time domain.     

Hydrodynamic analysis and shape optimization for vertical axisymmetric wave energy converters

The absorber is known to be vertical axisymmetric for a single-point wave energy converter (WEC). The shape of the wetted surface usually has a great influence on the absorber’s hydrodynamic characteristics which are closely linked with the wave power conversion ability. For complex wetted surface, the hydrodynamic coefficients have been predicted traditionally by hydrodynamic software based on the BEM. However, for a systematic study of various parameters and geometries, they are too multifarious to generate so many models and data grids. This paper examines a semi-analytical method of decomposing the complex axisymmetric boundary into several ring-shaped and stepped surfaces based on the boundary discretization method (BDM) which overcomes the previous difficulties. In such case, by using the linear wave theory based on eigenfunction expansion matching method, the expressions of velocity potential in each domain, the added mass, radiation damping and wave excitation forces of the oscillating absorbers are obtained. The good astringency of the hydrodynamic coefficients and wave forces are obtained for various geometries when the discrete number reaches a certain value. The captured wave power for a same given draught and displacement for various geometries are calculated and compared. Numerical results show that the geometrical shape has great effect on the wave conversion performance of the absorber. For absorbers with the same outer radius and draught or displacement, the cylindrical type shows fantastic wave energy conversion ability at some given frequencies, while in the random sea wave, the parabolic and conical ones have better stabilization and applicability in wave power conversion.     

Wave Propagation in A Converging Channel of Arbitrary Configuration

A numerical solution was derived to determine wave field in a converging channel bounded by rubble-mound jetties. The solution was achieved by applying boundary element method. The model was applied to analyze the effect of channel convergence, the cross-section of the jetties and their physical and damping properties on wave field in the channel. The study reveals numerous non-intuitive results specific for jetted and convergent channels. The analysis shows that wave reflection is usually low and is of secondary practical importance. Wave transmission strongly depends on the channel geometry and transmitted waves may be higher than incident waves, despite reflection and damping processes. Moreover, wave transmission depends on physical and damping properties of rubble jetties and the results show that wave transmission may increase with the increasing damping properties of jetties, which is a non-intuitive feature of wave fields in jetted channels. The analysis reveals several novel results of practical importance. It is shown that the rubble-mound jetties should be constructed from the material of high porosity, which ensures low transmission. More attention should be devoted to hydraulic properties of porous materials. It is recommended to use the material of moderate damping properties. The material of high damping properties often increases the wave transmission. It is possible, by a selection of rubble-mound material, to obtain lower transmission level for steep waves than for waves of moderate steepness. A series of laboratory experiments were conducted in the wave flume to verify the theoretical results. The comparisons show that theoretical results are in fairly good agreement with experimental data.     

Tests and Applications of An Approach to Absorbing Reflected Waves Towards Incident Boundary

If the upstream boundary conditions are prescribed based on the incident wave only, the time-dependent numerical models cannot effectively simulate the wave field when the physical or spurious reflected waves become significant. This paper describes carefully an approach to specifying the incident wave boundary conditions combined with a set sponge layer to absorb the reflected waves towards the incident boundary. Incorporated into a time-dependent numerical model, whose governing equations are the Boussinesq-type ones, the effectiveness of the approach is studied in detail. The general boundary conditions, describing the down-wave boundary conditions are also generalized to the case of random waves. The numerical model is in detail examined. The test cases include both the normal one-dimensional incident regular or random waves and the two-dimensional oblique incident regular waves. The calculated results show that the present approach is effective on damping the reflected waves towards the incident wave boundary.     

Single Mode Simulation Calculation of Oscillating Buoy Wave Energy Converter with A Slider

This paper presents an oscillating slider wave energy device which is based on a seabed anchoring and uses eagle beak as the absorber. The self-compiled program uses the boundary element theory based on the simple Green's function to solve the wave forces and hydrodynamic parameters. And the equation of motion, the oscillation of the float and the capture width ratio are obtained by the modal method. The influences of the shape of the eagle beak, the angle of the slider and the wave heading on the capture ability of the device are investigated. According to the calculation results and the wave resources in the sea area, the optimal shape of the eagle beak and external damping can be selected to maximize the wave energy capture capability.     

Hydrodynamic Analysis and Power Conversion for Point Absorber WEC with Two Degrees of Freedom Using CFD

Point absorber wave energy device with multiple degrees of freedom (DOF) is assumed to have a better absorption ability of mechanical energy from ocean waves. In this paper, a coaxial symmetric articulated point absorber wave energy converter with two degrees of freedom is presented. The mechanical equations of the oscillation buoy with power take-off mechanism (PTO) in regular waves are established. The three-dimensional numerical wave tank is built in consideration of the buoy motion based upon the CFD method. The appropriate simulation elements are selected for the buoy and wave parameters. The feasibility of the CFD method is verified through the contrast between the numerical simulation results of typical wave conditions and test results. In such case, the buoy with single DOF of heave, pitch and their coupling motion considering free (no PTO damping) and damped oscillations in regular waves are simulated by using the verified CFD method respectively. The hydrodynamic and wave energy conversion characteristics with typical wave conditions are analyzed. The numerical results show that the heave and pitch can affect each other in the buoy coupling motion, hydrodynamic loads, wave energy absorption and flow field. The total capture width ratio with two coupled DOF motion is higher than that with a single DOF motion. The wave energy conversion of a certain DOF motion may be higher than that of the single certain DOF motion even though the wave is at the resonance period. When the wave periods are high enough, the interaction between the coupled DOF motions can be neglected.     

Numerical simulation of non-linear regular and focused waves in an infinite water-depth

Inviscid three-dimensional free surface wave motions are simulated using a novel quadratic higher order boundary element model (HOBEM) based on potential theory for irrotational, incompressible fluid flow in an infinite water-depth. The free surface boundary conditions are fully non-linear. Based on the use of images, a channel Green function is developed and applied to the present model so that two lateral surfaces of an infinite-depth wave tank can be excluded from the calculation domain. In order to generate incident waves and dissipate outgoing waves, a non-reflective wave generator, composed of a series of vertically aligned point sources in the computational domain, is used in conjunction with upstream and downstream damping layers. Numerical experiments are carried out, with linear and fully non-linear, regular and focused waves. It can be seen from the results that the present approach is effective in generating a specified wave profile in an infinite water-depth without reflection at the open boundaries, and fully non-linear numerical simulations compare well with theoretical solutions. The present numerical technique is aimed at efficient modelling of the non-linear wave interactions with ocean structures in deep water.     

The hydrodynamic coefficients for an oscillating rectangular structure on a free surface with sidewall

Based on a two dimensional linear water wave theory, the boundary element method (BEM) is developed and applied to study the heave and the sway problem of a floating rectangular structure in water to finite depth with one side of the boundary is a vertical sidewall and the other boundary is an open boundary. Numerical results for the added mass and radiation damping coefficients are presented. These coefficients are not only depend on the submergence and the width of the structure, but also depend on the clearance between structure and sidewall. Negative added mass and sharp peaks in the damping and added mass coefficients have been found when the clearance with a value close to integral times of half wave length of wave generated by oscillation structure. The important effect of the clearance on the added mass and radiation damping coefficients are discussed in detail. An analytical solution method is also presented. The BEM solution is compared with the analytical solution, and the comparison shows good agreement.     

Underwater spreading and surface drifting of oil spilled from a submarine pipeline under the combined action of wave and current

Tremendous economic loss and environmental damages are caused by oil-spilling accidents in sea. Accurate prediction of the underwater spreading and surface drifting of oil spills is important for the emergency response. In the present study, numerical investigation on the underwater spread and surface drift of oil spilled from a submarine pipeline under the combined action of wave and current was carried out to examine the effects of physical ocean environment, leaking flux and spilled oil density and viscosity. Reynolds-Averaged-Navier-Stokes (RANS) equations, realizable k-ε turbulence model and volume of fluid (VOF) model are employed to describe the multiphase flow, and velocity-boundary wave-making technique combined with the sponge layer damping absorber technique realizes the numerical wave flume. Oil spill experiments were conducted to validate the numerical model. The calculation results indicate that compared with the environmental conditions of still water, only current and only wave, a larger scope of underwater spreading and relatively slower rising rate and relatively faster drifting rate of oil droplets are observed under the combined action of wave and current. The leaking flux affects the floating time and dispersion concentration, while the ocean environment affects the horizontal migration and surface drifting. Under the specific conditions of present work, oil density has obvious effect on the underwater spread but limited effect on the surface drifting, while oil viscosity has little effect on both the two processes.     

Laboratory and numerical studies of wave damping by emergent and near-emergent wetland vegetation

Wetlands protect mainland areas from erosion and damage by damping waves. Yet, this critical role of wetland is not fully understood at present, and a means for reliably determining wave damping by vegetation in engineering practice is not yet available. Laboratory experiments were conducted to measure wave attenuation resulting from synthetic emergent and nearly emergent wetland vegetation under a range of wave conditions and plant stem densities. The laboratory data were analyzed using linear wave theory to quantify bulk drag coefficients and with a nonlinear Boussinesq model to determine numerical friction factors to better represent wetland vegetation in engineering analysis.     

Dynamics of elastically moored floating objects

The two-dimensional problem of wave transformation by, and motions of, moored floating objects is solved numerically as a boundary value problem by direct use of Green's identity formula for a potential function. The cross-sectional shape of the floating object, the bottom configuration and the mooring arrangements may be all arbitrary. For a given incident wave, the three modes of body motion, the wave system and mooring forces are all solved at the same time. A laboratory experiment is conducted to verify the theory. Generally good agreements between the theory and experiments are obtained as long as the viscous damping due to flow separation is small. A numerical experiment indicates that a conventional sluck mooring is to worsen the wave attenuation by a floating breakwater and that a properly arranged elastic mooring can considerably improve the wave attenuation by a floating breakwater.