三维湍流高速进气道内外流场的高效高分辨率解

本文提出一种高效、高分辨率求解三维复杂流场的隐式算法，并成功地用于求解高速进气道内外流场的三维、可压缩、雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，获得了三个工况五种流场的三维数值解；与有关实验数据比较，结果令人满意。这种算法由Ｊａｍｅｓｏｎ－Ｔｕｒｋｅｌ的ＬＵ分解法和张涵信ＮＮＤ格式的数值通量耦合而成，湍流采用ＢａｌｄｗｉｎＬｏｍａｘ模型。大量数值计算表明，该算法能捕捉到复杂的激波系，有较好的激波     

等截面矩形隔离段内流场的三维数值模拟

针对三维矩形等截面隔离段中由激波／附面层干扰诱导的复杂流场，采用有限体积方法，OC－TVD格式结合B－L代数湍流模型隐式求解了三维雷诺平均N－S方程，数值模拟了三维流场内复杂的激波串和旋涡结构，计算结果表明所采用的数值方法能够很好地捕捉隔离段内流场的复杂流动现象。     

压力恢复系统扩压器激波串现象的数值模拟

通过求解由S-A一方程湍流模式封闭的三维雷诺平均Navier-Stokes方程,采用有限体积方法对激光器压力恢复系统变截面扩压器常温状态下的三维流场进行了数值模拟.数值计算和实验符合良好,较好模拟了扩压器中由激波/边界层干扰诱导的复杂流场的流场特性并再现了流场中的"激波串"和"伪激波"现象.本文对激波串的形成机理进行了分析,发现了几种变截面扩压器特有的激波串结构,并研究了扩压器扩张角对扩压效率的影响.     

超声速进气道湍流流场的数值模拟

采用有限体积法数值模拟了超声速进气道的湍流流场,改进了离散纳维尔－斯托克斯方程中粘性项的方法,使之具有较好的守恒性,研究使用了变系数的隐式残差光顺等加速收敛技术,并讨论了超音来流、收敛型进气道复杂流场初场的给定方法,所得到的数值解正确地展示了流场的复杂激波结构。     

固体火箭冲压发动机进气道三维数值模拟

基于RANS平均的N-S方程、标准k-ε双方程湍流模型,利用FLUENT软件对某型固体火箭冲压发动机进气道内外流场进行了三维数值模拟.研究了不同马赫数下进气道前体、入口、进气道出口流场和流动特性;考虑了弹体,缩比,网格划分对进气道流场出口参数的影响.模拟结果表明:随着飞行马赫数的增加,由前缘点产生的激波逐渐靠近弹体,进气道内部总压恢复系数降低,在进气道出口处压力比较高,出口流场基本均匀.      

运用GAO-YONG湍流模型对扩压器内跨声速流动的数值模拟

运用GAO-YONG可压缩湍流方程组,采用同位网格SIMPLE算法,对扩压器跨声速流动中的二维激波/湍流边界层干扰现象进行了数值模拟。将计算得到的流场的时均参数与实验值进行比较,数值模拟结果在激波强度、壁面压力分布以及分离点和再附点位置等方面,与实验值吻合较好,表明GAO-YONG可压缩湍流方程组能够比较准确的模拟较强激波/湍流边界层干扰流动,从而进一步为GAO-YONG湍流模型的正确性及其在可压缩流场模拟方面的适用性提供了佐证。     

高超声速内转式进气道流动的壁面丝线显示

针对风洞实验观测高超声速内转式进气道内部流动困难,不易获得内流道三维激波/边界层干扰主导的复杂流场结构的问题,通过拓展壁面丝线流动显示技术的潜力,在Ma∞5.9激波风洞中,借助高速摄影实时拍摄丝线流谱,并结合纹影、壁面压强测量以及数值模拟分析,验证了丝线的动态响应特性,丰富了内转式进气道的观测技术,获得了进气道的流场结构。采用预设堵块实验方法,在激波风洞中考核了内转式进气道模型的自起动能力。结果表明,直径约为0.1 mm,长度约为15 mm的402号涤纶/棉缝纫线的跟随性较好,能够直观、动态地显示出壁面流动分离区的范围,为判断内转式进气道是否起动提供了依据。内转式进气道模型在实验条件下能够自起动,起动状态下进气道唇口附近的波系结构以及前体压缩面的丝线流谱和压力分布与数值模拟符合较好。     

旋转冲压压缩转子三维进气流道数值研究

借鉴典型三维超声速进气道的设计方法,设计了一种旋转冲压压缩转子三维进气流道,并采用三维雷诺平均N-S方程和Spalart-Allmaras湍流模型对其流场进行了数值研究.计算结果表明:与三维超声速进气道中的平直激波系不同,三维进气流道中的激波为曲线激波;与二维进气流道中的激波系相比,三维进气流道中的激波系要相对模糊和复杂;转速和背压对三维进气流道的性能有较大影响.      

气动塞式喷管的流场数值模拟和设计参数探讨

针对气动塞式喷管流场中存在的激波、激波与附面层相互作用、大尺度分离流动等复杂的物理现象 ,本文采用扩展的压强校正法求解用紊流模型封闭的可压缩湍流平均 N-S方程组 ,利用交错网格系统抹平迭代过程中的数值振荡 ,并通过高分辨率的 TVD格式进一步改善压强校正法的激波捕获效果 ,对气动塞式喷管的流场进行了详细的数值模拟。并进一步探讨了喷管倾角、塞锥长度和二次流流量比对喷管流场结构和性能的影响。     

超声速三维进气道内激波干扰的研究

为了为三维进气道的设计提供有用的分析数据,对高超声速和超声速来流下三维进气道内激波干扰进行了理论和数值研究。进气道模型选取"箱式"以及三面侧压进气道作为研究对象。理论分析采用了"空间降维"方法,即将进气道各个角落处的三维双楔定常激波干扰问题转换为二维非定常激波干扰问题,并利用激波动力学对其进行求解。数值验证方法利用2阶NND差分格式求解三维无粘欧拉方程,网格数量为1200多万,并采用MPI并行进行计算。该理论分析方法很好地对进气道各个角落处的激波干扰波系结构进行了判别,并得到了干扰区马赫构型三波点附近以及规则构型反射点附近的解析解,理论分析结果与数值模拟结果吻合较好。此外,针对进气道截面内各个流场区域的总压恢复系数以及压力、密度和温度进行了研究,并考虑了箱式进气道和三面侧壁压缩进气道内的流场区域的非均匀性,干扰区马赫杆后的总压损失要比其他区域高10%左右。通过研究表明,"空间降维"方法适用于进气道压缩部分,将为进气道的设计和性能评估提供一种理论分析手段。     

数值模拟纵向隔热屏加力室热态流场

数值研究纵向隔热屏加力燃烧室的在三维贴体坐标系下燃烧整体流场,所用的数学模型有:k-ε模型、EBU-Arrehenius和二阶矩紊流燃烧模型及六通量模型等.数值分析两种不同几何形状、进口气流参数分布及紊流燃烧模型等对加力室内各气流参数、隔热屏和加力室筒体壁面温度分布的影响,计算结果与试验数据比较表明:不同几何形状对加力室紊流燃烧流场的影响要比进口气流参数分布大些,二阶矩模型更适用模拟紊流燃烧流动.      

考虑附面层影响的二元混压式进气道设计

采用等激波强度的方法,考虑附面层修正,设计了一种飞行马赫数Ma=3.0的二元混压式进气道.通过数值仿真,模拟了激波-边界层的相互影响,研究了附面层抽吸对内流场的影响,获得了进气道内部复杂的流场分布,以及不同背压下进气道的起动特性.计算表明所设计的进气道性能较好,附面层抽吸对稳定正激波有明显的作用,提高了进气道抗反压能力.给出的方法可用于二元混压式进气道的初步设计和验证.      

扩压器内二维湍流流动的时间推进数值解

ＡＦ有限差分理论。指出薄层Ｎ－Ｓ方程抛物化后右端粘性项是主元素变量矢的一个简单线性函数,这对实际计算很有用。应用隐式ＡＦ方法计算了为发动机进气道设计的典型扩压器内的复杂湍流流场。研究表明,计算结果与实验结果基本一致,这为三维非定常湍流流动进一步研究提供了基础。     

冲压发动机进气道/燃烧室/尾喷管耦合流场计算

采用耦合求解的方法进行了冲压发动机进气道/燃烧室/尾喷管流场一体化计算,进气道采用单一气体、变比热比模型,燃烧室和尾喷管采用多组分化学动力学模型,湍流模型为k-g湍流模型。进气道与燃烧室界面通过流量、静压等参数实现耦合。结果表明,通过合理设计参数的传递方式,可以快速实现进气道/燃烧室/尾喷管流场的耦合求解,流场分布合理。     

任意边界扩压器内紊流流动的时间推进数值计算

研究了隐式ＡＦ有限差分理论。指出薄层抛物化控制方程中的右端粘性项是主元素变量矢的一个简单线性函数。给出了任意边界扩压器内紊流流动的数值计算结果，并和实验数据进行了对照。两者吻合得相当好。研究结果对ＡＦ有限差分理论以及推广至三维复杂管道内的非定常流动计算具有重要的参考意义。     

双侧进气突扩燃烧室冷态两相流场的LDV实验研究

本文用二维 L DV研究了两种双侧进气突扩燃烧室中三维冷态两相流场。结果表明 ,这两种燃烧室可在头部产生较强的气相及颗粒相逆流运动 ,且有一种结构还可使头部气相及颗粒有旋流运动。本研究还表明两相间有速度滑移 ,气相及颗粒相均有较高的湍流强度。这些特征都有利于点火 ,火焰稳定及强化燃烧。     

大涡模拟进口合成湍流方法

进口瞬时湍流脉动条件直接影响大涡模拟的精度。对比研究了更适合工程应用的合成涡和合成谱方法,并通过考虑进口速度时间相关和随机数的概率密度分布,对合成谱方法进行改进。对零压力梯度平板湍流模拟结果的对比表明:不同湍流合成方法得到的合成湍流均存在发展为真实湍流的过渡过程;考虑进口湍流脉动的时间相关有助于合成湍流更快地发展为真实湍流,需要的流向距离更短;合成谱方法采用正态分布系数优于均匀分布,改进的合成谱方法优于常规合成谱方法。     

斜切径向旋流器环形燃烧室数值模拟

采用三维贴体坐标系统,对包括突扩扩压器、帽罩、旋流杯、火焰筒以及内、外环通道的环形燃烧室全流程2相燃烧流场进行数值研究。采用k-ε湍流模型、2阶矩-EBU湍流燃烧模型、6通量辐射模型和颗粒轨道模型等模拟2相湍流燃烧流动,分析了进口工况对全流程燃烧流场的影响。计算结果表明:随着进口工况的改变,燃烧室出口温度场也发生相应变化;计算结果与试验数据比较表明:所用的数学模型合理、计算方法可行,其结果可为某型燃烧室优化设计提供可靠依据。     

An integration method based on a novel combined flow for aerodynamic configuration of strutjet engine

In this paper a novel design method of aerodynamic configuration is proposed to integrate forebody, strut and inlet for strutjet engine, and a model at design point of Mach number 6 is generated to investigate the aerodynamic performance by both simulations and experiments. The basic flow field employed by proposed method is a combined flow named IBB, which is combined by Internal Conical Flow A (ICFA), truncated Busemann flow I (BI) for external section, and truncated Busemann flow II (BII) for internal section. The model configuration is generated by streamline tracing method from basic flow field, in which the forebody section is traced from ICFA and BI flows, and the inlet as well as strut section is traced from BII flow. The simulations in Mach number 4, 5, and 6 demonstrate uniform starting flow fields with relatively high total pressure recovery, which agree well with experiments in wind tunnel. Additionally, in low Mach number cases, this inlet could start at Mach number 3 while it is unstarted at Mach number 2.7; in high Mach number cases, a uniform flow could still exist in Mach number 6.5 while a relatively strong shock wave boundary layer interaction is found in cowl area of Mach number 7 case, indicating the inlet designed by proposed method works in a relatively wide Mach number range.     

辅助动力装置环形回流燃烧室数值研究

利用Fluent软件计算辅助动力装置(auxiliary power unit,简称APU)环形回流燃烧室三维两相喷雾燃烧流场,研究不同进口温度和油气比对两相喷雾燃烧流场的影响,采用标准 k-ε 模型模拟湍流黏性,离散相模型(DPM)追踪油珠运动轨迹,燃烧模型采用非预混平衡化学反应模型.计算结果表明:随着进口温度和油气比的增加,燃烧室出口温度相应增加,但温度分布规律基本保持不变;计算结果与试验测量结果比较吻合,说明采用的数学模型和计算方法可用于预估实际APU回流燃烧室三维两相喷雾燃烧流场.