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

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

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

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

绕大攻角翼型低速非定常流动的SUPG有限元数值模拟

采用既有良好稳定性又有较高精度的ＳＵＰＧ有限元方法,通过数值求解涡量一流函数形式的不可压缩流ＮＳ方程组,模拟大攻角翼型在高Ｒｅ数下几种典型情况的非定常绕流流场,同时验证了ＳＵＰＧ有限元方法在求解对流占优的二维非定常对流的－扩散方程时的有铲性。     

粘性流绕俯仰振动翼型流动的数值模拟

本文直接从N-S方程出发,利用LU-ADI格式和Baldwin-Lomax代数湍流模型研究了粘性流绕俯仰振动翼型的流动。数值实验表明,本文的数值结果同实验吻合较好。通过数值模拟手段研究了振动对流场中激波和分离这两大主要特征的影响,结果发现;(1)激波滞后于攻角的变化,如当NACA0012翼型在负攻角状态下上翼面存在激波而下翼面无激波。(2)由于翼型的振动,分离被减轻升力相对提高。(3)随翼型攻角和振幅的变化,翼型振动的升力回线走向可不同。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

细长体大攻角绕流的数值模拟

用数值模拟方法研究了细长体大攻角的绕流特性（求解在周向３６０°范围内进行）。数值模拟的出发方程是三维可压流的全Ｎ－Ｓ方程,格式为二阶迎风ＴＶＤ格式。通过数值模拟研究了附近,攻角１９°、３２．５°和５０°状态下,细长体绕流场中的激波、分离流动及表面压力分布等重要特征。     

绕翼型分离流结构的数值研究

利用大涡模拟技术,对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟,详细地给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数。在数值模拟中,采用了弱压缩流的控制方程,用贴体坐标技术进行了网格生成。     

基于LBM-LES方法的多段翼型流场数值模拟

运用Lattice Boltzmann(LB)方法与Large Eddy Simulation(LES)结合的方法数值模拟研究了二维翼型绕流流场,其中,粒子速度模型采用D2Q9模型,松弛时间中对应的湍流粘性系数由Smagorinsky模型计算,网格生成方法采用无网格技术。对NACA0012翼型进行数值模拟,验证其有效性;为进一步研究其模拟复杂物体绕流的能力,以三段翼型30P30N为研究对象进行数值研究,计算出了给定攻角条件下的气动特性,并与实验值进行比较。计算结果表明,LBM-LES方法的计算结果能够与实验值一致,是一种可行、有效且较新的数值模拟方法。     

三角翼大攻角绕流场数值模拟研究

本文基于N-S方程采用Jameson的显式中心有限体积格式和Roc的上风格式对三角翼大攻角绕流场进行初步数值模拟，空间网格为H-0型网格。计算结果表明两种计算格式均可较好的模拟大攻角绕流场。本文计算结果为今后开展全机大攻角绕流场数值模拟提供了计算经验。     

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

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

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

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

大厚度翼型分离流场的稳定性探讨

基于k-ω的SST两方程湍流模型,求解了雷诺平均Navier-Stokes方程获得定常和非定常气动力,耦合翼型弹性振动方程,在时间域内模拟了弹性支撑下的大厚度对称翼型流场及其翼型动态过程,研究结果发现：在所给定的雷诺数条件下,弹性支撑下的大厚度翼型在一定大迎角范围内存在分离流场稳态的转化现象,此结果的意义在于：1）同一基本迎角下的翼型分离涡存在多个稳态流场;2）弹性支撑的翼型在大迎角分离流条件下气动力在物理上本身就存在非唯一性,从一个新的角度解释了大迎角气动力测量中的分散性。     

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

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

Computations of unsteady transonic flow over airfoil at low Reynolds number

Transonic flow over a thin airfoil at low Reynolds number was studied numerically by directly solving two-dimensional full Navier-Stokes equations through 5th order weighted essentially non-oscillatory(WENO) scheme without using any turbulence model.A series of distinguished unsteady phenomena for a thin 2-D transonic airfoil flow were presented.Due to continuous adverse pressure gradient in the subsonic flow downstream of the sonic line, the unsteady separated boundary layer with main vortex and secondary vortex was developed at the rear of the airfoil.At the trailing edge,the vortex-shedding was characterized by periodical connection of the main vortex and secondary vortex on the other side of the airfoil.The unsteady separation and vortex-shedding occurred with the same period.On the airfoil surface,the average pulse pressure related to the unsteady supersonic region was obviously smaller than that related to the vortex-shedding at the trailing edge.With the attack angle increasing from 0° to 2°, the frequency of vortex-shedding decreases about 4.2%.At last, the turbulence intensity and many second-order statistics in the wake region were investigated.      

基于CFD方法对微型扑翼翼型设计的研究

基于CFD方法从扑翼设计的角度出发对扑翼的外形进行深入的研究.使用基于广义无限插值理论的代数网格生成方法开发适合扑翼运动的三维结构N-S动网格,并应用有限体积法结合双时间推进技术求解预处理后的低速三维非定常N-S方程,完成了对低速扑翼绕流的数值模拟.通过与文献算例的比较验证了本文方法的可靠性和正确性.在此基础上,对文献[8]中介绍的一种混合翼型的相关参数变化对扑翼气动特性的影响开展了深入的研究.     

复合振动翼型跨声速非定常流计算

在文献[1]、[2]的基础上,本文把绕翼弦中点作旋转振动的翼型跨声速流计算推广到翼型作水平、垂直方向的平移振动以及平移与旋转振动的叠加。本方法主要特点是网格与翼型一起运动,对无穷远及翼面边界条件作了精确的处理以及计算网格数很少。本文结果对进一步研究旋翼、叶轮转子等的非定常气动力将有重要的意义。     

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。     

旋转坐标系N-S 方程无网格/网格混合算法

为了模拟涡桨发动机等绕固定轴旋转的桨叶流场,发展了一种求解旋转体非定常黏性绕流的无网格/网格混合算法。算法基于求解旋转坐标系纳维尔-斯托克斯方程展开,避免了旋转角速度特征物理量的插值运算;计算域采用整体网格和物面附近局部无网格离散,通过引入无网格和网格对偶点,实现了混合算法绕流信息的跨区交换,并给出了一种对偶点调整选取的方法;基于无网格点云空间导数逼近方法离散控制方程,并结合隐式LU-SGS算法,给出了混合算法求解非定常问题的双时间步隐式推进格式。结果表明:所提算法通过2维振荡翼型非定常绕流、旋转圆柱黏性绕流和模拟发动机桨叶旋转运动的3维悬停旋翼绕流进行了考核,所得升力系数等重要气动数据的变化趋势与文献试验或计算值一致,典型截面处所示的激波强度和位置亦与文献值吻合,展现出算法在模拟旋转部件绕流问题方面具有广阔的工程应用前景。     

弹性振动对翼型失速迎角附近流场的影响

通过求解雷诺平均Navier-Stokes方程得到气动力,结合翼型振动方程,计算了翼型不同迎角下的动态过程,分别通过层流和湍流情况的计算,重点研究了弹性振动对翼型失速迎角附近流场的影响,研究结果表明,在中低雷诺数、翼型具有弹性振动的情况下,翼型的失速迎角会比传统定常意义上的失速迎角提前出现,为长期以来数值计算得到的失速迎角与风洞实验、飞行试验结果的不同给出了一种物理解释。