脉冲爆震发动机管外复杂波系的数值计算

为了研究脉冲爆震发动机外流场复杂波系的变化特点,用时空守恒元与求解元(CE/SE)方法对单爆震管和三爆震管脉冲爆震发动机外流场进行了计算。三个爆震管并排排列,管内填充按化学当量比混合的乙烯氧气混和物。通过计算,获得了脉冲爆震发动机外流场的压力分布、温度分布、激波与膨胀波的变化特点。计算结果表明该时空守恒元与求解元方法是一种较好的爆震波模拟方法,能有效地捕捉激波等强间断;三爆震管脉冲爆震发动机相对于单爆震管脉冲爆震发动机外流场复杂,存在多道相互作用、相互影响的激波和膨胀波。     

脉冲爆震燃烧室内缓燃向爆震转捩数值模拟

为了研究脉冲爆震发动机燃烧室内火焰加速及缓燃向爆震转捩过程，利用7组分8反应的氢气详细化学反应机理进行 了2维数值模拟。结果表明：在火焰传播的初始阶段，障碍物、旋涡与火焰的相互作用是主导火焰加速的主要因素。在燃烧区域 产生的膨胀波向前传播并诱导未燃混气流动；障碍物后的回流区及障碍物顶端脱落的旋涡使火焰面拉伸、卷曲，增加火焰面面积， 同时提高燃烧放热强度，使火焰传播速度加快。在火焰发展的后期阶段，激波与障碍物的相互作用会使脱落在已燃混气中的未燃 混气微团发生局部爆炸，产生新的激波进一步推动火焰加速，缩短火焰锋面与前导激波的距离。由于交错型障碍物增加了火焰传 播的距离，其缓燃向爆震的转捩时间比采用对称型障碍物时的更长。     

一种修正的Osher通量在高阶WCNS中的特性

将一种基于多维修正的Osher通量运用于高阶加权紧致非线性格式(WCNS)中,该修正通量主要在垂直于激波面的界面上增加耗散,能够改善Osher通量的激波捕捉稳定性,同时对边界层和接触间断的分辨率影响非常小。对修正的Osher通量在高阶WCNS中的特性进行研究,通过数值模拟测试了基于Osher通量的WCNS的激波稳定性、热流预测精度、边界层模拟能力、激波边界层干扰模拟能力,并与Steger-Warming通量和Roe通量进行了对比。结果表明修正后的Osher通量比Harten修正的Roe通量具有更好的激波捕捉鲁棒性,而边界层、驻点热流值和激波边界层干扰的模拟则明显优于Steger-Warming通量。上述结果说明了基于修正的Osher通量的高阶WCNS具有较好的激波捕捉特性、热流预测精度和边界层计算能力。     

非定常同向激波边界层干扰的实验和数值模拟

运动激波在波前同向气流中传播时与壁面附近力界层发生的相互作用，既不同于定常的激波边界层干扰，也与普通激波管中端壁波与边界层的作用有很大差异。通过实验和数值模拟对这一问题进行研究表明，激波与边界层作用后，波面在壁面附近发生了明显弯曲，且触及壁面；弯曲的激波与壁面进一步作用后发生反射，由于波前气流速度不同，这种反射存在规则反射和马赫反射两种类型；由于激波反射，波后壁面附近形成了一个高压区；激波与边界层     

环形爆震室中火焰加速数值模拟及验证

对带不同数量孔板的环形爆震室进行数值模拟,并通过试验对数值计算进行验证,来研究火焰加速现象、爆燃向爆震转变过程和不同当量比下起爆距离.数值计算采用二维轴对称非定常Navier-Stokes方程来模拟流体动力学过程.研究发现用较低的点火能量对火焰混合区的可燃气点火产生低速火焰,低速火焰向环形爆震室射流并改变方向向出口传播,火焰在孔板的阻碍作用以及火焰诱导激波和反射波的加速作用下,由层流变为湍流,湍流火焰与其诱导激波相互加强,最终引爆未燃混气;还对爆震波在孔板区的传播过程进行了分析,对不同当量比下的火焰速度和起爆距离进行了模拟研究.      

电弧放电等离子体激励控制超声速压气机叶栅激波/边界层干扰仿真研究

为研究电弧放电等离子体激励对超声速压气机叶栅激波/边界层干扰的控制作用,建立了模拟等离子体激励作用效果的唯象学模型,进一步以ARL-SL19超声速叶栅为对象,通过数值仿真研究了电弧放电等离子体与叶栅通道内部流动的相互作用及其对叶栅流动损失的影响。结果表明:等离子体唯象学模型能够较好模拟电弧放电等离子体诱导产生冲击波的气动特性。电弧放电等离子体激励对叶栅通道内部流动主要具有三种作用效果:在放电区,注入的热量会产生阻塞效应,增加近壁面气流的流动损失;在激波/边界层相互作用区,能够改变激波系结构,减小激波损失;在尾迹区,冲击波会诱导产生脱落涡。     

突跃型与平滑型斜爆震波起爆机制数值模拟

采用带有化学反应的Euler方程,对突跃型与平滑型2种形态的斜爆震波(ODW)在起爆机制和结构特点等方面的差异进行了研究.结果表明:斜激波(OSW)后是否存在亚声速区是判定斜爆震形态的依据.当波后火焰抬升斜激波面使得亚声速区存在时,横向激波与三波点结构就会出现并形成突跃型斜爆震.当波后所有区域皆为超声速区时,波后火焰燃烧无法影响到上游激波面,则会形成平滑型斜爆震.通过斜激波关系式给出了横向激波形成所需的临界斜激波角度,只有当斜激波角度大于此临界角度时才能形成突跃型斜爆震波,反之则形成平滑型斜爆震波.      

非结构网格通量重构算法下三种紧致WENO限制器对比研究

在通量重构算法框架下对比研究了3种最新提出的用于间断伽辽金算法的紧致型WENO(weighted essentially non-oscillatory)限制器，即简单WENO限制器、p阶加权WENO限制器和多精度WENO限制器。这3种WENO限制器能够较高精度地模拟流场并捕捉激波，其紧致性在于问题网格单元中的重构只涉及问题网格单元及其相邻网格单元。同时，使用保精度的保正限制器用以避免可能出现的非物理的密度和压力负值。最后，采用非结构四边形网格在双马赫反射、激波与涡相互作用、激波反射、黏性弓形激波和激波与层流边界层相互作用等多个二维算例中，对这3种WENO限制器进行分析比较。结果显示：多精度WENO限制器与p阶加权WENO限制器能够高精度模拟流场并捕捉激波，同时一定程度抑制通量重构方法本身在激波附近的数值伪振荡；p阶加权WENO限制器与多精度WENO限制器相比，其稳态收敛性相对更好；简单WENO限制器则性能较差。根据以上研究，提出亟须发展能够使高阶WENO限制器在稳态问题中收敛的间断探测器。     

爆震管内波与火焰相互作用机理的试验

为了解波与火焰相互作用的机理,在60 mm×60 mm×2000 mm方爆震管内,用乙炔(C₂H₂)和空气混合物进行了单爆震性能研究.利用压力传感器与离子探针同时测得爆震管内的压力和火焰传播速度.根据波和火焰触发的不同时刻来分析爆震管内波与火焰相互作用的过程.结果表明:光滑爆震管内没有产生爆震波,压力波始终在火焰前面;加入阻塞比为0.4的扰流器后,爆震管内产生了爆震波.在扰流器内部,爆震波比较复杂,弱压缩波在火焰前面,但火焰在激波前面;在扰流器出口位置以后,激波在火焰前面.      

带尖针杆的钝体粘性绕流的数值模拟

用无波动、无自由参数耗散的差分（ＮＮＤ）格式数值求解了Ｎ－Ｓ方程,模拟了尖锥柱与球头组合体的高超音速粘性流动;清晰地给出了激波、膨胀波、剪切层和旋涡等流动图象。计算给出的定常流场与实验图象甚为一致。通过数值模拟研究了激波、旋涡的形成和发展过程以及相互作用情况。发现分离区中的旋涡存在分叉、合并及二次分离等演变历程。钝体主激波和分离激波从开始形成到相互撞击;发展到在与分离区相互干扰下向前、向后运动和变形;最后趋于定常状态。     

高负荷低压风扇三维气动设计与损失分析

在给定的限制条件下,通过从S2反问题计算出发,依靠在粗、细网格上求解三维黏性Navier-Stokes(N-S)方程的设计方法完成了某高负荷、跨声速低压风扇的气动设计.该设计方法抛弃了准三元设计对损失模型的依赖,较快的实现了跨声速风扇的气动设计.在设计结果的基础上,对风扇流场结构的分析表明,动叶顶部激波,以及激波、泄漏流和附面层的相互作用是风扇损失的主要来源.全三维的弯扭叶片设计实现了静叶的高亚声大折转角扩压流动.      

环形火焰引发爆震的数值研究

采用改进的适于多组分反应流的波传播算法,对直管道中甲烷．空气预混气的爆震引发过程进行数值研究,考虑了两种不同构型的初始环形火焰对引发爆震的影响。对这两种构型下爆震引发过程中的压力场、温度场、反应放热率以及组分浓度进行了计算,由此探讨了“热点”形成机制及其对爆震引发的作用。结果表明,火焰发展产生的燃烧压缩波在管壁或管轴心发生汇聚、碰撞所形成的复杂波系与火焰相互作用,导致火焰前锋与压缩波之间的局部区域反应剧烈,反应放热率明显增加,从而促使热点形成。形成的热点迅速放大成为过驱爆震波,并逐渐衰减为稳定爆震波。稳定爆震波阵面参数与CJ理论结果进行了对比,两者符合得较好。     

旋转爆震发动机燃烧室壁面烧蚀热防护研究

为了研究碳化烧蚀材料对旋转爆震燃烧室的被动热防护作用效果,建立了旋转爆震燃烧室热环境计算模型和多层复合壁一维烧蚀热响应分析模型。基于旋转爆震燃烧过程的数值模拟,确定了沿燃烧室轴向的壁面平均温度与热流的分布特征;针对典型的壁面热负荷输入条件,采用动边界隐式差分格式,对碳化型烧蚀材料的烧蚀过程和传热过程进行了耦合计算,分析了碳化层剥蚀、热解热、热解气体质量流率、材料厚度等因素对燃烧室壁面热防护效果的影响。研究结果表明,旋转爆震燃烧室温度和热流密度分布不均匀,斜激波扫掠区域热负荷更为剧烈;碳化层剥蚀后会使烧蚀材料的烧蚀率增加,壁面温度迅速上升;增大热解热及热解气体质量流率都有利于延长工作时间、降低壁面温度; 旋转爆震燃烧室不同区域的烧蚀有一定差异,增加烧蚀材料厚度无疑有利于壁面热防护,但是在实际应用中应综合考虑壁面重量和空间尺寸。     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

Steady and unsteady numerical investigation of transitional shock-boundary-layer-interactions on a fan blade

The present numerical study was performed to investigate the impact of both the Reynolds number variation and the used turbulence model to capture the boundary layer development on the characteristic of a BR710 fan blade. A one-equation model of Spalart Allmaras with an optional semi-empirical transition model of Abu-Ghanam Shaw has been applied. The transition model allows the boundary layer development from a laminar to a turbulent behaviour to be taken into account. This is of particular importance at low Reynolds numbers and thus high operating altitudes of airplanes when the transition location moves further downstream to the trailing edge. Therefore the interaction between the shock waves of the transonic fan blade row and the boundary layer leads to a significant change of the characteristic. The steady numerical studies of the 3-D blade passage have been carried out with the commercial 3-D Navier–Stokes solver NUMECA. Additionally an unsteady calculation has been applied to explain the penalty in efficiency on high operating altitudes.     

三维激波／附面层相互作用的起始分离预测

给出了一种预测三维激盘／附面层相互作用诱导的流动起始分离的方法。研究结果表明：激波与附面层相互作用所诱导的二次流动是影响起始分离的重要因素。当来流相对马赫数大于１．５时，强的激波／附面层相互作用可能导致跨音风扇转子叶尖区域的流动分离。     

一类关联转子分离尺度与损失的工程模型

将分离与损失相关联，提出了利用转子进出口参数来确定跨音风扇转子内部分离流动的方法。应用该方法对几个转子流场计算的结果表明，在转子的叶尖区域，由于强的激波与附面层相互作用，使得流动在激波后不久便发生分离。叶尖漏流涡可能对激波与附面层的相互作用诱导流动分离有较大影响。     

Conjugate heat transfer investigations of turbine vane based on transition models

The accurate simulation of boundary layer transition process plays a very important role in the prediction of turbine blade temperature field. Based on the Abu-Ghannam and Shaw (AGS) and c-Re h transition models, a 3D conjugate heat transfer solver is developed, where the fluid domain is discretized by multi-block structured grids, and the solid domain is discretized by unstructured grids. At the unmatched fluid/solid interface, the shape function interpolation method is adopted to ensure the conservation of the interfacial heat flux. Then the shear stress transport (SST) model, SST & AGS model and SST & c-Re h model are used to investigate the flow and heat transfer characteristics of Mark II turbine vane. The results indicate that compared with the full turbulence model (SST model), the transition models could improve the prediction accuracy of temperature and heat transfer coefficient at the laminar zone near the blade leading edge. Compared with the AGS transition model, the c-Re h model could predict the transition onset location induced by shock/boundary layer interaction more accurately, and the prediction accuracy of temperature field could be greatly improved.     

高温区碰并诱导气相爆轰二次起爆的数值研究

针对心形管内气相爆轰二次起爆过程进行了数值研究。采用二阶附加半隐的龙格-库塔法和5阶WENO格式求解二维欧拉方程,采用基元反应来描述爆轰化学反应过程。结果表明:气相爆轰在心形管中传播,由于壁面几何结构的作用,爆轰波熄灭后会发生二次起爆。该二次起爆由高温区碰并所诱导,且在管道不同位置诱导了两处二次起爆。爆轰流场中存在未反应气囊、射流、爆轰波马赫反射、激波与火焰作用及其所诱导的不稳定性等现象。计算得到的二次起爆过程胞格结构演变与实验基本一致。