高雷诺数时串列双圆柱脉动压力的实验研究

本文通过风洞实验给出了二维串列双圆柱在高雷诺数时(分别为Re=3.25×10~5和Re=6.5×10~5)的脉动压力分布。实验在低湍流度的均匀流中进行,模型表面光滑。实验结果表明:在这二种分别属于高亚临界和低超临界雷诺数下,串列双圆柱不同间距比时的脉动压力分布差别很大。Re=3.25×10~5时前后圆柱间的相互干扰对圆柱表面脉动压力分布的影响要比Re=6.5×10~5时大得多,而且二者干扰的规律和特性也截然不同。     

圆柱-平板角区湍流流动控制的风洞试验研究

角区流动中马蹄涡系的存在通常会造成不良影响。对圆柱-平板角区流动,在圆柱上游放置一倾斜的小圆棒能够改变角区流动结构。利用油流法和平板表面压力测量方法探讨了湍流流态下不同的小圆棒对平板表面的摩擦力线和压力分布的影响。油流实验揭示了倾斜棒能够改变角区的三维分离,新的分离线由倾斜棒和圆柱共同作用引起;倾斜棒对角区的作用可归类为两种不同的物理现象;倾斜棒能够引起圆柱侧面的分离线向下游发生极大的迁移,导致圆柱底部区域尾迹变窄。平板表面压力测量实验揭示附加的倾斜棒能够极大地改变压力分布情况,角区的逆压梯度相应减小;由此,逆压梯度引起的三维分离必然被削弱。     

利用分区ALE算法数值模拟圆柱湍流涡致振动

用分区算法结合任意拉格朗日-欧拉法(ALE)数值模拟了圆柱湍流涡致振动.求解基于非交错网格系统,压力求解采用压力泊松方程提法,湍流模型采用标准k-ε模型和重整化群(RNG)k-ε模型.计算取中等雷诺数Re=5000、10000、15000、25000、50000等,质量系数M=10.阻尼系数ζ=0.00331,自振频率fn=0.18315、0.1628.计算结果表明:湍流涡致振动下圆柱时均阻力系数大于孤立圆柱绕流,而升力系数(振幅)值都小于孤立圆柱绕流.随着雷诺数增大,湍流粘性系数随之增大,但湍动能和湍流耗散率变化趋势不明显.对孤立圆柱绕流,研究结果与前人实验结果基本相符.     

三圆柱绕流的实验研究

通过风洞实验研究了在亚临界雷诺数下成等边三角形布置的三个相同直径的二维圆柱在不同风向角下的表面压力分布的变化。研究结果表明,在较小间距比时,随风向角的变化会产生由邻近、剪切层和尾流三种不同机理引起的相互干扰。在一定的风向角下,由于上游柱分离剪切层的直接作用,处于下游位置圆柱的内侧或外侧会出现大面积的负压分布,因而产生很大的横向力。这种在圆柱表面单侧形成的大面积负压分布与某种同期性的旋涡脱落或分离剪     

高雷诺数时串列双圆柱平均压力的实验研究

通过风洞实验给出了二维串列双圆柱在高雷诺数时的平均压力分布，实验在低湍流度的均匀流中进行，模型表面光滑。实验结果发现；在亚临界雷诺数时串列双圆柱临是距时有一些新的现象。在超临界雷诺数时串列双圆柱的平均夺力分布显示出与亚临界雷诺数时截然不同的特性。就总体而言，超临界雷诺数时前后柱间的相互干扰要比亚临界时小得多。     

采用上游扰动体被动控制减阻技术的实验研究

通过风洞测压和水洞流动显示实验,研究了不同直径的扰动体尾迹下圆柱的流动结构和气动特性。结果表明当扰动体和圆柱串列布置时,流动呈现出两种流动模式:空穴流动和尾流撞击模式。当圆柱处于空穴流动模式时,其阻力可以大幅度降低。当扰动体直径为圆柱直径1/2时,圆柱阻力最大降低98%。综合绕流结构和圆柱表面压力分布的变化,指出阻力减小的原因主要是扰动体对圆柱前体的保护作用以及对后体流动的干扰作用,其中后者占主导地位。随着扰动体直径的增加,扰动体对圆柱前体的减阻效果增加。     

民用飞机静压孔布局规律

对民用飞机静压孔布局规律进行了研究。确定相关条例中静压孔的安装使用要求,将其转换成具体的设计分析技术指标;利用数值模拟方法对翼身组合体构型进行气动计算,得到不同状态下机身表面压力分布;对机身表面压力分布均方根进行分析,得到机身表面静压随马赫数、攻角变化较小的区域,定义为稳压线;将机身等直段等效为圆柱体,将来流分解为轴向平直流动和圆柱周向绕流,通过圆柱绕流Cp分布理论分析,获得与计算结果一致的稳压线分布规律,表明稳压线分布规律的普适性。开展高低速测压风洞试验,对稳压线规律进行验证,结果表明稳压线上压力分布变化范围趋近于0,试验结果与CFD计算结果以及理论分析结果高度一致,证明了稳压线计算及分析方法合理、可行,静压孔布局规律的正确、普适。所得稳压线分布规律可为常规布局民用飞机静压孔布局提供直接参考。     

正三角形排列三圆柱绕流与涡致振动数值模拟

基于有限体积法求解二维粘性不可压缩N-S方程,对低雷诺数下正三角形排列等直径三圆柱的绕流与涡致振动进行数值模拟。圆柱振动简化为两自由度的质量-弹簧-阻尼模型,圆柱的动力响应通过Newmark-β方法求解。间距比在1.5~6.0之间变化,重点研究各圆柱的气动力和响应及相关的频率特性、尾流流动模式随间距比的变化。研究结果表明:三圆柱发生振荡时的气动力要远大于绕流结果,当间距比在3.5~4.0之间气动干扰最为强烈。当间距比大于3.0时,下游圆柱的横向和流向位移相位差出现周期性切换,导致升力系数均值出现波动,并且下游圆柱的质心运动轨迹为斜向上和斜向下椭圆形的叠加。三圆柱的最大横向位移幅值均能达到0.9D,下游圆柱的最大流向振幅达到1.1D,这说明组成多圆柱振荡系统单柱的横向位移远大于孤立圆柱发生涡激共振的最大位移,并且下游圆柱的流向振荡不可忽视。     

用切比雪夫多项式对实验数据的曲线拟合

本文研究利用切比雪夫多项式对实验数据的曲线拟合。插值点的选择利用了代数插值余项的小化原理。利用切比雪夫多项式的性质，导出多项式系数的计算公式，它具有简洁特点，同时也易于计算机的实现。本文选择了两个算例，其一是二元机翼绕流，其二是圆柱绕流型阻。对这两个算例分别对表面压力分布及圆柱尾流区相对速度分布进行曲线拟合。在此基础上，对切比雪夫多项式拟合方法进行了讨论。     

基于微型压力传感器阵列的边界层分离点检测方法

为了检测边界层的分离点,探讨了根据表面压力分布变化判定分离点的检测方法,提出了基于微机电系统(MEMS)技术和柔性衬底的微型压力传感器阵列的结构方案,可实现在线实时测量,并满足非平面表面的流体测量要求;同时提出了1种分离点检测的判定方法,并通过二维圆柱绕流仿真验证了该方法的正确性和有效性。微型压力传感器阵列的引入,拓宽了边界层分离点检测的解决途径。     

前置圆柱列涡轮静叶栅流场特性的数值研究

数值模拟了前置圆柱列涡轮静叶栅的三维流动,详细分析了圆柱列动静状态下叶栅流向的气动特性变化,发现不同时刻上游尾迹扫过下游叶片排不同位置,是造成多级涡轮非定常特性的主要原因。通过圆柱列动静两种状态下叶片表面压力脉动、流场变化特性的比较,发现静叶表面的压力分布存在差异,静叶前缘受上游圆柱列影响较大;随着流动的继续至叶栅流道下游,差异逐渐变小,直至消散;而叶栅前缘的差异主要通过叶片自身具有的大转折角逐步被抹平。     

谱单元法及其在多圆柱绕流分析中的应用

系统阐明谱单元方法,基于谱单元方法对低雷诺数Re=200时不同间距下的顺排两圆柱和Re=150正方形排列的四圆柱绕流及其阻力系数、升力系数等进行数值模拟。研究比较不同间距比L/D(两圆柱圆心距离与圆柱直径之比)对两圆柱和四圆柱绕流的影响,计算分析了涡量图分布、平均阻力系数和斯托罗哈数随间距比的变化。研究表明,间距比对顺排两圆柱和正方形四圆柱绕流影响显著;顺排两圆柱绕流存在临界间距比,在Re=200时临界间距比约为3.6。正方形排列四圆柱存在三种流态。当流场从一种流动形态变成另外一种流动形态时,力学参数发生显著变化,在某些间距比区间内出现骤升或骤降现象。     

圆柱振荡绕流对淤泥床的冲蚀

本文实验研究不同圆柱直径和不同初始圆柱与床面距离时圆柱振荡绕流对其下方淤泥介质床面的冲蚀作用.实验得到了在产生冲蚀的临界条件下圆柱与床面初始无量纲距离随振荡流Kc数的变化规律.由于不同的无量纲间距下圆柱振荡绕流分离流动的旋涡运动类型的不同,形成了冲蚀床表面的不同类型的表面冲蚀图象.     

气膜冷却对涡轮叶片表面压力分布影响的试验研究

本文在中温，中压条件下，对有气膜冷却的复合式涡轮叶片的表面压力分布进行了试验研究。探讨叶片冷气流量比，气膜孔开设位置，燃气与冷气的绝对温度比对表面压力分布的影响。介绍了在中温，中压条件下，各冷却参数对涡轮叶片表面压力分布综合影响的试验研究结果。     

基于表面压力信息的空间流向涡识别方法

旋涡与表面的相互作用广泛存在于各类飞行器的绕流中，旋涡影响飞行器表面的压力分布，引起其气动特性的改变。而表面压力分布同样反映飞行器绕流中复杂涡系的空间流动特征，结合迎角、侧滑角等来流参数，可以判断飞行器受力状态和运动趋势。以平面点涡和“镜像涡”理论为基础，建立基于表面展向压力分布曲线的空间涡识别方法。根据截面展向压力分布曲线是否存在因主涡诱导下洗气流产生的正压区，定义流向旋涡的近物面流动和远物面流动；根据二次涡相对于主涡的位置，定义压力分布曲线的“近二次涡侧”和“远二次涡侧”。利用展向压力分布曲线“远二次涡侧”的1/4峰值、峰值及其展向位置，识别流向旋涡空间位置特征和强度特征。搭建涡-面相互作用试验平台，比较空间流场测量结果和表面压力信息识别结果，验证该方法有效性。研究结果表明：通过表面压力分布曲线可以辨识流向旋涡的空间位置和强度特征，空间流场测量与基于表面压力信息的旋涡识别结果的关联性分析验证了该方法的有效性。为重构飞行器周围旋涡流动结构，以及实现飞行器气动力的预测奠定了重要的技术基础。     

多孔介质空心圆柱绕流的摄动解

本文主要讨论了绕多孔介质空心圆柱的有势流动问题。运用摄动方法,给出了速度势的四阶近似,然后得到了圆柱表面的无量纲的压力系数和阻力系数。结果发现在势流的情况下多孔介质圆柱受到附加阻力作用,并且在圆柱的喷射或吸收流体的特征速度与来流特征速度的比值是一个小量的时候,阻力系数与特征速度比成简单的线性关系。     

DPIV对圆柱绕流的研究及运动参考坐标系的转换

采用以氢气泡为示踪粒子的数字粒子图像测速技术(digital particle image velocimetry，DPIV)对三种不同雷诺数下的圆柱绕流进行了研究。一般DPIV拍摄时相机固定，可得出一种速度场分布。笔者采用静止坐标，用简单的方法可以达到相机随流场移动的效果，DPIV处理结果可以看到一些更深入的情况。     

随机分布介尺度油滴群的燃烧特性

对液雾燃烧中随机分布油滴群的蒸发/燃烧特性进行了数值模拟研究,详细考虑了油滴尺寸、坐标分布、数密度等相关参数的随机性。采用一步总包反应机理对其进行大涡模拟,在不同的来流速度、环境温度下,剖析了不同数密度油滴群的阻力特性及蒸发率,并结合油滴群燃烧模态分析油滴间的相互作用规律。研究表明,油滴群燃烧模态同时包含单油滴三种油滴燃烧火焰模态;各油滴表面燃料蒸气的浓度分布呈现非对称性,反应核心区发生角度偏转;在强对流环境下,蒸发/燃烧油滴群的平均表面蒸发率与阻力系数的变化与蒸发/燃烧单油滴相似,油滴间的强烈相互作用加速了油滴群蒸发/燃烧的传热传质。      

垂直交叉双圆柱绕流强迫振动数值模拟

采用虚拟体法数值模拟了平面单圆柱绕流强迫振动,以及空间垂直交叉双圆柱绕流下游圆柱强迫振动流场。验证了单圆柱强迫振动中的锁定状态以及相位突变现象,从而证实了该数值方法模拟振动流场的可靠性。研究了在雷诺数Re=150、间距为5倍圆柱直径、下游圆柱按正弦曲线振动时,对尾流场的影响。下游圆柱两端尾流场在振动的作用下,涡的横向间距增大,而中心尾流场由于受上游圆柱尾流的影响而保持原先的状态。下游圆柱在锁定区的振动使尾流场变得稳定,表现在流向二次涡结构的减少,以及尾流场的速度场频谱趋向单一化。     

一种结合壁面模型模拟湍流流动的间断界面浸入边界法

提出了一种具有鲁棒性的间断界面浸入边界法，用来模拟不可压湍流的物体绕流问题。开发了一种新的OpenFOAM求解器，该求解器结合了压力隐式分裂算法和k-ωSST湍流模型求解Navier-Stokes方程。为了减少对近壁面网格的要求，在实施边界条件时，采用壁面模型计算得到的壁面剪切力将速度和湍流变量修正联系起来。壁面切应力采用远离壁面的镜像点求得，镜像点上流场变量通过反距离线性插值由周围流场信息得到。使用该求解器模拟了平板绕流、NACA0012表面压强和切应力分布、Buice扩散管、圆柱绕流、方柱绕流5个湍流流动问题，本文数值模拟结果与文献结果有较好的一致性，验证了该方法模拟湍流流动的有效性和准确性。