平面波压力扰动下压气机逐级动态响应分析

应用轴流压缩系统的逐级动态响应模型对平面波型压力扰动下的多级压气机瞬态响应及气动稳定性进行了较详细的数值分析。对－ ８级跨声压气机的数值模拟结果表明, 基于流动阻塞机理建立的压缩系统稳定性判别准则是合理和可靠的。结果还进一步建立进口扰动频率与压气机瞬态响应间的关联关系     

增升装置微型涡流发生器数值模拟方法研究

 针对数值模拟带微型涡流发生器(VG)的增升装置绕流过程中出现的典型问题,研究相应解决途径.引入面搭接网格技术,避免传统结构网格方法造成的网格数量过大问题;在微型涡流发生器网格上采用特殊边界条件处理以分别对应有/无涡流发生器的状态,排除网格变动对计算结果造成的干扰;研发出了面搭接网格生成技术与特殊边界处理方法结合的求解加装微型涡流发生器的增升装置流动数值模拟方法.研究结果表明:该方法在处理增升装置这类复杂外形、复杂流动问题以及加装涡流发生器的流动控制问题时,具有良好的数值计算精度,能够满足研究需要,为开展涡流发生器精细设计和参数优化的数值模拟研究提供了可靠的技术途径.     

提升机构微型燃气发生器的探索研究

介绍了一种宇航用的利用产生燃气做功并通过提升机构完成提升动作的微型燃气发生器的设计。建立了描述燃气发生器内部流场的数学物理方程。通过数值方法模拟了燃气发生器内部的燃烧和流动以及提升机构的运动过程。由模拟结果及理论分析来指导药柱的设计, 使得燃气发生器能满足高性能和高可靠性的要求     

等离子发生器燃烧流场数值模拟

采用漩涡破碎（EBU）燃烧模型、k-ε双方程湍流模型及SIMPLEC算法,对等离子发生器内部的燃烧流场进行了数值模拟,得到了关于等离子燃烧流场相关特性参数分布。该数值模拟对等离子发生器的研究具有指导作用。     

高超声速进气道强制转捩流动的大涡模拟

吸气式高超声速飞行器常在进气道边界层内布置粗糙颗粒或涡流发生器强制流动转捩为湍流以确保发动机正常启动。为了清晰认识强制转捩过程,采用隐式大涡模拟方法,对强制转捩问题开展了数值模拟研究。针对钻石形和斜坡形涡流发生器,计算得到涡流发生器诱导的流动结构,显示出强制转捩流动由涡流发生器产生的反向旋转流向涡对主导。扰动沿流向增长和发展,导致流向涡对以偶模式或奇模式失稳,偶模式失稳产生对称形式的涡对破碎,而奇模式失稳则导致非对称（弯曲）形式的涡对破碎。流向涡对破碎后产生一系列发卡涡并最终促使边界层转捩为湍流。最后就计算网格和数值耗散对隐式大涡模拟结果的影响以及计算的收敛性进行了讨论。     

平面波导器件封装工艺软件的研究与开发

平面波导器件的耦合封装是制约其发展的一个瓶颈,通过研究平面波导器件的封装工艺特点,分析了封装过程对控制系统的要求.基于Visual Basic.Net开发了封装工艺软件,研究了手动调整、视觉引导粗对准、自动化对准以及点胶固化四个模块的设计方法,并对编程过程中的关键技术进行了讨论.     

旋转开槽圆筒式阵风发生器流场特性数值模拟

在OpenFOAM开源代码框架下开发求解器,对旋转开槽圆筒式阵风发生器的动态等效升力系数进行数值模拟,与试验结果吻合较好,验证了求解器的正确性。在此基础上,对该型阵风发生器的下游流场扰动风速进行了数值模拟和流场分析,分析了开槽圆筒旋转方向、下游流场位置、减缩频率等因素对下游流场扰动阵风的影响规律。结果表明,开槽圆筒的旋转方向对中心线上下两侧流场产生不同影响,对中心线附近的流场干扰较大。阵风幅值的变化与减缩频率有关,当减缩频率值在0.2～0.25之间时,产生阵风的正弦特性最佳。该数值模拟的结果对阵风发生器试验装置的设计具有一定的指导意义。     

横截面形状对超声速磁流体发生器性能的影响

通过三维数值模拟研究了不同横截面形状的磁流体(MHD)发生器的热、电磁流动特性,发生器入口采用超声速来流条件,旨在应用于磁流体旁路超燃冲压发动机推进系统。数值模拟采用低磁雷诺数下磁流体五方程模型,通过熵条件格式和超松弛迭代(SOR)法联合算法分别求解Navier-Stokes方程组和电势方程。研究表明在电极壁面上电场、电流以及洛伦兹力等电磁场参数呈周期性变化,由于电极端点效应,电磁场参数在电极端点出现周期性的极值,而在管道中心电极端点效应对电流影响并不大,沿着流向电流保持连续变化。对不同管道横截面形状的发生器的数值模拟表明:由于涡电流和二次流效应的影响,形状参数较小的发生器焦耳热较为严重,二次流较弱;而形状参数较大的发生器焦耳热较小,二次流增强。因此,当形状参数α=0.8时可获得较好的焓提取率和电效率。     

数值涡流发生器技术在平板绕流数值模拟中的应用研究

本文将数值源项模型应用于叶片式涡流发生器的数值模拟中,以叶片在当地流场中产生的升力作为源项,该源项由叶片的坐标位置、安装方向、高度等参数根据Plandd升力理论确定。数值模拟中,叶片位于平板上,其高度约为当地边界层厚度的1／6。在相同计算条件下,对比分析采用该源项模型数值模拟结果和相关资料的对应数据,证明采用本文源项模型可以得到与资料采用全网格以及采用其涡流发生器模型基本一致的诱导速度场。     

环形爆震波聚焦起爆数值模拟

对三种带不同结构形式的环形预爆管的爆震发生器进行了数值模拟,来研究环形爆震波聚焦起爆现象及其气动特性.数值计算采用多组分理想气体详细的化学反应机理、二维轴对称非定常流动的NavierStokes方程来模拟化学动力学和流体动力学过程.研究发现用低的点火能量对环形预爆管中的燃料和氧化剂点火产生层流火焰,层流火焰在狭窄管壁的...     

动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

Hysteresis processes in the regular reflection↔Mach reflection transition in steady flows

Ernst Mach recorded experimentally, in the late 1870s, two different shock-wave reflection configurations and laid the foundations for one of the most exciting and active research field in an area that is generally known as Shock Wave Reflection Phenomena. The first wave reflection, a two-shock wave configuration, is known nowadays as regular reflection, RR, and the second wave reflection, a three-shock wave configuration, was named after Ernst Mach and is called nowadays Mach reflection, MR.A monograph entitled Shock Wave Reflection Phenomena, which was published by Ben-Dor in 1990, summarized the state-of-the-art of the reflection phenomena of shock waves in steady, pseudo-steady and unsteady flows.Intensive analytical, experimental and numerical investigations in the last decade, which were led mainly by Ben-Dor's research group and his collaboration with Chpoun's, Zeitoun's and Ivanov's research groups, shattered the state-of-the-knowledge, as it was presented in Ben-Dor (Shock Wave Reflection Phenomena, Springer, New York, 1991), for the case of steady flows. Skews's and Hornung's research groups joined in later and also contributed to the establishment of the new state-of-the-knowledge of the reflection of shock waves in steady flows.The new state-of-the-knowledge will be presented in this review. Specifically, the hysteresis phenomenon in the RR↔MR transition process, which until the early 1990s was believed not to exist, will be presented and described in detail, in a variety of experimental set-ups and geometries.Analytical, experimental and numerical investigations of the various hysteresis processes will be presented.     

进口斜激波、膨胀波干扰下等直隔离段内的激波串特性

 超燃冲压发动机的隔离段在实际工作中会受到进气道唇罩激波及肩部膨胀扇的显著干扰,本文针对这一特定问题进行了专门研究。提出了唇罩入射激波及肩部膨胀扇的模拟方法,并利用德国Achen的风洞试验对其进行了检验,而后以此研究了入射激波及肩部膨胀扇干扰下隔离段内激波串的基本形态,并分析了出口反压和激波入射位置的影响。仿真结果表明：当激波串在隔离段内不断前移时,受唇罩入射激波及其反射激波的干扰,其高速核心区交替地偏向上下壁面;与无激波入射的情况相比,此时激波串的耐反压能力显著降低,且入射点位置越高,降低幅度越大,管道内的沿程静压分布规律与Waltrup经验公式偏离程度也越来越大。该文结果可为进气道/隔离段的一体化设计提供依据。     

关于平面激波在凹两段楔面上的传播

众所周知,当一平面激波同任意形状的障碍物相互作用时,将产生十分复杂的波系,其中激波在曲面壁上的绕射和反射,将产生非定常波系。关于平面激波在凹圆柱壁面上的传播,我们可以用一连串、长度为无限小的楔面所构成的多边形来近似和代替凹圆柱面,因此分析和研究激波在两段楔面上反射和绕射的目的在于得出激波在圆柱壁面以及其它较复杂形状壁面上的传播规律。本文采用马赫-曾德干涉仪在激波管上进行光测模型实验。根据所拍摄的流场干涉图像,以分析其波系的演变过程及传播规律,得到有关流场的密度分布和压力分布,更为重要的是给出流场参数的定量结果。     

超音速多圆弧(MCA)叶型叶珊的研究

本文提出了含脱体激波的叶栅有旋超音速流场的特征线起始线计算方法。采用激波嵌入法解完全的欧拉方程,计算了含脱体激波的多圆弧叶栅的超音速进口流场。本文发展了一种将超音速流场计算与亚音速扩散损失计算的尾迹法相关联的计算超音速叶栅损失的方法,可分别确定激波损失和亚音速扩散损失。      

一种基于聚类分析的二维激波模式识别算法

在激波捕捉求解器计算的可压缩无黏流场基础上，提出了一种探测并识别二维激波干扰模式的新算法，从3个层面详细介绍了该算法的实施流程。首先，采用基于当地流场参数设计的传统激波探测方法，辨识出激波附近的一系列网格单元；其次，通过经典的K-means聚类算法将这些激波单元划分成许多簇，并根据簇的相邻信息定义每个簇的类别；最后，设定相关准则对某些紧邻的簇进行合并，进而确定各个激波干扰点的位置，记录各条激波分支所对应的簇，采用Bézier曲线拟合算法分别对其聚类中心进行拟合以获取更加光滑的激波线。数值试验表明，该算法不受网格类型的限制，不仅可以保证最终拟合的激波线具有较高的位置精度，还可以清晰地识别出流场中多激波干扰的模式，同时对分析非定常流场中激波的运动与演化过程也提供了一种有效的可视化手段。     

非等强度多道冲击波作用下空泡溃灭机制分析

航空航天、生物医学等领域均存在冲击波与空泡作用问题,冲击波尤其是多道冲击波作用下空泡溃灭包含着复杂的多相瞬变行为与物理现象。基于自主搭建的可压缩多相流并行数值平台,对液体中单一/多道冲击波与空泡的作用过程进行了数值模拟。通过对冲击波作用下空泡内波系性质,以及多道冲击波与空泡作用后的系列反射波系在液体中的相互作用过程进行详细解析,分析了不同冲击波作用下空泡的形变演化过程,探究了空泡的溃灭机制与特性。研究发现,相较于单一冲击波的作用,多道冲击波作用下空泡内部及周围液体流场中的波系结构更为复杂,然而无论是单一冲击波或是多道冲击波工况,其最终诱发的空泡溃灭形态均十分相似。而且,当空泡发生溃灭前所受到多道冲击波的总强度与单一冲击波强度相等时,空泡的溃灭波强度也十分相近,为多道冲击波与空泡作用效果的定量化评估提供了理论基础和依据。     

双弯曲入射激波的可控中心体内收缩基准流场设计

采用有旋特征线法设计了一种双弯曲入射激波的可控中心体内收缩基准流场,两道入射激波交于中心体起始点,入射激波和反射激波通过给定激波径向总压恢复系数分布进行反设计,壁面通过给定轴向马赫数分布规律进行反设计.该基准流场分为“三波四区”且压缩效率较高.基于该基准流场设计了圆形进口内收缩进气道并进行了黏性修正,数值计算结果表明:内收缩进气道设计点核心区的流场特征和激波形状与基准流场基本一致;在来流马赫数为4.0~7.0时进气道具有较高的压缩效率和良好的流量捕获能力,设计点喉道截面增压比和总压恢复系数分别为17.7和0.729;来流马赫数为5.0~7.0时内部总阻力系数变化平缓,从0.23下降为0.22.      

马赫数4二维直管内拟似冲击波实验和数值研究

通过计算流体力学数值仿真和高速彩色纹影照片显示技术对马赫数4拟似冲击波进行研究。实验是在日本室兰工业大学的压力—真空型超声速风洞中进行的。数值计算是采用三阶精度的QUICK格式和Spalart-Allmaras湍流模型进行的。用实验验证CFD模型,计算结果与实验显示出较好的一致性。可以得出结论,从目前的CFD模型获得的计算结果是精确的。因而,实验很难获得的拟似冲击波流场内部一些流动量可以通过数值仿真结果进行分析。另外,推测了马赫数4拟似冲击波非对称的原因。     

激光引致水诱导空气击穿所获无约束流场演化过程的实验研究

流场演化过程是揭示激光推进机理的重要研究内容之一。利用纹影系统和PCO-HSFC 高速相机,首次拍摄了自来水诱导CO₂激光击穿空气产生的激波向空气和水中的演化过程。阐述了产生两次气蚀空穴的原因：初始时刻产生的空穴为激波作用于水面所致,第二次空穴出现在相机被触发后约8 μs 时,是激光烧蚀水蒸气作用于水面所形成的低压区。实验结果表明：击穿后流场向激光入射方向（空气中）演化较快,激波初速度较大,约为6 km/s;流场向水中演化相对较缓,激波初速度约为3.33 km/s。两者都服从指数衰减,分别在约73.368 μs 和41.649 μs 时衰减到声速。研究结果对于把水作为工质应用于激光推进有一定意义。