固壁夹道中多体Stokes绕流的数值分析

本文利用格林函数方法,数值求解了固壁夹道中多体的Stokes绕流问题,计算了两平板间的二圆柱体的绕流,三圆柱体的绕流,三正方体的绕流以及五个不规则形状物体的绕流,对其流动特性进行了分析。表明格林函数方法是目前对这类流动问题进行数值分析的最有效方法。     

跨流域三维复杂绕流问题的气体运动论并行计算

通过研究求解描述跨流域三维绕流问题的Boltzmann模型方程气体运动论耦合迭代数值格式,分析气体运动论数值计算方法的内在并行性;从变量依赖关系、数据通信与并行可扩展性三方面开展基于离散速度空间区域分解计算的研究,发展求解稀薄流到连续流跨流域三维复杂绕流问题的并行算法。通过对不同Knudsen数、不同马赫数、不同攻角三维球体及返回舱绕流的并行计算,计算结果与实验数据和理论分析吻合较好。研究表明,该并行算法负载平衡和并行可扩展性较好,对不同并行计算机系统具有很好适应性,显示直接求解分子速度分布函数的气体运动论计算方法有良好的并行计算特性。     

欧拉方程组的隐式近似因子分解差分方法的研究

本文对求解二维欧拉方程组的隐式近似因子分解差分方法进行了简要的分析,并利用它计算了简单直扩张通道内的亚音速、超音速和跨音速管流、绕无限长拐角和对称棱形翼型的超音速绕流以及绕NACA 0012翼型的亚音速、跨音速和超音速绕流等典型定常流动。计算结果表明,它具有良好的适应性和精度,可以正确地捕获强激波和自动满足Kutta条件,Courant数可以取大于1的数值,从而加速收敛。分析和计算都是在贴体坐标下进行的。方法本身可以推广到三维流的计算。     

平面叶栅跨音绕流时间推进法中特征相容条件在边界点上的应用

在平面叶栅定常跨音绕流计算中,采用时间推进法,由于基本方程组是非定常欧拉方程组,属于双曲型,可以应用特征理论。在边界点上,特征相容方程精确地刻画了流动参量间的关系,使得数值处理更精确地反映无粘流的物理模型。本文从基本方程组出发,推导出适合于边界点应用的特征相容方程组的通用形式;并应用特征理论,讨论了平面叶栅绕流问题中边界条件的提法;对于边界条件的几种不同的近似处理方法,进行了数值试验,比较了它们对数值结果的影响。     

SLAU格式在大涡模拟方法中的应用

采用格心格式的有限体积法求解三维大涡模拟控制方程,成功地将最新改进的SLAU格式推广到大涡数值模拟方法中。SLAU格式的空间重构方式采用五阶WENO格式,时间推进采用LU-SGS隐式时间推进。湍流模式为最新提出的基于随机分析的动态亚格子模型。计算了圆柱绕流、NACA4412翼型绕流、NACA0012翼型绕流三个算例,分别从可压缩较大迎角翼型绕流、可压缩跨音速翼型绕流以及不可压缩低马赫数圆柱绕流三个方面验证了将SLAU格式应用到湍流大涡数值模拟计算中的可行性。结果表明,结合SLAU格式的大涡模拟方法可以达到较好的数值模拟效果,同时间接检验了基于随机分析的动态亚格子模型的数值模拟效果。     

模拟二维周期性非定常流动的谐波分析方法

介绍了一种从频域出发模拟二维周期性非定常翼型绕流的数值方法.考虑到很多工程中遇到的非定常流动问题都是周期性的,因此非定常流场中的守恒变量可以展开为傅立叶级数,其系数随空间位置的不同而改变.将这样的谐波假设引入欧拉方程,经过处理的方程不显式地含有时间导数项,实际求解的方程将具有定常方程的形式.在使用多重网格技术后,计算可以很快达到收敛.对NACA64A010翼型典型的周期性非定常绕流进行的数值模拟表明,这种方法在能够给出较为准确的计算结果的同时,与传统的双时间方法相比,可以大幅提高计算效率.     

基于混合RANS/LES方法与FW-H方程的气动声学计算研究

在气动噪声的数值计算过程中,非定常流动的求解精度对声学计算结果有着重要的影响.以机体噪声计算的标准算例双圆柱绕流气动噪声问题为研究对象,采用基于非线性k-ε湍流模型的限制数值尺度(LNS)方法对双圆柱绕流进行了数值模拟,将计算得到的气动特性和流动特征与相应的试验结果进行了对比分析.为了求解远场观测点处的气动噪声,在精确求解双圆柱绕流流动的基础上结合基于FW-H(Ffowcs Williams-Hawkings)方程的声类比方法进行数值计算,并通过圆柱体的展向相关性对计算结果进行了修正,将得到的最终结果与相应的声学试验结果进行了对比,两者吻合良好,表明该数值方法是准确、可靠的.     

电磁力控制翼型体绕流的数值研究

为考察包覆电磁激活板的翼型体绕流特性,采用数值模拟的方法研究了低雷诺数下电磁力对翼型体流场结构和绕流过程的影响,探讨了不同迎角和不同电磁力大小对翼型体尾流涡街形态和流体边界层结构的影响。计算结果表明,作用于翼型体表面的电磁力可以明显改变流体边界层的结构,抑制边界层分离,并有效提高翼型体的升力。     

壁面喷流可控环量机翼的数值计算

本文提出一种适用于可控环量机翼研究的势流/粘性流的数值解法,包括钝体分离绕流的一种势流解,以及改进的Spalding统一边界层理论的解。通过对具有壁面切向喷流可控环量圆柱绕流和类椭圆翼绕流所作出的数值计算表明,在其计算模型确定后,就不必依赖实验数据来模拟实际问题,计算结果与实验数据符合较好。     

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

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

NND格式在航天飞机头部段N—S方程求解中的应用

本文采用实质上为二阶精度的无波动、无自由参数的耗散差分格式(NND格式)求解了非定常的二维和三维完全N-S方程,对类似于航天飞机头部外形的超声速粘性气体绕流流场进行了模拟。在计算中,使用了代数网格生成技术,对脱体激波采用了装配法。作为验证算例的二维圆柱绕流计算结果同文献[4]的数据进行了比较。二维凹陷外形绕流的计算结果同无粘流结果进行了对比,由于粘性流动出现了分离,二者的壁面压力是有差异的。最后,对三维外形的绕流也进行了模拟,计算结果提供了流场的细节。     

壁判据用于计算流体力学（CFD）可信度评估

本文把作者提出的近壁干扰剪切流动(ISF)全域理论与流体运动方程组及流体在壁面上无滑移条件相结合导出一组壁面判据.壁判据为计算流体力学(CFD)仿真的可信度评估提供了基于流体理论的一条直接验证途径.对不可压缩流动的十一个熟知的NS方程组精确解,包括二维驻点和斜入射再附点流,二维分离点和背风驻点流,轴对称驻点和背风驻点流,旋转圆盘附近的三维Von Karman 流、收缩和扩张渠道流和非定常斜入射三维驻点流;以及经典边界层及其无粘外流和相似性边界层及其粘外流的NS方程组解,提出可用于验证近壁流动计算的几个壁相关函数.证实它们准确满足所有壁判据;说明壁判据可用来检验NS方程组数值解近壁计算结果的计算精度并验证其可信度.     

回流环形燃烧室反应流场数值模拟

采用自编计算程序,在任意曲线坐标系下对带扩压器、涡流杯、火焰筒、内外环冷却通道与弯管的回流环形燃烧室三维两相化学反应整体流场进行计算.数值预测火焰筒头部结构参数对燃烧流场的影响.所得计算结果与实验数据符合较好,结果表明:采用的数学模型和计算方法合理,该软件可供回流燃烧室优化设计与研制使用.      

轴对称矢量喷管外流对内流干扰研究

应用有限体积法在同位网格中对控制方程进行离散化,对某一轴对称矢量喷管的内外流场进行了数值模拟,数值模拟和试验结果进行了比较,分析了该喷管在矢量偏角20°下,不同的落压比、不同的外流速度对轴对称矢量喷管内流的影响。研究结果表明,喷管内的分离激波是影响轴对称矢量喷管性能的重要因素。      

三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。     

封严间隙对盘缘密封性能影响的数值研究

通过对真实工况下高压涡轮的盘缘封严进行三维数值模拟,研究了3种不同轴向间隙和径向间隙对封严效率和盘腔内流场变化特性的影响规律,结果表明:针对复合封严结构及研究参数范围,由于静盘封严环将盘腔划分为封严外容腔和封严内容腔,而封严外容腔可以有效地将主流燃气阻隔在内,因此封严内容腔具有较高的封严效率。当轴向间隙增大时,封严环轴向重合度会逐渐减小,同时冷气出流阻力降低,盘腔的封严效率因而提高。当轴向间隙一定径向间隙逐渐减小时,封严外容腔径向距离减小,此时盘腔低半径处封严效率升高。     

高超声速粘性激波层高雷诺数流场的数值求解方法

本文将拉伸变换引入Hosny的高超声速粘性激波层的数值求解方法中,数值计算结果与Little所给出的实验结果相符合。从而不论在雷诺数范围上,或在后身流场的推进计算上,本方法均比Hosny原方法有明显改进。     

跨声速条件下轴对称收扩喷管内外流场的数值研究

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用不同湍流模型、壁面函数对轴对称喷管内外流场进行三维数值模拟.计算结果与实验数据对比表明,采用涡粘修正的RNG模型,内流场使用增强壁面函数,外流场使用标准壁面函数可以得到准确的计算结果.研究表明:在亚声速条件下,喷管外壁面压力分布随落压比的变化不大;超声速条件下,喷管外壁面激波随着落压比的增大前移.      

在电离层飞行的航天器近尾迹内电子数密度渐近解

本文讨论了航天器近尾迹流场电子数密度分布,给出无碰撞假设下等离子体绕带电圆盘流动的近尾迹流场渐近解。用匹配渐近展开法求解实际可能现出的三种情形的电子数密度分布。所得的摄动解对近尾迹流动的理论和数值研究有指导意义。     

数值研究涡流器对环形燃烧室燃烧性能的影响

数值研究不同类型涡流器对环形燃烧室流场与性能的影响.采用分区耦合方法分别生成两级轴向、斜切径向及三级轴向涡流器的环形燃烧室的结构化网格,在任意曲线坐标系下对带有扩压器、涡流器、火焰筒和内外环冷却通道的三种环形燃烧室三维两相燃烧整体流场进行计算,气相和液相分别采用Euler与Lagrange法处理,采用多维/经验分析法估算燃烧室性能.计算与实验比较相符表明,可利用本文计算方法预测不同涡流器对燃烧室流场和燃烧室性能的影响.