模拟厚度翼型非定常运动的涡方法

本文研究二维厚度翼型的大攻角非定常运动。用沿翼型表面分布的源和面涡来模拟翼型的厚度效应和升力效应,应用离散涡方法模拟前缘分离涡层和尾涡层。计算了大攻角翼型作俯仰运动的情形,与实验结果符合较好。     

翼型—扰流片的分离气动特性计算

本文介绍了用涡面元法模拟带扰流片的翼型低速无粘分离绕流。在翼型和扰流片的面元上分布线性变化涡。在扰流片后的上下分离流线的面元上分布等强度的涡。上分离流线始自扰流片的梢部，下分离流线自翼型的后缘引出。分离所泡由两离散涡结尾。气泡内总压为常值，它与涡强大小一同求解。分离气泡的形状在迭代求解过程中确定。压强分布和升力系数的计算值与现存文献的数值结果和实验数据是一致的。     

绕翼型低速湍流的数值模拟

采用求解低速流动的SIMPLE算法,用标准K-ε两方程模型、非线性K-ε两方程模型和一种改进的K-ε两方程模型来对绕两种翼型（A翼型和GAW-1翼型）在接近失速攻角情况下的低速湍流流动进行了数值模拟。计算结果表明,非线性K-ε两方程模型和改进的K-ε两方程模型较好地模拟出了翼型表面上的分离流动特性,较准确地预测出了分离点位置以及翼型的气动力系数,对翼型与机翼的工程设计具有一定的参考价值。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

涡对与物体的相互作用研究

本文从非定常形式的Ｎ－Ｓ方程出发，采用有限差分法数值模拟了入射涡对与圆柱及翼型相互作用的流场结构，旨在分析讨论其复杂的非定常流动现象，如二次涡及三次涡的形成，多涡配对，非定常分离及剪切层的演化等，以及不同的物形及涡对与物体间不同的相对位置对影响。此外，文中还通过分析物面压力分布随时间的变化来探讨非常分离的流动特性。     

圆柱／翼型干涉流场的试验研究

在风洞出口低速流中,以NACA0012尾缘钝化翼型为模型,利用粒子成像测速系统,研究了圆柱/翼型结构干涉流动时翼型前缘、近壁和尾缘区域的流场.试验结果表明,由于上游圆柱引起的卡门涡街和翼型相互干涉,在翼型前缘存在大尺度涡的变形、拉伸和破裂,在翼型表面近壁区域和尾迹流场中仅存在小尺度湍流涡,由此可以推断翼型前缘可能是干涉噪声的主要声源区.     

不同翼型厚度和来流马赫数下的桨涡干扰噪声分析

桨涡干扰噪声是直升机气动噪声主要组成之一，为了正确预测和降低直升机噪声，必须开展气动噪声相关物理参数研究。在对声场进行计算流体力学（CFD）直接数值模拟的基础上，分析了不同厚度和来流马赫数下二维平行桨涡干扰噪声传播特性和声源位置，分析了翼型厚度和来流马赫数对桨涡干扰噪声的影响，并得到了可压缩情况下远场声压预测公式。研究表明，低马赫数下，翼型厚度对噪声指向性影响不大，高马赫数下，翼型厚度对噪声指向性影响程度增大；噪声强弱主要随来流马赫数变化，翼型厚度对其影响较小；翼型厚度和来流马赫数变化不会改变声源点位置。开展不同翼型厚度和来流马赫数下的桨涡干扰噪声分析可以为进一步了解并控制直升机桨涡干扰噪声提供一定的参考。     

充气式翼型绕流大涡模拟分析

针对多梁式充气机翼阿翼型,采用大涡模拟方法对其表面的复杂流动进行了数值仿真研究,得到了其迎角在0°-20°时的流动情况和气动力特性,并与其他仿真结果和风洞试验结果进行了比较分析。结果表明：采用大涡模拟方法得到的计算结果与风洞试验结果较为吻合,可以更好地模拟充气式翼型的表面流动情况。     

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

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

翼型近尾迹流中大尺度涡结构的理论模型

根据流动稳定性理论相干结构的相关理论 ,提出了一种二维波模型来描述翼型近尾迹流动中的大尺度涡结构 ,通过采用对近尾迹流动区域分区处理的方法计算了近尾迹区域的大尺度涡结构。计算所得的流线图 ,涡量分布以及雷诺应力分布具有与实验结果相似的分布规律。说明该二维模型能够较好的描述近尾迹流动中这种大尺度涡结构的运动学和动力学特性      

考虑地面效应的翼型动态特性数值模拟

为了研究考虑地面干扰的飞行器非定常气动特性,基于滑移网格技术,通过求解二维非定常Euler方程,对不同近地高度下的NACA0012翼型的动态振动特性进行了数值模拟。分析了存在地面干扰时的升力系数以及俯仰力矩系数迟滞环变化,并结合强迫振动法,计算了翼型各个近地高度下的纵向组合动导数。结果表明,地面效应不仅对定常流场产生影响,更显著地影响了非定常气动力及力矩,近地高度越小,这种效应越明显,使得升力系数的迟滞环面积变小,幅值减小,而俯仰力矩系数的迟滞环变化不规律,力矩系数的动导数随着近地高度减小而增大(绝对值减小),纵向的动态稳定性产生损失。地面效应干扰在飞行器非定常气动研究中应该引起重视。     

多段翼型绕流的数值模拟研究

采用C—H型混合的结构化网格,选用Wilcox提出的k-ω模型求解可压的N-S方程,模拟了多段翼型的流场。结果表明：在4°和6°迎角下的非定常计算结果与试验结果的对比,计算结果与试验值吻合的较好。本文应用了亚迭代时间法,正确地描述了30P30N多段翼型前缘结构区域的流动特征和流动机理。     

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

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

地面对带后缘分离的翼型绕流影响的计算

本文提出了用在地面上分布奇点来模拟地面的方法。计算结果与精确解和镜像法结果很吻合。     

翼伞弧面下反角、翼型和前缘切口对翼伞气动性能的影响

为了研究翼伞弧面下反角、翼型和前缘切口对翼伞气动性能的影响,对带气室的、展弦比为3的不同特征几何参数翼伞模型的流场进行了三维、定常数值模拟。运用有限体积法对三维坐标系下不可压雷诺时均Navier-Stokes (RANS)方程进行了直接求解,采用剪切应力输运(SST)k-ω两方程湍流模型对湍流进行模拟。数值模拟得出的原始翼伞的气动性能参数与试验数据在总趋势上符合很好,不同几何参数翼伞模型计算结果表明:翼伞弧面下反角越大,升力及诱导阻力越小,升阻比变化不大;前缘半径、厚度小的翼伞翼型,阻力更小,升阻比大;前缘切口对翼伞影响区域限于前缘附近,压力分布同干净翼类似,降低了其失速迎角,对升力影响不大,但明显增大阻力。该数值方法可为进一步研究更多不同几何参数的翼伞模型提供参考。     

纳秒等离子体激励控制翼型流动分离机理研究

为研究纳秒介质阻挡放电(NSDBD)等离子体控制翼型流动分离的物理机理,采用已建立的NSDBD唯象学模型耦合非定常Navier-Stokes方程模拟纳秒等离子体对流场的作用。使用非定常雷诺平均NavierStokes方程(URANS)和大涡模拟(LES)两种求解方法,研究纳秒等离子体激励对NACA0015翼型流动分离控制。结果表明:NSDBD等离子体激励促使边界层提前转捩,转捩对控制流动分离起重要作用;NSDBD激励开始时在翼型前缘形成展向涡,展向涡促使分离剪切层失稳并最终进入尾迹,展向涡贴近壁面运动,将外区的高能气流带入近壁区,使上翼面流场结构发生变化,然后翼型前缘流动提前转捩促使流动经过一个小层流分离泡后发生湍流再附,最终在上翼面形成稳定的附着流动。     

考虑间隙内流的二元机翼跨声速气动力分析

基于Navior-Stokes (N-S)方程和Spalart-Allmaras湍流模型,建立了含间隙二元机翼的跨声速绕流场计算模型和间隙内流计算模型.通过FLUENT数值仿真,分析了二元机翼的激波压力分布和马赫数分布,并与风洞试验 结果进行了对比验证.数值结果表明,连通波后高压区和波前低压区的间隙内流会提高机翼绕流场中激波前后的速度梯度和压力梯度.局部间隙内流会对机翼结构带来“向外吸”和“向前顶”的气动载荷作用.     

低速两段翼型的数值模拟

在多段翼型风洞实验中,很难观察翼型不同迎角下主翼、襟翼上的绕流及缝道流动。使用了O-H型混合的结构化网格,采用k-ε二方程湍流模型求解可压的N—S方程来模拟二元风洞中两段翼型的流场。计算结果与实验进行了比较,结果表明：计算与实验结果吻合良好,说明本方法可以较好地模拟两段翼型的绕流,并能很好地显示流场的变化情况。另外,从结果中可以看出：随着迎角的变化,主翼尾流和缝道间的流动对襟翼附面层分离有很大的影响作用。     

基于Lattice Bohzmann方法的翼型绕流数值模拟

将通过物理平面的结构化C型贴体网格变换到计算平面的均匀正交网格,引入非均匀网格插值方法,采用基于Lattice-Boltzmann方程的D2Q9模型,对低雷诺数下的NACA0012翼型绕流进行了数值模拟,对五个不同位置的边界层速度型与基于有限体积法的CFL3D程序的结果进行了比较.结果表明,两种方法的计算结果一致.改进后的Lattice-Boltzmann方法适用于曲边边界,非均匀网格,同时对计算速度影响不大.Lattice-Boltzmann方法计算较简单,因而用LB方法来模拟翼型绕流是一种新的选择,并可应用于更复杂的低雷诺数流动中.     

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。