民用飞机起落架气动噪声数值仿真研究

对某型民用飞机前起落架进行气动噪声进行数值仿真分析。采用分离涡模拟方法,计算起落架周围非定常湍流流场,提取声源信息,利用FW-H方程积分外推法求解模型辐射的声场,获得噪声的频谱特性及远场指向特性曲线。通过频谱对比,研究粗网格和细网格对仿真精度的影响以及不同积分面对噪声预测结果的影响。结果表明,分离涡模拟方法能够较好地模拟流场中存在的大量复杂的涡结构;当选取模型固体表面为积分面时,基于粗网格比基于细网格计算的声压级低约1dB;基于细网格定义空间可穿透面为积分面比定义模型固体表面为积分面计算的声压级高约5～7dB;起落架气动噪声源与其后部有规律的涡脱落相关,并且噪声辐射具有明显的指向特性。     

带偏转角气膜出流发散壁冷却特性的数值分析

运用数值计算方法对带偏转角气膜出流发散孔壁的冷却特性进行了研究,分析了气膜孔偏转角对气膜冷却效率和流量系数的影响。研究结果表明,带偏转角气膜孔喷出的气膜出流具有侧向涡旋特征,增强了气膜射流的横向能力,对于改善发散壁前端的冷却效率具有一定的作用,在高吹风比下更为显著;随着偏转角的增加,发散壁前端气膜孔的流量系数则呈降低趋势。     

改善横向波纹隔热屏综合冷却效率的发散冷却结构优化

为改善横向波纹隔热屏综合冷却效果，采用三维数值模拟和基于支持向量机的代理优化模型，在给定的单位面积冷却空气流量下对发散冷却结构参数进行了优化研究。设计变量选取为气膜孔直径、气膜孔排布的展向间距和流向间距，以面积平均综合冷却效率作为目标函数，通过遗传算法搜索获得了设计变量区间内的优化设计点。在隔热屏单位表面积冷却空气流量G_f = 2.647 kg/（m²?s）的工况下，优化后的d、P和S分别为0.8、4和5 mm。研究结果表明，优化的隔热屏冷却结构应具有较小的气膜孔直径d和展向间距P，以及适中的流向间距S。相对于参考结构，优化后的发散冷却结构能够改善沿流向的气膜覆盖，缩减发散冷却起始段局部高温段，起到增强隔热屏发散冷却综合冷却效率的作用。     

基于ANSYS的齿根应力和齿面变形的合理分析

根据齿轮啮合原理,准确求出复杂的齿廓曲线;采用MATALB—APDL混合建模方法,创建精确的齿轮有限元模型。在不同网格密度等级下,分别对齿轮的齿面变形、齿根应力分布进行分析。对照经验公式的计算结果,得出如下结论：在齿轮数值模拟中,由于有限元软件计算中舍入误差的存在,网格划分合理与否将直接影响仿真程度的高低;网格密度等级对齿面最大变形的仿真结果有较大影响,而对齿根应力影响甚微。     

基于混合网格的直升机旋翼/机身非定常干扰流场数值模拟方法

基于混合网格建立了一个适合于直升机旋翼/机身非定常干扰流场的数值模拟方法。其中,采用混合网格系统生成旋翼/机身干扰流场网格以降低网格生成难度;将二阶MUSCL格式和Roe格式相结合进行空间离散以提高流场的计算精度;同时引入隐式LU-SGS时间格式和OpenMP并行模式以提高干扰流场的计算效率。应用建立的计算方法,对Georgia和Robin旋翼/机身模型进行了数值模拟和分析,与试验结果良好的相关性表明了该方法的有效性。     

陶瓷基复合材料旋转爆震燃烧室壁面冷却研究

为了降低旋转爆震发动机燃烧室壁面温度，设计了陶瓷基复合材料燃烧室主动冷却结构。对燃烧室主动冷却结构的传热特性进行数值模拟，获得主动冷却燃烧室壁面温度响应和温度分布规律。对燃烧室主动冷却结构进行了模型简化，将模拟旋转爆震波获得的不同壁面温度下的热流密度参数加载在冷却模型上，提高了壁面温度模拟的计算效率。结果表明：燃烧室内壁面热流密度随着壁面温度的升高而降低，扩散区的平均热流密度最大；陶瓷基复合材料燃烧室主动冷却结构可以有效降低燃烧室壁面温度，在相同冷却流量下，矩形冷却截面的冷却效果优于圆形冷却截面，可以将燃烧室壁面的温度降到1 200 K以下；燃烧室壁面最高温度在燃烧室中段区域。     

CO2中激波脱体距离的弹道靶实验测量和数值计算

在中国空气动力研究与发展中心超高速所超高速弹道靶进行了CO₂条件下圆球和火星着陆巡视器模型的激波脱体距离测量实验，为数值模拟提供验证依据。实验模型为φ10mm圆球和头部半径12.5mm的着陆巡视器模型。圆球模型的飞行速度为2.122~4.220km/s，靶室压力为2.42~12.30kPa；着陆巡视器模型的飞行速度为2.802km/s，对应靶室压力为1.836kPa。实验数据与采用双温度非平衡模型计算的结果进行了对比。得到以下结论:采用双温度非平衡模型能够较准确地再现模型头部激波脱体距离；根据计算结果推测绕模型流动主要为非平衡流动；需补充更高模型飞行速度（>5km/s）的实验数据，验证CO₂中更高流速状态下双温度非平衡模型的适用性与准确性，并进一步研究多温度模型和不同化学反应动力模型对CO₂下非平衡流数值计算准确性的影响。     

大涡模拟在加力燃烧室冷态流场中的应用

对带V型稳定器模型加力燃烧室二维冷态流动进行了大涡模拟。采用了代数和k方程两种亚网格尺度模型，并用SIMPLE算法和混合差分格式求解离散方程，用壁面函数处理固壁边界条件，得到了模型加力燃烧室二维瞬态流场的模拟结果。通过对不同亚网格尺度模型和不同入口速度下的瞬态流场的比较，揭示了不同亚网格尺度模型和不同入口速度分布对带V型稳定器燃烧室二维冷态瞬态流场的影响。     

二维非结构网格DSMC方法的实现及其应用

研究了非结构网格DSMC方法实现的过程，将Bird位置元方案中的子网格思想引入到非结构网格上来，只存储子网格的总体标识号，做到了利用较少的计算网格提高分子的分辨率和计算精度；提出了将面积元坐标搜索算法与交替数字二叉树搜索算法(ADT)相结合来跟踪模拟分子在网格之间的迁移；使用ADT方法判别分子与物面是否相交，节约了计算时间；最后对过渡流域高超声速绕流进行了数值试验．结果表明这种方法的可行性、高效性及算法的通用性。     

SSTk—ω模型用于冲击射流冷却的可靠性

文章基于FLUENT软件,首先采用SST k-ω湍流模型,对单孔冲击射流在有横流和无横流两种情况下的换热特性进行了数值模拟,获得了被冷却面的Nu分布规律;然后将计算结果与Goldstein&Behbahani公开的实验结果作了比较.通过比较发现,数值计算结果与实验结果吻合良好.     

收敛形进气道数值研究的验证

介绍一种新型的、具有最小喉道面积的三维高超声速进气道 (称之为收敛形进气道 )的数值和实验研究结果。表明使用这种形式的进气道 ,在整个飞行速度范围内可以降低阻力和高超声速发动机表面的热防护要求 ,通过降低外压缩表面的倾斜度和减少进气道及燃烧室壁的面积就可以做到这一点。在采用低维次流动的气体动力设计方法的基础上设计成这种形式的进气道。计算是在无粘气体模型构架内用有限体积法进行的。同时用边界层方程计算出计及粘性的气流特性和进气道特性。数值算法是通过收敛形进气道的有限宽楔形外压缩表面的计算和实验数据来验证的。进行实验研究的马赫数M=2～ 1 0 7,基于模型进气道高度的雷诺数Re=( 1～ 5) × 1 0 6。数值计算与实验结果一致性很好。这些结果也和通常的二维进气道的数据作了比较。     

基于五孔探针的大S弯进气道总压畸变测量与评估

进气道总压畸变的测量与评定是进气道/发动机一体化的重要内容。大S弯进气道具备出色的隐身性能，但其出口流场非常复杂，传统总压测量方法造成的误差显著增大，进而引起总压畸变评估误差，阻碍进/发一体化设计。为了提高大S弯进气道的总压畸变测量与评估的准确性，本文提出了一套基于五孔探针的测量方法。分析测量结果表明:分区拟合方法更能适应大S弯进气道强旋流场的总压数据处理；随着马赫数从0.2增加到0.6，周向总压畸变指数从0.005左右递增到0.09左右，径向总压畸变指数最大不超过0.055，马赫数越大，总压畸变越剧烈，周向总压畸变占据主导；出口截面主要总压畸变区的总压恢复系数最低不到0.85；相比数值计算、总压耙测量，五孔探针测得结果更加全面、合理。     

密切曲面内锥乘波前体进气道设计和试验研究

介绍了密切曲面内锥乘波前体进气道(Osculating Inward turning Cone Waverider Inlet,OICWI)的一体化设计方法,对该型乘波前体进气道的性能进行了数值分析,针对该型一体化乘波前体进气道完成了风洞试验研究。理论设计结果和设计状态无粘模拟结果一致,设计状态下的计算结果表明,前体进气道具有较高的总压恢复、较好出口流场均匀度及较高的流量捕获率。试验研究结果表明,改型一体化前体进气道在马赫数5~7条件下顺利启动,流场波系及压力分布同数值分析结果吻合。     

轴对称进气道外罩跨音速非守恒型全位势流AF3算法

本文针对轴对称进气道外罩跨音速绕流问题的特点,采用了贴体坐标网格;并对这类网格拓扑构造了三种Baker的隐式近似因式分解AF3有限差分迭代算法。通过理论分析与数值实验,找到了既能计算亚临界情形,又能计算超临界情形的AF3格式;并就此格式讨论了适合于进气道绕流特性的AF3迭代格式中间变量的提法,具体提出了格式实施的相应方法,研制了可供分析任意轴对称进气道外罩跨音速特性的计算程序。由于采用了任意曲线坐标系上的轴对称非守恒型全位势方程适合于用不同方法生成的计算网格,因此,本文的工作还有可能进一步推广到计算更复杂的轴对称进气道跨音速绕流问题。计算表明,本方法收敛快,结果好。     

一种可利用部分冲压的直升机进气道流场特性实验研究(英文)

针对可利用部分冲压、带前输出轴直升机进气道结构特点 ,实验研究了在侧滑角从 0～ 1 35°状态下的直升机进气道流场特性 ,分析了沿程静压分布、进气道出口截面流场畸变指数、总压恢复系数等进气道性能参数。研究表明 ,该类进气道在各种侧滑状态下总压恢复系数较高。当侧滑角大于 90°时总压恢复系数随着来流速度的增加而减小 ;进气道内气流分离的区域和出口截面流场畸变指数与侧滑角和来流速度的大小有关。其中在侧滑角小于 90°时 ,进气道出口截面流场品质较好。当侧滑角大于90°时 ,随着来流速度或侧滑角的增加出口流场迅速恶化     

平面扩压叶栅流场犘犐犞与三孔尾迹探针对比测试研究

针对PIV技术在暂冲式高亚声速平面叶栅流场中遇到的示踪粒子投放问题,通过采用高压雾化式粒子发生器以及安装在稳定段前的撒播器,有效地使示踪粒子均匀地与主流混合,并成功开展了某扩压平面叶栅叶片槽道及出口尾迹可视化测量,获得了零迎角、进口马赫数0.2~0.8状态下的二维速度矢量场。为了验证PIV在叶栅流场测试结果的可靠性,在相同工况下,将PIV测量结果分别与数值计算结果和三孔尾迹探针测量结果进行了对比分析。结果表明:采用PIV技术测得的叶栅中截面二维速度矢量场合理地反映了叶片槽道及尾迹的流动结构,与数值模拟结果较为接近;PIV与三孔楔形尾迹探针在叶栅出口尾迹的测量所获得的气流速度和主流区的出口气流角重合性较好;尾迹分离区的出口气流角重合性略差,主要原因是尾迹区气流角超出了探针校准范围,这也说明了PIV测试技术优势。本文提出的PIV测量技术也可用于连续式叶栅风洞中。     

气膜-发散冷却结构冷却效果的实验研究

为了研究气膜-发散组合冷却结构的冷却特征,设计了5种组合冷却结构形式,采用实验的方法对气动参数和几何参数对冷却效率的影响规律开展了研究,结果表明:(1)由于气膜的存在使得发散冷却的起始段冷却效率有很大的提高,极大提高了发散冷却结构的整体冷却效率;(2)随着发散孔纵向间距的增大,发散段整体冷却效率逐渐降低,但对发散冷却起始段的影响不大;(3)发散孔复合角对冷却效率对冷却效率的影响随吹风比的变化而变化,在小吹风比时,45°复合角的冷却效率最高;在吹风比较大时,在发散冷却段的中后部,0°复合角的冷却效率最高。     

离心泵叶轮内含颗粒的幂律流体的密相两相湍流研究

阐述了幂律流体控制方程的特点，首先对单相的幂律流体与密相液固两相流进行了计算，计算结果与实验数据吻合较好。其次对泵离心叶轮内的含颗粒的幂律流体的两相流动进行了计算，从进口到出口，幂律流体的速度是减小的，压力随半径增大逐渐升高，而且压力面上的压力大于吸力面上的压力。吸力面上的颗粒拟温度小于压力面上的颗粒拟温度，而且靠近壁面处的大于叶道中的。最后比较了液固两相流与含颗粒的幂律流体的两相流的流动。流体相的紊动能有一定的差别，而且颗粒拟温度也有较大差别。     

囊腔液体受挤压和抽吸时内流场 PIV 实验研究

主动脉旁反搏装置作为治疗心力衰竭的心脏辅助装置,装置腔内的流场结构对装置的机理分析及优化设计具有重要作用。将主动脉旁反搏装置简化为液囊模型,利用粒子成像测速仪(PIV)对液囊内流场进行测量,选取3种液囊截面,每种截面取3种抽吸频率和3种吸气量,共27种工况,采用整体平均法对实验数据进行处理,得到了平均速度矢量图和湍流切应力图。结果表明:装置腔内流场主要是以旋转的形式存在,在装置进出口附近存在二次流动,而且进出口的位置对内部流场影响较大,同时进出口附近流体的动量交换较为强烈,呈现出较大的雷诺应力,另外进出口附近流体流动能够同步响应抽吸频率,且湍涡中心点随时间在不断地移动。     

冲击多斜孔壁中冲击孔与多斜孔相对开孔位置对换热特性的影响

用数值计算的方法研究了冲击多斜孔壁复合冷却中冲击孔与多斜孔相对开孔位置变化对换热特性的影响。研究发现,随着冲击孔与多斜孔偏距位置变化,两壁狭缝中冲击效果变化显著,多斜孔冷侧壁面冲击换热系数分布规律发生相应变化。偏距位置的变化会使得冲击孔两侧多斜孔抽吸作用不对称,导致位置较近的多斜孔侧涡漩难以形成,削弱了换热效果。当多斜孔在冲击孔两侧对称分布时,总体换热效果较好;而非对称分布时,有利于局部换热增强。