An experimental and numerical study of heave added mass and damping of horizontally submerged and perforated rectangular plates

Forced harmonic heave motions of horizontally submerged and perforated rectangular plates are studied experimentally and numerically at both a deep and shallow submergence. The steady-state vertical forces are expressed in terms of added mass and damping coefficients. The numerical results are partly obtained by combining potential flow with linear free-surface conditions and a nonlinear viscous pressure loss condition at the mean oscillatory plate position. A domain decomposition technique is applied with a boundary element method in the inner domain and an analytical representation of the velocity potential in the outer domain. A drag term accounts for the vortex shedding at the outer plate edges. The numerically predicted Keulegan–Carpenter number dependent heave added mass and damping coefficients agree reasonably with experimental values, in particular for the deeper submergence.     

Dynamics of dual pontoon floating structure for cage aquaculture in a two-dimensional numerical wave tank

The trend of using floating structures with cage aquaculture is becoming more popular in the open sea. The purpose of this paper is to investigate the dynamic properties of a dual pontoon floating structure (DPFS) when attached to a fish net by using physical and numerical models. A two-dimensional (2-D) fully nonlinear numerical wave tank (NWT), based on the boundary element method (BEM), is developed to calculate the wave forces on the DPFS. The wave forces on a fish net system are then evaluated using a modified Morison equation. The comparisons of dynamic behaviors between numerical simulations and experimental measurements on the DPFS show good agreement. Results also display that a fish net system causes the resonant response of body motions and mooring forces to be slightly lower due to the net's damping effect. Finally, for designing the rearing space of cage aquaculture, the influences which net depth and net width have on the DPFS dynamic responses are also presented in this paper.     

Friction Drag Variation via Spanwise Transversal Surface Waves

The introduction of spanwise velocity is a promising technique to effect the near-wall turbulent flow field to influence friction drag. However, the essential physical mechanism which significantly reduces friction drag has not been understood, yet. It is the objective of this numerical study to improve the fundamental knowledge on the drag reduction mechanism. The investigation is based on spanwise traveling transversal surface waves which are applied to modify the near-wall flow field and to influence friction drag. Two actuation configurations are analyzed in detail. Compared with an unactuated flat plate boundary layer simulation the first wave setup, which represents a low frequency wave at an amplitude larger than the viscous sublayer, leads to a reduced wall-shear stress resulting in friction drag reduction of up to 9%. The second wave setup, which possesses a higher frequency and an amplitude in the range of the viscous sublayer, yields an increase of friction drag of about 8%. Unlike previous investigations which focus on excitation setups to lower friction drag, the comparison of the two wave setups in this study allows to identify the effects which on the one hand, lead to drag reduction and on the other hand, result in drag increase. That is, due to the pronounced differences the major effects determining the friction distribution are more evident. The two key features for drag reduction are the damping of the wall-normal vorticity fluctuations above the entire surface and the decrease of turbulence production. Furthermore, the effect of rearranging streamwise vorticity, which has been stated to be responsible for drag reduction, is found to occur at increasing and decreasing drag, i.e., it is not the effect that lowers the friction drag.     

Capturing the experimental behaviour of a point-absorber WEC by simplified numerical models

The paper presents a wave basin experiment of a direct-driven point-absorber wave energy converter moving in six degrees of freedom. The goal of the work is to study the dynamics and energy absorption of the wave energy converter, and to verify under which conditions numerical models restricted to heave can capture the behaviour of a point-absorber moving in six degrees of freedom. Several regular and irregular long-crested waves and different damping values of the power take-off system have been tested. We collected data in terms of power output, device motion in six degrees of freedom and wave elevation at different points of the wave basin. A single-body numerical model in the frequency domain and a two-body model in the time domain are used in the study. Motion instabilities due to parametric resonance observed during the experiments are discussed and analysis of the buoy motion in terms of the Mathieu instability is also presented. Our results show that the simplified models can reproduce the body dynamics of the studied converter as long as the transverse non-linear instabilities are not excited, which typically is the case in irregular waves. The performance of the more complex time domain model is able to reproduce both the buoy and PTO dynamics, while the simpler frequency domain model can only reproduce the PTO dynamics for specific cases. Finally, we show that the two-body dynamics of the studied wave energy converter affects the power absorption significantly, and that common assumptions in the numerical models, such as stiff mooring line or that the float moves only in heave, may lead to incorrect predictions for certain sea states.     

Suppression of transients in wave-excited response testing

When a floating body in a wave tank has low hydrodynamic damping, for example in the heave mode, very long duration transient responses can arise if it is excited from a state of rest by sinusoidal waves. Such behaviour can be undesirable when steady state response characteristics are the object of investigation in a numerical tank, because of the consequential need for very long computations. The present paper develops a method for suppressing such transient behaviour in computational models. The success of the approach is demonstrated in the context of the heaving motion of a simple buoy. A linear model of such a buoy initially at rest in a wave tank, excited by propagating sinusoidal waves, is used here for a preliminary investigation of the removal of transients. The technique is then incorporated into a fully nonlinear potential flow simulation of the buoy, and the approach is shown to be effective.     

带环形隔板的圆柱储箱内液体晃动阻尼分析

根据现有的阻尼理论，在线性势流的假设下，分别对带刚性和弹性隔板的圆柱形储液箱内液体晃动做了 有限元仿真计算. 计算结果与实验测定结果相比较发现，在晃动频率上符合得很好，但阻尼 比相差较大. 同时, 通过ALE有限元数值仿真，对储箱内无隔板和有隔板情况下的 黏性流体晃动做了比较，初步推定隔板对液体晃动的``滞阻'的一个重要表现形式为涡耗散.     

Width effects on hydrodynamics of pendulum wave energy converter

Based on two- and three-dimensional potential flow theories, the width effects on the hydrodynamics of a bottom-hinged trapezoidal pendulum wave energy converter are discussed. The two-dimensional eigenfunction expansion method is used to obtain the diffraction and radiation solutions when the converter width tends to be infinity. The trapezoidal section of the converter is approximated by a rectangular section for simplification. The nonlinear viscous damping effects are accounted for by including a drag term in the two- and three-dimensional methods. It is found that the three- dimensional results are in good agreement with the two-dimensional results when the converter width becomes larger, especially when the converter width is infinity, which shows that both of the methods are reasonable. Meantime, it is also found that the peak value of the conversion efficiency decreases as the converter width increases in short wave periods while increases when the converter width increases in long wave periods.     

浮式垂直轴风机的动力学建模、仿真与实验研究

考虑气动力和水动力的耦合研究浮式垂直轴风机系统的运动响应,将固定式垂直轴风机的气动载荷计算方法进一步推广到海上浮式垂直轴风机的气动载荷计算.考虑阻尼力、波浪力、风载荷、系泊力等,建立了浮式垂直轴风机系统的纵荡-垂荡-纵摇运动方程.考虑动态失速和浮式基础运动,基于双致动盘多流管理论,推导了风机叶片气动载荷计算公式,编制了数值计算程序.以Sandia 17 m风机为例,验证了气动载荷计算程序的正确性.最后进行了模型实验,其中模型的风机为Φ型达里厄垂直轴风机,支撑基础为桁架式Spar型浮式基础,将模型实验结果与数值计算结果进行了对比,验证了耦合计算程序.结果表明,数值计算得到的风机系统的垂荡、纵摇运动的RAO(幅值响应算子)曲线与模型实验结果吻合较好,验证了耦合程序的正确性.然而,由于数值计算与模型实验在运动自由度、阻尼、风载荷等方面存在差别,数值计算结果与模型实验结果仍有一定的差异.     

一种考虑薄壁散射效应的声学计算模型

采用薄壁边界元/FW-H理论混合方法建立了考虑薄壁声学散射效应的数值计算模型.这种声学计算模型可以预测存在薄壁如风扇机匣、蜗壳等条件下的声波的传播及散射问题.计算模型的建立主要包含噪声源的计算和声源的传播两方面:首先建立FW-H的频域方程,并采用计算流体力学方法计算流场,通过流场数据计算气动噪声源;然后采用薄壁面边界元法计算固壁对声波的散射,并计算声波在固壁散射后的声场分布.数值计算结果和实验结果及经典的叶轮机管道风扇噪声理论进行了对比,结果表明,这种计算模型与理论计算结果及实验结果吻合较好,可以准确的预测机匣壁的散射效应对声源传播的影响.     

Time-dependent hydroelastic analysis of cavitating propulsors

A 3-D potential-based boundary element method (BEM) is coupled with a 3-D finite element method (FEM) for the time-dependent hydroelastic analysis of cavitating propulsors. The BEM is applied to evaluate the moving cavity boundaries and fluctuating pressures, as well as the added mass and hydrodynamic damping matrices. The FEM is applied to analyze the dynamic blade deformations and stresses due to pressure fluctuations and centrifugal forces. The added mass and hydrodynamic damping matrices are superimposed onto the structural mass and damping matrices, respectively, to account for the effect of fluid–structure interaction. The problem is solved in the time-domain using an implicit time integration scheme. An overview of the formulation for both the BEM and FEM is presented, as well as the BEM/FEM coupling algorithm. Numerical and experiment validation studies are shown. The effects of fluid–structure interaction on the propeller performance are discussed.     

New analytical solutions to water wave radiation by vertical truncated cylinders through multi-term Galerkin method

New analytical solutions to water wave radiation by vertical truncated circular cylinders are developed based on linear potential flow theory. Two typical cylinder configurations of a surface-piercing cylinder and a submerged floating cylinder are considered. The multi-term Galerkin method is employed in the solution procedure, in which the fluid velocity on the interface between different regions is expanded into a set of basis function involving the Gegenbauer polynomials, and the cube-root singularity of fluid velocity at the side edges of the truncated cylinders is correctly modeled. The present solutions have the merits of very rapid convergence. The results with six-figure accuracy for added mass and radiation damping can be obtained using a few truncated numbers in the basis function for three motions (surge, heave and roll). The calculated results of the present solutions agree well with that by a higher-order boundary element method solution. Calculation examples are presented to investigate the influence of the motion frequency on the added mass and the radiation damping of the truncated cylinders with different geometric parameters. The present solutions can be used as a reliable benchmark for numerical solutions to water wave radiation by complicated structures.

     

Response of a dynamic system to flow-induced load

An approximate analytical method is used to study in-line vibrations of a linear system induced by oscillatory flow. The hydrodynamic drag force is accounted for by an equivalent viscous dashpot. The obtained equivalent linear system is used to determine the amplitude and the phase of the oscillatory component, and the offset component of the steady-slate periodic response of the linear system. Several parametric studies are presented and discussed in detail. Particular attention is given to the magnitude of the effective viscous damping.     

Fully nonlinear modeling of radiated waves generated by floating flared structures

The nonlinear radiated waves generated by a structure in forced motion, are simulated numerically based on the potential theory. A fully nonlinear numerical model is developed by using a higher-order boundary element method （HOBEM）. In this model, the instantaneous body position and the transient free surface are updated at each time step. A Lagrangian technique is employed as the time marching scheme on the free surface. The mesh regridding and interpolation methods are adopted to deal with the possible numerical instability. Several auxiliary functions are proposed to calculate the wave loads indirectly, instead of directly predicting the temporal derivative of the velocity potential. Numerical experiments are carried out to simulate the heave motions of a submerged sphere in infinite water depth, the heave and pitch motions of a truncated flared cylinder in finite depth. The results are verified against the published numerical results to ensure the effectiveness of the proposed model. Moreover, a series of higher harmonic waves and force components are obtained by the Fourier transformation to investigate the nonlinear effect of oscillation frequency. The difference among fully nonlinear, body-nonlinear and linear results is analyzed. It is found that the nonlinearity due to free surface and body surface has significant influences on the numerical results of the radiated waves and forces.     

Steady flow around,and drag on a circular cylinder moving at low speeds in a viscous liquid between two parallel planes 1

The present paper deals with the two-dimensional viscous flow around a circular cylinder moving at low speeds along the central plane parallel to two parallel side planes. It is semi-analytically solved by applying the collocation method for the stream function. The present results for velocity, pressure and viscous stress agree well with other numerical results. The pressure drag and the friction drag on a cylinder are obtained from the pressure distribution and the viscous stress distribution at the surface of a circular cylinder, respectively. It is shown that the pressure drag is larger than the friction drag and the ratio of the former to the latter decreases with increase of the distance between the two parallel side planes. An empirical formula for the drag which fits the present results is obtained by modifying Faxen's analytical formula. The new formula covers a wider range of the distance between the two side planes than Faxen's one.     

Water wave radiation by a sphere submerged in water with an ice-cover

Using the multipoles method, we formulate the problems of radiation (both heave and sway) of water waves by a submerged sphere in deep as well as in uniform finite depth water with an ice-cover, with the ice-cover being modelled as an elastic plate of very small thickness. In each case this leads to an infinite system of linear equations which are solved numerically by standard techniques. The added-mass and damping coefficients for a heaving and swaying sphere are obtained and depicted graphically against the wave number for various values of the radius of the submerged sphere and flexural rigidity of the ice-cover to show the effect of the presence of ice-cover on these quantities. When the flexural rigidity is taken to be zero, the numerical results for the added-mass and damping coefficient for water with a free surface are recovered.     

基于卷积非粘滞阻尼模型的超高层结构风效应分析

有限元方法中相对于对结构质量与刚度特性的描述,结构阻尼的描述仍具有较大的模糊性。随着新型建筑材料与复杂结构体系的发展,以及对计算机模拟要求的提高,阻尼作用的机理与相应阻尼模型的研究成为值得关注的问题。基于一种阻尼力与质点速度历程相关的卷积非粘滞阻尼模型,采用微分求积求解算法,对一个大型复杂超高层建筑结构的风振响应进行了分析,并与常用的比例粘滞阻尼模型进行了对比。对卷积非粘滞阻尼力模型系统的响应特征进行了分析,特别是该模型的松弛效应对结构响应的影响。另外,作为将这种新阻尼模型应用于实际工程的一次探索,本文采用微分求积算法,建立了一套可将该阻尼模型及其求解算法嵌入通用有限元软件的求解系统,可用于复杂结构的动力响应分析。     

Solution of the Navier–Stokes equations in velocity–vorticity form using a Eulerian–Lagrangian boundary element method

This paper describes the Eulerian–Lagrangian boundary element model for the solution of incompressible viscous flow problems using velocity–vorticity variables. A Eulerian–Lagrangian boundary element method (ELBEM) is proposed by the combination of the Eulerian–Lagrangian method and the boundary element method (BEM). ELBEM overcomes the limitation of the traditional BEM, which is incapable of dealing with the arbitrary velocity field in advection‐dominated flow problems. The present ELBEM model involves the solution of the vorticity transport equation for vorticity whose solenoidal vorticity components are obtained iteratively by solving velocity Poisson equations involving the velocity and vorticity components. The velocity Poisson equations are solved using a boundary integral scheme and the vorticity transport equation is solved using the ELBEM. Here the results of two‐dimensional Navier–Stokes problems with low–medium Reynolds numbers in a typical cavity flow are presented and compared with a series solution and other numerical models. The ELBEM model has been found to be feasible and satisfactory. Copyright © 2000 John Wiley & Sons, Ltd.     

一种求解瑞利波散射问题的修正边界元方法

边界元法的一大优势是用于求解半空间等无限域问题,然而对于弹性波的传播问题,传统边界元法在采用全平面或全空间格林函数时,在截断边界处仍会产生虚假的反射回波,直接影响到散射场的求解准确性。因此,本文在传统边界元法基础上提出一种修正边界元法,用于计算无限大半平面中的弹性波场问题。该方法以瑞利波形式的远端散射场代替原本因截断而舍去的部分,通过互易定理建立单位瑞利波和全平面格林函数的积分方程,求得修正系数,并代入修正边界元矩阵,计算出瑞利波的散射场。为验证本文所提方法,文中将多个算例的结果与解析解对比,并用该方法计算了不同缺陷的散射场。这些对比结果表明,本文所提修正边界元法可准确求解瑞利波散射场,为基于表面波的缺陷反演问题研究提供了有效的正演途径。     

NONLINEAR VIBRATION RESPONSE AND BIFURCATION OF CIRCULAR CYLINDRICAL SHELLS UNDER TRAVELING CONCENTRATED HARMONIC EXCITATION

The nonlinear vibration of a cantilever cylindrical shell under a concentrated harmonic excitation moving in a concentric circular path is proposed. Nonlinearities due to large-amplitude shell motion are considered, with account taken of the effect of viscous structure damping. The system is discretized by Galerkin's method. The method of averaging is developed to study the nonlinear traveling wave responses of the multi-degrees-of-freedom system. The bifurcation phenomenon of the model is investigated by means of the averaged system in detail. The results reveal the change process and nonlinear dynamic characteristics of the periodic solutions of averaged equations.     

Comparison of linear spring and nonlinear FEM methods in dynamic coupled analysis of floating structure and mooring system

This paper compares the dynamic coupled behavior of floating structure and mooring system in time domain using two numerical methods for the mooring lines such as the linear spring method and the nonlinear FEM (Finite Element Method). In the linear spring method, hydrodynamic coefficients and forces on the floating body are calculated using BEM (Boundary Element Method) and the time domain equation is derived using convolution. The coupled solution is obtained by simply adding the pre-determined spring constants of the mooring lines into the floating body equation. In FEM, the minimum energy principle is applied to formulate the nonlinear dynamic equation of the mooring system with a discrete numerical model. The ground contact model and Morison formula for drag forces are also included in the formulation. The coupled solution is obtained by iteratively solving the floating body equation and the FEM equation of the mooring system. Two example structures such as weathervane ship and semi-submersible structure are analyzed using linear spring and nonlinear FEM methods and the difference of those two methods are presented. By analyzing the cases with or without surge-pitch or sway-roll coupling stiffness of mooring lines in the linear spring method, the effect of coupling stiffness of the mooring system is also discussed.