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

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

方孔横向射流的大涡模拟研究

采用大涡模拟方法,对射流雷诺数ReD=4700、吹风比R=1.5的方孔横向射流进行了计算,为处理完全三维的复杂流动,亚格子模型使用了局部平均的动力模型.通过对时间平均流场的分析,了解了射流的穿透、扩散性质以及肾形涡的形成原理;而瞬态流场的时间演化序列则揭示了大尺度相干结构在流场演化中的主导作用.计算结果表明:使用局部平均动力模型的大涡模拟方法健壮性好,具有较高的数值精度,能够准确给出方孔横向射流的复杂流动图画,适合进行这类复杂流动问题的机理研究.     

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

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

突然起动圆柱的非定常绕流计算

本文用分区方法定量模拟高Re数情况下突然起动圆柱的初期对称流动。绕圆柱的二维不可压流场可分为无粘区、附着粘性区和分离区。无粘区有解析解，用非定常层流边界层方程解附着粘性区，用离散涡方法处理随时间变化的分离区。试用了混合的Rankine和Lamb涡模型及二次离散涡。给出了Re=300、9500两个算例，基本符合有关实验结果。     

冲击射流流动换热超大涡模拟研究

为了准确预测发动机热端部件中广泛采用的冲击射流冷却复杂的流动和换热特性,发展了基于BSL k-ω模型的超大涡模拟(VLES)高精度模拟方法,并对高雷诺数Re=4×104,两种不同射流距离2和6的单孔冲击射流及三孔冲击射流这一经典的流动传热问题进行三维非稳态高精度数值计算。同时,将分离涡方法 (DDES)和k-ωSST,RNG,Transition SST三种RANS方法的数值模拟和开发的超大涡模拟(VLES)方法进行对比。研究表明,VLES方法均能够准确捕捉冲击射流流场的复杂非稳态流动及传热特征,包括自由射流区、壁面射流区小尺度涡系和大尺度湍流结构的演化和破碎,同时冲击壁面的换热系数计算结果与实验值吻合较好。DDES方法未能准确捕捉流场复杂的小尺度湍流结构,壁面换热计算结果与实验值差异较大。RANS方法计算的换热结果与实验数据差异最大,基本未能预测到壁面换热特性。在相同的计算网格和计算方法下,VLES方法计算结果优于DDES方法,DDES方法一般好于RANS方法。这表明新开发的VLES方法能够准确地计算冲击射流相关的流动及换热问题。     

基于各向异性非结构网格的超声速流动自适应计算

基于各向异性非结构网格,实现了自适应求解技术,将其应用于超声速楔形体绕流及超声速横向喷流问题的计算研究,对各向异性网格自适应计算过程的收敛性和求解精度进行了对比分析,展示了各向异性网格自适应算法在降低问题求解规模。提高各向异性流场分辨率等方面的优势及存在的问题。超声速楔形体无粘绕流场的自适应计算,清晰地捕捉到了激波附近的流场信息,激波前后的马赫数、压力、密度。温度等物理量与理论分析解吻合很好。超声速横向喷流流场存在激波、分离涡、边界层等流场结构的相互干扰。计算研究表明,单纯基于Hes-sian度量张量的各向异性网格生成及自适应算法不能有效模拟边界层内的流动情况,是将来需要进一步开展的研究挑战。     

对称驱动方腔定常流动的多解现象研究

对方腔内顶盖正弦速度驱动的定常流动进行了数值研究.采用非线性动力系统数值分岔分析方法,讨论了Re<2000时Navier-Stokes方程定常解的多解共存现象,计算了对称和非对称定常解的分岔曲线和相应的流动模式,给出了非对称解出现的临界参数Re=320,稳定的对称和非对称解同已有文献的结果一致,并获得了一组新的不稳定的非对称解.随着Re数的增加,侧壁出现流动分离并有小涡出现,随之发生流场拓扑结构分岔,但并不伴随N-S方程解的分岔.     

非稳态单喷口射流控制扩压器内边界层的分离流动

对扩压器内单喷口射流产生的流向涡控制扩压器内部边界层分离流动行为进行了研究。研究的对象是一个扩张角为10°的二维扩压器,扩压器进出口面积比为1∶4.7,入口高度为15 mm。采用数值计算方法研究了当射流出口速度为非稳态时对分离流动的影响。给出了不同时刻,不同截面上的旋涡涡量分布。结果表明,采用非稳态喷射流向涡方法可以在不改变主流流量特性的情况下,有效地减小分离流动区域。     

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

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

三个丁坝流态研究

利用流动显示和片光源技术对绕丁坝的流动结构在北航水槽中进行了研究。研究发现在丁坝的上游会形成马蹄涡，在丁坝的头部会形成卡门涡，按Re数的大小可以把绕丁坝的流态分成5类。实验表明丁坝的交角及丁坝的个数对丁坝的流态有重要影响。     

在激波管中粒子阻力系数的实验研究

在激波管中用消光法测量了气固两相流动激波松弛区的粒子浓度。基于一维静止激波理论和理想“等效气体”模型,实验确定了粒子阻力系数。拟合实验结果,得到了阻力系数与粒子雷诺数之间的实验关系式:C_D=4300/Re~(1.57)。     

双尾撑倒V尾无人机气动特性及静稳定性计算分析

在双尾撑无人机上配置倒V 形尾翼，其夹角会影响全机的气动特性和静稳定性。提出一种双尾撑倒V形尾翼布局，采用数值方法与U 形尾翼布局进行对比分析；对倒V 形及U 形尾翼布局的计算模型进行校核和验证，分析两种布局在不同迎角和侧滑角下的纵向、横向和航向气动特性及静稳定性。结果表明：倒V 形尾翼可减小浸润面积，使全机升阻比提高2.2%，增大了尾翼失速迎角，但横向静稳定度明显下降；当尾翼夹角小于90°时，倒V 形尾翼能够增强全机横向稳定性，夹角增大也会使航向静稳定度略有降低，而纵向静稳定度有一定增强。     

趋于临界马赫数的圆柱跨声速绕流特性分析

通过深化认识趋于临界马赫数Ma_cr的圆柱跨声速绕流特性,明确新型飞行器增升减阻设计的空气动力学理论依据。采用大涡模拟方法数值研究了来流马赫数Ma_∞为0.75和0.85、雷诺数Re为2×10⁵的圆柱跨声速绕流。结果表明:当Ma_∞趋于临界马赫数(Ma_cr≈0.9)时,圆柱的阻力下降且升力系数振荡被抑制;通过力的分解,得知圆柱的阻力减小来自旋涡力的影响,而非可压缩性;圆柱的阻力减小与其背压上升有关;剪切层初始阶段的对流马赫数Ma_c随Ma_∞的增加而增大,而增长率相反,这使得剪切层更为稳定、柱体背压更高。此外,由于Ma_∞=0.85时边界层分离点处的激波和尾迹处的激波向下游推移,使得近尾迹处的湍流脉动减弱,进而导致柱体的表面压力振荡和升力系数振荡被抑制。     

后掠翼边界层横流不稳定性的实验研究(英文)

在西北工业大学的低湍流度风洞中,采用升华法研究有无粗糙带情况下,45°后掠翼三维边界层内的横流驻波不稳定性及其转捩模式。在未引入人工粗糙带,雷诺数为5.50×105~1.65×106范围内,模型的转捩分界始终为位于最小压力点之后的一条直线,转捩由T-S波触发。当Re≥1.38×106时,对应最不稳定横流驻波的3.5~4.0mm条纹出现在转捩的上游区域,条纹间距与线化稳定性理论的结果吻合。当Re=1.65×106时,实验证实了横流驻波扰动对前缘粗糙度的极度敏感。考虑到抑制最不稳定横流驻波就很有可能抑制后掠翼飞行器上由其主导的转捩,在机翼前缘布置不同间距的粗糙带,研究其对边界层内横流驻波的影响。当Re=1.38×106时,2.5mm间距的粗糙带有效的抑制了3.5~4.0mm最不稳定驻波,该现象为后掠翼上的转捩控制技术提供了一条新思路。此外,当6.0mm、7.0mm和8.0mm的粗糙带被引入时,条纹间隔表现为3.0mm、3.5mm和4.0mm的谐波波长。     

Numerical investigation of the shock interaction effect on the lateral jet controlled missile

A computational study on the supersonic flow around the lateral jet controlled missile has been performed. A three-dimensional Navier–Stokes computer code (AADL 3D) has been developed and case studies have been performed by comparing the normal force coefficient and the moment coefficient of a missile body. Different jet flow conditions including jet pressures and jet Mach numbers, and the circumferential jet positions have been incorporated into the case studies. The missile surface is divided into four regions with respect to the center of gravity, and the normal force and moment distribution at each region are compared. The results show conspicuously different normal force and moment variations according to each parameter variation. From the detailed flow field analyses, it has been verified that most of the normal force loss and the pitching moment generation are taking place at the low-pressure region behind the jet nozzle. Furthermore, it is shown that the pitching moment can be efficiently reduced by the lateral thrust obtained through higher jet Mach number rather than high jet pressure. Thus, an angle of yaw is more effective for missile control by side jet than an angle of attack.     

均匀来流中旋转振荡圆柱绕流的数值研究

通过求解二维非定常不可压缩 Ｎ-Ｓ方程,研究了均匀来流中旋转振荡圆柱的流动特性,系统地改变流动的振荡频率、振幅,分析了圆柱的受力特性及尾迹结构,发现了与横向、流向振荡的圆柱相似,旋转振荡也是流动控制的一个有效方法,有其自身的流动规律及特点。     

八字形出口合成射流激励器机翼分离流控制

设计研制了一种适于机翼分离流动控制的八字形出口合成射流激励器，对其出口射流与主流的相互作用特性进行了研究，粒子图像测速仪（PIV）流场测试和边界层速度型测试结果揭示了其控制机制为促进边界层与主流的诱导掺混，提升边界层底层能量。利用该激励器阵列对NACA633-421三维直机翼模型开展了针对射流能量比Cμ和阵列位置两个参数的分离流控制研究，天平测力及翼型表面测压结果显示该激励器可有效抑制翼面流动分离、推迟失速迎角。在设计范围内，射流能量比Cμ值越大，控制效果越好，当Cμ＝0．00168时，机翼最大升力系数提升了5．92％，失速迎角推迟了2．5°（激励器阵列位于0．3c处）。激励器阵列的弦向布置位置是一个重要控制参数，阵列位于0．3c处时最大升力系数提升量大于位于0．55c时。     

圆锥脉冲等离子体诱导流场特性分析

对等离子体诱导流场特性进行研究,有利于解决双稳态非对称分离涡带来的连续比例控制困难的问题。在封闭光学玻璃箱体内,应用介质阻挡放电等离子体对20°顶角圆锥附近静止大气进行了定常和脉冲循环控制,对等离子体诱导的圆锥截面绕流速度场进行了二维PIV测量,对定常控制和脉冲循环控制下最大绕流速度及最大轴向涡量进行了比较分析。实验结果表明：相对于定常控制模式,脉冲循环控制下沿垂直于圆锥截面对称面径线分布的时间平均切向速度和轴向涡量范围较广;在脉冲循环控制下,动量传递的主要表现在离散涡的形成而不是气流的加速。     

战斗机垂尾脉动压力数值模拟

在亚跨超计算流体力学（CFD）软件平台（TRIP）上开发了基于RANS/LES混合思路的IDDES流动模拟技术,并通过NACA0021翼型60°大迎角分离流动与串列圆柱绕流模拟对RANS/LES混合方法的精确度进行了验证,针对某战斗机外形的垂尾脉动压力开展了数值模拟研究。战斗机来流马赫数为0.1,基于全机长度的雷诺数为2×10⁶,模型迎角为20°、30°和40°。分别通过脉动压力系数、脉动压力功率谱密度、空间流动结构以及侧向力响应曲线等对战斗机的垂尾脉动压力进行了分析。脉动压力模拟结果表明：当垂尾完全沉浸在边条翼脱体涡破碎后的宽频湍流脉动气流中时,垂尾翼梢位置的脉动压力会发生明显的增大。     

典型气动布局高超声速飞行的气动力数值评估

采用自主开发的数值软件SPACER对典型气动布局在Ma=6.0(Re=1×107)高超声速巡航时的气动力性能进行了数值评估。评估的气动布局包括类乘波体(仿X-51A)、翼身融合体(仿ISR)、传统升力体(仿X-33)和轴对称锥形体(仿Fasthawk导弹)。由于升力体范围较广,还对一种相对扁平的升力体模型进行了数值评估。评估结果表明:在考察条件下,类乘波体具有较大的升阻比和较小的阻力,是很有潜力的高超声速巡航气动布局;翼身融合体和传统升力体需一定的改进和进一步的研究,如扁平升力体的气动性能可得到大幅提升;轴对称锥形体以其总阻力小,也具有一定优势。