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

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

基于卷积神经网络和多项式混沌方法的翼型鲁棒性优化

在常规的翼型优化设计方法中,设计点处最优翼型的气动性能会在非设计点处有所恶化,因此有必要对翼型鲁棒性优化方法进行研究。提出一种基于卷积神经网络和多项式混沌方法的翼型鲁棒性设计方法,首先搭建基于卷积神经网络的气动力预测模型;其次采用多项式混沌方法对马赫数和攻角进行不确定度量化,构建翼型鲁棒性气动优化设计系统;最后对以 RAE2822 翼型为基准翼型的气动优化设计问题进行优化设计验证。结果表明：本文提出的翼型鲁棒性设计方法可行,优化后翼型的气动性能和鲁棒性气动优化设计效率在较宽的设计范围内都有所提升。     

激波控制鼓包提高翼型跨声速抖振边界

翼型抖振边界是仅次于升阻比的一项重要气动指标.采用定常雷诺平均Navier-Stokes方程,以升力线斜率平缓及激波位置振荡作为基本判据确定了RAE2822翼型在指定跨声速来流条件下的抖振边界.通过大量计算流体力学(CFD)验证,针对RAE2822翼型设计了一种特定外形的激波控制鼓包并确定了其具体安装位置.该激波控制鼓...     

超临界翼型低雷诺数流动分析及优化设计

以高空长航时无人机(UAV)翼型研究为背景,对超临界RAE2822翼型在高空高亚声速下的低雷诺数气动特性进行了数值模拟及优化设计研究。采用求解雷诺平均N-S方程的有限体积法,对典型低雷诺数下RAE2822翼型绕流进行数值模拟,验证了SST k-ω湍流模型的可靠性和准确性;基于不同高度不同雷诺数下RAE2822翼型的计算气动力对比分析,研究了高度增大所带来的低雷诺数效应;通过对低雷诺数下超临界翼型表面流场结构及流动机理的详细分析,提出了一种弱化激波的翼型设计思想,并通过优化算例验证了该思想的可行性。     

翼型不确定性量化中正交匹配追踪的应用

不确定性在实际系统中广泛存在，为了研究不确定性因素影响下系统输出的随机响应特性，传统的不确定性量化方法如蒙特卡洛采样、混沌多项式展开等需要大量的样本，制约了在飞机机翼等复杂系统中的应用。近年来，在信号处理领域发展迅速的压缩感知技术，利用原始信号的稀疏性可以用少量的样本精确重构信号。这一特性促使研究人员探索将压缩感知技术应用于不确定性量化研究中。以RAE2822实际翼型为研究对象，使用类函数/形函数变换将原始翼型参数化，考虑加工、装配过程和实际飞行工况下的几何不确定性，将压缩感知技术与混沌多项式展开相结合，利用正交匹配追踪算法实现多项式系数的稀疏重构，获得翼型气动力系数和流场参数在考虑几何不确定性影响下的均值和标准差，并与蒙特卡洛采样和满秩概率配点法获得的结果进行对比。通过对收敛性能、样本数需求和重构精度等方面的对比分析表明，正交匹配追踪算法能够利用相对较少的样本获得与传统不确定性量化方法相近的精度。考虑到实际系统的随机响应在混沌多项式基底上大多具有稀疏的展开形式，因此将压缩感知技术应用到不确定性量化中可以显著降低样本数需求，从而降低时间成本，提高计算效率。     

无网格法耦合RNG k-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

基于移动最小二乘无网格方法,耦合RNG(Re-Normalisation Group)k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM(Advection Upstream Splitting Method)+-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP(Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNGk-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

无网格法耦合RNGk-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

 基于移动最小二乘无网格方法,耦合RNG (Re-Normalisation Group) k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM (Advection Upstream Splitting Method)⁺-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP (Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNG k-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

翼型模型几何误差对气动性能影响的自动微分分析方法

基于自动微分原理和NS方程有限体积方法建立了一套翼型敏感性导数计算方法和程序,可以一次性获得翼型不同气动力系数、压力分布对模型几何外形误差的敏感性导数和不确定度。计算结果表明,在跨声速范围内,即使翼型的外形误差只有63μm(弦长1m),也可以给翼型压力分布带来0.312(以来流动压为参考)量级的不确定度,而激波处的流动最为敏感。这种自动微分方法对于分析数值模拟结果的分散性、风洞试验结果的分散性或不确定性具有很好的指导意义。     

无尾布局翼型的DISC设计研究

采用DISC方法耦合N-S方程求解程序研究了无尾布局翼型的设计。通过NLF1015,RAE2822翼型的反设计表明,该方法可用于亚、跨声速翼型设计,具有简单、计算量小的优点。完成了翼型改进和新翼型设计,对某BWB飞机内翼段翼型进行了改进设计,采用人工干预对现有翼型压力分布进行局部修正,达到快速削弱激波、改善翼型气动性能的目的;针对小展弦比翼身融合布局内翼段新翼型设计,采用基于多目标多约束压力分布优化的设计思想,得到满足气动和几何约束的目标压力分布。设计结果表明,数值优化目标压力分布方法较易满足总体与气动的设计目标和约束。     

固定转捩和自由转捩流动预测方法验证

采用固定转捩位置和基于SST k-ω湍流模型的γ-Reθt转捩模型自由转捩模拟两种方法研究了RAE2822翼型流动转捩问题。先采用γ-Reθt转捩模型确定转捩位置,再采用所获得的转捩位置实施固定转捩位置计算,并与全湍流计算结果和其他软件计算结果进行了对比。表明两种方法均能有效捕捉超临界翼型流动转捩现象,固定转捩位置模拟结果基本上重现了自由转捩预测结果。     

基于分布估计算法的翼型稳健设计(英文)

A transonic airfoil designed by means of classical point-optimization may result in its dramatically inferior performance under off-design conditions. To overcome this shortcoming, robust design is proposed to find out the optimal profile of an airfoil to maintain its performance in an uncertain environment. The robust airfoil optimization is aimed to minimize mean values and variances of drag coefficients while satisfying the lift and thickness constraints over a range of Mach numbers. A multi-objective estimation of distribution algorithm is applied to the robust airfoil optimization on the base of the RAE2822 benchmark airfoil. The shape of the airfoil is obtained through superposing ten Hick-Henne shape functions upon the benchmark airfoil. A set of design points is selected according to a uniform design table for aerodynamic evaluation. A Kriging model of drag coefficient is constructed with those points to reduce computing costs. Over the Mach range from 0.7 to 0.8, the airfoil generated by the robust optimization has a configuration characterized by supercritical airfoil with low drag coefficients. The small fluctuation in its drag coefficients means that the performance of the robust airfoil is insensitive to variation of Mach number.     

基于特征正交分解的跨声速流场重构和翼型反设计方法研究

在二维流场的重构问题中应用特征正交分解(Proper Orthogonal Decomposition,POD)数据处理方法。利用CFD技术计算得到的流场快照对气动力模型进行降阶,然后利用基于POD的降阶模型(Reduced Order Model,ROM)对所需的流场参数进行重构,在快照范围内可以得到高精度的结果,且具有一定的外插能力。在翼型反设计问题中该方法仍然是成功的,通过修正快照向量,利用基于POD降阶模型的数据重构方法,由已知的翼型表面压力分布通过反设计就能够高效精确地得到对应的最优翼型形状,这极大地简化了翼型反设计问题。本文分别在跨声速范围对RAE 2822翼型的流场重构和Korn翼型及NACA 63212翼型的反设计进行了验证,证明了基于POD的流场重构和翼型反设计方法的有效性和高效性。     

Transonic buffet control research with two types of shock control bump based on RAE2822 airfoil

Current research shows that the traditional shock control bump (SCB) can weaken the intensity of shock and better the transonic buffet performance.The author finds that when SCB is placed downstream of the shock,it can decrease the adverse pressure gradient.This may prevent the shock foot separation bubble to merge with the trailing edge separation and finally improve the buffet performance.Based on RAE2822 airfoil,two types of SCB are designed according to the two different mechanisms.By using Reynolds-averaged Navier-Stokes (RANS) and unsteady Reynolds-averaged Navier-Stokes (URANS) methods to analyze the properties of RAE2822 airfoil with and without SCB,the results show that the downstream SCB can better the buffet performance under a wide range of freestream Mach number and the steady aerodynamics characteristic is similar to that of RAE2822 airfoil.The traditional SCB can only weaken the intensity of the shock under the design condition.Under the off-design conditions,the SCB does not do much to or even worsen the buffet performance.Indeed,the use of backward bump can flatten the leeward side of the airfoil,and this is similar to the mechanism that supercritical airfoil can weaken the recompression of shock wave.     

基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。     

一种新的边界处理方法在笛卡尔网格中的应用

采用自适应笛卡尔网格方法求解Euler方程,固壁边界条件通过一种离散强迫浸入边界方法引入,边界切割网格与流场网格保持一致,从而有效解决了小切割网格单元的时间步长限制问题.针对自适应笛卡尔网格的特点提出了新的虚拟点反射方法和插值方法,解决了数值吸力峰值问题.对RAE2822翼型、双NACA0012翼型和Suddhoo三元翼型的流动状况进行了无黏数值模拟,并与现有的理论解、单域以及分区结构网格解进行了对比.结果表明:该离散强迫浸入边界方法有效地引进了固壁边界条件,结合自适应笛卡尔网格技术能够准确模拟复杂几何外形下的流动状况.      

基于天线安装的菱形翼无人机翼型优化设计

以菱形翼布局高空长航时(HALE)传感器无人机研究为背景,分析了翼型内部安装平面相控阵雷达天线的基本原理。通过推导雷达探测性能的估算方程,建立了在翼型中安装平面相控阵雷达天线的数学模型。采用数值模拟的方法,对典型雷诺数下的RAE2822翼型的气动特性进行了研究并与试验数据进行了对比,验证了使用有限体积法k-ωSST(Shear Stress Transport)湍流模型求解雷诺平均Navier-Stokes方程在这一状态下的可靠性和适用性;同时对不同厚度的NACA 64A系列翼型的流场结构和流动机理进行了分析,在此基础上提出了一种满足平面相控阵雷达天线安装情况下翼型的优化设计思想。优化结果验证了这一设计思想的可行性。