含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

动力对全机水滴收集率的影响计算

针对飞机在飞行中遭遇过冷水滴撞击并结冰现象,建立了适合于发动机带动力情况下结冰过程水滴收集率计算的三维数值方法和计算程序.其基本思路为:采用多块技术与SIMPLE方法计算空气流场,以流场分布的计算结果为基础,求解水滴相的控制方程,进而获得物体表面的水滴收集率.空气相控制方程和水滴相控制方程均写成典型输运方程的形式,采用一致的有限体积法离散求解,方便了计算程序的编制.对某型运输机巡航构型有/无动力条件的水滴收集率进行了比较计算,获得了不同直径水滴在飞机表面的撞击特征以及水滴收集率在飞机机翼、平尾、垂尾和发动机进气道唇口上的分布规律.研究表明:(1)发动机是否带动力对机翼、平尾、垂尾的水滴收集率基本无影响;(2)飞机带动力主要影响发动机进气道唇口处的水滴收集率,带动力后唇口的收集率比无动力情况高,水滴撞击范围增大,在进行防除冰研究和设计时需引起重视.     

火花放电合成射流与超声速来流相互干扰特性数值模拟研究

为了指导火花放电式合成射流激励器在超声速流动控制中的应用,数值模拟研究了火花放电合成射流与超声速来流的相互干扰特性。研究表明火花放电式合成射流在超声速流场中产生强烈扰动,产生较强的激波结构;随着射流的喷出,激励器上游分离区和流场中激波呈先增强后减弱的趋势,激波由弓形激波逐渐弱化为斜激波,并且随着放电能量的增加射流与主流的动量通量比不断增大,射流的干扰和控制能力显著增强。由于超声速流的较大惯性及其对腔内气体的引射作用,激励器的腔体回填速率大幅下降、回填时间明显增长,使得激励器的工作频率受到很大限制。     

基于星载可见光相机的空间碎片探测

空间碎片在轨识别与参数辨识为空间轨道预警、航天器规避空间碎片提供了重要依据。文章在轨道相对动力学的基础上,模拟了空间碎片在光学探测过程中的拖尾成像特性;随后,采用 Hough变换对空间碎片尾迹特征.r进行提取与识别,获取碎片方位角信息,通过匹配多帧图像,获取碎片的方位角速度信息。同时,结合激光测距仪的测距信息,获取空间碎片的位置和速度信息;通过数值仿真验证,该方法能够实现对空间碎片探测、识别和定位,因而具有一定的工程应用价值。     

某航机燃烧室煤油改柴油燃烧特性研究

为了给某型航空发动机改为地面用柴油型燃气轮机的设计提供重要的技术支持,本文借助数值计算的方法,采用FLUENT稳态压力求解器、P1辐射模型和涡耗散破碎(EDU)燃烧模型对某航空发动机燃烧室在巡航工况和最大工况下煤油与柴油两种燃料的燃烧特性进行了计算及对比研究。得到了该燃烧室使用航空煤油(RP-3)和0号柴油的热态流场、空气流量分配、温度场、出口温度分布、污染物排放及头部燃油蒸发量。研究结果表明：当该燃烧室的燃料由航空煤油改为0号柴油后,燃烧室的热态温度场分布基本一致,流量分配最大差异在0.45%之内;燃烧效率降低约4.3%和NO、Soot排放量相当;出口温度分布和总压损失差异分别在1%和4.1%之内。     

某大型地效飞行器及其驮运巡航状态气动特性数值模拟

研究地效飞行器及其驮运负载飞行器巡航状态的气动特性,对大型地效飞行器的研发具有重要意义。采用求解N-S方程的有限体积法,分别针对地效飞行器巡航状态及其驮运负载飞行器巡航状态进行数值模拟分析。结果表明:综合考虑升阻比与阻力的大小,地效飞行器巡航状态机翼的最佳安装角为5.5°;综合考虑升阻比与巡航飞行的安全性、操纵性,地效飞行器巡航状态的最佳飞行高度为8 m;综合考虑巡航经济性以及从巡航状态过渡到分离状态的安全性与可靠性,两机驮运巡航状态下负载飞行器的最佳迎角为3°。     

高负荷压气机静叶根部叶型气动性能实验研究

为研究高负荷压气机静叶根部流动状态,抽取该静叶根部叶型并模化成为平面叶栅进行吹风实验,冲角变化范围为-6°~+6°、进口马赫数变化范围为0.5~0.7。对叶栅出口截面气动参数进行了详细的测量,结果表明:非正冲角时叶展中部节距平均总压损失系数小于等于0.036,尾迹宽度和总压损失峰值随马赫数变化均不明显,正冲角时尾迹宽度和总压损失峰值急剧增加;端壁处的压力梯度随马赫数和冲角的增加而增加;0°冲角下随马赫数的增大,出口叶展中部截面二次流动能系数增加,二次流动得到加强,高能量损失区域增大并且尾迹略有变宽。     

高度欠膨胀声速射流的自相似特性

构建了高压气体射流模型和高分辨率的计算网格,对喷压比（NPR）为5.60、流动雷诺数Re为105量级的高度欠膨胀射流进行了三维大涡模拟（LES）计算.讨论了时均的射流近场结构,发现大涡模拟成功捕捉到了高度欠膨胀射流近场的典型波系结构,并与文献结果吻合较好.研究着重定量考察流向速度和氮气质量分数的分布规律,以揭示高度欠膨胀射流的自相似特性.结果表明:在射流核心区之后,不同流向位置上流向速度和氮气质量分数沿径向的分布呈现出与亚声速射流类似的自相似特性.但射流流场开始呈现自相似的位置仍然为超声速,这是高度欠膨胀射流自相似特性的独有特征.提出了高斯拟合公式,在靠近和远离射流中心线的区域内均可较好地表征射流的自相似特性.此外,还考察了射流剪切层的发展特征,比较了计算得到的射流剪切层增长率与先前实验测量结果的差异,并分析了其中可能的原因.      

扰流片式推力矢量喷管气动特性数值模拟研究

为了满足工程应用中扰流片式推力矢量喷管控制率建模以及优化设计的需求,本文选取了三种几何外形的扰流片,并通过数值模拟手段,研究了扰流片几何形状对于轴对称喷管推力矢量气动特性的影响规律,提出了减小扰流片推力损失的设计方法。数值模拟结果表明,推力矢量角与推力损失系数都随着扰流片插入高度的增加而增加;对于矩形扰流片,可以通过增加扰流片宽度的方式减小推力损失,对于扇形扰流片,可以通过减小上圆弧圆心与扇形顶点距离的方式减小推力损失;在插入高度及面积一定时,对比不同形状的扰流片,弧顶矩形扰流片的推力矢量角及推力损失系数均为最大,圆形扰流片均为最小。     

第一届航空CFD可信度研讨会总结

总结了第一届航空可信度研讨会(AeCW-1)的数值模拟结果。AeCW-1会议由中国空气动力学会计算空气动力学专委会主办,并在第一届空气动力学大会上做了专题研讨。本次研讨会吸引了14家单位参与,共提交了26组计算结果。活动的主要目的是针对典型运输机巡航构型,评估国内CFD当前技术状态,探索CFD进一步发展的研究方向,促进CFD验证和确认工作的开展,为大型客机的研制提供技术支撑。AeCW-1选择的研究模型为中国空气动力研究与发展中心(CARDC)设计的单通道运输机模型(CHN-T1),CHN-T1包括机身/机翼/平尾/立尾等部件,分别在CRADC的2.4m×2.4m风洞(FL-26)、欧洲的2.4m×2.0m风洞中(European Transonic Wind Tunnel,ETW)开展了风洞试验。AeCW-1组委会提供了系列结构网格和非结构网格,研究工况主要包括网格收敛性研究、抖振特性研究、雷诺数影响研究,风洞模型支撑装置和静气动弹性变形对数值模拟结果的影响是重要的研究内容。本文对提交组委会的数值模拟结果进行统计分析,并与FL-26及部分ETW的风洞试验结果进行了对比,对未来的研究工作提出了意见和建议。