用N-S方程计算三维叶栅流

发展了基于微元中心有限体积方法,采用雷诺平均N-S方程求解三维叶栅流动的程序,特别注意了格式中人工粘性的加法;程序采用Baldwin-Lomax模型。格式经过平板层流及湍流流动的校验,层流与Blasius解、湍流与实验经验公式及Wilcox的计算值符合很好,最后以Goldman涡轮导向器为例进行了叶栅流动计算校核,结果符合于实验以及Chima等人的计算。      

RAE2822翼型跨音速流动CFD计算的可信度分析

采用格点格式的有限体积法对RAE2822翼型的气动特性进行了粘流数值模拟,流动模型为雷诺平均NS（RANS）方程。在此基础上系统研究了不同网格类型以及不同湍流模型的计算精准度。通过与经典的试验数据进行确认研究,得出了适宜二维跨音速粘性流动计算的网格拓扑分布规律,并总结出了适宜的湍流模型。     

高阶精度方法下的湍流生成项对低速流动数值模拟的影响研究

基于雷诺平均的Navier-Stokes方程和结构网格技术,采用五阶空间离散精度的WCNS格式,开展了SST两方程模型不同湍流生成项对低速流动数值模拟的计算分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR7301两段翼型和Trap Wing高升力构型,研究内容主要包括不同湍流生成项对收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布、气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于低速二维流动,不同湍流生成项对收敛历程有比较明显的影响,对附面层湍流粘性系数分布和速度型影响不明显,不同湍流生成项主要影响主翼前缘的吸力峰值,进而影响升力系数和压差阻力系统;对于低速三维流动,不同湍流生成项对低速流动的收敛特性影响不明显,对翼梢涡的模拟精度有比较明显的影响,进而影响翼梢站位的压力分布和总体气动特性。     

一种新型多维通量函数在求解平面叶栅跨音速流动的应用

采用一种改进的Ｒｏｅ格式对平面叶栅的无粘流动进行了数值模拟,并与Ｒｏｅ格式的计算结果进行了比较。采用四阶Ｒｕｎｇｅ－Ｋｕｔａ时间差分格式,空间三阶有限体积法。网格生成采用贴体网格近壁面修正的方法。计算结果表明,该格式在减少格式对网格的依赖性,抑制数值粘性方面具有明显的作用。     

高阶精度方法下的湍流生成项对跨声速流动数值模拟的影响研究

利用五阶空间离散精度的WCNS格式和多块结构网格技术,通过求解雷诺平均NS方程,开展了SST两方程模型不同湍流生成项组合方式对跨声速流动数值模拟影响的计算分析。研究的主要目的是为高阶精度格式在复杂外形上的工程应用提供技术支撑。计算模型采用了RAE2822超临界翼型和DLR-F6翼身组合体构型。研究内容主要包括不同湍流生成项对残差收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布以及模型整体气动力特性的影响。不同湍流生成项组合方式的流场计算结果还与风洞试验数据进行了对比。研究结果表明:对于小迎角不存在明显分离的跨声速流动,不同湍流生成项对流场的高精度计算结果的影响很小,可以不用考虑。     

带有分离的二维可压粘性喷管流动计算

采用二维、时变、可压纳维-斯托克斯(N-S)方程计算喷管流动。流动可为无粘、层流或湍流。湍流模型可采用代数混合长模型,一方程和二方程模型。喷管形状为二维(平面或轴对称)。计算网格内点采用MacCormack显式格式。采用局部细网格,快速解法以及时间间隔光滑方法,以加速高雷诺数下收敛。边界点采用参考平面特征线格式计算,粘性项作为源项处理。给出了初步算例结果,与其它计算和实验数据比较,符合良好。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

复杂回旋流场的数值模拟和实验研究

本文通过四角喷射炉膛内回旋流场的数值计算和实验测试的对比研究，比较了几种计算方法的预估结果与实验结果的符合程序。结果表明：为克服当流动方向与网格线交角过大时所产和珠过分数值扩散，引入流向平衡粘性方法与κ－ε湍流模型相结合可获得与实验结果相一致的计算结果，而采用混合分格式加κ－ε模型所得到的计算结果比流向平衡粘性方法加等效常湍流粘性系数这样的简化处理的计算结果还要差。     

飞船高超声速粘性绕流的数值模拟

本文采用稳式ＮＮＤ格式，通过求解三维薄层近似Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了飞船高超声速粘性流场。在计算中充分考虑了网格雷诺数和网格分布对计算结果的影响，不但得到了满意的表面压力、热流分布和流场的物理量分布，而且得到与实验观察一致的表面油流图像和分离结构。     

应用不同湍流模型计算叶栅流场的结果与比较

采用有限体积TVD-LW格式求解N-S方程,使用了B-L模型、高雷诺数k-ε模型和低雷诺数k-ε模型等多种湍流模型。并且用具体算例说明本文所用程序的可靠性,比较了3种的湍流模型的特点。B-L模型最大的优点是简单、使用方便。k-ε方程适应性强,分为低雷诺数k-ε方程和高雷诺数k-ε方程两种。高雷诺数k-ε方程比低雷诺数k-ε方程精度略低,但所需计算时间和存储量都比后者小,网格也可以粗些。k-ε模型对粘性层、分离区有比较高的分辨能力,高雷诺数k-ε模型比较适于一般数值计算。      

旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

喷管化学反应粘性流场的数值模拟

用时间相关法的显格式求解氢氧发动机喷管的粘性化学反应流动,采用６种组分,８个基元反应有限速率的化学反应模型,湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流动参数和产物的质量分数在喷管流场中的分布。计算结果表明,采用显格式方法求解喷管化学反应粘性流动是可行的。     

离心叶轮内稀释固液两相流动的数值模拟

 对固液两相流动进行了数值模拟。液相和固相分别采用 k- ε湍流模型和代数颗粒流( AP)湍流模型并考虑了液固两相间的相互作用。用此模型结合 SIMPLE- C算法对离心叶轮内的稀释固液两相流动进行了计算并研究了颗粒对流场的影响。计算结果同试验结果基本一致     

喷管化学反应流的对角化点隐式法求解

采用弱耦合对角化点隐式方法的ＭａｃＣｏｒｍａｃｋ格式,求解喷管的粘性化学反应流动,解决了控制方程组的刚性问题,又大大减少了常规点隐式方法中求逆矩阵的大量工作。对液氢液氧火箭发动机,采用６种组分、８个反应有限速率的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流场参数在喷管中的分布。计算结果表明,采用对角化点隐式方法求解化学反应喷管粘性流动能提高计算效率。     

翼型风洞侧壁干扰的数值模拟研究

运用Navier-Stokes数值模拟对翼型模拟试验对风洞侧壁干扰进行模拟，将简单代数湍流模型扩展用于机翼／风洞侧壁结合区流动，分析了风洞实验侧壁干扰问题的形成机理和影响翼型实验侧壁干扰的各种因素，如翼型展长、风洞侧壁边界层厚度及侧壁边界层抽吸等，对实验结果的影响，得出了一些有用的结论。计算格式空间采用中心有限体积离散，时间采用多步Runge-Kutta时间步长格式进行积分。计算结果证明了该方法的可行性和优越性。     

跨音速压气机级三维紊流流场数值模拟

在对跨音速压气机级的动／静双叶排中的三维定常紊流流场数值的模拟中，利用三阶高分辨率ＮＮＤ格式和ＬＵ－ＳＧＳ隐式推进迭代法，既保证了流场中激波的模拟质量和粘性流动特征的正确预估，也实现了求解过程的高效率。对某单级跨音速压气机的计算及其与实验数据的对比证实了本文方法的有效性。将计算推广至多级轴流叶轮机械的情况将会是十分简便的     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

超声速开式空腔流动的数值模拟

研究了不同长深比、不同雷诺数条件下,超声速开式空腔流动特性的变化。算法采用空间上5阶的WENO格式,时间上具有TVD性质的3阶Runge-Kutta方法,湍流模型选用B－L模型。数值模拟获得了不同的时刻的注场图像,压力波、涡以及空腔结构之间的相互耦合作用,使得空腔内存在典型的周期性压力脉动,但不同长深比空腔的流动性存在较大差别。     

用粘性体积力方法计算高速离心压气机叶轮内部流场

高速离心压气机是用途非常广泛的一种动力机械,充分了解其内部复杂的三维流场对设计高压比、高效率离心压气机具有很重大的意义。本文采用一种简化的Navier-Stokes方程的计算方法-粘性体积力方法,计算了Krain叶轮内部流场。使用H型网格离散微分方程,并采用多重网格方法和当地时间步长方法加速计算收敛。粘性应力通过Baldwin-Lomax模型计算。为了证明计算结果的有效性,文中也给出了相应的激光测量结果。      

梯度加权局部分析解差分格式及其应用

提出了一种新的针对Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程符合含有梯度项的对流扩散方程的梯度加权局部分析解的有限差分格式，且给出了通用的系数表达式。ＧＷ－ＬＡＤ格式是基于积分型的有限体积法并且通过附加源项的方法来实现。能够适用于粘性和无粘性及不可压缩和可压缩的压力驱动流动。通过对叶栅流道中二维湍流流的计算表明ＧＷＬＡＤ格式与ＥＤ格式相比，具有较高的精度和很快的收敛速度并且得到更好的压力场和速度场。