翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构，层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟，研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型，基于结构化拼接网格，对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性，雷诺数对翼型气动特性具有显著影响。随雷诺数降低，时均分离泡外形增大、位置后移，平均阻力系数增大，特别是在较低雷诺数下，翼型升阻力系数随时间出现振荡现象。进一步研究表明，造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低，流动黏性增大，导致分离剪切层速度梯度减小，流动发生转捩及再附位置后移，直至翼型表面不再发生转捩和再附。     

简化微小液滴在超声速气流横向喷射中的运动探索

从研究超声速气流中简化液滴(刚性小球)的气动力着手,比较分析使用CFD方法与使用两相流理论中已有的颗粒阻力系数模型得到的简化液滴气动力结果,得出Charles B.Henderson给出的阻力系数经验关系式适用于计算简化液滴在超声速气流中的运动。进一步,对不同直径简化液滴的运动开展工程方法的计算。在来流Ma=2.7的二维平板超声速流场中选取一个截面,作为气相流场,结果显示:(1)简化液滴与主气流存在相对超声速运动,当简化液滴直径dk≤0.12mm时,纵向相对超声速运动区域约为0.15m~0.4m,当dk0.12mm时,作用区域明显增大;(2)以10m/s喷射出的简化液滴,其横向穿透深度与纵向运动距离比约为0.004m/1m~0.021m/1m;(3)以100m/s的速度喷射出的30~120μm直径简化液滴,其横向穿透深度与纵向运动距离比约为0.02m/1m~0.055m/1m,实际过程中,小尺寸简化液滴的汽化很快,其穿透深度很小。     

锯齿形转捩片触发高超声速进气道边界层转捩的大涡模拟

为促进锯齿形转捩片在高超声速进气道中的应用，以地面风洞条件下的二元进气道为研究对象，采用高精度大涡模拟方法对锯齿形转捩片在三级压缩楔面上触发的边界层转捩现象开展了研究。数值方法基于隐式亚格子模型，空间离散采用高精度通量限制型紧致格式，时间推进采用显式Runge-Kutta方法。数值模拟清晰捕捉到了边界层转捩的空间发展演化过程，并获得了统计平均流场以及流场脉动特征。数值模拟结果表明转捩片能够有效触发进气道压缩面边界层转捩；通过与等熵压缩面及单楔面数值模拟结果的对比分析，获得了转捩片触发边界层转捩的内在机理，为后续研究工作奠定了基础。     

超声速气流中简化微小液滴横向喷射的汽化过程探索

为研究超声速气流中简化液滴的汽化过程问题,本文分析了两相流计算中已有的两相传热模型,并对简化液滴绕流开展数值计算.在来流Ma≤0.6的条件下,数值计算得到的简化液滴-气流之间的传热速率与已有模型得到的结果相一致,而在来流Ma≥0.9的条件下,数值计算得到的简化液滴-气流之间的传热速率与已有模型得到的结果存在很大偏差.由此建立了考虑简化液滴与气流相对超声速相互作用的两相传热模型.进一步,采用Charles B.Henderson阻力系数关系式与新建立的传热模型,对不同直径简化液滴的运动与汽化开展工程计算.在来流2.7Ma的二维平板超声速流场中选取一个截面,作为气相流场,结果显示,(1)简化液滴与主气流存在相对超声速作用.当简化液滴直径dk≤0.12mm时,作用区域约为0.1m～0.4m,当dk>0.12mm时,作用区域明显增大,(2)简化液滴的穿透尺度不超过0.011m/m(深度/长度),时间尺度约为0.28ms～3ms,(3)简化液滴完成汽化的空间尺度约为0.1m(dk>=0.03mm)、0.45m(dk>=0.05mm)与1.24m(dk>=0.075mm),而当dk>0.09mm时,简化液滴完成汽化的空间尺度则大于1.9m.使用考虑简化液滴与气流相对超声速相互作用的两相传热模型与使用传统的传热模型对简化液滴的运动轨迹没有影响,而对简化液滴的汽化过程有较大影响.     

高超声速钝双楔绕流流动转捩与分离流动的壁温影响

为研究壁温对吸气式高超声速飞行器进气道转捩流动的影响,选取钝双楔这一典型外形,基于德国Aachen工业大学Thomas与Herbert所开展的双楔高超声速风洞试验,分析了一些已有的计算流体力学(CFD)研究内容,并结合本文不同方法的CFD数值模拟结果,讨论了不同壁面温度对该双楔模型高超声速绕流流动转捩与流动分离的影响。对于双楔模型,流动分离一般发生在拐角附近,由于流动分离旋涡的剪切作用会诱发流动转捩,转捩又会改变流动分离强度、分离涡尺寸,若分离流动存在非定常特征则将导致非定常旋涡运动与流动转捩的复杂相互作用。通过比较已有文献的CFD数值模拟结果与本文计算结果,表明只有按照转捩思路开展的数值模拟才能够反映该风洞试验情况。计算结果与试验数据的比较显示,文献中按照第一压缩面层流与第二压缩面湍流状态计算得到的结果能够在一定程度上与风洞试验数据相符,本文使用MUSCL格式、剪切应力输运(SST)湍流模型与γ-Reθ关联转捩模型这种计算方法,得到的结果与试验数据符合较好,正确地反映了风洞试验情况。分析还表明,在分离流动之前的区域,随着壁面温度的升高,壁面热流会下降,近壁区域黏性系数变大,边界层内速度剖面不饱满,速度边界层较厚,厚的速度边界层容易发生流动分离现象。