激波运动与反射流场的ENO算法

在ＴＶＤ格式的基础上，以修正限制函数的方法，将通常的ＴＶＤ格式改造为性能更为优越的ＥＮＯ格式，并利用这种有限差分格式求解欧拉方程和薄层Ｎ－Ｓ方程。分别数值模拟了激波管中运动激波与三种楔形物体相互作用流场；给出了激波在不同物体上的反射及波涡干扰过程中非定常流场的变化情况；算出的波系结构在与ＣＡＲＤＣ小型激波管中做的相关实验结果一致。     

一种用于TSTO级间分离CFD计算的网格动态优化技术

两级入轨空天运输系统具有难度较低和经济性好的优点，级间分离是该运输系统设计的一项关键技术。相比前后串联级间分离，两级入轨上下两级间存在严重的激波相互干扰，直接影响分离物体的气动力矩。为了更好地模拟分离过程中两级之间的运动激波多次反射问题，本文结合非结构混合网格分布优化技术和压力比值激波识别技术，建立了用于级间分离运动过程的并行网格动态优化技术。基于计算流体力学软件NNW-FlowStar，完成了两级入轨运输系统并联分离过程的数值模拟。结果表明，采用网格动态优化技术，可实现级间复杂激波的精细捕捉，更精准地获得了强干扰下飞行器的气动特性，显著提高了并联分离过程中运动轨迹和姿态的预测精度。     

带凹腔钝头体第IV类激波干扰特性研究

采用激波风洞实验与数值模拟相结合的方法,对高超声速带凹腔钝头体的第IV类激波干扰非定常振荡特性进行了考察和分析。结果表明,随着入射激波与弓形激波相对位置的不同,第IV类激波干扰既可能出现稳定的波系干扰结构,也可能出现较为剧烈的非定常振荡现象。在实验条件下,观察到了高频前后振荡和低频上下振荡两种不同的振荡模式。数值模拟结果与实验吻合较好,超声速射流与凹腔的耦合作用在振荡过程中起到了关键作用。     

电控双驱动激波管和激波风洞波系和状态计算方法

本文给出了一种计算电控双驱动激波管和激波风洞内波系干扰及实验状态的简单方法,并提出了一种计算运动激波在喷管内运动规律的公式。     

高超声速风洞多体干扰与分离试验技术

在FL-31风洞中进行了某高超声速飞行器的多体干扰与分离试验技术研究,成功建立了多体干扰与分离试验技术。试验模型是某典型构型的可重复使用航天飞行器,由助推器以及再入体两部分组成。利用风洞上下投放机构实现两模型间的相对运动,采用两台天平对模型的气动力进行测量,同时利用纹影仪记录模型分离过程中的激波干扰情况。结果表明：试验系统设计合理,能准确模拟物体间分离过程,并能精确测量多体干扰的气动力特性,激波干扰清晰可见。     

在激波管中实现“固定模型”双波斜相互作用的研究

本文在电控双驱动激波管实验技术的基础上,探讨了国外双波斜相互作用的各种实验方案和装置,提出并初步实现了在电控双驱动激波管中进行平面运动激波与固定模型上的头激波斜相互作用的实验方案。在第一激波名义马赫数M_s1=3.5,4.0,4.5下,对双驱动激波管进行了波系运动状态的计算以及双波斜相互作用流场干扰角(或称相互作用角)等的计算。荻得了在这些状态下的双波斜相互作用的纹影照片。实验结果表明,实验所测得的平面运动激波与尖劈模型的干扰角与计算结果基本一致。     

减小翼型激波阻力的鼓包流动控制技术

针对2020年使用的N+2代民用飞机的翼身融合(BWB)布局发展需要,以减小激波阻力为目标,采用计算流体力学(CFD)方法,开展弱化激波、减小激波阻力的鼓包流动控制技术研究.提出了λ形激波结构“强干扰”和等熵压缩“弱干扰”两种鼓包激波减阻流动控制原理,给出了两种鼓包基本形状设计方法和工程应用的可行性分析,指出λ形激波结...     

连续双扫掠激波/湍流边界层干扰流动特性研究

为了探究连续双扫掠激波/湍流边界层干扰的流动特性,采用仿真方法对一双尖鳍/平板物理模型进行研究。结果表明:双扫掠激波/湍流边界层干扰形成的两个干扰区存在明显的相干现象,虽然第一道扫掠激波/边界层干扰流动仍具有典型的准锥形相似特性,但受其干扰所形成的非均匀流的影响,第二道扫掠激波/边界层干扰却不再具有准锥形相似特性,同时第二个干扰区将影响其上游临近气流的运动甚至影响第一个干扰区的再附线和分离线等结构。两个干扰区形成各自的λ波结构,并且沿着流向两个干扰区内的激波结构相互汇聚,最终合并为单个更强的λ波结构;不仅如此,两个干扰区内还形成了复杂的旋涡结构,包括一级主旋涡和二级主旋涡,这些旋涡向下游运动,最终融合成一个尺度更大的锥形主旋涡。     

斜激波与弓形激波相互作用的彩虹纹影实验研究

斜激波与弓形激波相互作用形成的波系结构可以分为六类,其中IV型激波干扰产生了超声速"喷流",产生了极高的热载和压载。基于彩虹纹影显示技术,实验获得了六类激波干扰彩虹纹影照片,分析了各类激波干扰的波系结构。通过数值计算,获得了有无斜激波干扰的条件下,壁面附近温度与压强的分布,讨论了IV型激波干扰产生极高热载和压载的原因。研究表明,激波干扰发生以后,钝头体表明压载与热载急剧增加,尤其是IV型激波干扰,形成的超声速"喷流"射向钝头体壁面的超声速射流,会导致峰值压力剧烈升高,其数量相当于来流静压的近40倍,驻点附近的热流升高将近8倍。     

激波与旋涡对干扰研究:激波动力学特性及声波产生机理

采用五阶WENO格式,通过数值求解二维非定常Navier-Stokes方程,系统模拟了激波与旋涡对的干扰过程.我们研究了强度为1.05和1.2两个典型的激波同不同强度的旋涡对的干扰过程.根据激波动力学特性,我们把激波与掠过型旋涡对干扰分成四类,而碰撞型干扰分成七种不同的类型.研究表明,激波与旋涡对干扰过程中,声波产生有三种机制,激波与涡对的干扰,涡对之间的耦合以及激波与声波的干扰.     

弹体穿越平面激波三维非定常流场的数值模拟

用数值方法求解三维非定常可压缩Navier-Stokes方程，计算了复杂弹体在相同马赫数、不同迎角下的流场和气动力特征，计算结果与实验结果吻合，证明本文采用的方法能够准确计算此类复杂流场；然后采用该方法模拟了弹体穿越正激波的非定常过程，得到了弹体周围流场及气动力特征随穿越过程的非定常变化，阻力特征穿越后急剧下降。     

不同攻角条件下偏转头弹箭流场与气动特性

基于DES（detached eddy simulation,分离涡模拟）方法及AUSM（advection upstream splitting method,对流迎风分裂方法）格式,对不同偏角的偏转头弹箭在非零攻角下的超声速流场和气动特性进行了数值模拟.结果表明:在非零攻角情况下,与常规弹箭流场相比,偏转头弹箭头部迎风面激波更强;而当有效入射角相等时,即使弹箭的攻角不一样,其头部流场结构仍相同,且尾翼流场基本不受头部偏角影响.另外,由气动力系数变化曲线可知,偏转头弹箭具有良好的纵向静稳定性,且升力系数随攻角的增长率大于阻力系数,即升阻比大于常规弹箭,因此适当增大攻角可以提高偏转头弹箭的气动特性.      

弹头外形不同尺度模型对二维射流流场的影响

超声速射流装置能很好地模拟榴弹引信发电机飞行时大气环境的工况，并成为地面考核和验收必不可少的设备。从二维轴对称N-S方程出发，利用二阶精度有限体积法和S-A湍流模型，对超声速射流装置的自由射流流场以及不同模型尺寸对超声速流场结构的影响进行了模拟计算。给出了激波与射流边界之间复杂干扰流场结构图，结果证明不同模型的尺寸对超声速射流流场有较大的影响。     

旋转弹体马格努斯效应数值模拟方法研究

旋转的弹体通常会存在马格努斯效应,严重影响飞行稳定性和弹道轨迹。提出基于旋转壁面法的多对称面弹体旋转效应数值模拟方法,无需借助动网格和多参考系模型,仅需要将旋转速度附加壁面,即可使用定常方法模拟马格努斯效应;同时,使用SOCBT标准弹体模型对该方法进行验证分析。结果表明:计算值与试验值吻合较好;由于旋转引起的附面层堆积以及大攻角分离效应、激波的干扰,计算时附面层的网格节点要足够密,并且应该针对不同的攻角选择精度与效率兼顾的计算模型。     

变截面管道中激波诱导的两相流动

本文研究扩张型与收缩型变截面管道中激波诱导的非定常两相流动。对于稀相气固悬浮体采用了双连续介质模型,对于两相流动按照准一维近似处理。控制方程则是利用算子分裂技术和二阶GRP方法数值求解。文中给出了波后气固两相的流动结构,并讨论了粒子对两相流场的影响。     

含侧喷流弹丸三维绕流干扰流场数值分析

从全Navier-Stokes方程出发，应用三维隐式有限体积TVD格式，在超声速和零攻角条件下，分别对无侧喷流和含侧喷流的低阻远程弹丸三维粘性绕流流场进行了数值模拟，为避免分区计算而用代数方法生成了弹丸绕流场O型网络，计算结果得到了清晰的流场波系结构，并从流场波系结构、弹体和弹底表面压力分布等角度将两者进行了对比与分析，表明侧喷流不仅提供增程推力，而且改变了弹丸阻力特性，其结果可用于研究弹丸的火箭增程技术。     

超燃冲压发动机内外流场三维燃烧数值模拟

以某型煤油燃料超燃冲压发动机整机模型为研究对象,对发动机外流场的流动与发动机内的流动燃烧过程开展了耦合三维数值模拟研究.燃烧室壁面静压计算结果与实验数据吻合良好,平均相对误差为4.9%,验证了数值方法的有效性.结果表明:内流场的流动燃烧状况既受到外流场捕获空气流量的限制,又会反过来影响着外流场的激波结构;泄压孔对保持发动机中各部件流量匹配起到关键作用;由于进气道内激波与边界层的相互作用以及泄压孔诱发的斜激波的影响,流场在第1处燃油喷射孔附近呈现出燃烧流场不对称特征.      

非均匀来流下三维激波反问题的微元密切轴对称解法

传统的密切轴对称理论被广泛应用于均匀来流下的三维密切曲面激波反设计，为解决非均匀来流条件下的三维曲面激波反问题，提出了一种微元密切轴对称流场（MOA）求解方法。该方法沿激波面的周向和流向构建一系列微元密切面，在每个微元面内进行三维向二维流动的等效转换，从而突破了传统密切方法中不能有横向波后流动的限制。利用该方法编写设计程序，分别基于带攻角来流条件和外锥型流来流条件重构了标准内锥曲面激波，并与数值仿真结果进行了比较。结果表明，非均匀来流下激波曲面的三维形状均与预设形状完全一致，实现了非均匀来流下曲面激波形状可控。MOA方法在吸气式高超声速推进领域中前体/进气道一体化设计方面有重要应用前景。     

密切内锥乘波前体进气道一体化设计和性能分析

采用特征线方法设计了具有直线初始激波、内收缩段消除激波反射、出口参数均匀可控的基准内锥流场。基于密切内锥(Osculating Inward turning Cone,OIC)乘波体设计方法,发展了密切内锥乘波前体进气道(Os-culating Inward turning Cone Waverider Inlet,OICWI)一体化设计技术。基于基准内锥流场和前体进气道一体化设计技术,设计了密切内锥乘波前体进气道。采用数值方法对设计的密切内锥乘波前体进气道进行了计算分析,结果表明无粘流场结构和基准内锥流场吻合,无粘模拟结果和理论设计结果吻合。粘性数值模拟结果显示一体化进气道具有较高的流量捕获率及总压恢复特性,进气道出口流场分布均匀。     

动态微涡流发生器对激波/边界层干扰的影响研究

本文基于OPENFOAM数值仿真平台,采用动态网格技术和湍流离散涡(DES)模型,研究了微涡流发生器以一定速度向下游移动时,激波/边界层干扰(SWBLI)流场特性的变化,重点关注干扰区域内的流向和展向的流场特性。来流马赫数为4,微涡流发生器向下游移动速度为0m/s,20m/s和40m/s。研究表明：当MVG向下游移动时,SWBLI区域的“弓”形高压区会演化成“双弓”形；入射激波形成高压区的压力明显降低,同时,入射激波和反射激波形成高压区的峰值位置均会向下游移动；流场下游 “双圆弧”状高压区的高度逐渐降低；SWBLI区域边界层的高度逐渐降低,同时边界层底部的速度也有所降低；随着MVG移动速度的增加,对SWBLI流场的控制效果更加明显；动态MVG对流场的控制是通过尾迹涡和波系结构实现的。