使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

超声速细长梯形翼背风面流型

给出了前缘后掠角65°，双弧形剖面的细长梯形翼背风面流动显示结果。实验马赫数为1．10、1．53、2．53、3．01和4．01，攻角范围5°～25°。借助于蒸汽屏、纹影和油流技术拍摄了脱体和表面流型照片。蒸汽屏显示表明：在机翼背风面三角形区域的脱体流型可在垂直于前缘的法向攻角和法向马赫数构成的坐标平面上，区分出七种不同的流型；在切尖区域，有侧缘分离涡形成，后缘拖出尾涡。从摄取的纹影照片与横截面上的蒸汽屏照片一起，可获得机翼弓形激波位置随马赫数变化，以及激波－诱导分离线位置随马赫数和攻角变化曲线。在机翼上表面通过油流显示出主再附线、二次分离线、二次再附线和侧缘涡区。显示出的流型与其它有关实验和数值计算结果比较，符合得很好。     

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

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

低雷诺数下50°后掠三角翼的旋涡流动

采用数值模拟和流动显示的方法研究了50°后掠角三角翼在低雷诺数下的旋涡流动,结果表明:低雷诺数下,非细长三角翼在5°攻角时就形成了稳定的前缘涡,较小攻角时前缘主涡就开始破裂,并观察到泡型和螺旋型两种旋涡破裂方式。另外,在一定的攻角范围内,前缘主涡的外侧又生成一对新的集中涡,构成双涡结构;随着攻角的增大,前缘涡涡核不断升高,主再附线向中心移动,二次分离区扩大。     

跨音速压气机转子中三维湍流流场计算及涡系分析

从叶轮机械转动坐标系中的三维、可压、湍流、Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组出发,用Ｂａｌｄ－ｗｉｎ－Ｌｏｍａｘ湍流模式使方程组闭合;结合有限体积离散并采用隐式矢通量分裂技术;求解了西德宇航院（ＤＦＶＬＲ）的单级跨音速压气机转子内的三维湍流流场。计算结果可分辨出该转子内三维流场的细微结构、激波的空间曲面;得到马蹄涡、通道涡、尾迹涡和角涡的形态及其发展,流场的三维性非常明显。计算得到的壁面极限流线图可反映出压气机中分离流动的拓扑结构     

涡对与物体的相互作用研究

本文从非定常形式的Ｎ－Ｓ方程出发，采用有限差分法数值模拟了入射涡对与圆柱及翼型相互作用的流场结构，旨在分析讨论其复杂的非定常流动现象，如二次涡及三次涡的形成，多涡配对，非定常分离及剪切层的演化等，以及不同的物形及涡对与物体间不同的相对位置对影响。此外，文中还通过分析物面压力分布随时间的变化来探讨非常分离的流动特性。     

计算带有层流分离的粘性/无粘流相互作用的一个方法

基于粘性/无粘流相互作用模型和边界层积分方程,推导出计算二维和三维不可压层流粘性/无粘流相互作用的方程组。该方程组可以用于带有中等程度的分离和再附情况下。文中给出了计算方法和算例。所得的二维计算结果与其他计算和实验结果相符。本法具有简便迅速的优点。     

圆柱尾流中的二次分离及环量衰减问题

在圆柱绕流的非定常流动中,相应于旋涡的周期性脱落柱后的二次分离是周期性间断发生。二次涡的影响不能忽略。但对长时间流动的模拟还必须考虑尾涡的环量衰减。模拟Bluff Body流动有效的离散涡模型应包括二次分离及尾涡衰减两种机因。本文数值计算结果证明了上述结论。     

超音速细长梯形翼背风面流型

给出了前缘后掠６５°、双弧形剖面的细长梯形翼背风面流动显示结果。实验Ｍａｃｈ数为１．１０,１．５３,２．５３,３．０１和４．０１,攻角范围为５°～２５°。应用蒸汽屏、纹影和油流技术拍摄了空间和表面流型照片。蒸汽屏显示表明：在机翼背风面三角形区域的空间流型随法向攻角αＮ（在垂直于前缘的平面内流速与弦线间的夹角）和法向Ｍａｃｈ数ＭａＮ（来流Ｍａｃｈ数在垂直于前缘平面内的分量）变化,并可在αＮ和ＭａＮ为坐标的平面上划分出７种流型存在的区域。侧缘区有侧缘分离涡形成;后缘有尾涡拖出。从纹影照片与横截面上的蒸汽屏照片对照可获得机翼锥面激波位置随Ｍａｃｈ数的变化;以及激波－诱导分离线位置随Ｍａｃｈ数和攻角变化曲线。机翼表面油流谱显示出了主再附线、二次分离线、二次再附线和侧缘涡区。显示出的流型与其他有关实验和数值计算结果比较符合得很好     

大迎角粘性复杂流场的数值模拟

本文使用于存在较大分离区的Ｊｏｈｎｓｏｎ－Ｋｉｎｇ湍流模型对多段翼型流场及某型歼击机外形全机流场的雷诺平均ＮＳ方程进行数值计算，研究了多段翼型及有三角翼的全机外形的最大升力和三角翼翼尖涡，二次涡等粘性特性。     

强粘性层流理论及在粘性流计算中的应用

对高雷诺数流动计算，为了解决部分粘性项甚至全部惯性粘性项落入误差与流动物理尺度的关系问题，本文提出强粘性流动理论。强粘性流中至少有一个粘性项与惯性项同量阶，理论包含了物理尺度各向极限、经典边界层和多层边界层理论为其特例，给出了从经典边界层向物理尺度各向相同极限演化的尺度规律和粘性惯性诸项变化的量阶关系，阐明了粘性与惯性力强相互作用将在剪切层的法向以流向同时“激发”涉尺度结构。对粘性流计算，利用强粘性剪切流尺度律重新标度NS格式的修正微分方程，给出临界网格尺度与流动物理尺度和差分格式精度的关系，得到部分粘性项落入误差和计算结果为非物理数值粘性解的二个判据。并以流场中的边界层、驻点和分离点领域计算为例说明理论的应用，对强粘性剪切流计算、证实部分粘性项甚至全部粘性惯性项落入误差的问题值得重视。     

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

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

基于N-S方程的跨声速翼型多目标多约束优化设计

本文将粘性流场分析和数值优化方法耦合起来，发展了一种跨声速翼型设计方法，用以提高翼型在一个或多个设计点、在多种约束条件下的气动性能。由粘性流场分析程序计算得到的升力、阻力等气动参数构成目标函数，数值优化程序对其进行最小化。粘性流场分析采用了雷数平均Ｎ－Ｓ方程，这比过去翼型设计中使用的全速势方程或Ｅｕｌｅｒ方程更能模拟流动的本质，因而设计结果的可靠性大大提高了。优化算法采用传统的拟牛顿法（Ｑｕａｓｉ     

球心非对称的钝圆头体近区绕流场气动特性分析研究

本文采用数值模拟对球心非对称的钝圆头体附近流场气动特性进行分析研究。数值求解三维可压缩NS方程，采用NND格式，以修正的Baldwin-Lomax代数模型对不同飞行状态进行数值模拟，分析了球心非对称钝圆头体对流场的影响，与部分试验比较结果满意。     

高分辨率格式与钝舵绕流数值模拟

 采用Beam-Warming/TVD格式，Roe二阶、三阶以及四阶通量差分分裂格式数值计算雷诺平均Navier-Stokes方程，湍流附加粘性采用Baldwin-Lomax模型计算，针对超声速绕钝圆舵流动进行了数值模拟比较研究，计算给出了流场的详细结构，物面压力分布与实验值符合良好，尤其是Roe三阶格式给出了较好的计算结果，两个超声速回流区以及湍流边界层分离也得到了很好的数值模拟。     

冲压加速器燃烧流场的数值模拟

本文利用ＥＮＮ格式,通过求解包含非平衡效应的完全Ｎ－Ｓ方程,数值模拟了超声速燃烧流动中很有兴趣的Ｈ２／Ａｉｒ斜爆轰波冲压加速器的绕流。     

超音速气流绕局部透气凹角流动的数值模拟

以具有压力分裂形式的简化Ｎ Ｓ方程为控制方程,数值模拟了超音速来流条件下的激波 边界层干扰被动控制（ｐａｓｓｉｖｅｃｏｎｔｒｏｌｏｆｓｈｏｃｋ ｂｏｕｎｄａｒｙｌａｙｅｒｉｎｔｅｒａｃｔｉｏｎ）。模拟是以预先给定激波前吹气和激波后吸气的流量来实现的。为了定性地确定吹气或吸气对激波 边界层干扰的影响,首先计算了单独吹气和单独吸气两种情况。数值计算时采用了多重扫描法对控制方程差分离散,以反映亚音速区压力对流场的椭圆性影响。     

一种高效的多重网格三维N-S方程计算方法及其应用

介绍一种基于 Jameson中心差分和 Runge-Kutta时间推进的三维 N-S方程计算的有限体积方法 ,计算点选取体积单元中心 ,为了提高收敛效率和程序稳定性 ,采用了当地时间步长和多重网格方法 ,特别是针对粘性计算中计算网格形状比高的特点 ,隐式残差光顺采用当地变系数的方法 ,并专门设计了考虑计算网格形状比的人工粘性模型。应用本文方法对收扩喷管的三维流场和NACA0 0 1 2翼型二维流场进行了计算 ,结果表明程序具有良好的收敛能力和较高的稳定性。     

Numerical simulation of multi-dimensional inviscid compressible flows by using TV-HLL scheme

This paper deals with the numerical solution of inviscid compressible flows. The threedimensional Euler equations describing the mentioned problem are presented and solved numerically with the finite volume method. The evaluation of the numerical flux at the interfaces is performed by using the Toro Vazquez-Harten Lax Leer(TV-HLL) scheme. An essential feature of the proposed scheme is to associate two systems of differential equations, called the advection system and the pressure system. It can be implemented with a very simple manner in the standard finite volume Euler and Navier–Stokes codes as extremely simple task. The scheme is applied to some test problems covering a wide spectrum of Mach numbers, including hypersonic, low speed flow and three-dimensional aerodynamics applications.     

边界层转换阶段直接数值模拟的影响矩阵法

提出一种求解时间相关、三维不可压完全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的影响矩阵法。考虑流动区域为上半空间的平极边界层流动,在水平方向,利用周期性条件,采用Ｆｏｕｒｉｅｒ离散,在半无穷方向,定义了一种特殊映像并利用Ｃｈｅｂｙｓｈｅｖ多项式展开。由此将Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程转化为以水平波数（ｋ１,ｋ２）为参数的一组法向Ｈｅｌｍｈｏｌｔｚ方程组,同时可更有效地求解压力。本文讨论了连续问题的影响矩阵法与原问题的等价性。影响矩阵法的离散形式,以及影响矩阵法的算法实现等。最后给出了一维Ｈｅｌｍｈｏｌｔｙ问题的数值实验结果。