错列角度对双圆柱涡激振动影响的数值模拟研究

为研究错列角度\alpha对双圆柱涡激振动问题的影响,采用自主研发的基于CIP (constrained interpolation profile)方法的数值模型,对雷诺数\mathit{Re}=100、错列角度\alpha =0^{°}\mathrm{~}90^{°} (间隔15 #x00B0;)的等直径双圆柱涡激振动问题进行数值模拟. 模型在笛卡尔网格系统下建立,采用具有三阶精度的 CIP 方法求解 N-S (Navier--Stokes)方程,采用浸入边界法处理流--固耦合问题,避免了任意拉欧方法下的网格畸变和重叠动网格技术中的大量信息交换问题,保证了模型的计算效率. 重点分析不同错列角度\alpha上下游圆柱的升阻力系数、位移响应、涡脱频率和尾涡模态等. 结果表明:折合速度U_r=2.0\mathrm{~}3.0时,上下游圆柱升阻力随错列角度的增大基本呈单调增大的趋势;U_{\mathrm{r}}=5.0\mathrm{~}8.0时,随错列角度的增大,上下游圆柱阻力变化较小,升力呈“上凸”趋势,在\alpha =15^{°}\mathrm{~}30^{°}取得最大值;U_r=10.0\mathrm{~}13.0时,随错列角度的增大,上下游圆柱阻力变化较小,升力呈“下凹”趋势,在\alpha =30^{°}\mathrm{~}45^{°}取得最小值,且柱体横流向振幅和升力没有明显的对应关系. 最后,结合尾涡模态对以上规律的成因进行分析. 研究结果可为相关海洋工程设计提供参考.      

涡激诱导并列双圆柱碰撞数值模拟研究

圆柱类结构物的涡激振动是工程中较为常见的一种现象,如果圆柱结构物之间的距离较小,就会产生涡激诱导碰撞现象,而涡激碰撞会比涡激振动对结构物疲劳破坏产生更严重的威胁.采用浸入边界法模拟流体中的动边界问题,避免了传统贴体网格方法在求解流体中存在固体间碰撞问题时出现数值求解不稳定问题,采用有限元方法对圆柱的运动和碰撞进行求解,通过数据回归方法建立了流体流动条件下的润滑模型,对不同间隙比下涡激诱导并列双圆柱振动及碰撞过程进行了数值模拟,数值结果表明,如果两圆柱产生了碰撞将会有连续的碰撞发生,碰撞时出现了多阶频率,振动主频率要比无碰撞时大,两圆柱碰撞时的相对速度比自由来流速度小;当两圆柱相互接近时,随着涡环分离角度的逐渐倾斜,横向流体力先逐渐减小,当两圆柱间涡环开始相互影响发生挤压时,横向流体力开始逐渐增大;当两圆柱开始反弹时,两圆柱间形成了低压区,改变了横向流体阻力的方向,使两圆柱又产生了接近运动,如此反复从而产生了碰撞后横向流体力和圆柱速度的振荡现象.     

斑海豹胡须涡激振动及其尾流循迹机理直接数值模拟

在视力和听力受到限制的条件下,斑海豹可以通过具有特殊外形的胡须识别和追踪水中游鱼的尾迹.从仿生学的角度,对斑海豹胡须在尾迹流场中的振动特性和循迹机理进行研究,有助于开发研制新型水下探测器.本文采用嵌入式迭代浸入边界法,在雷诺数Re=300和折合流速Ur=6.0的条件下,对均匀流场和尾迹流场中斑海豹胡须模型的涡激振动进行...     

不同控制角下附加圆柱对圆柱涡激振动影响

在弹性支撑的圆柱周围布置直径更小圆柱会影响剪切层发展以及旋涡脱落,进而改变其涡激振动状态.通过不同的布置形式和附加小圆柱个数可以实现对圆柱涡激振动的促进或抑制.激励更大幅值的振动可以更好地将水流动能转化为可利用的机械能或电能,抑制其振动则可以实现对海洋平台等结构物的保护.采用基于迭代的嵌入式浸入边界法对前侧对称布置两个小圆柱的圆柱涡激振动进行数值模拟研究,系统仅做横向振动,其中基于主圆柱直径的雷诺数为100,质量比为2.0,折合流速为3~11.小圆柱与主圆柱的直径比为0.125,间隙比为0.125.结果表明,在研究的控制角范围内(30°~90°),附加小圆柱可以很大程度上改变圆柱涡激振动的状态.当控制角较小(30°)时,附加小圆柱对主圆柱的振动起抑制作用;当控制角为45°~60°时,圆柱的振动分为涡振和弛振两个阶段,在弛振阶段,圆柱振幅随折合流速增加而持续增加;当控制角较大(75°~90°)时,附加小圆柱的促进作用随着控制角增加而减小.进一步地,结合一个周期内不同时刻旋涡脱落以及圆周压强分布,解释了附加小圆柱对主圆柱涡激振动的作用机制.应用能量系数对圆柱系统的进一步分析发现,弛振阶段由流体传递到主圆柱的能量系数随折合流速的增加逐渐下降,旋涡结构的改变是产生这种变化的直接原因.      

绒毛控制圆柱涡激振动的风洞试验研究

利用绒毛对圆柱涡激振动抑制进行了风洞试验研究。通过改变附属绒毛无因次长度L/D比(L为绒毛长度,D为圆柱外径),研究L/D分别为0.6、1.2和1.8的模型在约化速度2~40的范围内对弹性支撑大质量阻尼系数圆柱涡激振动的抑制作用。试验采用激光位移传感器采集圆柱的横向(Y)和顺流向(X)位移,并用烟线测流场以揭示流动控制机理。结果表明,三种无因次长度的绒毛对大质量阻尼系数圆柱的涡激振动都有显著的抑制作用,随着L/D的增加,圆柱Y向无因次位移及功率谱密度幅值减弱,多达73.5%的无因次位移被抑制;且随着L/D的增加,圆柱附属绒毛频率比远离原始圆柱频率比。绒毛改变了圆柱的边界层分离点位置、抑制了边界层的相互作用并改变尾涡结构,从而抑制振动。     

不同剪切率来流作用下柔性圆柱涡激振动数值模拟

采用浸入边界法对细长柔性圆柱在线性剪切流条件下的涡激振动进行三维数值模拟。细长柔性圆柱振动采用三维索模型模拟,其两端铰接,质量比为6,长细比为50,无量纲顶张力为496。来流为线性剪切流,剪切率从0到0.024变化,最大雷诺数为250。研究发现:剪切流作用下柔性立管横流向振动表现为驻波模式,而顺流向振动表现为行波-驻波混合模式。随着剪切率增大,振动频谱呈现多频响应,振动能量逐渐向低频转移。阻力系数平均值随着展向变化,脉动阻力系数和升力系数的均方根值均表现为“双峰”模式。流固能量传递系数沿立管轴向的分布表明,振动激励区集中于高流速区,而振动阻尼区多位于低流速区。剪切率较小时,圆柱的泻涡为平行交叉模式;剪切率较大时,圆柱的泻涡为倾斜泻涡模式,且由于泻涡频率沿立管轴向变化,尾流发生涡裂现象,形成泻涡频率不同的胞格结构。      

层流条件下并列三圆柱涡激振动响应与尾流形态

应用基于嵌入式压强-力迭代的高精度浸入边界法研究等间距并列三圆柱涡激振动。其中,雷诺数Re=100,间距比T/D=2.0～5.0,圆柱质量比m*=2.0,折合流速Ur=2.0～10.0,忽略振动系统的阻尼且三圆柱仅横向振动。研究发现,圆柱的振动响应随折合流速的增加呈现初始响应分支和下端响应分支两种模式;振幅响应出现不连续现象,且随着间距比的增加,该不连续现象对应的折合流速增加;尾流模式与间距比和折合流速密切相关。共发现六种尾流形态,分别为窄宽窄尾流、不规律尾流、反相同步尾流、调制尾流、同相同步尾流和偏斜尾流。总结并绘制了尾流形态在参数空间[Ur,T/D]内的分区图。     

倾斜圆柱结构涡激振动研究进展

圆柱结构涡激振动现象在生活中十分常见,如海洋工程中的管道、土木工程中的高耸建筑、桥梁斜拉索,核工程中的热交换器等频繁受到涡激振动影响,诱发结构的疲劳损伤,甚至破坏失效.现阶段,人们对垂直来流作用下圆柱结构涡激振动机理已有较为全面的认识.然而,当圆柱倾斜置于流场中,结构后缘的尾流形态与垂直放置差异显著,结构与流体的耦合作用机理更为复杂.为简化倾斜圆柱涡激振动问题,提出了不相关原则,来流速度被分解为垂直圆柱结构轴向和平行圆柱结构轴向的两个速度分量,仅考虑垂直结构轴向速度分量的影响,忽略平行结构轴向速度分量的影响.近年来,针对倾斜圆柱涡激振动及不相关原则的适用性,出现了大量实验和数值模拟研究成果.为了深化对倾斜圆柱结构涡激振动相关机理的认知,本文全面阐述了倾斜圆柱结构涡激振动响应规律、尾迹流场模式和流体力特性等方面的研究进展,分析了不相关原则的适用范围,探讨了倾斜圆柱结构涡激振动抑制措施,并对今后该领域的研究进行了力所能及的展望.     

防风网透流风空气动力学特性大涡数值模拟研究

基于有限体积法建立不可压缩粘性流体运动的大涡模拟模型,采用Smagorinsky-Lilly亚格子模型,并引入浸入边界法(IBM)实现无滑移固壁边界条件,对雷诺数30~30000之间防风网透流风进行模拟研究。基于模拟结果,提出蝶型防风网透流风存在4个典型分区结构,流场中存在由蝶型形态引起的大尺度分层剪切流动,加强流体动能耗散。透流风在雷诺数300时发生层流至湍流的转捩,而在雷诺数增长至3000以上时,湍流充分发展,纵向流速脉动强度可达70%。防风网整体空气阻力远大于单个孔口射流阻力的线性叠加,射流间的相互作用以及大尺度的分层剪切结构大大增加流体阻力损失,这为通过优化孔口布置和网板形态来节省材料提供了科学依据。     

低雷诺数下串列布置双圆柱涡激振动特性研究

基于四步半隐式特征线分裂算子有限元方法,对串列布置双圆柱双自由度涡激振动问题进行了数值模拟计算,并分析了间距比、剪切率、频率比以及折减速度4个参数对圆柱结构动力响应的影响.研究发现:不同固有频率比与剪切率对下游圆柱振动幅值影响较大,然而对上游圆柱振动幅值影响较小.上游圆柱在两个自由度方向达到最大值的折减速度不同,然而下...     

并列双圆柱流致振动的不对称振动和对称性迟滞研究

对雷诺数Re = 100 间距比s/D = 2.5 和5.0 的并列双圆柱流致振动进行了数值模拟研究, 其中圆柱质量比m = 2.0, 折合流速U_r 在2.0~10.0 之间, 两圆柱仅能做横流向振动. 研究发现, 当间距比s/D = 2.5 时, 在折合流速4.4 < U_r< 4.8区间内, 两圆柱流致振动响应出现不对称振动现象, 在折合流速4.4 < U_r< 4.8 区间内, 两圆柱流致振动响应出现对称性迟滞现象; 而当间距比s/D = 2.5时, 圆柱流致振动响应与单圆柱涡激振动响应相似, 没有出现不对称振动和对称性迟滞现象. 在不对称振动区间内, 两圆柱的升、阻力参数也出现了不相等的情况. 此外, 当两圆柱不对称振动时, 圆柱间隙流稳定地偏斜向其中的一个圆柱; 相应地, 尾涡也出现了宽窄不等的模式. 窄尾流圆柱的振幅和升、阻力均较宽尾流圆柱的大. 通过对比不对称振动现象发生前后的尾涡模式, 对新现象的产生机制进行了阐述.      

海洋热塑性增强管(RTP)涡激振动数值计算

基于Van der Pol尾流振子模型和多体系统传递矩阵法(transfer matrix method for multibody systems,MSTMM),建立了可以快速预测海洋热塑性增强管(reinforced thermoplastic pipe, RTP)振动特性和涡激振动响应的动力学模型.仿真结果与ANSYS软件仿真结果以及文献实验数据对比,验证了本文模型的准确性.研究了考虑RTP立管刚性接头,不同顶张力,不同来流分布等情况对RTP立管涡激振动响应的影响.计算结果表明:流速越大,立管涡激振动激发出的模态越高;立管涡激振动主要受低阶模态控制;立管的刚性接头对立管的湿模态影响较小,但是对较高阶模态为主所激发出的涡激振动振幅分布影响较大;剪切流对沿立管轴向的涡激振动振幅分布影响较大,低流速能量小所引起的涡激振动幅值较小,但是当剪切流流速达到能激发出较高阶模态时,相比同等流速的均匀流所引起的涡激振动振幅要大.     

随机压力脉动对涡致振动的影响

用标准随机平均法研究两个线性振子与两个Ｖａｎ ｄｅｒ Ｐｏｌ振子耦合系统在宽带随机激励下的平稳响应，在一定的参数条件及假定下，由各种共振状态得到了以各振子能量及相位差为基本变量的移居记概率密度和概率意义上的分叉条件，用数字模拟证实了解析解的正确性。     

风力机气动特性的浸入边界法模拟1)

风力发电机的空气动力学性能是决定风力机安全与效率的最重要因素之一。但由于影响风力机气动性能的参数众多,更加高效精确地模拟风力机气动特性一直是风力机的重要发展方向。本文提出了基于浸入边界法的风力机建模,网格离散,以及数值模拟的统一性框架。利用同伦变形来生成光滑的叶片模型,并且使用仿射变换来处理叶片的渐缩与扭转问题。首先,针对二维翼型的升阻力,检验了算法的数值精度。表明此方法对于阻力的模拟具有非常严格的一阶精度,进而提出采用理查森外推法来精确高效修正升阻力模拟结果。同时,模拟研究了拱曲度以及厚度对二维翼型升阻力的影响。随后,模拟研究了单风力机(包含塔架)在不同尖速比下的功率系数,并对塔架与叶片间的相互气动作用进行了初步分析。最后,模拟研究了双风力机在风场中不同前后间隔距离下的气动干涉问题。本文主要意义在于验证建模,离散,与数值模拟的一体化框架的有效可行性,进而为后续研究(给定约束下风力机自动优化选型)提供坚实基础。     

基于浸入边界法的扑翼飞行数值模拟

本文提出了基于浸入边界法的扑翼鸟建模与仿真,首先检验了数值方法的精确性,而后对NACA 0012翼型的升沉运动与俯仰运动进行了研究,最后对三维扑翼的翅膀拍动时间非对称性进行了研究。结果表明：浸入边界法对拍动翼型的模拟能够很好地和文献结果吻合。升沉运动的推进能力由翼型前缘涡的大小和位置决定,升沉运动推进效率的峰值主要集中在0.3≤St≤0.4时。升沉运动耦合俯仰运动时,在俯仰角25°及相位差85°时,推进效率达到峰值。在三维模拟中,适当增加翅膀下拍速度,能提供更大的升力,同时耗能也更高。研究结果可以为微型扑翼飞行器的扑动参数设置提供参考。

     

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法. 鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究. 先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解. 对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析. 分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成. 当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导. 当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导. 随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.      

NUMERICAL SIMULATION OF THE PHENOMENA DUE TO THE PASSING-BY OF TWO BODIES USING THE UNSTEADY BOUNDARY ELEMENT METHOD

In this paper a numerical analysis was made to investigate the aerodynamic forces surrounding two bodies in relative motion in a fluid at rest in three dimensions. The unsteady boundary element method was employed in the numerical calculations. This method is very convenient for obtaining an approximate expression of the velocity potential, especially for practical use. The passing-by of two spheres in an incompressible perfect fluid which extends to infinity is treated by the present method. The resultant pressure coefficients on two spheres passing each other in opposite directions are calculated and discussed numerically. Numerical examples are presented to show the validity of the present method. The method is also applied to the calculation of the passing-by of two trains in an open area in order to investigate its applicability.