复杂前机身—进气道组合体的流场计算

本文提出了一种能够对复杂前机身-进气道组合体进行流场计算与气动分析的方法。通过求解Euler方程来模拟前机身-进气道组合体在亚、跨、超声速,不同攻角和发动机流量条件下的各类流场。     

直升机旋翼对机身气动干扰的计算

基于旋翼自由尾迹模型和三维机身源面元模型，建立了一个全耦合的旋翼／机身气动干扰迭代计算方法。为正确模拟旋翼桨尖涡与机身表面间的大定常贴近干扰，采用了一个“分析数值匹配法”的“贴近涡／面干扰模型”。应用该计算方法，分别计算了旋翼干扰状态下的机身表面点的非定常压强分布和机身的非定常气动升力，俯仰力矩随前进比的变化，以及干扰前后的机身定常压强分布。计算表明，机身上由于旋翼的干扰引起的非定常气动载荷呈现出与旋翼相同的周期性。然后，对比了本“贴近涡面干扰模型”和先前“截断涡线模型”对干扰计算的影响，表明了前对计算结果的重要改进。     

一种新的自由涡尾迹计算方法

在使用自由涡理论进行旋翼气动计算时,尾迹几何结构的确定是计算旋翼气动载荷以及空间流场计算的关键。发展了一种新的自由涡尾迹结构计算方法,该方法能快速得到收敛的尾迹结构,能比较准确地计算桨叶载荷分布。算例与实验数据的比较验证了该方法的正确性。     

直升机机身外形对气动特性的影响

基于FLUENT流体力学软件,进行了直升机机身流场的数值模拟和分析。描述了网格生成方法,给出了计算迭代步骤和流程图。针对ROBIN机身绕流场进行了计算,并与可得到的实验结果进行了对比,验证了计算结果的可靠性。在此基础上,以NUAA模型机身为算例,改变机身尾部、截面和头部形状对其流场进行了对比计算,给出了表面压强系数分布、机身阻力等计算结果,讨论了直升机机身外形参数对气动特性的影响,得出了一些有意义的结论。     

亚、跨音速导弹气动力的有限差分计算

本文介绍亚、跨音速导弹纵向气动力特性计算的有限差分方法。采用园柱座标系,用变步长的混合差分格式进行线超松弛改进迭代计算。速势方程考虑了x向大扰动,并把精确边界条件嵌入到头部的速势方程中,以适用于计算钝头外形的要求。对尾部收缩的弹体,翼涡沿其平面伸展到收缩部分,以模拟涡面对尾翼和尾段的影响。 所建立的程序适用于计算任意形状旋成体弹身,翼身组合体及翼身尾组合体的表面压力分布和纵向气动力系数。 为鉴定方法的有效性和程序的正确性,对三个不同的气动外形进行了计算,与风洞实验结果比较,基本符合。     

小展弦比扭转机翼三维流场的数值模拟

以高速战斗机的气动特性分析为研究背景,对小展弦比带扭转机翼的亚跨超声速气动特性进行了数值分析,首先对机翼关键数据进行参数化,选取几何扭转角作为设计变量对机翼进行建模,然后利用数值模拟方法,对参数化建模后的机翼在不同扭转角下的气动力和流场结构进行了计算与分析,揭示了机翼不同扭转角对失速特性,机翼展向不同截面的压力系数分布的影响,同时,对最佳扭转角的机翼在不同马赫数下的三维流场进行了数值模拟,验证了小展弦比机翼在亚跨声速下的优越性能。     

不同压力下微秒脉冲表面介质阻挡放电流场实验

采用粒子图像测速（PIV）技术,在不同空气压力条件下,测量了微秒脉冲等离子体气动激励诱导流场的演化过程,分析了不同压力下的流场启动涡、流场结构和壁面射流.根据实验数据计算研究了诱导力随压力变化的空间分布趋势.实验结果表明:常压下和5500Pa压力下产生一个启动涡,19000Pa和11700Pa压力下产生两个启动涡.稳定流场结构随压力减少分别为L型、∽型和V型.压力减小,诱导流场对等离子体气动激励的响应时间减少,射流切向距离变短,距壁面法向距离增加.最大诱导力随压力降低减小,x坐标逐渐向表面介质阻挡放电（SDBD）激励器靠近.      

超声速主流中横向喷流干扰流场数值模拟

本文用二阶迎风 TVD格式求解层流 N- S方程 ,数值模拟了超声速主流与横向喷流相互作用复杂干扰流场 ,得到了合理的涡系结构 ,计算波系结构及沿物面压力分布也与实验结果一致。计算还表明 ,TVD格式不但具有较强的捕捉激波的能力 ,而且用来模拟复杂涡系结构也是有效的     

基于火焰面模型的超声速燃烧混合LES/RANS模拟

为了明晰超燃冲压发动机燃烧室内部燃烧过程的细节,建立了超声速湍流燃烧稳态火焰面亚格子模型,并采用混合LES/RANS方法对氢燃料超燃冲压发动机进行算例验证.控制方程对流项用五阶精度WENO格式离散,时间方向采用二阶Runge-Kutta方法.研究表明:(1)冷流流场中燃料分布与大尺度结构分布相似,说明混合过程受大涡控制;(2)燃烧流场中涡的尺寸明显变大,且仅存在于火焰面上,另外温度分布和主要生成物分布与涡量云图基本相同,说明燃烧过程也由大涡控制;(3)时均计算结果与实验阴影基本符合,速度剖面和温度剖面与实验测量值定性一致,说明本文的数值模拟方法和燃烧模型可以较好地描述和预测超声速流动燃烧过程.     

超声速自由旋涡气动窗口的气动光学特性计算与分析

笔者讨论了超声速自由旋涡气动窗口的气动结构，对设计的超声速自由旋涡气动窗口射流流场及超声速自由旋涡气动窗口的光学性能进行了分析研究。研究了自由涡射流对透射激光产生的气动透镜效应，给出了计算结果。     

马赫数2.5超/亚声流流向涡掺混超燃室冷态实验

结合对应的数值计算,对一种带亚声速预燃室和流向涡掺混器的超声速燃烧模型燃烧室实验台,在其进口马赫数为2.5的来流条件下,进行了冷态流场的实验研究.实验测得其壁面静压分布和激波系结构与流场的CFD计算结果基本一致.实验结果表明,模型燃烧室全流场超声速,达到设计状态.马赫数2.5下的冷态实验数据和CFD计算数据为进行点火实验提供了依据.      

马赫数2.5超燃室中超/亚声流流向涡掺混冷态流场数值研究

结合实验,对一种带亚声速预燃室和流向涡掺混器的超声速燃烧模型燃烧室,在其进口马赫数为2.5的来流条件下,进行了冷态流场的数值研究.计算与实验得到的燃烧室沿程压力分布相当一致,计算与实验得到的激波结构也基本吻合.计算结果表明:首先,在马赫数2.5的来流条件下,亚燃预燃室易于达到启动状态;其次,流向涡掺混器增强超/亚声速流之间的掺混的效果明显,但其掺混深度尚有限;最后,流向涡掺混的超燃室掺混段有着复杂的激波膨胀波波系,波涡干涉和激波附面层干涉结构.     

基于WENO-分段线性格式的直升机流场数值模拟

发展了一套适合于格心格式求解器的基于加权本质无振荡(WENO)-分段线性格式的旋翼/机身气动干扰高精度CFD计算方法。应用该方法对多种旋翼/机身组合模型前飞算例进行了数值模拟,计算得到的机身表面压力系数分布与实验结果吻合良好,说明了该方法对旋翼/机身气动干扰研究的有效性。之后进一步将该方法应用到悬停状态的X3构型复合式高速直升机旋翼/机翼/螺旋桨组合模型的复杂流场模拟中,并与悬停状态孤立旋翼和旋翼/机翼组合模型的流场进行了对比。研究发现:机翼对旋翼下洗流起阻滞作用,引起机翼下方不规则流动,且螺旋桨滑流与旋翼下洗流会相互干扰产生一定的偏折,旋翼下洗流速度更大,螺旋桨滑流会产生明显的向下偏折。      

三维超声速流动的压力反问题

为了进一步探索三维超声速流道的设计方法,采用一种预设壁面压力分布计算壁面型线的思想,并结合双特征线方法提出一种全三维超声速流动压力反问题的求解方法。在三维超声速流场设计中,可直接根据来流条件和壁面压力分布求解壁面的三维坐标,通过空间步进的方式,使得解在一系列解平面上推进,从而使得所设计的型面与预设的壁面压力分布相容。通过Prandtl-Meyer膨胀波的理论解验证了该格式的设计精度。根据预设的压力分布,设计了圆形和椭圆形入口的三维超声速喷管,并将设计方法与数值模拟进行对比验证。验证结果表明:所设计的流场与CFD计算得到的等值线符合得较好,因此基于双特征线的压力反问题求解方法具备三维超声速气动设计的能力,并具有纯三维、高精度、壁面压力分布可控的优势,对未来高超声速气动设计应用将起到重要的支撑作用。     

前飞状态直升机旋翼／机身非定常气动干扰的分析

一个计算旋翼／机身／尾迹间非定常气功干扰的分析方法是建立在二阶升力线／全展自由涡模型和机射面元模型基础之上的。通过迭代机身在桨盘平面、阡迹定位点的诱导速度和旋翼／尾迹在机身表面的诱导速度,形成一个耦合的分析模型。在分析中计入了非定常项。作为算例,对两种孤立机身表面的平均压强系数分布以及旋翼机身组合体中机身表面的非定常压强系数分布进行了计算,其结果与实验值相吻合。     

考虑激波串的超声速燃烧流场分析模型

发展了一种考虑激波串结构的超声速燃烧流场分析模型,以应用于吸气式高超声速飞行器的设计优化过程。模型通过求解耦合有限速率化学反应的刚性常微分控制方程组来描述燃烧室内点火、燃烧等气动热力现象,采用Billig激波串模型模拟燃烧高压前传过程。采用该模型对"等直段+扩张段"、"后向台阶"和"支板喷射"三种典型构型燃烧室流场进行了计算。结果表明:所建模型正确地反映了超声速燃烧室流动中的物理化学过程;与实验数据比较,模型计算出的壁面压力分布比没有考虑激波串结构的模型符合的更好。     

双三角翼及其翼身组合的滚转运动特性比较研究

采用Navier-Stokes方程与滚转运动方程耦合计算方法,比较研究了不同后掠角的双三角翼和翼身组合体的滚转运动特性,分析了机翼前缘后掠角及细长机身对非定常滚转力矩时滞环、动态流场结构和物面瞬时压力分布的影响。研究结果表明：主翼迎风面上的融合涡能量在80°/60°双三角翼上耗散较小,而在76°/40°双三角翼上耗散严重,这是造成两模型滚转力矩稳定性与时滞特性差异的主要因素;机身对气流的扰动作用,大幅增强了滚转力矩的线性分量;机身对气流的上洗作用,增强了边条涡与融合涡吸力及其时滞性,同时加剧了主翼背风面的两涡干扰;大滚转角时机身对横流流动的干扰,使得主翼背风面压力分布的时滞差异显著增加。该研究结果有助于认识后掠角与细长机身影响双三角翼滚转运动特性的物理机理。     

跨音速组合体纵向空气动力的差分计算

 本文建立机翼-机身-平尾-立尾组合体跨音速纵向空气动力的混合差分线松弛迭代计算方案。在物理空间采用直角座标网格,列出组合体的远场边界条件,在流场内部和物面、涡面网格点上,满足速势方程,利用物面边界条件,将速势延拓到机身内部。做了两个算例,得到了收敛或接近收敛的结果,与实验尚符合。对收敛过程和超松弛加快收敛的方法作了一些研究,发现超音速点的松弛参数应取大于1。     

TRIP2.0软件的确认：DPWⅡ复杂组合体的数值模拟

 采用“亚跨超声速计算流体力学软件平台”（TRIP2.0）数值模拟了阻力预测小组（AIAA CFD Drag Prediction Workshop Ⅱ,DPWⅡ）翼/身/架/舱复杂组合体运输机构型,数值模拟采用的多块对接网格、测压和测力的试验结果均来自DPWⅡ,对比计算采用了CFL3D的结果。重点针对DLR-F6翼/身/架/舱复杂组合体构型,详细研究了网格密度和湍流模型对总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致。采用SA一方程和SST两方程模型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布有一定的影响。     

直升机机身、平尾、垂尾的气动干扰

运用面元法模拟直升机机身流场,计算前飞状态的ROBIN(rotor body interaction)模型孤立机身顶部中线的压力系数分布,并与参考数据、CFD计算结果对比,验证了本方法的准确性.采用离散涡系模拟平尾、垂尾和短翼等升力面,以带短翼的UH-60直升机为例,研究了升力面参数变化对机身、平尾、垂尾气动干扰的影响.结果表明:改变升力面的安装角显著改变了气流对其周围的绕流情况,在参数中影响最大;平尾参数尤其是安装角会对垂尾和短翼的压力系数产生较大影响;减小垂尾展长和增大短翼安装角会提升各自的压力系数.