聚焦波浪与直立圆柱作用的数值模拟

基于Fluent平台建立了可用于波浪传播模拟的数值计算水槽,利用不同波浪参数的线性规则波模拟分析了网格尺度L/Δx对产生波浪质量的影响,确定数值水槽的网格划分尺度,并验证了阻尼消波的有效性。在此基础上,用模拟频谱的方法对聚焦波浪的传播进行了数值模拟,通过数值模拟与物理试验不同位置处的聚焦波波面和聚焦点下水质点速度时间过程线对比,验证了模拟频谱方法数值计算的有效性。进一步对聚焦波浪与直立圆柱的作用进行了数值模拟,数值模拟聚焦波面过程、聚焦波浪水质点速度、圆柱周围波浪波动过程、圆柱上所受的波压力等,并与试验结果进行了对比,数值计算结果与试验结果吻合很好,验证了所建立的数值水槽的有效性。对聚焦波波面和压力特征值进行了误差分析,并分析了误差的产生原因,为进一步研究聚焦波浪与工程结构物的作用奠定基础。     

THE RADIATION AND DIFFRACTION OF WATER WAVES BY A BOTTOM-MOUNTED CIRCULAR CYLINDER IN A TWO-LAYER FLUID

An analytical method was proposed to analyze the radiation and diffraction of water waves by a bottom-mounted circular cylinder in a two-layer fluid. Analytical expressions for added mass and damping coefficients, as well as the wave excitation forces of the circular cylinder were obtained by an eigenfunction expansion method. The hydrodynamic forces on the bottom-mounted circular cylinder in a two-layer fluid include not only the added mass and damping coefficients, but also the wave forces of the surface and internal-wave modes. This is different from the case of a homogenous fluid. Some examples were given, showing that density stratification can have a relative large effect on these hydrodynamic forces over a wide range of frequencies.     

SECOND-ORDER DRIFT FORCES AND MOMENTS BETWEEN TWO SUBMERGED BODIES IN WAVES ON 3D METHOD

Two near field methods, namely the integral method and differential method, were presented for giving second-order mean drift forces and moments between two fixed submerged bodies in regular waves. For the integral method, with a series of mathematical manipulations, second order drift forces and moments could be easily expressed by distributed sources which could be calculated by source distribution techniques with the assumption that the amplitude of ship motions are small on the basis of the linear 3D frequency theory. For the differential method, drift forces and moments could be expressed by the derivative of velocity potential with respect to space coordinate. Because two bodies would behaveas a single body while the clearance is very large, the numerical results of one sphere in such case were given and compared with analytical results of a single sphere which does not involve the effect of free surface. When submerged depth becomes enough large, a good agreement can be reached. Then the integral method was used to predict the second order drift forces and moments of two submerged spheres and spheroids with a small lateral separation distance in waves compared with the numerical results obtained by the differential method and they agree well. By comparison, it indicates the interaction effects between two submerged bodies have a profound influence on the drift forces and moments. In this paper, the forward speed effect on submerged spheres was also considered.     

REFLECTION OF OBLIQUE INCIDENT WAVES BY PERFORATED CAISSONS WITH TRAVERSE WALL

The interaction of oblique incident waves with infinite number of perforated caissons is investigated. The fluid domain is divided into infinite sub-domains by the caissons, and eigen-function expansion is applied to expand velocity potentials in each domain. A phase relation is introduced for wave oscillation in each caisson, and the structure geometry is considered in constructing the models of reflection waves. The reflected waves with the present analysis include all of the waves traveling in different directions when incident wave period is short. Numerical examinations show that velocities at the inner and outer sides of the front walls of caissons are close to each other, and reflection coefficients satisfy the energy conservation relation very well when porous effect parameter is infinite. Numerical results show that the reflection coefficients of oblique incident waves are smaller for shorter caissons at low frequency, and decrease with the increase of wave incident angle.     

有限深两层流体中水波与水平圆柱潜体的相互作用

研究了有限深两层流体中表面波模态和内波模态的水波与水平圆柱潜体的相互作用问题。在线性势流理论框架内,建立了两种模态入射波作用下,水平圆柱体绕射势的多极展开求解理论。基于所建立的求解模型,对作用在圆柱体上的波浪力进行了数值计算分析,并且与均匀流中情况进行了比较。结果表明了,在入射波的某个频率范围内,流体的分层效应对圆柱体波浪力的影响是不可忽略的。同时,对流域内的绕射速度矢量场进行了数值分析。     

用时域法求解二维二阶非线性水波的绕射问题

本文用二阶理论在时期围内研究非线性水波的绕射问题。一、二阶问题分别满足各自的自由表面条件和物面条件。用一阶Ａｄａｍｓ格式并配合时步处理来计算自由表面上的波高和速度势，而流场内的速度势通过求解有限元方法得到。采用人工阻尼区来吸收反射波，对水平半圆柱在水面的绕射问题进行了计算，得到了自由表面波高时间历程和圆柱所受到的波浪力，计算结果和有关文献进行了比较，其符合程度是令人满意的。     

OSCILLATION AMPLITUDE OF THE FREE END OF A FLEXIBLY MOUNTED RIGID CYLINDER IN TRIANGLE ARRANGEMENT

In this paper,the oscillation amplitude of the free end of a group of threecylinders are studied.The three cylinders are erected at the vertices of an equilateral trian-gle,and three arrangements are made in this study.One of the three cylinders is served asa working cylinder which is mounted by a flexible cantilever and two accelerometers areattached on the free end in two directions.The results of the three-cylinder arrangementare compared with those of two cylinders in tandem/side-by-side arrangement in the samecondition.The experimental results show that the rear cylinder of a group of three cylin-ders is affected more seriously than that of the corresponding two cylinders in the contrastarrangement.The main direction of oscillation is changed.The transversal oscillation isstronger than the longitudinal oscillation in the three-cylinder arrangement.On the con-trary,the transversal oscillation is much weaker than the longitudinal oscillation in two-cylinder arrangement.The largest transversal oscillat     

多层水平板衰减波浪的影响因素分析

基于线性波浪理论,进行多层水平板对波浪的反射和透射作用研究。通过分解速度势为对称势和反对称势,结合特征函数展开法,在考虑水平板厚度的条件下由匹配条件确定速度势展开式中的待定系数,进而得到结构对波浪的反射系数和透射系数。所得结果与现有相关的理论解、数值解以及试验结果吻合较好。对相对水深、相对宽度、相对厚度以及水平板层数等对结构衰减波浪的影响进行了深入探讨和分析,给出了反射系数和透射系数随上述参数的变化规律。     

基于OpenFOAM的挡板透空式防波堤水动力特性数值分析

基于有限体积法理论,应用开源计算流体力学软件OpenFOAM求解黏性不可压缩流体的连续性方程和动量方程,建立了二维数值波浪水槽,并通过规则波作用下挡板式透空结构的物理模型试验数据(波面、压力和流速)对数值波浪水槽进行了验证;计算并比较分析了不同入射波浪条件下挡板相对入水深度对透射系数和挡板底部流速的影响.结果表明,当挡...     

Finite element analysis of second order wave radiation by a group of cylinders in the time domain

A finite element based numerical method is employed to analyze the wave radiation by multiple or a group of cylinders in the time domain. The nonlinear free surface and body surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a Finite Element Method (FEM). The matrix equation of the FEM is solved through iteration and the initial solution is obtained from the result at the previous time step. The three-dimensional (3-D) mesh required is generated based on a two-dimensional (2-D) hybrid mesh on a horizontal plane and its extension in the vertical direction. The hybrid mesh is generated by combining an unstructured grid away from cylinders and two structured grids near the cylinder and the artificial boundary. The fluid velocity on the free surface and the cylinder surface are calculated by using a differential method. Results for various configurations including the cases of two cylinders and four cylinders and a group of eighteen cylinders are obtained to show the joint influences of cylinders on the first- and second-order waves and forces, including the effects of spacing ratios and wave frequency on the second order waves and the mean force, in particular.     

水轮机及其调速系统对电力系统低频振荡的影响分析

基于水轮机及其调速系统的典型模型,推导了转子的机械阻尼力矩系数;采用非线性规划的方法,对典型参数范围内的机械阻尼力矩系数进行了极值分析。分析发现,水轮机及其调速系统提供的负阻尼转矩,可能导致电力系统呈现弱阻尼甚至负阻尼状态,进而可能发生强迫共振和负阻尼两种形式的低频振荡。该分析结果在频域和时域仿真中都得到了验证,对低频振荡的分析及抑制有现实意义。     

受控振动柱体对波浪影响的数值模拟研究

以不可压缩黏性流体的Navier-Stokes方程为控制方程,采用RNG k-ε湍流模型,利用VOF方法追踪运动的自由表面的方法,建立了三维数值波浪与振动柱体相互作用的模型.基于经典绕射理论验证了波浪绕固定柱体传播的数值模型的可靠性,再采用动边界技术赋予柱体振动规律,研究受控振动柱体对波浪传播的影响,并通过控制柱体初始振动相位变化,系统地给出了在柱体振动的不同初始相位下柱体前后波浪变化的规律.     

电力系统稳定器在喜河水电站的应用研究

针对喜河水电站水轮发电机励磁控制系统中的励磁调节器产生负阻尼转矩可能导致电力系统低频振荡的问题,采用电力系统稳定器（PSS）以抑制低频振荡。建立了喜河水电站3号水轮机组的PSS模型,通过模型参数测试结果与仿真结果的比较,验证了参数的准确性,并进行了现场试验研究。结果表明,该机组PSS能够有效地抑制系统低频振荡,对本机振荡有良好的正阻尼作用。     

Vortex shedding in the flow around two side-by-side circular cylinders of different diameters

In this paper,the 3-D turbulent flow around two side-by-side circular cylinders of different diameters,at sub-critical Reynolds number(Re=3 900),is numerically simulated by the large eddy simulation(LES).The spacing ratios(T/D)between the two cylinders are considered in four cases(T/D=1.2,1.5,1.8 and 2.7)to study the vortex shedding and turbulent properties in the flow field.The main results are focused on the drag and lift coefficients,the vortex shedding frequency,the coherent structure,and the scale properties.It is shown that when T/D is equal to 1.2,the vortex shedding of the main cylinder is strongly suppressed by the small cylinder,the drag and lift coefficients of the main cylinder are smaller than those in other three cases.While T/D is equal to 1.5,the vortex shedding of the main cylinder can be improved,the drag and lift coefficients of the main cylinder are larger than those in other three cases.The empirical mode decomposition(EMD)method is applied to decompose the velocity signals traced by the LES.It is shown that there is a linear relationship between the mean period and the mode in the semi-log coordinates.The vortex shedding period of the main cylinder is consistent with the period of the restructured coherent structures quantitatively.     

Numerical simulation of rolling for 3-D ship with forward speed and nonlinear damping analysis

Strongly nonlinear characteristics of ship roll owing to viscous effect can be usually observed. To describe the nonlinear roll behavior, the CFD method has been frequently employed with obvious advantages compared with the traditional semi-empirical formula method in estimating the roll damping. Numerical simulations of free decay and forced rolling at various forward speeds and amplitudes for a 3-D ship hull are conducted in the present research to predict ship roll damping, in which a RANS solver is employed and a dynamic mesh technique is adopted and discussed in detail. Numerical results, including nonlinear flow characters around ships, rolling decay curves and damping coefficients, show that they are all in good agreement with available experimental data. The linear and nonlinear damping coefficients are estimated and analyzed by fitting with exponential functions for various rolling amplitudes, frequencies and speeds in the free decay simulations, and the damping coefficients are obtained by a polynomial fitting in the forced roll simulations. It is indicated that the damping coefficients increase with increasing rolling angle amplitude and velocity. It is also emphasized that the effect of forward speed is significant to roll damping and the nonlinear damping decreases with increasing velocity.     

Numerical simulation and comparison of potential flow and viscous fluid models in near trapping of narrow gaps

Numerical simulation results of fluid oscillation in narrow gaps subjected to incident water waves are presented in this paper. Both viscous fluid model and potential flow model with artificial viscous damping are employed to investigate wave motions in a single narrow gap separated by twin bodies and double narrow gaps formed by three identical bodies. Variation of wave heights in the narrow gaps with incident wave frequency is compared with experimental data available in literatures. The numerical results demonstrate that both the viscous fluid model and the potential flow model are able to predict the resonant frequency reasonably well. However the conventional potential flow model (without introducing the artificial viscous damping) significantly over-predicts the wave height in narrow gaps around the resonant frequency. The potential flow model with artificial damping predicts resonant wave heights well with a properly tuned (or calibrated) damping coefficient. The horizontal wave force on the individual body is also examined in this work. It is found that a properly calibrated damping coefficient is also very important for estimating the horizontal wave forces on the structures. A damping coefficient of 0.4 in the potential flow model is found to result in satisfactory predictions of wave heights in the narrow gaps and the horizontal hydrodynamic force for the cases considered in this work.

     

寒区流冰撞击输水隧洞的数值模拟分析

位于西北寒区的长距离输水工程在流冰期输水过程中,不同流速和不同碰撞角度的流冰对输水隧洞衬砌会形成不同程度的撞击作用,导致输水隧洞衬砌产生破损,因此为了研究该地区流冰对输水隧洞的撞击作用,采用了理论分析结合数值模拟的方法。研究结果表明:当流冰速度一定时,随着流冰与隧洞衬砌碰撞角度的增大,其对隧洞衬砌的影响呈现出近似的非线性(多项式)关系;当碰撞角度一定时,随着流冰速度的增大,其对隧洞衬砌的影响呈现出近似的线性关系;综合不同组合工况还可以发现流冰速度对隧洞衬砌的影响作用要强于碰撞角度对其的影响作用,且碰撞角度对隧洞衬砌所受的撞击力及隧洞碰撞区位移的影响较为不明显。     

三体船在斜浪规则波中运动响应预报方法研究

对于三体船在波浪中的水动力预报,航速效应和主体与片体间水动力干扰效应是水动力载荷合理预报的关键之处,需合理考虑.二维半势流理论可较为准确反映上述效应,已在三体船顶浪中垂向运动响应理论预报上获得了较为满意的结果.在此基础上,该文探讨了该方法在三体船斜浪中运动特别是横摇预报上的应用.为较准确获得三体船横摇运动时粘性阻尼贡献,采用三体船模静水中自由横摇衰减试验获得了自由横摇运动衰减曲线,并经能量法处理获得了船舶横摇阻尼系数.然后,结合二维半理论,获得了三体船运动数值解.数值解与斜浪规则波中运动响应模型试验结果符合较好,反映本文所采用的数值方法可较好地预测三体船斜浪中水动力载荷和运动,可用于三体船在波浪中耐波性评估.     

高速冲头有界入水初始阶段的数值模拟

冲头以高速撞入比它最大外径稍微大一点的盛水容器，在特殊工程中由撞水阻力起到制动作用。本文采用实用的方法处理自由面条件及狭缝两侧的物面条件计算了高速冲头有界入水初始阶段的运动过程。分别对不同质量、不同初始速度、以及不同外形的冲头入水过程进行了计算，得到了冲头阻力系数、以及速度和加速度的变化规律，分析了各因素对高速冲头有界入水效果的影响。数值模拟结果表明：这样处理边界条件对解决高速冲头有界入水问题是有效的     

Stability of a tumblehome hull under the dead ship condition

Some methods for direct stability assessment under the dead ship condition were currently developed by the international maritime organization(IMO)under the Second Generation Intact Stability Criteria.Model tests and simulations are carried out to validate the numerical methods used in assessing the stability under the dead ship condition.This is done in three stages.Firstly,the uncoupled roll mathematical model(1 DOF)is adopted to calculate the roll motion based on the irregular beam waves and the steady wind.Secondly,a drift free experiment is conducted to measure the roll motion under irregular beam waves with zero speed,and then two restrained experiments with counter weights and four springs are performed under the same condition.Finally,the effects of the drift and sway motions on stability under the dead ship condition are then verified by experimental results,and the results of the numerical methods are compared to the results of the model experiments.It is concluded that more accurate numerical methods could be developed for assessing the direct stability under the dead ship condition.