高精度ENO格式在喷管流动模拟中的应用

本文分析了ＥＮＯ格式的特点,并应用于Ｅｕｌｅｒ方程和全ＮＳ方程的迁移项和压力项,模拟了各种喷管流动。首先计算了ＪＰＬ喷管流动,取得了和实验一致的结果,分析了喷管内产生弱激波的原因;其次计算了非定常二维喷管流动,给出了相应的计算结果,并分析了壁面附近流场的拓扑结构;最后计算了摆动喷管的流动及推力矢量计算,取得了推力随摆角变化的规律,和实验结果表明基本一致。     

收敛-扩张喷管中运用激波矢量控制法的计算研究

采用二阶TVD格式的有限体积法耦合RNG k-ε湍流模型和非平衡壁面函数求解二维守恒型雷诺平均N-S方程，数值模拟了运用激波矢量控制法的二维收敛一扩张喷管中的流动现象，很好的模拟出了由于次流喷射所产生的激波及上下壁面的压力分布，并将计算结果与国外的实验结果进行了对比。     

激波诱导矢量喷管流场的数值模拟及试验

为研究固体火箭发动机激波诱导矢量控制效率的影响因素,及验证数值模拟方法的准确性,对激波诱导轴对称推力矢量喷管模型进行了壁面测压试验,采用二阶精度Roe格式和k-ωSST两方程湍流模型求解强守恒型Navier-Stokes方程对矢量喷管复杂干扰内流场进行数值模拟。根据试验和数值模拟结果分析了喷管内主流和次流相互作用产生的复杂流场结构,比较了在不同喷管落压比NPR和次主流压力比SPR下喷管壁面静压的分布情况。结果表明,数值计算和试验结果基本吻合,验证了计算方法的准确性;在一定范围内减小喷管落压比,增大次主流压力比可以增大喷管周向壁面静压差,提高喷管的推力矢量偏角。     

三维数值模拟再生冷却喷管的换热

为了解液体火箭发动机喷管再生冷却的换热特点，采用数值模拟的方法，对内喷管燃气、壁面和冷却液建立不同的三维控制方程，进行流动和传热的耦合计算。在计算中，假定喷管流动为冻结流动，考虑燃气向壁面的对流换热和辐射换热；采用二阶迎风格式离散控制方程，采用DO模型离散求解辐射换热方程，水蒸气的吸收系数根据Leckner公式计算。计算模型采用缩比热试车发动机，数值计算结果与实验结果吻合较好，较准确地模拟出了喷管的壁面热流密度，得到了喷管燃气和冷却液的流场和温度场，对高压再生冷却喷管的设计具有指导意义。     

不同矢量偏角下二元喷管内流特性的数值研究

应用多重网格MacCormack显式格式对不同矢量角下二元喷管的内部流动进行了模拟分析。计算时, 在细网格上采用了雷诺平均N-S方程与k-ε湍流模型方程;在粗网格上则采用Euler方程。数值计算给出了较清晰的流谱, 壁面压力分布的计算也与实验结果相吻合。通过计算结果分析了矢量喷管几何参数对流场的影响, 并对矢量喷管的推力特性进行了数值研究。      

非对称超声速喷管内流动分离非定常特性

综合采用流场观测和动态压力测量技术,对非对称超声速喷管分离流状态下喷管内激波结构和壁面动态压力进行了试验测量,通过壁面压力时频特性分析获得非对称喷管内不同流动分离模态的流动非定常特性。结果表明：在喷管落压比(NPR)为1.8～12.7范围内,喷管内流场结构由下偏向上偏转换;上壁面流动分离模态经历了受限激波分离(RSS)、末端振动状态和自由激波分离(FSS);下壁面流动分离模态主要为FSS;流动分离模态为RSS时,Schmucker分离预测标准不再适用。RSS和末端振动状态模态下,尽管分离间歇区壁面动态压力具有相似的低频特征,激波运动呈显著低频特征,但末端振动状态模态下频率值略高,主要是由于流动再附点近喷管唇口,分离剪切层撞击喷管出口位置,剪切层的不稳定性影响了分离激波的振荡特性。     

管内正激波/湍流边界层干扰的数值模拟

本文采用时均Ｎ－Ｓ方程和Ｂａｌｄｗｉｎ／Ｌｏｍａｘ代数湍流模型计算了典型拉伐尔喷管内正激波与湍流边界层干扰流场。计算与实验结果的比较表明，方法可准确地预测激波结构，激波与边界层干扰区流动特征，激波位置、激波前马赫数和壁面压力分布等。     

轴对称矢量喷管壁面静压分布的试验研究

在模拟实际发动机喷管落压比的情况下, 利用轴对称矢量喷管的缩比模型, 对其壁面静压分布进行了试验研究。研究结果表明, 当矢量偏转后, 轴对称矢量喷管的收敛段与扩张段内流场均呈现出复杂的三维流动的特征; 当喷管落压比小于完全膨胀落压比时, 喷管内部存在周向不对称的斜激波和分离区, 且激波位置随几何矢量角的增大而前移。本文的研究结果可用于检验喷管三维流场的计算精度。      

双钟型喷管高度补偿特性的数值分析

针对双钟型喷管流场中存在的激波、激波与附面层相互作用、大尺度分离流动等复杂的物理现象，应用RNGκ-ε湍流模型封闭二维粘性可压缩N-S方程组，采用二阶迎风格式及隐式Gauss-Seidel迭代方法进行耦合求解，对双钟型喷管流场进行数值模拟研究。结果表明，双钟型喷管在低空出现壁面可控流动分离，避免了过膨胀现象，性能优于大面积比钟型喷管；在高空下为完全附着流动，性能优于小面积比钟型喷管，实现了高度补偿特性。     

二元矢量喷管内流特性的数值模拟

应用多重网格ＭａｃＣｏｒｍａｃｋ显示格式对不同矢量角下二元喷管的流动进行了模拟分析。计算时,在细网络上采用了雷诺平均Ｎ－Ｓ方程和Ｋ－ε湍流模型方程;在粗网格上则采用Ｅｕｌｅｒ方程。数值计算给出了较清晰的流谱。壁面压力分布的计算值也与实验结果吻合。通过计算结果分析了喷管几何参数对流场的影响。所得结果显示出多重网格算法高效可靠,可用于开展参数研究     

激波诱导圆形矢量喷管数值研究

利用时间推进的有限体积法求解三维雷诺平均Navier-Stokes方程,对激波诱导矢量控制方案的圆形收敛-扩张喷管全流场进行了数值模拟.从射流压强比及落压比的变化计算了次流对主流的控制效果.结果表明,利用激波诱导方式可有效的迫使主流偏转,在落压比为4、射流压强比为1.499时,矢量偏转角最大可以达到15.35°.通过分析,圆形喷管中形成的激波为三维曲面激波,结构比较复杂,数值研究可观察到清晰的流场结构,为试验研究提供参考.      

跨声速风洞高速段一体化数值模拟研究

追求高亚声速经济巡航的民机、跨声速高机动特性的战斗机对高性能跨声速风洞的需求日趋紧迫,开 展跨声速风洞高速段一体化数值模拟研究,对跨声速风洞设计具有一定的参考意义。通过非对称平板扩压器 算例,初步验证计算方法的可行性,并对跨声速风洞高速段进行计算收敛评判方法、不同初始条件和槽壁扩张 角等因素研究。结果表明:采用模型区前后两个监测点马赫数变化作为收敛判据,方法可行且模型区流场均 匀;不同初始化条件对收敛结果总体影响较小,特别是各截面流场分布和槽道流动方向上,两者结果基本相同; 跨声速状态槽壁扩张角0.3°得到的试验段模型区域流场品质较槽壁扩张角0.0°更均匀。     

二元喷管跨音速流场的数值计算

本文采用有限体积法进行空间离散,Runge-Kutta法进行时间推进的办法,通过求解Euler方程来模拟二元收敛—扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡,在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解。计算结果与试验数据和其他数值方法吻合良好。      

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

机械式轴对称矢量喷管动态偏转的矢量特性研究

采用限体积法在同位网格中对满足几何守恒定律的控制方程进行离散化,对某一轴对称矢量喷管在几何偏角速度大小为40°/s下从0°偏转到14.2°过程中的内外流场采用压力隐式分离算法(PISO)进行了数值模拟,分析了该喷管在不同的落压比、不同的外流速度下偏转时,喷管内外流场的动态变化特性,几何偏角和气动矢量角的关系,当喷管偏转后在内壁面产生分离时,喷管的气动矢量角振荡。研究结果为矢量喷管的设计和发动机控制技术的研究提供了参考。      

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

Flow characteristics around airfoils near transonic buffet onset conditions

In transonic flow, buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation. The phenomenon is common in transonic flow, and it has serious impact on the structural strength and fatigue life of aircraft. In this paper, three typical airfoils: the supercritical OAT15A, the high-speed symmetrical NACA64A010, and the thin, transonic/supersonic NACA64A204 are selected as the research objects. The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method, and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD). Qualitative and quantitative analysis of the shock wave motion, shock wave intensity, shock foot bubble and trailing edge separation, and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions. The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack, while the growth rate of the primary mode increases inversely. It was found that at the onset of buffet, there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle, and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.     

喷管跨声速流场计算

本文用时间相关法求解定常跨声速喷管流场。对守恒型和非守恒型方程采用全场统一时间步长、当地时间步长和修正超阻尼法进行了计算。计算表明:采用当地时间步长可大大提高收敛速度,采用全场统一时间步长时,数值解非守恒型方程收敛较快;采用修正超阻尼法也能提高收敛速度,但与阻尼系数和时间间隔有关。     

喷管激波贴口压比计算

本文根据试验数据,利用多元线性回归分析方法建立了轴对称收敛—扩散喷管激波贴口时落压比与面积比,收敛半角,扩张半角,无量纲喉道曲率半径之间关系的数学模型。计算结果为喷管试验研究和跨音速流计算提供了依据。      

马赫数连续可变跨声速湿蒸汽风洞的研制

用于跨声速气动测量的探针须从亚声速到超声速范围进行标定。变质量槽式喷管通过扩张段壁面上槽缝流出部分气流的自适应特性可在不同背压下得到不同出口马赫数,从而使标定气动探针的风洞实现马赫数从0到超声速的连续变化。为了研究采用湿蒸汽为工质的变质量槽式喷管的性能及优化其结构,采用三维犖-犛方程以及可实现犽-ε湍流模型对其进行了详细的数值仿真。结果表明收缩段型线、扩张段长度及壁槽尺寸等对喷管流场特性有重要影响,喷管进出口压比在一定范围内,槽式喷管有最优的收缩段型线、扩张段长度和开槽尺寸。根据数值仿真结果研制了马赫数从0到1．6连续可变的跨声速湿蒸汽风洞,对此风洞性能进行验证,表明该风洞在马赫数从零到超声速范围内可获得均匀、稳定的出口气流,满足跨声速湿蒸汽气动探针的标定要求。