Wave characteristics past a flexible fishnet

The scattering of surface waves by a flexible fishnet is studied analytically. The fishnet is modelled as a porous flexible barrier displaced solely by hydrodynamic force like a catenary. The objective is to investigate how a flexible permeable barrier affects the passing waves in the way they are transmitted and reflected, as observed by the fact that the water inside a fishfarm surrounded by fishnets is significantly calmer than that outside. The boundary value problems are solved by defining the reflection coefficient in terms of velocity potential and then the full solutions are obtained by suitable application of the eigenfunction expansion method and the least squares approximation method. The variations of the reflection coefficient, hydrodynamic pressure, barrier deformation and surface wave elevation are determined with respect to the barrier length, porosity and stiffness. It is observed that as the fishnet gets more flexible, its deformation increases and the reflection coefficient decreases, whereas as the fishnet gets more porous, more water can pass through it and thus the reflection coefficient, barrier deformation and the hydrodynamic force are reduced. The flexibility of the barrier behaves like its porosity by allowing more wave energy to act on it through its deformation and hence reduce the reflection and hydrodynamic force of the incident waves acting on the barrier.     

柔性仿金枪鱼月牙形尾鳍水动力分析

为研究鱼类高速游动机理,以金枪鱼月牙形尾鳍为对象,采用面元法计算分析了其非定常水动力性能。假设尾鳍在横移和摇摆的同时,以某一匀速向前运动,并假设其在弦向和展向以某一给定规律发生变形,以模仿柔性变形。探讨了前进速度、横移和摇摆的幅度、频率及其相位差对柔性尾鳍推进性能的影响,并与做相应运动的刚性尾鳍进行了对比。     

基于桁架模型和多孔介质模型的柔性网衣和周围流场单向耦合方法研究

柔性网衣作为渔业养殖的重要组成部分,其水动力的计算分析对网箱的安全使用至关重要。基于桁架模型和多孔介质模型对柔性网衣和流场单向耦合计算的数值模拟方法进行了研究。通过code_aster软件基于桁架模型计算出网衣变形后节点坐标和网衣受力,然后将节点坐标导入OpenFOAM识别出网衣变形后的多孔介质区域,并基于多孔介质模型模拟水流作用下网衣周围的流场和计算网衣的受力。为了验证模型的正确性,将不同工况对圆形网衣的数值模拟结果与物理试验数据进行比较,包括code_aster中计算的网衣变形形状、OpenFOAM模拟的网衣周围流场以及在code_aster和OpenFOAM中计算的网衣受力。计算结果表明,圆形网衣在单纯来流下变形明显,与试验结果保持一致,证明了网衣变形和水动力单向耦合方法的准确性。     

A practical surface panel method to predict velocity distribution around a three-dimensional hydrofoil including boundary layer effects

A practical, low order and potential-based surface panel method is presented to predict the flow around a three-dimensional rectangular foil section including the effect of boundary layer. The method is based on a boundary-integral formulation, known as the “Morino formulation” and the boundary layer effect is taken into account through a complementary thin boundary layer model. The numerical approach used in the method presents a strongly convergent solution based on the iterative wake roll-up and contraction model including the boundary layer effect. The method is applied to a three-dimensional foil section for which the velocity distribution around the foil was measured using a 2D Laser Doppler Velocimetry system in a large cavitation tunnel. Comparison of the predicted velocity distributions both inside and outside of the boundary layer of the foil as well as the boundary layer shapes obtained from the numerical model show fairly good correlation with the measurements, indicating the robustness and practical worthiness of the proposed method.     

Numerical simulation of the motion of an ice keel in a stratified flow

The results of the three-dimensional numerical simulation for the study of the stratification effect and wave processes associated with it on the drag of the underwater part of the hummocked ice are considered. The numerical model is based on the sampling of equations on a rectangular grid using the immersed boundary method that makes it possible to explicitly describe the interaction of moving ice with a stratified flow. The dependence of the drag force on the Froude number was established based on these calculations. This dependence has expressed points of maximum and minimum. The form of this dependence is common for the considered models of ice keels. The obtained estimations of drag force consistent with the known results of laboratory experiments show the need for the construction of parametrizations of the drag coefficient on the ice–ocean boundary, taking into account wave effects.     

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

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

Hydrodynamic Characteristics of Square Heaving Plates with Opening Under Forced Oscillation

The existence of the heaving plates can improve the heaving motion performance of an offshore structure significantly by providing both extra added mass and damping. In the current research, numerical investigation is carried out on the hydrodynamic characteristics of both isolated square heaving plate and double square heaving plates with opening by an immersed boundary-lattice Boltzmann method. The effects on hydrodynamic performance of plates due to Keulegan-Carpenter(KC) number, frequency number, opening ratio, opening distribution and spacing of plates are examined. It is found that the heaving plates with optimized opening ratio can provide additional damping compared with the plates without opening. Better hydrodynamic characteristics of double plates can be obtained with the increase of plate spacing.     

波浪与外壁透空双方形沉箱相互作用

在线性波浪理论下,利用复合边界元素法(composite BEM)数值解析在等水深、规则波浪入射二种外壁透空双方形沉箱的无因次波力及Kd绕射分布图,并分别与其他研究者所作双圆筒内、外圆柱无因次波力及绕射分布图的计算结果进行比较,都说明本模式的合理性和可行性。在考虑不同透水参数下,分别计算波浪作用在内、外结构物的波力及沉箱四周绕射系数大小分布图。计算结果显示:波浪作用于外壁透空全透水双方形沉箱,外方柱无因次最大波力值会随着透水参数增加而降低;内方柱无因次最大波力值则随着透水参数增加而增加。本研究结果可供设计外海透水方形沉箱结构交互作用参考。     

细长输流管内外流耦合振动特性研究

管内流动会影响输流管的振动响应,目前关于输流弹性管涡激振动方面的研究较少。基于计算流体力学(CFD)方法,开展内外流对细长输流弹性管振动特性影响的研究。首先在不考虑内流的情况下将弹性管涡激振动数值预报结果与模型试验数据进行对比,验证了数值方法的可靠性。再者考虑内外流耦合作用情况下,对不同内流流速下细长输流弹性管振动位移时—空分布、顺流向最大平均偏移、振动轨迹、内部横向涡的形成与分布等进行了对比分析。结果发现,与外流流速相比,内流流速的增加虽然难以改变弹性管的主振模态,但对沿管体的振动强度影响显著。顺流向最大偏移处管体运动轨迹发生明显的变形和跳跃。在剪切外流和均匀内流对弹性管的联合作用下,沿管跨方向模态间能量转换频繁,伴随着间歇性出现或消失的沿弹性管传播的行波组分,这主要归因于复杂的双重流固耦合系统(外流—管体,内流—管体)。在内流以附加质量力、离心力和科氏力形式的激励下,弹性管内二次流现象明显。在振动过程中,内部横向涡沿管壁生成、脱落并逐渐散布于整个横截面。     

Numerical simulation of the effects of fish behavior on flow dynamics around net cage

The computational fluid dynamics study is performed to analyze the impact of the cultured fish on the flow field through net cage and the deformation of net cage. The shear stress turbulent k-omega model is applied to simulate the flow field through the net cage, and the large deformation nonlinear structure model is adopted to conduct the structural analysis of the flexible net cage. To validate the net-fluid interaction model of the net cage in current, a series of physical model tests are conducted, which indicate that the numerical model can accurately simulate the flow field around the net cage and the deformation of the net cage. A fish model is used to simulate the effect of fish behavior on the flow pattern around the net cage and the deformation of the net cage. In addition, the flow fields around the net cage in current are investigated considering different fish group structures, fish swimming speeds, fish distributions and fish stocking densities. The results indicate that the circular movement of fish in the still water leads to a low pressure zone at the center of net cage, which causes a strong vertical flow along the center line of the net cage. The drag force on the net cage is significantly decreased with the increasing fish stocking density, but the most severe deformation of net cage occurred in the case of medium fish stocking density.     

An Iteration Method for Computation of Flexible Fender Piles

—In the designing of the flexible fender pile,M-method is usually adopted.But M-method isusable only in the case that the value of the force(moment M and shear force Q)exerted on the pile abovethe soil is known.In practical cases,only the percussive energy of the ship can be found,while the exactvalue of the percussive force cannot be computed directly.In this paper,a method is introduced which canbe used to calculate the inner force of flexible fender piles by a computer in these cases.First,according tothe known docking dynamic energy of a ship,the absorbed energy can be found from the deformation ofthe pile and the rubber fender respectively,at the same time the value of the force acting on the fender pilesand the reaction of fenders can be calculated.Then the inner force and displacement of the fender piles canalso be calculated further by M-method.The whole computing process is realized by a computer programand the results obtained are reasonable.In this paper,the computation,the iteration method a     

直驱式波浪能转换系统的研究和性能优化

A direct-drive wave energy conversion system based on a three-phase permanent magnet tubular linear generator （PMTLG） and a heaving buoy is proposed to convert wave energy into electrical energy. Sufficient experimental methods are adopted to compare the computer simulations, the validity of which is verified by the experiment results from a wave tank laboratory. In the experiment, the motion curves of heaving buoy are with small fluctuations, mainly caused by the PMTLG＇s detent force. For the reduction of these small fluctuations and a maximum operational efficiency of the direct-drive wave energy conversion system, the PMTLG＇s detent force minimization technique and the heaving buoy optimization will be discussed. It is discovered that the operational efficiency of the direct-drive wave energy conversion system increases dramatically after optimization. The experiment and optimization results will provide useful reference for the future research on ocean wave energy conversion system.     

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

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

Numerical simulation of oblique water entry of an asymmetrical wedge

The hydrodynamic problem of a two-dimensional asymmetrical wedge entering calm water obliquely at constant speed is analyzed based on the velocity potential theory. The gravity effect on the flow is ignored based on the assumption that the ratio of the entry speed to the acceleration due to gravity is much larger than the time scale of interest. The problem of this similarity flow is solved by a boundary element method together with an analytical solution for the jet based on the shallow water approximation. Various results are provided for the wave elevation, pressure distribution and force at different deadrise angles and at different oblique entry. The effects of asymmetry and horizontal speed on these results are investigated.     

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

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

Hydrodynamic characteristics of semi-immersed breakwater with an attached porous plate

In the present study, wave interaction with a fixed, partially immersed breakwater of box type with a plate attached (impermeable-permeable) at the front part of the structure is investigated numerically and experimentally. The large scale laboratory experiments on the interaction of regular waves with the special breakwater were conducted in the wave flume of Laboratori d’Enginyeria Marνtima (LIM) at Universitat Politecnica de Catalunya (UPC) in Barcelona. Experimental results are compared with numerical results obtained with the use of the Cornell breaking Wave and Structures (COBRAS) wave model. The effects of an impermeable as well as a permeable plate attached to the bottom of the breakwater on its hydrodynamic characteristics (wave transmission, reflection, dissipation, velocity and turbulence kinetic energy) are investigated. Computed velocities and turbulence kinetic energy in the vicinity of the structure indicate the effects of the breakwater with the attached (impermeable/permeable) plate on the flow pattern and the turbulence structure. The attached impermeable plate at the front part of the breakwater enhances significantly the efficiency of the structure in attenuating the incident waves. The permeable plate reduces the efficiency of the structure since wave energy is transmitted through the porous body of the plate. Based on the hydrodynamic characteristics it is inferred that the breakwater with an impermeable plate attached to its bottom is more efficient. The comparison of horizontal and vertical forces acting on the breakwater for all cases examined reveals that plate porosity influences slightly vertical force and severely horizontal force acting on the structure, reducing maximum values in both cases.     

Time stepping solutions of the two-dimensional nonlinear wave radiation problem

The two-dimensional nonlinear time domain free surface flow problem is analysed using potential flow theory. The problem is solved by a time marching method. At each time step two numerical approaches are used. One is based on the boundary element method in the complex plane. The complex potential is assumed to vary linearly within each element and the solution is obtained by imposing the boundary conditions at the nodes of the elements. The other approach is based on the finite element formulation. Triangular elements and linear shape functions are used. The solution is obtained by the Galerkin method. Numerical results are obtained for the wave elevation generated by a vertical wave maker. Results are also provided for a circular cylinder oscillating below the free surface. For these cases the finite element method is found to provide substantially more efficient computations than the boundary element method using equivalent discretizations.     

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.     

An Approximate Method for Calculation of Fluid Force and Response of A Circular Cylinder at Lock-in

In this paper, equations calculating lift force of a rigid circular cylinder at lock-in in uniform flow are deduced in detail. Besides, equations calculating the lift force on a long flexible circular cylinder at lock-in are deduced based on mode analysis of a multi-degree freedom system. The simplified forms of these equations are also given. Furthemore, an approximate method to predict the forces and response of rigid circular cylinders and long flexible circular cylinders at lock-in is introduced in the case of low mass-damping ratio. A method to eliminate one deficiency of these equations is introduced. Comparison with experimental results shows the effectiveness of this approximate method.     

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.