斜切反喷管粘性流场数值模拟

为了在级间分离期间提供反推力，许多固体火箭发动机前端都装有一组斜切反管，由于反管气动型面具有尖点，并在超音速区有台阶，喷管内存在一系列激波，并伴有流动分离现象。从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，利用时间相关法，采用ＭａｃＣｏｒｍａｃｋ两步显格式，结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，数值模拟了斜切反喷管流场。     

真空羽流场的 N-S 和 DSMC 耦合数值模拟

根据调姿发动机真空羽流场的特点,采用差分求解Ｎ－Ｓ方程数值模拟喷管内流场,用ＤＳＭＣ方法数值模拟喷管外羽流场,得到了符合流动规律的数值模拟结果。为姿控发动机喷管设计提供了参考,并为真空羽流污染研究打下了良好的基础。     

斜切反喷管粘性流场数值模

为了在级间分离期间提供反推力，许多固体火箭发动机前端都装有一组斜切反喷管，由于反喷管气动型面具有尖点，并在超音速区有台阶，喷管内存在一系列激波，并伴有流动分离现象。从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，利用时间相关法，采用ＭａｃＣｏｒｍａｃｋ两步显格式，结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，数值模拟了斜切反喷管流场。计算得到的壁面压强分布与风洞吹风实验测得的压强分市相当一致。     

塞式喷管流场数值模拟

采用ＮＮＤ２Ｍ差分格式,从考虑к－ε湍流模型的理想气体三维薄层近似Ｎ－Ｓ方程出发,对塞式喷管内外流场进行了数值模拟。计算结果与国外文献十分一致。通过流场分析,对塞式喷管流动机理有了清楚的认识。计算表明ＮＮＤ２Ｍ差分格式及相应计算软件适合于模拟复杂喷流流场,可用于塞式喷管发动机的研究。     

分离流动对矢量喷管性能的影响

采用有限体积法计算了轴对称矢量喷管三维粘性流场，通过改进人工粘性及其边界条件提高了计算的收敛性和数值模拟的精度。利用数值模拟研究了矢量喷管的分离流动对矢量推进性能的影响，通过同实验结果对比，验证了数值模拟的精度。研究表明，分离流动对矢量喷管的性能参数尤其是有效矢量角具有较大影响。     

低密度小喷管流场的DSMC直接模拟

用ＤＳＭＣ方法数值模拟了低密度轴对小喷管的流场流动，获得了相应的结果，并将所得结果与用差分求解Ｎ－Ｓ方程方法获得的结果进行了对比，二者吻合较好。     

固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。     

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

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

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

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

轴对称超临界粘性喷流的数值模拟

从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，采用隐式的Ｂｅａｍ－Ｗａｒｍｉｎｇ格式中的矢通量分裂技术，结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，数值模拟了收敛喷管的内流场和喷射流场，通过对趋临界工况下喷管跨音速流场的数值计算，表明本文的结果与试验结果是吻合的。     

二元收-扩喷管与球形收敛调节片喷管的性能对比研究

为了对比二元收-扩喷管(2DCD)与球形收敛调节片喷管(SCFN)的性能,基于CFD数值计算软件对2种喷管在相同工作点下进行了内外流场计算与分析。数值模拟结果表明:在相同高度及飞行马赫数条件下,SCFN的流量系数相较于2DCD呈现出较大不同,推力系数则略差于2DCD;在低落压比条件下,2种喷管内流对后体阻力特性影响不大;在高落压比条件下,2种喷管内流对后体阻力特性影响较大。     

Experiences in accurately predicting time-dependent flows

As computational fluid dynamics matures, researchers attempt to perform numerical simulations on increasingly complex aerodynamic flows. One type of flow that has become feasible to simulate is massively separated flow fields, which exhibit high levels of flow unsteadiness. While traditional computational fluid dynamic approaches may be able to simulate these flows, it is not obvious what restrictions should be followed in order to insure that the numerical simulations are accurate and trustworthy. Our research group has considerable experience in computing massively separated flow fields about various aircraft configurations, which has led us to examine the factors necessary for making high-quality time-dependent flow computations. The factors we have identified include: grid density and local refinement, the numerical approach, performing a time-step study, the use of sub-iterations for temporal accuracy, the appropriate use of temporal damping, and the use of appropriate turbulence models. We have a variety of cases from which to draw results, including delta wings and the F-18C, F-16C, and F-16XL aircraft. Results show that while it is possible to obtain accurate unsteady aerodynamic computations, there is a high computational cost associated with performing the calculations. Rules of thumb and possible shortcuts for accurate prediction of massively separated flows are also discussed.     

翼型分离流动的等离子体控制数值分析方法研究

利用等离子体进行流动控制是当前的研究热点之一.本文研究通过数值方法模拟等离子体对流场作用的实现方法,及利用数值模拟方法研究翼型大迎角分离流动的等离子体控制.利用CFD软件Fluent中的自定义函数接口,通过C语言编程在软件中引入DBD等离子体激励模型外加体积力源项,对NACA0015翼型大迎角下的等离子体控制进行数值模拟.验证DBD等离子体激励在抑制流动分离与增升减阻方面的作用.结果表明：流动控制效果与DBD激励器布置位置有直接关系;激励器的数量与激励强度均会影响流动控制作用.     

平面压气机叶栅通道分离流动高精度直接数值模拟

采用高精度有限差分格式直接求解二维Navier-Stokes方程组,数值模拟V103平面压气机叶栅分离流动,数值结果表明:在瞬时流场中,吸力面后部发生流动分离,在分离区前端存在大尺度分离涡,分离涡下游是由二次涡和脱落涡交替形成的涡串,直至叶片尾缘,形成以脱落涡为主结构的尾迹;在时均流场中,吸力面后部存在短分离泡,分离区压力分布存在明显压力平台。与逆压力梯度下平板边界层分离流动相比,瞬时和时均流场结构相似;叶栅通道内无量纲涡脱落频率是前者的两倍。与文献计算结果对比表明:叶片表面压力分布除分离区外吻合很好;非定常计算所得分离区轴向长度比定常计算大41%。在分离区内三个二阶统计量均达到最大值,表明流场强非定常性集中在分离区。     

球形收敛调节片推力矢量喷管的发展

球形收敛调节片喷管(SCFN)是一种正在研制中的先进的轴对称/非轴对称混合式推力矢量喷管。简述了它的发展现状并总结了其结构特点,并指出,它是一种具有发展潜力的推力矢量喷管。     

绕钝体分离流动的非定常数值模拟研究

利用非定常RNGk-ε模型对三种典型的钝体绕流问题进行了研究,这三种流动分别为:方柱绕流、圆柱绕流和绕立方柱流动。研究表明了采用非定常雷诺平均方法可以研究钝体绕流的非定常流动,但所获得的结果各不相同。对于分离角固定的方柱绕流和绕立方柱流动,数值模拟结果和实验以及大涡模拟结果吻合较好,对于分离角不固定的圆柱绕流,随着雷诺数的增加,数值模拟结果和实验结果相差越大,分析认为这是由于壁面模型无法准确预测流动中不固定的分离角所致。改进了绕立方柱流动的模拟方法,采用非定常模拟方法后获得的结果要比定常模拟方法得到的结果有了根本改善。     

双三角翼俯仰振荡运动的流场特性数值模拟

 建立了适用于双三角翼大迎角非定常分离流场模拟的数值方法,研究双三角翼俯仰振荡时的动态流场特性,给出动态流场结构和气动力性能随迎角的变化规律,重点考察了减缩频率、转轴位置、平均迎角和振幅等参数对动态流场迟滞效应和气动力曲线迟滞环的影响。研究结果表明:俯仰振荡到相同大迎角时上仰和下俯的流场存在明显差异;减缩频率对气动力迟滞效应的影响相对大于转轴位置;平均迎角的变化导致双三角翼背风区流场结构呈现不同流态,而振幅的大小决定这些流态的数目,事实上俯仰运动时如果跨越的流态数目越多则流场结构的动态响应滞后现象就越显著。通过数值分析,有利于提高对双三角翼在俯仰振荡运动条件下的非定常特性和流场滞后效应等非线性现象的认识。     

蒸汽轮机长叶片颤振预估方法研究

本文根据叶栅实验和流场计算结果对蒸汽轮机末级流动特征进行了分析。在小容积流量工况,末级流场可分为根部脱流区、脱流区之上的主流区叶栅槽道存在附体流区和分离区。三区大小随容积流量变化。由于末级流动复杂,发展工程上实用的算法很有吸引力。实验结果和理论分析表明,在小容积流量工况容易诱发叶片自激振动。为了能预估叶片颤振,本文发展了系列变形激盘法(机时少,适于工程应用)、数值方法 (能给出叶片表面压力分布和激波振荡,有助于了解叶片颤振发作机理)。经实验证明,可以用于叶片设计阶段颤振预估。      

球面收敛二元喷管电磁散射特性

采用迭代物理光学法和等效边缘电磁流法,对球面收敛二元喷管电磁散射特性进行研究.计算分析了该喷管电磁散射的构成及散射机理、散射场的空间分布特性、散射的频率特性,分别计算和对比分析了电磁波在俯仰平面和偏航平面入射时喷管在中间和最大两种状态下的电磁散射情况.研究表明该喷管在垂直极化下边缘绕射场对散射总场的贡献在大角度下体现.电磁波入射频率越高,该喷管后向散射越强烈,雷达散射截面分布对姿态角度越敏感.      

介质涂覆位置对球面收敛喷管电磁散射特性影响

外形设计和隐身材料使用是缩减目标雷达散射截面积(RCS)的两种常用方法。为了研究雷达吸波材料(RAM)对于球面收敛矢量喷管(SCFN)的RCS减缩效果,采用引入阻抗边界条件后的迭代物理光学(IPO)法,研究了球面收敛二元喷管及8种吸波材料涂覆方案的电磁散射特性,并获得了X波段下9种模型的后向RCS随探测角度的变化规律。研究结果表明:吸波材料涂覆可以有效地缩减球面收敛喷管的RCS;合理的涂覆方案可以在保证RCS缩减效果的基础上,降低吸波材料的使用量;相比全涂覆方案,仅在喷管出口和球面段进行涂覆,可以在吸波材料减少30%使用量的情况下,达到全涂覆方案80%的缩减效果。