影响冷气掺混数值模拟精度的若干问题分析

对冷气掺混数值模拟中的计算网格、湍流模型、射流边界条件等影响精度的若干问题进行分析。在吹风比0.5下,对不同网格分布、湍流模型、射流边界条件进行数值模拟,所得结果与Ajersch的实验结果进行对比。结果表明,出口下游的网格分布可适当稀疏以减少计算量;在所研究的几种湍流模型中,k-ε模型所得结果与实验结果吻合得最好;考虑射流通道内流动能提高精度,在射流出口给定流场分布也可保证计算结果的精度。     

三维槽道湍流的大涡模拟

利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。为了比较网格和亚格子模型对计算结果的影响,采用了三种不同的网格数目,以及Smagorinsky模型和动力模型等两种亚格子模型,并与直接数值模拟结果进行了比较。结果表明,大涡模拟可以比较准确地给出速度分布、脉动速度均方根分布等湍流统计参数。在此基础上,对近壁区小尺度湍流结构进行了研究。     

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

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

液体火箭发动机推力室冷却通道流动与传热数值研究

采用气固耦合算法对液体火箭发动机推力室再生冷却通道的流动与传热过程进行了三维湍流流动与传热数值模拟，冷却工质为氢气，其密度、导热系数、动力粘度随着温度和压力而变化。应用大涡模拟及标准k-ε双方程模型两种湍流模型分别进行数值模拟，详细揭示了再生冷却通道固体区和流体区内的速度场和温度场，并在不同的计算网格数目下对两种湍流模型的计算结果进行了对比。结果表明，在相同的网格条件下，标准k-ε双方程模型与实验数据的吻合精度比大涡模拟模型更好，且满足工程计算精度。随着网格数的增加，大涡模拟的计算精度逐渐得到改善。     

应用混合网格分析复杂积冰翼型气动性能

通过发展一种点-点对接的结构/非结构混合网格生成方法,避免了复杂积冰外形难以划分结构化网格的问题,并对复杂积冰翼型的气动性能进行了分析计算.在靠近积冰边界的内层采用了非结构网格以拟合复杂的边界,在非结构网格外采用结构化网格以节省存储空间和计算时间.为了检验网格对选用湍流模型的影响,整个流动区域分别采用SST和SA湍流模型,求解了雷诺平均N-S方程.数值计算结果表明SST模型更适于模拟复杂分离流动,积冰对翼型的气动性能造成了严重的影响.     

不同矢量偏角下二元喷管内流特性的数值研究

应用多重网格MacCormack显式格式对不同矢量角下二元喷管的内部流动进行了模拟分析。计算时, 在细网格上采用了雷诺平均N-S方程与k-ε湍流模型方程;在粗网格上则采用Euler方程。数值计算给出了较清晰的流谱, 壁面压力分布的计算也与实验结果相吻合。通过计算结果分析了矢量喷管几何参数对流场的影响, 并对矢量喷管的推力特性进行了数值研究。      

湍流模型离散精度对数值模拟影响的计算分析

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS),通过改变物面法向第一层网格间距,开展了剪切应力输运(SST)两方程模型不同离散精度的数值分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR 7301两段翼型和高速RAE2822翼型,研究内容主要包括湍流模型的二阶精度离散和五阶精度离散两种方式对收敛历程、边界层湍流黏性系数分布、边界层速度分布、压力系数分布以及气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于不同的第一层物面法向网格间距,湍流模型离散精度对低速绕流计算结果有比较明显的影响,对于高速小迎角附着流动计算结果影响不明显;相对于湍流模型二阶精度离散,湍流模型高阶精度离散网格敏感性较弱,具有更高的数值模拟精度,但收敛性略差。     

基于S-A与SSG/LRR-ω两种湍流模型的CHN-T1标模计算与分析

通过求解雷诺平均Navier-Stokes方程,对第一届AeCW-1研讨会上要求的所有工况进行了计算,着重研究S-A湍流模型和雷诺应力模型在CHN-T1标模上的表现,评估两种湍流模型的网格收敛性,对抖振特性的模拟能力以及对雷诺数的敏感性。根据与实验数据的对比分析,得出主要结论如下:两种湍流模型都能趋向于网格收敛,雷诺应力模型的网格收敛更快。在较大迎角的工况,两种湍流模型对机翼上激波诱导的分离流有较大的差别。CHN-T1标模的尾撑,改变了飞机平尾的压力分布,并且在飞机机尾上产生复杂的三维分离流动,对俯仰力矩系数的准确计算产生了较大的影响。考虑静气弹效应,两种湍流模型计算的力矩系数比不考虑静气弹效应更大,雷诺应力模型的偏差更大。在高雷诺数工况,两种湍流模型的计算结果趋于一致。     

燃烧室内有旋射流的数值计算

本文用N-S方程描述燃烧室内单相有旋气体喷射的轴对称非定常湍流流动。计算了不同Swirl数的有旋射流进入突放管的速度场和流线分布。湍流模型采用SGS方法。计算结果与其它模型所得结果相比,具有所需网格少、精度高和计算方便的优点。     

高阶精度格式WCNS在三角翼大攻角模拟中的应用研究

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟了65°三角翼的大攻角绕流流场,主要目的是考核高阶精度格式WCNS在大攻角旋涡流动方面以及跨声速流场的激波附面层干扰、涡破裂位置的模拟能力,重点研究不同网格规模和湍流模型对尖前缘三角翼涡系之间的相互作用的影响。通过求解任意坐标系下的雷诺平均N-S方程,采用5阶精度的加权紧致非线性格式(WCNS-E-5)和多块对接结构网格技术,两种湍流模型分别是一方程SA和两方程SST湍流模型,在与相应试验结果对比的基础上,详细研究了WCNS-E-5格式在跨声速大攻角旋涡流动中的表现,以及不同网格规模、两种湍流模型对主涡二次涡相互作用、涡破裂位置和表面压力分布的影响。本文的研究结果表明,高阶精度格式WCNS-E-5能成功应用于三角翼的跨声速大攻角流动,网格规模的增加进一步提高流场分辨率,SST湍流模型相对SA湍流模型在三角翼大攻角流动中具有更好的适用性。     

风力机标模非定常数值模拟中的影响因素研究

基于自主研发的“亚跨超 CFD 软件平台”(TRIP3.0),采用刚性运动网格技术和动态拼接网格技术,开发了针对旋转类机械的非定常求解模块。本文开展了 NREL Phase VI 风力机标模非定常数值模拟中的影响因素研究,影响因素主要包括不同时间步长的影响、不同湍流模型的影响、刚性动网格技术和动态拼接网格技术的影响三个方面。本文数值模拟采用的控制方程为雷诺平均 N-S 方程,采用有限体积法离散控制方程,离散方程的时间方向采用“双时间步”方法求解,空间方向无粘项离散采用 MUSCL-Roe 格式,湍流模型采用 SA 和 SST 湍流模型,并引入多重网格和并行计算技术加速求解。数值模拟结果表明：时间步长、湍流模型、网格方案等因素主要影响风力机叶片吸力面的流动结构,进而影响吸力面的压力分布,而对压力面的流动结构和压力分布基本没有影响;采用刚性运动网格结合 SA 湍流模型所得到的压力分布更接近实验值。     

基于隐式嵌套重叠网格技术的阻力预测

 采用一种多层多块隐式嵌套重叠网格技术,对美国国家航空航天局通用化研究模型(NASA-CRM)翼身平尾(WBT)组合体进行了数值模拟与分析。多层多块隐式嵌套重叠网格技术是结合多层多块嵌套重叠网格处理策略和隐式切割方法,在建立重叠网格之间的流场信息传递关系时,基于网格单元切割准则选择"最优"重叠单元而无需人工设定插值边界。对美国AIAA委员会召开的第4届阻力预测研讨会(DPW-4)提供的CRM WBT组合体生成4种不同密度的结构化多层多块嵌套重叠网格,并采用计算流体力学(CFD)方法进行数值计算和阻力预测,计算结果与CFL3D和OVERFLOW的结果进行了对比。数值模拟结果表明:计算得到的压力分布和极曲线与CFL3D和OVERFLOW的结果几乎相同,说明了隐式嵌套重叠网格技术的有效性,同时也验证了流场求解方法与程序的可靠性。当迎角增大到3°左右时,在机身与机翼、尾翼连接处出现明显的分离涡,影响CRM WBT组合体的气动特性。在阻力预测方面,增加网格密度能够提高阻力预测的精度。采用不同的湍流模型会导致升、阻力系数的计算结果存在一定的差异,因此,湍流模型的选择也是阻力预测需要考虑的因素。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切.使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格.应用特征恢复技术恢复凹面特征并改善凸面处的网格质量.流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量.加入Spalart-Allmaras一方程湍流模型以获得湍流解.上下对称高斯-赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性.算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算.     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

k-kL两方程湍流模型的改进及验证研究

提高湍流数值模拟的准确性,从而明确湍流模型对数值模拟结果的影响具有重要的意义,应用K.S.Abdol-Hamid给出的尺度自适应k-kL两方程模型封闭RANS方程,并修改von Karman长度尺度的限制方法,通过平板、翼型、后台阶等流动的模拟,考察k-kL模型在湍流模拟中的准确性,及其反映主要流动特征的能力和网格收敛性,并对影响流动模拟准确性的因素进行讨论.结果表明:改进长度尺度限制之后的k-kL两方程模型无论是对附着流动还是分离流动都可以给出比较准确的结果.     

多段翼型大迎角下主翼、襟翼上的分离流及缝道流动

使用雷诺平均N-S方程、采用可用于较大分离区的Johnson-King紊流模型、嵌套网格和有限体积法研究大迎角下的多段翼型绕流,特别是主翼、襟翼上的分离流动及缝道流动。利用嵌合体技术对组合体每一部分生成高质量并适于高效求解的贴体网格。以具有17%相对厚度的GAW-1翼型带30%襟翼翼型为例进行了计算,计算结果与实验结果吻合很好,证实该方法可以较好地预计多段翼型上的分离流、缝道流动与最大升力。     

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

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

Robust and efficient computation of turbulent flows around civil transport aircraft at flight Reynolds numbers

This work aims at providing numerical methods that enable the robust and efficient simulation of turbulent flows around civil transport aircraft configurations at flight Reynolds numbers. Robustness problems related to the multigrid treatment of advanced transport equation turbulence models have been identified preventing convergence of simulations at high or flight Reynolds numbers. Therefore the application of multigrid to the turbulence equations is omitted while the multigrid treatment of the mean flow equations stays unchanged. Moreover, a fully implicit time integration scheme – a DDADI approach – is applied to the turbulence equations. The new approach raised robustness to a sufficient level to converge viscous computations at flight Reynolds numbers. Furthermore, the convergence speed of aerodynamic coefficients for three dimensional applications has also been improved significantly.     

非结构网格下含冷却孔的涡轮转子三维流场计算

在转动坐标系里,讨论非结构网格下,含多排冷却孔的涡轮转子三维流场雷诺平均Navier-Stokes方程组的求解。计算表明:含冷却孔的转子流场计算要比静子流场复杂的多。为了得到较为满意的流场数值解,本文对叶片压力面与吸力面采用了不同的湍流模式、并且还对流场迭代中初始场选取方式以及转子叶尖间隙计算站的设置等问题进行了一系列的改进。几个典型流场算例表明:初始场选取、湍流模式的合理选用对转动坐标系下获得一个收敛的三维流场数值解至关重要;另外,转子叶尖间隙站的设置直接影响着间隙附近的流场,尤其对现代高负荷涡轮转子,间隙流动的研究就更为重要。