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

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

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

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

非定常低雷诺数流数值模拟方法研究

低雷诺数翼型流动的一个重要特征是分离泡的产生,并且分离泡通常是非定常的。基于有限体积法,N-S方程和S-A模型的控制方程非耦合求解,时间均采用牛顿型LU-SGS方法推进,无粘项和粘性项用中心差分方法离散计算,在计算非定常流场时,采用伪时间子迭代方法,湍流模型计算时采用在固定点转捩,研究了二维翼型低雷诺数流粘性流动的数值模拟方法,并利用所开发的非定常可压缩粘性流计算方法,数值模拟了在低雷诺数下两种翼型的非定常流动。计算结果表明：气动力系数及分离泡呈规则的周期性变化,证明方法是成功的、可行的。     

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

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

翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。     

振动翼型跨音速绕流的数值模拟

从积分形式的非定常Euler方程出发,在固联于振动冀型的贴体坐标系下,采用有限体积法进行离散,并根据对称型TVD格式构造相应的数值通量,数值模拟了振动翼型的非定常跨音速绕流的变化过程。计算结果表明,该方法能较好地捕捉流场中的运动激波,准确地定出激波位置与强度;并且易于推广到求解任意振动方式的跨音速非定常绕流问题。     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

跨声速无粘流绕振动翼型流动的数值模拟

本文从Euler方程出发,利用LU-ADI格式对绕二元振动翼型跨声速流场进行了数值模拟。文章中给出了绕NACA0012和NACA 64A 010翼型不同状态的几组数值结果,在同某些实验和其它文献结果基本吻合的前提下,研究了振动对翼型表面压力、激波位置的非定常效应以及其它振动参数对流场非定常效应的影响。     

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

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

微型飞行器低雷诺数矩形扑翼非定常气动特性的数值模拟

采用双时间步方法求解三维可压缩非定常N-S方程,数值模拟了微型飞行器低雷诺数矩形扑翼的非定常绕流,首先将得到的结果与文献进行了对比,数据间具有较好的一致性.然后针对不同的展弦比、减缩频率及初始攻角,计算了矩形扑翼的非定常气动特性及表面流态和动态压力分布,并分析了翼尖涡对扑翼非定常气动特性的影响.     

二维机翼非定常气动力学反命题的变分理论

非定常流动的反命题迄今在世界文献中尚属空白。本文首次对它进行探讨，就二维振荡机翼的非定常气动力学反命的适当提法及其数学模型作了较深的研究，利用变域变分工具成功地建立了相应的变分原理及广义变分原理，为非定常机翼设计及有限元解的应用提供完整周密的理论基础。     

Application of Active Flow Control Technique for Gust Load Alleviation

A new gust load alleviation technique is presented in this paper based on active flow control. Numerical studies are conducted to investigate the beneficial effects on the aerodynamic characteristics of the quasi "Global Hawk" airfoil using arrays of jets during the gust process. Based on unsteady Navier-Stokes equations, the grid-velocity method is introduced to simulate the gust influence, and dynamic response in vertical gust flow perturbation is investigated for the airfoil as well. An unsteady surface transpiration boundary condition is enforced over a user specified portion of the airfoil’s surface to emulate the time dependent velocity boundary conditions. Firstly, after applying this method to simulate typical NACA0006 airfoil gust response to a step change in the angle of attack, it shows that the indicial responses of the airfoil make good agreement with the exact theoretical values and the calculated values in references. Furthermore, gust response characteristic for the quasi "Global Hawk" airfoil is analyzed. Five kinds of flow control techniques are introduced as steady blowing, steady suction, unsteady blowing, unsteady suction and synthetic jets. The physical analysis of the influence on the effects of gust load alleviation is proposed to provide some guidelines for practice. Numerical results have indicated that active flow control technique,as a new technology of gust load alleviation, can affect and suppress the fluid disturbances caused by gust so as to achieve the purpose of gust load alleviation.     

低雷诺数下微型飞行器主动变形机翼非定常气动特性数值分析

以微型飞行器的气动外形设计为研究背景,通过数值求解N-S方程分析了主动变形柔性机翼的气动特性。空间离散格式采用中心有限体积法,时间推进采用双时间推进方法,其中子迭代过程由隐式LU-SSOR方法完成。在动网格技术的基础上,模拟分析了薄翼面上作行波运动的流场。计算结果发现,该行波翼型的气动特性相比于刚性翼型有明显改善,有效地抑制了大攻角下大尺度流动分离,同时升阻比提高了约35%,起到了显著的增升减阻的作用。本文的工作对于进一步的柔性机翼优化具有良好的理论研究与数值模拟的参考价值。     

基于拉格朗日拟序结构的局部激励流动分离控制有效性分析方法

从流体输运角度提出了用于局部周期激励流动分离控制有效性研究的数值分析方法。基于有限时间不变流形理论建立用于非定常流动的流体输运分析方法,并采用数值方法从非定常流场中提取得到吸引拉格朗日拟序结构（LCSs）和排斥LCSs描述流体输运行为。通过对局部周期激励的流动分离控制规律的研究,结果表明存在三种与激励频率相关的影响翼型气动特性的流体输运模式,其中在锁频激励下吸引LCSs在前缘形成的尖楔结构有效增强主流与分离区流体的物质交换,减小翼面流动分离区的面积,显著提高翼型升力。     

时间谱方法中的高效GMRES算法

研究了时间谱方法求解周期性非定常流场的计算效率,并对时间谱方法应用于周期性非定常流动的隐式求解方法进行探讨。当采样点数增加或减缩频率增大时,时间谱方法对应的雅可比矩阵对角占优性质迅速恶化,导致很多传统的迭代方法失效。为了解决上述问题,论文采用带预处理的广义极小残差(GMRES)算法来提高雅可比系数矩阵的计算收敛性。使用时间谱方法对NACA0012翼型强迫振荡算例进行计算,并与时域差分方法的计算效率和精度进行对比。研究表明在保证计算精度的同时,时间谱方法普遍可将计算效率提高一个量级左右。对于跨声速周期性流动,广义极小残差算法不论是稳定性还是收敛性都优于对称SGS迭代算法。     

不同雷诺数下翼型气动特性及层流分离现象演化

低雷诺数下空气黏性效应突出,翼型表面普遍存在层流分离现象,相比常规雷诺数情况气动特性显著恶化。采用带预处理的Roe方法求解非定常可压缩Navier-Stokes方程的数值模拟技术和低雷诺数低湍流度风洞油流显示试验技术,对FX63-137翼型不同雷诺数下气动特性和流动结构展开深入研究。通过风洞油流显示试验可以清晰获得低雷诺数层流分离流动的两道油流汇集线。数值模拟结果表明其分别为时均化主分离线和二次分离线,两种结果定性定量均吻合较好,证明了本文的研究方法有效可靠;雷诺数从500 000降至20 000,翼型气动特性和层流分离流动结构均发生显著的变化,伴随阻力系数剧增和升力系数剧降,时均化流动结构从附体至出现经典的长层流分离泡,并最终演化为后缘层流分离泡,相应的两种分离泡的非定常流动结构也存在显著差异;对于阻力系数和升力系数而言,存在不同的临界雷诺数,因为导致阻力系数剧增的机理在于经典长层流分离泡的产生使翼型压差阻力大增,而造成升力系数剧降的主要原因在于后缘层流分离泡使得等效翼型后部弯度减小;非定常结果显示正是由于翼型表面漩涡周期性的生成与脱落,才造成了低雷诺数下升力系数的周期性波动。翼型上表面主分离涡即将脱落时,流线在后缘附近再附,升力系数达到峰值;而当流体从下表面向上卷起二次分离涡时,尾部流线大尺度分离,升力系数降至谷值。     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

基于适当控制翼型演化历程的气动设计方法

本文提出一种翼型设计的新方法,它基于适当控制翼型演化历程,采用求解非定常流动控制方程的方法,实现了翼型反问题的高效求解.同时在笛卡尔坐标系中应用Level-Set函数,避免了贴体动网格的采用,从而大大减少了所需计算工作量.     

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

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

Aerodynamic periodicity of transient aerodynamic forces of flexible plunging airfoils

For flapping micro-air vehicles, periodical aerodynamics in every plunging period is always desired in the design of aircraft controlling system. During the periodic plunges of rigid airfoil, either periodical or non-periodical aerodynamic forces can be generated. Real airfoils are usually flexible so that the effect of flexibility on unsteady characteristics of aerodynamic periodicity should be considered. In this study, a fluid structure interaction analysis is employed to systematically investigate the aerodynamic periodicity of a periodically plunging flexible airfoil. The influences of several dimensionless kinematic and structural parameters on aerodynamic periodicity are investigated. The results show that the aerodynamic periodicity of plunging flexible airfoil qualitatively resembles that of rigid airfoil as the mechanisms contributing to aerodynamic periodicity are kept for airfoil with various flexibilities. However, aerodynamic periodicity at small plunging amplitude and high Reynolds number can be improved as airfoil becomes medium flexible but worsened as airfoil is very flexible. The reason is that the deformation of the airfoil changes the strength and structure of the nearby vortices, thereby leading to weak vortex–vortex and vortex–airfoil interactions. The extent of improvement in aerodynamic periodicity is also found significantly influenced by incoming velocity.