翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

扰流片快速张开诱发的非定常升力过冲

通过数值求解二维非定常可压N-S方程,研究了二维翼型上扰流片快速张开诱发的复杂非定常流动问题.利用动态重叠网格方法处理扰流片和主翼相对移动,分析了扰流片快速张开诱发的升力过冲现象及其形成机理.着重探讨了来流攻角和扰流片高度对过冲量的影响.     

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

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

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

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

扰流片快速张开诱发的非定常流动研究

通过数值求解二维非定常可压N-S方程,研究了二维翼型上扰流片快速张开诱发的复杂非定常流动问题。利用动态重叠网格方法处理扰流片和主翼相对移动,分析了扰流片快速张开诱发的升力过冲现象及其形成机理。着重探讨了来流迎角、扰流片高度以及快速张开速率对过冲量的影响。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

旋翼三维非定常跨音速绕流的数值模拟

 利用小扰动位势方程，用有限差分方法求解了旋翼跨音速流场。导出了严格的三维非定常位势方程远场无反射边界条件，并成功地用于数值计算中，缩小了物理求解域的空间尺度。采用了贴体坐标和网格加密技术，以得到合理的非均匀网格分布。应用较少的网格点数成功地计算了无升力情况下旋翼的三维非定常亚、跨音速绕流流场，得到了与实验结果符合良好的数值解     

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

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

双时间法求解大迎角翼型绕流非定常N-S方程

本文采用双时间推进法、顶点格式有限体积法求解非定常ＮＳ方程、应用多层网格技术计算了ＮＡＣＡ００１２翼型大迎角振荡绕流,并对伪时间推进收敛精度,物理时间步长,多层网格迭代方式等对计算结果与计算ＣＰＵ时间的影响进行了讨论,另外还讨论了采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ湍流模型时用固定转捩或全湍流模式对计算结果的影响,并与实验结果进行了比较。     

跨声速粘性流绕振荡翼型的非定常计算

本文采用时间和穴是二阶精度的Ｂｅａｍ－Ｗａｒｍｉｎｇ格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式及ｑ－ω二方程微分模式，结合网格自适应技术，数值模拟Ｎ－Ｓ方程，计算了跨声速下的翼型非定常运动，包括俯仰，浮沉和前后平移振荡。结果表明，压力分布和气动系数与实验基本符合，微分模式和自适应网格能够显著提高激波和边界层计算精度。     

旋翼翼型非定常动态失速特性的CFD模拟及参数分析

构建了一套基于运动嵌套网格技术和可压缩RANS方程的旋翼翼型非定常流动特性模拟的高效、高精度的CFD方法。首先,发展了基于Poisson方程求解的围绕翼型的粘性贴体正交网格生成方法,并提出了基于最小距离法(MDM)改进策略的运动嵌套网格生成方法,克服了弹簧法可能导致网格畸变的不足;其次,为准确模拟由湍流分离和气流再附引起的气动力的迟滞效应,基于RANS方程、双时间方法和高阶插值格式,建立了旋翼翼型非定常气动特性分析的高精度数值方法,并采用能够较好捕捉气流分离现象的S-A湍流模型;再次,针对旋翼后行桨叶动态失速时桨叶剖面来流速度较低、迎角较大的特点,为解决低来流速度时L-B半经验模型在旋翼翼型非定常动态失速计算中的局限性,并克服可压缩方程对低速流场计算收敛困难和精度低的问题,建立了基于Pletcher-Chen低速预处理方法、FAS多重网格法和隐式LU-SGS方法相结合的高效数值方法。应用发展的方法,分别针对NACA0012、SC1095旋翼翼型静态和轻度、深度动态失速进行计算,精确捕捉了气动力迟滞效应以及翼型前缘脱体涡的产生、对流和脱落过程,验证了本文方法的有效性;最后,着重针对NACA0012动态失速状态,开展了振荡参数对旋翼翼型非定常动态失速特性影响的分析,研究结果表明翼型迎角平均值、振幅及减缩频率的变化均能引起迟滞效应的改变并使得气动力峰值发生有规律的前、后移现象等。     

带扰流片的二维翼型非定常计算及分析

通过求解二维非定常可压N-S方程,研究了带扰流片的NACA0012翼型绕流问题.首先计算了扰流片固定在三个不同张开角度时的流动,气动力系数变化规律和已有实验数据对比,结论一致.并着重描述了扰流片张开时,背风区涡形成、发展和脱落过程,分析了非定常气动力产生的机理.为了进一步探讨振动扰流片对主翼的影响,采用Chimera重叠网格技术处理多体相对移动问题,分析了振动频率和振动幅度两个因素对时均气动力的影响.     

振荡扰流片非定常流场的N-S方程数值模拟

 通过有限差分法求解不可压缩流体的二维非定常Navier-Stokes方程,研究扰流片作振荡运动引起的强迫非定常分离流特性。着重研究了雷诺数、缩减频率等因素的影响。得到涡在扰流片上生成、发展、脱落以及脱落后的演化过程;并给出扰流片和平板表面的压力分布。计算结果和实验结果符合良好。     

高超声速验证飞行器助推分离段流场数值研究

为了模拟有相对运动的多体非定常绕流，发展了基于弹簧近似的非结构动网格方法及耦合动网格的Euler解算器，作为方法验证，模拟了俯爷振动的NACA0012翼型绕流，计算结果与实验及文献结果非常接近。采用二维近似外形，对高超声速验证习行器助推分离段流场进行了数值研究，得到了不同时刻由于多体相对运动形成的干扰流场结构以及分离过程的气动力参数。     

用混合网格模拟绕X-38粘性流动

采用混合网格，利用DLR所发展的计算软件TAU对X-38高超声速粘性绕流进行了数值计算，网格自适应技术被用来生成有效的计算网格及提高对流场细节的分辨率。计算结果与现有的实验结果和用结构网格得到的结果吻合很好。     

CFX在军机微小部件优化设计中的应用

针对某型飞机在跨声速飞行中出现垂尾方向舵的铆钉被吸出的问题,本文采用CFX 5.7对全机加细微尺度的扰流条构型进行了一定状态的精细数值模拟。详细分析扰流条改进前后的流场差异,对原来的扰流条改进设计,通过流场的变化判断改进设计是否达到预期目的。计算前处理的网格采用ICEMCFD 5.1网格生成器生成HEX六面体结构网格,网格规模为450万。由于跨声速流场是风洞试验和数值计算模拟的难点,为了准确地反映全机复杂流动的物理现象,本文数值模拟选取M a=0.95,4°迎角,SST湍流模型,二阶差分格式,采用小的物理时间步长在600个时间步以内即可获得跨声速流场的收敛解。通过对流场计算结果的详细分析,证实了方案优化设计的合理性和有效性。本文的计算是CFX 5.7对复杂构型细微尺度变化流场模拟的一个典型算例,充分反映了ICEM 5.1中HEX六面体结构网格对流场细微变化的捕捉能力和CFX 5.7全隐式求解器的良好鲁棒性和快速的收敛性。对于这种特殊部件细微尺度变化的流场模拟,风洞试验显得无能为力,只有通过高精度流体仿真软件来模拟优化设计的效果,为优化设计提供参考的依据。     

基于笛卡尔网格的非定常气动力计算

在非结构自适应笛卡尔网格基础上构建了一种基于笛卡尔动网格的非定常气动力计算方法。为了适应非定常流场对网格品质要求较高的特点,网格生成采用了基于几何外形、模型表面曲率为基础的自适应生成及加密方法。采用显格式龙格库塔推进方法和中心格式空间离散的有限体积法对NACA0012翼型,三维机翼的跨音速刚性振动进行了数值模拟。对伪时间收敛精度控制,真实时间步长对结果的影响进行了讨论。将计算结果与试验数据进行了比较,得到了良好的效果。     

低雷诺数翼型层流分离现象大涡模拟方法

为准确预测翼型低雷诺数条件下出现的层流分离流动现象,发展了基于结构化拼接网格的大涡模拟方法。控制方程为Favre滤波Naver-Stokes方程,并选用多种亚格子模型。空间离散采用AUSM格式以及高阶WENO格式,时间推进采用显式方法和隐式方法。以SD7003翼型在雷诺数60 000及4°迎角下的层流分离流动为研究对象,对比分析了数值格式、亚格子模型、网格尺度对流场预测结果的影响。数值结果表明:划分的计算网格能够有效解析小尺度流动结构,基于隐式亚格子模型、采用AUSM格式和双时间步推进的大涡模拟方法能够准确预测流动分离、转捩、再附等复杂流动现象,计算得到的平均压强系数与雷诺应力分布与文献数据吻合较好。与转捩模式计算结果对比进一步表明:发展的大涡模拟方法能够准确预测翼型低雷诺数层流分离流动非定常演化过程,为下一步的研究工作建立了有效的数值模拟手段。     

可压缩混合层的涡结构演化与流质混合

为研究可压缩混合层的流动结构,采用七阶精度广义紧致格式离散对流项和用显式八阶精度的中心格式离散粘性项,数值求解了非定常三维可压缩Navier-Stokes方程。用约4亿规模的网格,直接数值模拟了对流马赫数为0.7的超声速可压缩混合层的空间发展流动,获得了自初始流动失稳直至充分发展湍流流动结构的精细演化历程,所得结果表明:大尺度涡结构的生成使得混合层的动量厚度快速增长,并主宰了由被动标量质量分数展示的可视混合厚度的量级,充分发展湍流的小尺度的结构主要使得该量级厚度内的流质混合趋于均匀,对可视厚度增长的贡献非常有限。     

基于伴随算子的翼尖小翼优化设计

采用基于伴随算子的优化设计方法,对跨音速翼尖小翼进行了优化设计。其中,采用基于多块结构网格的流场解算器计算气动力。解算器的空间离散采用Osher格式,湍流模型采用k-ω两方程模型。优化设计结果表明,采用优化设计后的翼尖小翼,能改善机翼的气动特性,在巡航状态下翼尖小翼能增加机翼的升阻比5.8%,从数值计算的角度表明了方法的有效性。