二维结构/非结构混合网格的生成及其流场模拟

采用分区求解方法在二维结构/非结构混合网格中对粘性绕流作了数值模拟。利用层面推进法生成绕物体的O形结构网格,利用推进阵面法在结构网格外生成非结构网格。在结构网格区内,采用B-L湍流模型,求解雷诺平均N-S方程。在非结构网格区内求解Euler方程。两区的解通过在交界面上满足通量守恒条件进行耦合。数值算例表明了方法的有效性和实用性。     

运输机组合体绕流非结构混合网格N-S方程数值模拟

本文利用通用网格生成软件生成运输类飞机翼身组合体（DLR-F6）非结构混合网格，在边界附面层区域内，生成三棱柱形网格，在其他流动区域采用四面体网格进行填充。粘性流场求解器采用基于N-S方程的CFX5计算流体软件求解，通过计算结果与风洞试验所得到的数据对比分析，验证了网格生成方法及网格数据的正确性和实用性。本文研究结果也为运输类飞机高速风洞试验提供了对比数据。     

非结构网格生成方法的改进及气膜冷却三维静子流场的求解

讨论五面体网格下,涡轮静子叶片含多排冷却孔时三维粘性流场雷诺平均Navier-Stokes方程组的数值求解问题。考虑到叶片含多排冷却孔的工程背景,本文发展了生成五面体非结构网格的快速有效方法即仅在一个拟S₁面上生成非结构网格,而沿叶高方向仍采用结构网格的生成办法。文中还采用了新的网格优化措施以改善生成网格的质量。流场的计算采用了非结构网格下的高分辨率迎风格式。几个典型算例表明:本文发展的网格生成方法快速有效,而且网格生成质量较高;发展的计算格式能够获得较为满意的三维流场数值解,并具有较高的激波分辨率。      

基于粘性非结构网格的三维增升装置N-S方程数值模拟

应用基于粘性非结构网格求解N—S方程方法数值模拟了某型飞机三维增升装置的绕流。在物面附近粘性主导作用区域内，采用三棱柱型非结构网格，在其他流场区域，采用传统四面体非结构网格的混合非结构网格生成技术。湍流模型采用Spalart—Allmaras一方程模型，运用格心格式有限体积法，对控制方程和湍流模型方程进行数值离散。用四步龙格-库塔方法作显式时间推进，采用当地时间步长、隐式残值光顺等加速收敛措施．求解某型飞机全机带增升装置构型的粘性绕流流场，其数值计算与风洞试验结果显示了良好的吻合度。     

基于Euler方程结合附面层修正方法的应用研究

采用欧拉加附面层修正的方法在结构网格上对民用飞机翼身组合体流场进行了数值模拟研究.采用一种适合翼身组合体模型的的网格生成技术来生成计算用的结构网格,该技术是先采用一种基于Higenstock源项修正的椭圆型偏微分方程生成二维网格,然后沿机身轴向合理的光滑性调整成三维网格.在求解Euler方程时采用中心离散格式和显式时间推进,并利用附面层修正的方法来考虑粘性的影响.为了验证此方法的正确性,本文对DLR-F4翼身组合体进行了数值模拟,结果表明此方法和试验结果能够很好的吻合.     

多段翼型粘性绕流计算研究

针对多段翼型外形结构复杂,结构网格生成难度大的缺点,运用割补法对重叠网格技术进行了研究,生成贴体及与边界正交的高质量的结构化网格,在识别贡献单元方面,介绍了一种基于kd树的找重方法.流场求解应用雷诺平均Navier-Stokes方程结合Spalart-Allmaras一方程湍流模型在所生成的重叠网格上对多段翼型的粘性绕流进行了模拟,为了兼顾流场计算效率问题,在重叠网格的各个子块上采用多重网格方法.计算结果表明,用该方法可以很好地预计压力分布,并且计算效率得到很大的提高.     

任意平面域的三角形网格和混合网格生成

本文以阵面推进法为基础,给出了一种可对任意平面域进行三角形网格和混合网格自动生成的方法。该方法包含了易于生成与调控的矩形背景网格技术,简便的外形输入和初始阵面划分技术,加速查寻的数据结构,以及节点松弛和基于Ｄｅｌａｕｎａｙ准则的对角线交换技术。网格生成的实例表明,方法能快速可靠地生成任意平面域的高质量三角形网格和适于粘性计算的三角形／四边形混合网格,具有较好的通用性与实用性。     

非结构网格上N—S方程解法研究及在多段翼型绕流中的应用

在非结构网格上求解N－S方程时，靠近物面的网格伸展比相当大，用有限体积法进行显式求解时不易收敛并且雷诺数有一定的限制。本文采用的人工耗散项中含有网络的方向尺度，因而在不同区域选取不同：在粘性区内该项比物理粘性项小得多，在粘性区外与Jameson的人工耗散项相同。生成非结构网格采用的是层推进技术与传统的阵面推进法相结合的方法，这样能保证在粘中性区内网格有一定的方向规则，而在粘性区外网格完全非结构化。     

三维非结构贴体直角网格的N-S方程数值模拟

基于八叉树数据结构,采用直接分层光顺的方法,解决了非结构直角网格在物面不正交贴体的难点,发展了任意复杂外形的三维非结构贴体直角网格,并且可以根据流场特征进行自适应处理。在此基础上采用最小二乘法对N-S方程进行了MUSCL有限体积法的非结构网格的粘性求解,提出了带约束边界条件的最小二乘法的处理办法。计算表明,三维非结构贴体直角网格的粘性求解,可以适应各种高超声速的复杂外形,计算结果与试验吻合较好。     

直升机旋翼悬停流场的粘性数值模拟

在旋转坐标系下求解可压缩NS方程,对Caradona-Tung旋翼的跨声速悬停流场进行了数值模拟,计算中采用嵌套网格方法,并且应用了多重网格技术加速流场收敛。计算结果表明,多重网格技术显著改善了流场的收敛效率,计算量仅为单重网格的四分之一。文中同时还对无粘和粘性流场进行了对比,NS方程考虑了物面的粘性作用,能正确模拟出激波和边界层的干扰现象,计算得到的桨叶表面压力分布与实验值更加吻合。     

指数信号模型中参数的稳健估计

 给出了阻尼指数信号模型中参数的M 估计和Bootstrap M 估计; 基于Prony 方法给出估计的重新定义和新算法; 并且在较弱的条件下证明了这些估计的相容性; 最后利用Huber 函数给出了两种估计的模拟结果, 结果表明这两种估计的稳健性和精度都比传统的最小二乘估计好。     

大后掠翼大迎角定常与非定常俯仰气动特性及其控制

 对装有“前端襟翼”和“前缘襟翼”的74°后掠三角翼, 在迎角0°～ 90°范围的定常和非定常 (俯仰振动, f = 0. 4, 0. 8, 1. 2Hz) 的空气动力特性以及襟翼的控制作用作了实验研究。定常流的实 验结果与解N 2S 方程的计算结果作了比较。非定常俯仰振动使气动特性的变化出现迟滞现象, 但 仍能保持襟翼对定常流控制所具有的优越性。不装襟翼的三角翼的部分校测结果与文献中相应数 据一致, 表明实验结果可信。     

飞机/外挂物组合体绕流的数值计算

采用嵌合体埋入方法计算带外挂翼-身组合体复杂外形的流动,将飞机/外挂物外形的流动区域分为机翼-机身和外挂物两个子区域,在机翼-机身子域上求解欧拉方程,在外挂物子域上求解全速势方程,通过两个流场的耦合迭代计算全流场的解。对C-5A飞机及某战斗机的计算表明,数值模拟与实验吻合。实现了不同子域流场使用不同方程的计算方法。     

超声速混合层液滴两相流动的大涡模拟

为研究超声速气流中液滴与气流的混合及液滴蒸发对混合的影响,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构,气相流场采用亚格子(SGS)模型和切应力输运(k-ωSST)湍流模型,液相模拟采用轨道模型和单液滴蒸发模型。在混合层前缘入口处均匀持续地投放液滴,并在液滴入口处下方添加非周期小扰动,并观察液滴蒸发过程对该小扰动造成的影响。分析了入口小扰动在流场中不同的发展情况,发现液滴的蒸发过程使混合层厚度增加并加速混合层的发展,对气相流场扰动较强,可能导致流动失稳,对混合过程有很大的促进作用。     

NACA0012翼型动态失速特性和测压方法的研究

全面介绍了应用测压方法,对ＮＡＣＡ００１２翼型,在俯仰正弦振动的情况下,进行广泛实验研究的情况和结果。实验时对模型的俯仰振动频率、雷诺数、平均迎角等的影响进行较为广泛的研究。此外,还应用两种动态数据采集方法——扫描法和采样保持法,在相同的实验条件下,同时进行数据采集,对动态测试系统截止频率也进行了广泛的实验研究,利用动态测试系统的软件包,对数据进行频谱分析,在大量实验的基础上,对测压方法和两种数据采集方法进行对比与分析     

跨声速翼型上激波／边界层干扰的自适应控制计算

在激波区使用自适应壁对跨声速翼型上的激波/边界层干扰进行控制,可改变机翼的气动性能。这种被动控制可通过在翼型的激波区开一凹腔,其上覆盖一弹性橡胶膜柔壁来实现。本文给出用N-S方程数值模拟这一自适应控制翼型的跨声速粘性绕流,提出了一个适用于本特殊情况(物面边界局部地区在求解过程中有变化)的处理办法。并探讨了自适应柔壁对当代跨声速翼型绕流的影响。     

绕三角翼纵向俯仰大迎角气动特性计算研究

 采用数值计算方法, 对三角翼从0°上仰至90°的动态流场结构进行了计算, 在此基础上, 对三角翼在上仰过程中受到横侧小扰动情况下的流场结构和气动力特性进行了计算研究。给出了三角翼纵向动态情况下的气动力系数变化, 特别是大迎角横侧力矩系数的变化特征, 并对受到横侧小扰动后横侧运动的稳定性进行了计算与分析。结果表明, 机翼的上仰运动延迟了机翼上翼面旋涡的破裂。同时, 随着机翼俯仰角速度的提高, 机翼抵抗旋涡非对称破裂的能力明显增强, 机翼运动的稳定性也明显提高。     

基于结构动网格的无人机地面效应研究

为了更好地分析地面效应对无人机气动特性的影响,采用一种新的结构动网格生成技术,运用CFD手段对某无人机的地面效应进行了研究。在一套静态网格CFD计算基础上,运用结构动网格生成技术,获得非定常的网格移动,得到一系列的不同攻角、地面效应高度的CFD网格及初始流场。在此基础上进行定常流场求解,计算得到地面效应的定常CFD结果。最后对某双尾撑布局无人机的地面效应进行了CFD计算研究,对相应的全机气动特性的影响进行了分析。计算结果表明：地面效应对平尾的影响大于对机翼的影响并使平尾的失速迎角提前;该方法是可行的,具有一定的工程应用价值。     

Linear stability analysis of interactions between mixing layer and boundary layer flows

The linear instabilities of incompressible confluent mixing layer and boundary layer were analyzed.The mixing layers include wake,shear layer and their combination.The mean velocity profile of confluent flow is taken as a superposition of a hyperbolic and exponential function to model a mixing layer and the Blasius similarity solution for a flat plate boundary layer.The stability equation of confluent flow was solved by using the global numerical method.The unstable modes associated with both the mixing and boundary layers were identified.They are the boundary layer mode,mixing layer mode 1 (nearly symmetrical mode) and mode 2 (nearly anti-symmetrical mode).The interactions between the mixing layer stability and the boundary layer stability were examined.As the mixing layer approaches the boundary layer,the neutral curves of the boundary layer mode move to the upper left,the resulting critical Reynolds number decreases,and the growth rate of the most unstable mode increases.The wall tends to stabilize the mixing layer modes at low frequency.In addition,the mode switching behavior of the relative level of the spatial growth rate between the mixing layer mode 1 and mode 2 with the velocity ratio is found to occur at low frequency.     

压气机转子内尖区近失速状态三维流场

在低速大尺寸压气机实验台上 ,借助于旋转四坐标探针位移机构 ,用锥形五孔探针测量了压气机近失速状态下 ,转子叶片通道后段尖区内的三维流场。测量结果表明 ,吸力面附面层径向潜移强烈 ,并出现气流分离 ,在尖区的近吸力面区域形成一个旋涡 ;压力面角区存在刮削涡 ;在叶尖槽道中部 ,吸力面角区附面层与压力面角区附面层气流掺混 ,造成高损失和高阻滞。所有这些构成了尖区的复杂流动