开口风洞剪切层流场特性

给出了开口风洞剪切层厚度简化计算公式,开展了数值计算,研究了扩散率和偏移量参数对剪切层速度剖面的影响;分别在两座开口风洞中开展了实验研究,得到了剪切层速度剖面、剪切层厚度、湍流强度和湍流功率谱的变化规律,并将理论结果与实验结果进行了对比分析,结果表明:离开喷口一定距离后,剪切层的轴向速度剖面强自相似,而横向速度剖面弱自相似;剪切层的厚度沿流向位置线性增长;湍流强度值在射流核心区较小且基本稳定,在剪切层区域快速增加;剪切层湍流功率谱为宽带谱,没有与周期性旋涡脱落相关的谱峰。     

可压缩剪切层与边界层流动相互作用的稳定性研究

为了探究可压缩剪切层与边界层组合流动稳定性规律,采用线性稳定性理论对这种组合流动的线性稳定性问题进行研究。平均流动模型取双曲正切函数描述的可压缩剪切层与可压缩相似边界层速度剖面的直接叠加。采用整体数值方法 (Global)求解组合流动线性稳定性方程,获得全部不稳定模态。结果表明,在这样的组合流动中,存在两个不稳定模态,一个与剪切层相关的不稳定模态,另一个是与边界层相关的不稳定模态。在此基础上,研究了剪切层与壁面之间的距离、扰动频率、雷诺数、马赫数、速度梯度、温度梯度对这种组合流动稳定性的影响。马赫数增大对不稳定模态起抑制作用,剪切层对边界层的不稳定模态的促进、抑制作用与马赫数、扰动频率有关,边界层对剪切层不稳定模态基本上起抑制作用。     

声学风洞内气动噪声源识别定位方法研究

声学风洞流场和试验环境会影响声音的传播,进而会影响声源识别。根据用于静态声源识别的经典Beamforming算法,利用声波在气流场和风洞射流剪切层中的传播规律,给出了一种能够用于声学风洞试验的基于麦克风阵列的气动噪声源识别分析方法,并进行了数值仿真和数值验算。数值验算结果表明风洞流场和剪切层对噪声源识别效果影响明显,采用本文所总结的计算方法能有效地修正这些影响。     

攻角与马赫数对进气道起动影响的可比拟性

针对高速风洞变来流马赫数的能力有限,不易获得高超声速进气道起动马赫数界限的问题,本文尝试通过改变进气道模型攻角的方法来模拟来流马赫数的变化,分析定来流马赫数变攻角与定攻角变来流马赫数两种途径下二元进气道起动性能之间的关系,以助于风洞实验获得进气道的起动马赫数界限。结果表明,在Ma∞=5.9来流条件下,进气道不起动和自起动的临界攻角分别在12°~13°和0°~1°。定来流马赫数变攻角与定攻角变来流马赫数两种途径下,进气道不起动和自起动临界所对应的内收缩段入口马赫数基本一致。对于同一内收缩段入口马赫数,当变攻角跨度小于4°时,两种路径下进气道的内部流场与壁面压力分布规律符合较好;当变攻角跨度较大时,两者的差别也增大。对于二元高超声速进气道模型,当常规风洞的来流马赫数在进气道起动边界附近时,在一定的攻角范围内,可以尝试通过攻角机构连续改变模型安装攻角的办法来模拟变马赫数引起的进气道起动特性的演变过程。     

树脂含量对F-8H3/602 芳纶复合材料性能的影响

探索研究了F-8H3/602芳纶织物增强树脂基复合材料层合板不同树脂含量对其冲击韧性、层间剪切强度及弯曲性能的影响,并采用显微镜及扫描电镜对破坏试样形貌进行观察。结果表明芳纶复合材料的树脂含量对其冲击韧性、层间剪切强度及弯曲性能影响较为明显,树脂含量为46.02wt%的芳纶复合材料层合板的0°冲击韧性为22 J/cm2、0°层间剪切强度为49.1 MPa、0°弯曲强度为506 MPa;在一定的树脂含量下,芳纶复合材料的0°冲击韧性、层间剪切强度及弯曲性能较90°的性能更优异;破坏试样的形貌表明纤维与树脂界面结合较为薄弱。     

碳纤维复合材料短梁剪切疲劳特性

研究了碳纤维增强双马来酰亚胺复合材料的层间剪切静态与疲劳加载性质,分析了不同基体、不同铺层界面对层剪性质的影响。给出了层剪疲劳加载的S-N曲线;讨论了层间剪切损伤和增强机理,观察了微观特征。与碳/环氧复合材料相比,双马来酰亚胺基体复合材料在碳纤维和基体之间粘接良好,从而改进了纤维抗疲劳裂纹扩展的能力。     

关于超声速剪切流动的数值模拟

本文通过采用混合型的四阶精度频谱关系保持格式和基于信号特征传播的边界处理方法,对对流马赫数分别为0.5、0.9和1.4的时间发展和超声速空间发展的剪切流动进行了数值模拟,目的是对剪切层不稳定结构的演化方式、相同对流马赫数下时间发展和超声速空间发展的剪切流动的异同作初步研究.研究表明,对时间发展的剪切流动来说(1)当Mc=1.4时,剪切层内流动演化的形态为首先压力等值线在其极值点处产生中心型结构,且极大值点和极小值点交替分布.然后中心型结构进一步发展并产生横向分裂,从而在剪切层内等压力线出现鞍点结构.随后剪切层内重新产生中心型的压力极值点结构,进而又出现鞍点结构.如此发展下去,剪切层内的压力等值线呈现中心和鞍点的组合.剪切层外的压缩波沿特征线方向传播并加强,形成类激波的结构,最后在流场中包含许多压缩和膨胀波.在计算时间内流向的周期性未受破坏.(2)当Mc=0.5时,初始的流向的周期结构经过一定时间后将丧失其稳定性.在计算域为两个扰动周期的情况下,扰动形成的旋涡的后期演化表现为涡的对并.(3)对Mc=0.9的流动,演化图象介于前两者之间.对具有相同对流马赫数的空间问题来说,由于粘性和非线性发展对周期性的破坏,随着空间距离的增长,流动演化与时同问题的差异增大.相比较而言,高对流马赫数的时间问题与空间问题的差异明显大于低Mc下的差异.     

二维可压缩自由剪切层的理论分析

本文从"二维可压缩自由剪切层存在相对稳定的大尺度拟序涡结构,流场应该满足什么样的条件"这样一个反问题出发,应用常微分方程奇点理论对NS方程进行了分析,维持大尺度拟序涡结构的条件是流场内有交替出现的能量不平衡点,鞍点是温度极大值点,结点是温度极小值点;压缩性增加抑制了结点的生成和持续,是导致可压缩自由剪切层厚度增长速度比不可压情况变小的机理;得到的结论可以合理解释现有的物理现象.其次,对于目前研究可压缩自由剪切层广泛应用的对流马赫数进行讨论,认为在建立过程中采用鞍点是压力极大值点假设条件不合理,得出的公式不能反映介质特性,本文采用鞍点是温度极大值点作为假设条件,提出新的求解对流马赫数的公式.根据流动机理分析,提出通过在流场内建立交替出现的鞍点和结点,触发大尺度拟序涡结构,解决超燃冲压发动机混合问题;理论上给出大尺度拟序涡结构频率条件.     

超声速/高超声速双拐点喷管设计

为实现直连式试验台、高温风洞等试验设备的多马赫数运行,提出了双拐点喷管设计方法.喷管分2段设计,第1段共用,采用3次B-Spline函数描述喷管轴线马赫数分布.首先采用特征线方法求解Euler方程,得到无黏的理想喷管型面.其次采用参考温度方法求解边界层位移厚度,对无黏壁面进行修正得到实际壁面.共用段喷管出口的平行均匀流作为第2段喷管设计的初值.为验证设计方法的可行性,设计了中间马赫数为3.0,出口马赫数分别为4.0,4.5和5.0的双拐点喷管,并采用雷诺平均的Navier-Stokes方程对设计的喷管流场进行数值模拟.计算结果表明:喷管出口流场均匀,试验菱形区的马赫数误差小于1.2%.该方法提高了喷管设计精度,保证消波干净,为直连式试验台、高温风洞等设备的多个喷管共用一套动力系统提供了基础.      

超声速变马赫数风洞流场参数线性变化验证

旋转喷管型面使超声速变马赫数风洞在单次运行过程中可连续调节实验区的马赫数,便于研究飞行器的机动过程、进气道起动过程中的气动问题。在控制喷管型面旋转过程中,流场参数能否线性变化是衡量超声速变马赫数风洞性能的一个重要指标。分析变马赫数风洞实验区流场参数的线性变化规律,利用弹簧光顺的动网格技术建立数值仿真模型,验证喷管位于马赫数3.041~3.215 范围所对应的位置时,实验区流场参数是否满足线性变化规律。结果表明：通过对喷管型面旋转的控制实现了风洞实验区流场参数的线性变化,动态计算结果与预期实验区流场参数线性变化规律吻合良好;在不同加速度的流场参数线性变化过程中,各时刻实验区的平均参数与预期参数之间的偏差均小于0.13%。     

边界层厚度对腔体气动声学特性影响数值模拟

为了研究来流边界层厚度对开式腔体气动声学特性的影响,基于分离涡模拟方法,计算了来流马赫数为2.0条件下,不同来流边界层厚度与腔体深度比时,长深比为5.88的腔体流动特性,得到了该腔体声压级的频谱特性.计算结果表明:随着来流边界层厚度增加,形成的剪切层稳定性增强,失稳后上下摆动幅度减少,失稳生成的大尺度涡与超声速主流的相互作用减弱,使得大尺度涡发展到腔体后缘时所具有的平动动能和转动动能降低.大尺度涡撞击腔体后缘在腔体内形成的气动噪声的声压级降低,最大减小幅度达7.5dB.同时各阶模态的频率也发生偏移,偏移值在100Hz左右.基于新的假设重新推导了Rossiter公式,明确了经验常数的物理意义,并以此解释了频率偏移现象.      

有/无尾迹作用下低压涡轮叶栅分离边界层转捩的大涡模拟

采用经过大量算例验证的可压缩大涡模拟求解器对雷诺数为60154、马赫数为0.402的低压涡轮叶栅T106D-EIZ进行了细致模拟,计算了定常来流和周期性尾迹来流两种工况.对计算结果的分析表明:定常来流工况下,叶片吸力面后部出现大尺寸的层流分离泡,分离剪切层的转捩过程受Kelvin-Helmholtz （K-H）不稳定性控制;尾迹来流工况下,由于来流尾迹的周期性扫掠,时均分离泡尺寸变小,叶栅总压损失降低.对相位平均和瞬态流场的分析表明,尾迹引起的逆射流使分离点后移,形成卷起涡结构,逆射流掠过卷起涡的过程中与其发生强烈的相互作用,产生大量气动损失,而后卷起涡破碎,流动转捩为湍流.      

Mach reflection in steady supersonic flow considering wedge boundary-layer correction

Mach reflection in steady supersonic flow is an important phenomenon having received extensive studies, among which simplified theoretical models to predict the size of Mach stem and other flow structure are of particular interest. Past efforts for such models were based on inviscid assumption while in real cases the flow is viscous. Here in this paper we consider the influence of wedge boundary layer on the Mach stem height. This is done by including a simplified boundary layer model into a recently published inviscid model. In this viscous model, the wedge angle and the trailing edge height, which control the Mach stem height, are replaced by their equivalent ones accounting for the displacement effect of the wedge boundary layer, with the boundary layer assumed to be laminar or fully turbulent. This viscous model is shown to compare well with numerical results by computational fluid dynamics and gives a Mach stem height as function of the Reynolds number and Mach number. It is shown that due to the viscous effect, the Mach stem height is increased, through increasing the effective wedge angle.     

计算实际机翼湍流边界层流动的积分法

本文采用一种拟坐标变换掺混方法,计算任意外形机翼的三维边界层流动,亦适用于边界层外流为位流或非位流情况。对C-5A机翼的计算结果表明,该方法具有同微分方法可以比拟的计算精度。并分析了在亚音速情况下粘性对机翼升力和力矩的干扰影响。     

用于可压缩自由剪切流动的湍流混合长度

抓住可压缩流动变密度特性,构造出基于有效涡量的三维von Karman混合长度。湍流模型采用仅依赖湍动能k的单方程KDO(Kinetic Dependent Only)模型,引入新构造的混合长度替换旧尺度得到CKDO模型。为了验证其描述可压缩自由剪切湍流的能力,选择无壁面束缚、密度梯度大和可压缩效应强的自由剪切混合层为算例,其对流马赫数Ma_c=0.8。计算结果表明,KDO模型对混合层的速度分布有着良好的控制和模拟,而经可压缩修正后的CKDO模型与原模型及其他可压缩修正模型相比,所计算的速度分布、主雷诺剪切力和混合层厚度与试验结果更加接近,说明了该混合长度对可压缩混合层这种自由剪切湍流有着良好的刻画能力。     

时间推进法求解叶轮机械内不可压流动

发展了一种用时间推进方法求解叶轮机械内部不可压流动的计算程序。在这个程序中, 使用分布体力方法简化湍流粘性的求解过程, 根据Ｂａｌｄｗ ｉｎ－Ｌｏａｍ ｘ 两层代数湍流模型确定粘性应力的大小。程序中还使用局部时间步长和多重网格数值技术加速收敛。火箭燃料主泵前诱导轮内部流场的实例计算表明, 采用这一计算程序能够正确预测叶轮机械内部不可压流场。     

"回收/调节"方法在混合LES/RANS模拟方法中的应用

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)模拟方法结合三阶加权基本无振荡(WENO)格式对马赫数为2.88的压缩斜坡流动进行了数值模拟,并采用"回收/调节"方法生成入口湍流脉动边界条件.当采用固定入口边界条件时,湍流边界层会缺乏合理的湍流能量,使之抵抗逆压梯度的能力减弱,会严重高估分离...     

高超声速进气道再起动特性分析

1引言未来高超声速飞行器飞行必然要经历低马赫数飞行过程,因此高超声速进气道同样要面临低马赫数下的不起动问题,那么进气道一旦进入不起动,如何才能再起动?再起动的特征是怎样的呢?常规内压式进气道再起动过程中存在迟滞回路现象,高超声速进气道的再起动过程是否也有相似的现     

二维超／高超声速进气道流场数值模拟

对超／高超声速三级压缩进气道流场进行了数值模拟,来流马赫数为4、6,进气道内流动为层流状态,根据二维Navier-Stokes方程,采用二阶精度Roe格式进行离散。按照流场特点,合理地设计网格分布及调整不同黏性范围的熵修正,防止了壁面附近过大的数值耗散,使计算结果更加合理。在进气道模型的各级压缩折转角处,获得了清晰的激波结构,在进气道内部的各种波系的相交、反射和激波诱导的边界层分离等现象都得到合理的描述。计算得到的压力分布,在各级压缩斜板上同简单波理论结果十分接近。用本文方法计算了另一个二级压缩进气道,沿上、下壁面的压力分布与试验比较符合得较好。     

冲压增程弹丸进气道特性分析

采用块结构网格与二阶精度流场分区求解技术，对固体火箭冲压发动机增程弹丸超声速进气道特性进行了深入研究。通过数值模拟得到了对应于不同来流马赫数和攻角情况下，临界工况时，超声速进气道内外粘性流场复杂的波系结构，详细分析了来流马赫数和攻角对进气道性能的影响。结果显示，随着来流马赫数的增大，总压恢复系数显著降低，流量系数增大，同时随着来流攻角的增大，总压恢复系数及流量系数逐渐降低，而流场畸变指数则明显增大。