一种风力机气动计算的全自由涡尾迹模型

采用全自由方式建立风力机尾流场的涡尾迹模型,引入“虚拟周期”的概念,并发展一种自适应松弛因子方法,从而改善了自由尾迹迭代的稳定性,提高了迭代收敛速度。利用建立的自由涡尾迹模型,计算了风力机叶片的尾流场结构、气动性能及叶片载荷,并与实验结果进行了对比分析。结果表明,尖速比越大,自适应松弛因子方法对缩小模型计算时间越有效;全自由涡尾迹模型能准确给出风力机尾流场的结构,包括尾迹的扩张以及叶尖涡和叶根涡的产生、发展和耗散的过程,风轮扭矩与实验数据吻合;叶片载荷分布的计算结果在低风速下与实验值基本一致,但是在大风速下差别较大,说明需要一个准确的失速模型。     

尾流激励的叶片气动力降阶模型

气动力降阶模型是研究叶片气动弹性振动快速高效的新方法。现有气动力降阶模型的研究主要集中在叶片颤振方面,没有涉及更为常见的上游尾流激励的叶片振动问题。本文提出基于Volterra级数的尾流激励叶片气动力降阶模型,为尾流激励下叶片振动和动静叶干涉振动研究提供了新的思路。采用行波法简化尾流的参数个数,用阶跃信号法识别降阶模型的核函数。二维叶片的算例结果表明,本文方法可以较准确地描述尾流激励引起的叶片气动力振荡,而且计算效率极高。     

尾流激励的叶片气动力降阶模型

气动力降阶模型是研究叶片气动弹性振动快速高效的新方法。现有气动力降阶模型的研究主要集中在叶片颤振方面,没有涉及更为常见的上游尾流激励的叶片振动问题。本文提出基于Volterra级数的尾流激励叶片气动力降阶模型,为尾流激励下叶片振动和动静叶干涉振动研究提供了新的思路。采用行波法简化尾流的参数个数,用阶跃信号法识别降阶模型的核函数。二维叶片的算例结果表明,本文方法可以较准确地描述尾流激励引起的叶片气动力振荡,而且计算效率极高。     

球栅的Oseen流动及阻力的尾流效应

本文通过大小相同且等距放置的球栅的轴对称Oseen流动研究了阻力的尾流效应。在每个球的球心放置Oseen流子得到了问题的级数解。截断无穷级数并采用配置法解线性代数方程组求出了球栅Oseen流动的近似解及每个球所遭受的阻力。 在不同的球的个数,不同的球的间距以及不同的雷诺数下计算了各个圆球的阻力系数,发现除球栅的遮蔽效应和端缘效应外还存在着尾流效应。研究了上述参数对这些效应的影响并与Stokes流动的结果进行了比较。文章还对方法的收敛性进行了数值研究。     

并列圆柱绕流的格子Boltzmann数值模拟

采用非均匀不可压格子Boltzmann模型对低雷诺数下并列圆柱绕流进行了数值模拟,给出了数值计算结果,分析了间距g对圆柱尾流及升力、阻力的影响,并在此基础上得到了4种尾迹模式.此外,研究了流场的初始扰动对流动分岔现象的影响,发现在适当的扰动下可以很快得到同步同相的尾流.对Re=160和200下圆柱的升、阻力进行了对比,结果表明升力和阻力受间距g的影响大于雷诺数.     

圆柱尾流雷诺应力与平均运动变形率的空间弛豫效应

针对圆柱尾流中沿流向存在的Karman涡街周期性涡旋结构,对湍流雷诺应力与平均运动变形率之间的空间弛豫效应进行了实验研究.在回流式水槽中,放入不同直径的圆柱模型,获得不同雷诺数下的圆柱尾流,利用二维高时间分辨率粒子图像测速(TRPIV)技术测量圆柱尾流二维瞬时速度空间分布图像的时间序列.经过数字图像处理,获得二维雷诺应...     

圆柱绕流流场结构的大涡模拟研究

为进一步揭示绕流现象的形成机理,本文分别对处于层流稳态区、尾流过渡区、剪切层转换区Re分别为26、200、1.4×105的三种典型流态下的单圆柱绕流进行了二维数值模拟研究。Re为26时应用层流模型直接求解N-S方程,而Re分别为200、1.4×105时使用大涡模拟的方法进行计算。数值模拟很好地再现了稳定的涡旋结构、周期性交替脱落的卡门涡街结构、不规则的涡旋结构,在此基础上分析了尾流结构的基本特征及其压强分布规律、平均的流场特性、积分参数(如升力系数、阻力系数、斯特劳哈尔等),并与有关研究成果进行了对比。研究发现,采用不同流动介质时流场特性有所差异,空气为介质时的计算结果更符合实验的成果,而水为介质时计算结果偏差较大,这主要是由尾流涡旋产生的不合理负压造成的。     

水面舰船尾流电导率信号分布规律的研究

在舰船尾流区与非尾流区之间存在显著的速度差和盐度(或密度)差,利用电导率探头可获得对应于这些差别的尾流电导率信号。在水槽中形成了与海洋环境类似的盐度分层流场,由双螺旋桨自航水面船模产生尾流,分别在盐度分层流场和非分层均匀流场中测量了尾流电导率信号的横向分布,并对盐度分层流场中水面舰船尾流的纵向速度分布进行了数值计算。结果表明:在盐度分层流场和非分层均匀流场中水面舰船尾流的电导率信号沿其横向近似呈高斯分布;尾流速度对其电导率信号的影响比盐度梯度的影响大得多;尾流的无量纲纵向速度亏损的数值计算结果与尾流的相对电导率信号横向分布的实测结果具有很好的一致性。     

风嘴对旋转振荡板旋涡脱落的控制

5∶1矩形柱体被认为是通用桥面几何形状的代表，其简化模型可以用来进行风振控制的研究。风嘴作为一种常用的流动控制装置，能起到减阻和增加矩形板气动稳定性的效果，但控制装置对强迫振荡柱体的控制机理仍缺乏研究。研究者通过对节段模型施加强迫振动，在实验模型前部施加边缘型对称型风嘴，研究控制装置对矩形板尾流旋涡脱落的影响。通过流动显示结果，总结了施加风嘴后的四种旋涡脱落模式。并通过快速傅里叶变化频域分析法，得到实验模型后缘X/D=10处的速度功率谱。施加风嘴控制能增大旋涡脱落频率，并抑制尾流旋涡脱落的能量。而数值模拟得到的不同实验工况下的升力均方根和力矩均方根显著减小，最大降幅分别为52%和23%，表明矩形板气动稳定性的提升和尾流不稳定性的减弱。     

通气空泡尾部气泡流仿真和特性分析

基于欧拉-欧拉双流体模型对通气空泡尾部气泡流进行数值仿真, 并采用基于多尺寸分组模型的总体平衡方法预估气泡尺寸分布. 应用改进后湍流耗散系数计算模型, 考虑了气泡体积含量对在湍流作用下气泡扩散现象的影响. 基于上述模型对两种试验工况下流场进行了数值仿真.结果表明模型对空泡尾部回流区特性进行了准确预示, 在回流区高湍流度作用下气泡迅速破碎成小气泡. 并进一步得到试验体尾流区空泡体积分数和速度分布. 尾流区水流速度分布保持了流体经过非流线型对称体时产生的尾流分布规律. 仿真结果与试验数据相一致, 模型适用性得到验证.      

A NONLINEAR FLUID–STRUCTURE INTERACTION ANALYSIS OF A NEAR-BED SUBMARINE PIPELINE IN A CURRENT

A negative lift force (attraction) can be induced on a near-bed pipeline in a horizontal current due to asymmetric flow. This negative lift force has a significant influence on the behaviour of the near-bed pipeline, causing two remarkable failure patterns. One failure pattern is due to stability loss, and the pipeline fully rests on the seabed. The other is due to the excessive stress or deformation even if the pipeline is stable in a position between the original equilibrium position and the seabed. A quantitative method to assess these two failure patterns by combining boundary element and finite element methods is proposed in this paper. This is a nonlinear fluid–structure interaction problem, and an iteration procedure is used herein to solve it. Numerical examples reveal that there exists a critical current velocity, above which the pipeline fails. The relationship between the critical velocity and the distance from the pipeline to the seabed is given.     

海底管道离地点的接触条件探讨

海底管道的铺设分析和提升分析都用到管道和海床之间的接触条件，这个同时含有挠度、转角和弯矩的接触条件，在一般的力学模型中是很难理解的．考虑海床变形，通过海床和管道的变形分析，给出这种接触条件的真实背景和使用条件．     

浅水区海底管道周围海床渗透压力计算研究

我国海上油田开采起步较晚,大部分油田处于浅水区,因此,在设计管道时,应充分考虑由浅水区波浪引起的管道周围海床渗流力。根据浅水波相关假设,考虑自由水面非线性影响,推导出椭圆余弦波的波面方程,在此基础上进一步得到一个关于速度势的表达式,并根据该表达式得出作用于海床表面的波压公式。考虑海床土的压缩性,推导出一阶近似椭圆余弦波作用下浅水区埋置管道周围海床的渗流压力解析解,最后将计算结果与大型水槽试验及以往研究成果作对比。结果表明,在椭圆余弦波的作用下,由一阶椭圆余弦波理论得到的计算结果与试验结果规律基本一致,与相似工况下的现有理论成果数值基本相同,具有一定的可行性和工程价值。     

On the interaction of water waves and seabed by the porous medium model

The interaction of water waves and seabed is studied by using Yamamoto's model, which takes into account the deformation of soil skeletal frame, compressibility of pore fluid flow as well as the Coulumb friction. When analyzing the propagation of three kinds of stress waves in seabed, a simplified dispersion relation and a specific damping formula are derived. The problem of seabed stability is further treated analytically based on the Mohr-Coulomb theory. The theory is finally applied to the coastal problems in the Lian-Yun Harbour and compared with observations and measurements in soil-wave tank with satisfactory results. The project supported by the National Science Foundation of China     

激励Stuart-Landau方程的研究--周期解、稳定性及流动控制

解析得出了有外部激励的Stuart-Landau（S-L）方程的频率锁定周期解,对这些解与外部激励振幅和频率的依赖关系做了详细研究,并用周期系统稳定性理论确定了解的稳定性边界.还对S-L方程所描述的流动控制效果进行了研究,发现由于外部激励的作用,稳定的锁频解可能比原来的饱和解能量减少了,外部的控制最多能使扰动能量减少为原来的一半.     

Absolute and convective instabilities of two-dimensional bluff body wakes in ground effect

Local and global instabilities are investigated of wakes of general two-dimensional bluff bodies placed near and parallel to a plane boundary or ground. A spatio-temporal linear stability analysis is first applied to a four-parameter family of local wake profiles to investigate the fundamental local stability characteristics of the wake in ground effect. The analysis shows significant dependencies of the stability characteristics of the wake on the distance from the wake centreline to the ground (normalised by the wake width), and also on the velocity ratio of the near- and far-ground sides of the wake. The analysis is then compared with earlier experiments on a circular cylinder to examine, according to the transition scenario of the steep global modes, the streamwise variation of the local stability characteristics of the wake in ground effect. The comparison indicates that the near wake region of the cylinder changes from being absolutely unstable to being convectively unstable when the cylinder comes down into the near-ground range in which the von Kármán-type vortex shedding from the cylinder is suppressed, being qualitatively consistent with the transition scenario for general wake-type flows. A possible explanation is also given for the counter-intuitive relation between the thickness of the boundary layer on the ground and the critical gap distance for the cessation of the von Kármán-type vortex shedding in ground effect.     

Analytical solutions of tsunamis generated by underwater earthquakes

Tsunamis induced by underwater earthquakes are theoretically analyzed by applying the linear potential theory. Special attention is placed on the initial state of tsunami. For instantaneous seabed deformations, analytical wave solutions induced by three fundamental seabed deformations at initial stage are derived rather than integral expressions in past studies. These analytical solutions constitute a fundamental base for analyzing waves generated by arbitrary seabed displacement with the help of Fourier analysis. Tsunamis induced by non-instantaneous seabed deformation are analyzed as well. For the sake of examining the contributions of all wave components involved in the tsunami waveform, the amplitude density is proposed to examine the effects of deformation width, water depth, harmonic mode and rising time on waveforms. Results show that a larger ratio of water depth to deformation width results in a greater difference between initial waveform and seabed deformation, and the effect of the rising time is significant in deeper-water configuration. For cosine and sine seabed liftings, the effects of higher harmonic modes might be ignored.     

Experimental investigation of the turbulent wake past real seabed elements for velocity variations characterization in the water column.

In high flow velocity areas like those suitable for marine energy application, bathymetry variations create strong velocity fluctuations in the water column. It is therefore essential to characterize the turbulence evolution in the wake of seabed elements which may impact the loads on tidal turbines. For that purpose, experiments are carried out in a flume tank with R_e as high as achievable in Froude similitude, with bathymetry variations experimentally represented with various wall-mounted square elements of height H: a cylinder or a cube as unitary obstacles and combinations of these elements followed by an inclined floor to resemble smooth bathymetry changes. The onset flow is a simple boundary layer profile with height 1.3 H and a low turbulence intensity. PIV and LDV measurements are used to investigate the wake past all test cases in order to distinguish high floor elevation cases (unitary obstacles) from mean roughness effect (obstacle combinations). Results show that the obstacle combinations produce a wake less extended than for a single wide cylinder that produces an extended wake and very energetic turbulent events. With a single cube, no downstream development of large turbulent events exist and the wake reduces by a factor of 3 compared to the wake cylinder case. An inclined floor downstream of a single wall-mounted obstacle reduces its wake length but does not alter the turbulent structures shed. Turbulent velocity profiles extracted from every wake topology investigated are also compared. The general conclusion is that: for small aspect ratio cases, the obstacle will not affect the water column. On the contrary, strong energetic turbulent events are emitted from large aspect ratio obstacles. Combinations cases stand in-between.     

固体颗粒对半球扰动尾迹影响的实验研究

采用粒子图像测速技术(particle image velocimetry, PIV)研究固体颗粒对放置在平板层流边界层中半球粗糙元尾迹的影响. 实验采集了清水和加入粒径为140 $\mu$m, 220 $\mu$m, 350$\mu$m聚苯乙烯固体颗粒4种工况下二维速度场信息, 基于半球半径的雷诺数为994 ($Re_R=RU/\upsilon$), 固体颗粒的体积浓度为$3.0\times10^{-5}$. 对比清水和两相工况下的平均速度剖面、湍流强度等宏观统计量, 分析固体颗粒对半球尾迹流动宏观特性的影响. 分别利用沿流向不同位置的流向脉动速度的二维空间相关系数和法向脉动速度的功率谱密度函数分析颗粒对尾迹结构演化过程及尾迹结构脱落频率的影响. 结果发现: 与清水相比, 回流区随颗粒粒径增大而逐渐增大; 颗粒使湍流强度增大, 回流区的存在导致在半球后流向位置$2R$前后区域湍流强度呈现不同变化趋势; 颗粒使尾迹结构的流向尺度减小并且随着颗粒粒径的增大先减小后增大; 在尾迹结构运动过程中颗粒的存在促进了尾迹结构的周期性加速和减速运动, 促进作用随着颗粒粒径的增大先增强后减弱; 颗粒的存在促进了尾迹结构的脱落, 脱落频率随颗粒粒径的增大先增大后减小.      

考虑床面边界层效应的波浪海床相互作用分析

在考虑海床表面附近黏性边界层的影响下, 考察了波浪与弹性多孔海床之间的相互作用．波浪场包括势流部分和边界层部分,海床域控制方程由比奥固结理论给出．波浪场和海床域通过交界面处应力连续和速度连续条件进行耦合．在简谐波和小变形的前提下,通过联立求解势流方程、波浪边界层方程和海床准静态比奥固结方程得到了波浪运动及相应的海床动力响应的解析解．通过与以往文献的实验结果进行对比,解析解的合理性得到了验证．通过参数分析讨论了实际问题中需要考虑波浪和海床相互作用机制的海床土质条件,以及流体黏性对波能衰减的影响规律．