翼型分离流动的等离子体控制数值分析方法研究

利用等离子体进行流动控制是当前的研究热点之一.本文研究通过数值方法模拟等离子体对流场作用的实现方法,及利用数值模拟方法研究翼型大迎角分离流动的等离子体控制.利用CFD软件Fluent中的自定义函数接口,通过C语言编程在软件中引入DBD等离子体激励模型外加体积力源项,对NACA0015翼型大迎角下的等离子体控制进行数值模拟.验证DBD等离子体激励在抑制流动分离与增升减阻方面的作用.结果表明：流动控制效果与DBD激励器布置位置有直接关系;激励器的数量与激励强度均会影响流动控制作用.     

翼型动态失速DBD等离子体流动控制的数值模拟研究

进行了翼型深度动态失速及基于DBD介质阻挡放电等离子体激励器的流动控制技术数值模拟研究。将激励器对流动空气的作用以彻体力源项形式加入N-S方程。通过数值求解此N-S方程,研究了DBD激励器对NACA0012翼型俯仰运动深度动态失速的控制作用;研究了DBD激励器工作方式对动态失速平均气动力、气动力迟滞曲线的影响,提出了控制效果较好的激励器工作方式。     

等离子体控制边界层法流动仿真研究

使用泊松方程、漂移-扩散模型和N-S方程对介质阻隔面放电等离子体控制边界层流动进行了一体化数值模拟.仿真结果表明介质阻隔面放电等离子体通过推-拉空气可以对低速气流进行控制,正向推力和反向拉力均可以在距壁面一定高度处形成射流,正向射流自身增加边界层动量,更重要的是通过引射主流将主流中的能量输运到边界层中,从而减小边界层厚度;反向射流与来流相撞后形成回流区将边界层向外推,能够增加边界层厚度.     

等离子体流动控制研究进展与展望

等离子体流动控制是基于等离子体气动激励的新型主动流动控制技术,具有响应时间短、激励频带宽等显著技术优势,在改善飞行器/发动机空气动力特性方面具有广阔的应用前景,已成为国际上等离子体动力学与空气动力学交叉领域的前沿研究热点。鉴于此,从介质阻挡放电(DBD)、电弧放电等离子体气动激励特性,等离子体气动激励抑制流动分离、控制附面层、控制激波与激波/附面层干扰、控制压气机与涡轮内部流动、控制管道流动和飞行控制等方面,综合评述了国际上等离子体流动控制的研究进展情况;从创新等离子体气动激励方式,揭示等离子体气动激励与复杂流动的非定常耦合机制,突破等离子体流动控制系统关键技术等方面,对未来的发展进行展望。     

常气压辉光放电等离子体控制翼型失速的数值模拟研究

基于Shyy提出的大气压下均匀辉光放电等离子体与空气干扰的物理模型,通过求解电位势方程得到电场分布,得到了作用于流体上的电场力.通过数值求解考虑等离子体作用的流体运动控制方程,以NACA0015翼型低速绕流为对象,研究了常压辉光放电等离子体位置和个数控制翼型绕流分离的影响.位于分离点上游的等离子体能够有效地抑止翼型分离,而在分离区的等离子体对流动影响很弱,同实验结果一致,并给出了等离子体对翼型气动力影响的规律.     

飞翼模型微秒脉冲等离子体控制低速风洞试验研究

为改善飞翼模型低速、大迎角气动特性,在试验段截面为4.5m×3.5m的低速生产型风洞中开展了大展弦比飞翼模型微秒脉冲等离子体流动控制的试验研究,所用的飞翼模型展长为2.4m,展弦比为5.79,试验研究采用了测力和PIV (Particle Image Velocimetry)两种试验方法。通过测力试验研究了等离子体激励位置和激励频率对飞翼模型失速特性的影响,通过PIV流动显示试验给出了等离子体对翼面流场结构的影响。试验研究表明:等离子体控制能显著改善大展弦比飞翼模型低速大迎角下的气动特性,激励位置和激励频率对流动控制效果具有较大影响;等离子体激励位置在机翼前缘驻点附近、激励频率为100Hz时控制效果最好;试验风速V=70m/s (Re=2.61×106),等离子体激励的峰峰值电压为10kV时飞翼模型的最大升力系数提高20.51%,失速迎角推迟6°。     

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

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

等离子体控制边界层流动仿真研究

使用泊松方程、漂移-扩散模型和N-S方程对介质阻隔面放电等离子体控制边界层流动进行了一体化数值模拟。仿真结果表明介质阻隔面放电等离子体通过推-拉空气可以对低速气流进行控制,正向推力和反向拉力均可以在距壁面一定高度处形成射流,正向射流自身增加边界层动量,更重要的是通过引射主流将主流中的能量输运到边界层中,从而减小边界层厚度;反向射流与来流相撞后形成回流区将边界层向外推,能够增加边界层厚度。     

新型等离子体激励器的数值模拟

为了对新型多电极对等离子体激励器进行可靠的数值模拟,通过求解带有体积力源项的Navier—Stokes方程,模拟单电极对等离子体激励器在二维不可压流中对周围流场的影响,计算结果表明所使用唯象模型的参数基本正确。在此基础上,改进原有模型,对新型双极性等离子体激励器诱导流动进行计算。通过与实验结果的对比,发现改进后的模型能较好地模拟新型双极性等离子体激励器的控制作用,可为今后新型等离子体激励器的应用提供可靠的数值模拟方法。     

空气动力学发展中值得关注的一些问题

本文阐述了空气动力学研究的特点；分析了21世纪初空气动力学面临的形势，提出了空气动力学发展中值得关注的一些问题，包括新流型的研究与先进气动布局、流动控制与增升减阻、流动的逢适应控制与智能化空气动力学、低雷诺数流动与微流体力学、等离子体减阻技术、等离子体隐身与电空气动力学、机载激光武器的发展与气动光学等。     

用粘性体积力方法计算高速离心压气机叶轮内部流场

高速离心压气机是用途非常广泛的一种动力机械,充分了解其内部复杂的三维流场对设计高压比、高效率离心压气机具有很重大的意义。本文采用一种简化的Navier-Stokes方程的计算方法-粘性体积力方法,计算了Krain叶轮内部流场。使用H型网格离散微分方程,并采用多重网格方法和当地时间步长方法加速计算收敛。粘性应力通过Baldwin-Lomax模型计算。为了证明计算结果的有效性,文中也给出了相应的激光测量结果。      

基于新型丝状电极等离子体激励器的前体涡控制

本文首次将新型丝状暴露电极DBD等离子激励器应用于大迎角下细长体非对称涡控制。丝状暴露电极的材料的选择对DBD推力以及推力效率至关重要,通过地面精细推力测量对丝状暴露电极等离子体激励器进行了优化,结果表明,本文研究材料中采用钨丝作为暴露电极,其推力效率最优;且随着电极直径从d=0.3 mm减小到d=0.08 mm,DBD推力效率显著提升。基于优化后的DBD激励器,将其应用于前体非对称涡控制:未施加等离子体控制时,压力测量以及PIV结果均表明细长体背风区流场为明显的非对称涡结构;在等离子体激励下,该非对称涡结构可变为对称甚至反向非对称,且非稳态激励控制能力明显优于稳态激励。研究发现,大迎角下细长体非对称涡控制与背风区原始涡系结构有关,其中包含对称涡系和非对称涡系。本文研究为大迎角下细长体非对称涡控制提供了一种新思路,同时也为丝状暴露电极DBD等离子体激励器的应用提供参考。     

基于雷达散射截面约束的气动优化设计

采用拉格朗日乘子法优化设计了雷达散射截面约束条件下的锥形融合气动外形.拉格朗日乘子法中的极小化问题采用动态演化的优化设计方法求解.该方法是一种基于非定常演化的优化设计方法,即在求解非定常流动支配方程的时候同时履行优化过程,较其它基于定常解的优化方法具有高得多的计算效率.其中的雷达散射截面通过求解非结构的笛卡儿网格上的时域麦克斯韦方程来得到,而升阻比则通过求解锥形流方程来计算.通过优化设计,得到了M∞=8.0时,升阻比为4.98,雷达散射截面只有1.66m2的锥形融合气动外形.     

喷管化学反应流的对角化点隐式法求解

采用弱耦合对角化点隐式方法的ＭａｃＣｏｒｍａｃｋ格式,求解喷管的粘性化学反应流动,解决了控制方程组的刚性问题,又大大减少了常规点隐式方法中求逆矩阵的大量工作。对液氢液氧火箭发动机,采用６种组分、８个反应有限速率的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流场参数在喷管中的分布。计算结果表明,采用对角化点隐式方法求解化学反应喷管粘性流动能提高计算效率。     

叶轮机械三维粘性动静叶干涉的数值模拟

本文发展了 1套模拟叶轮机械三维粘性动静叶干涉的数值方法。在成熟的空间离散方法基础上 ,提出了 1种高效的双时间步隐式时间推进方法 ,并构造了一种保证守恒的动静叶交界面处的数值通量计算方法。本文对某三维亚音速涡轮的动静叶干涉进行了数值模拟 ,结果与实验符合较好。     

求解叶轮机粘性流场的再生解析数值方法

应用再生核方法求解叶轮机粘性流场是一种新颖的解析数值方法, 它可以解决计算中的存贮、计算速度和精度问题。本文给出了用再生核方法求解叶轮机粘性流场的思路和方法, 以及如何在曲线坐标系中应用再生核方法, 并给出了流函数积分方程的求解过程。最后对平面扩压管道进行了数值试验, 结果与有关文献是吻合的。      

Numerical research on airfoil transition delay by alternative current dielectric barrier discharge actuation

A Dielectric Barrier Discharge (DBD) plasma actuator can create a body force which locally accelerates the base flow leading to an attenuation of broadband disturbance to delay the transition. In this study, numerical simulation on an NLF0416 airfoil is conducted to investigate transition delay and drag reduction by a DBD plasma actuator. To simulate plasma’s effect more accurately, boundary-layer data is acquired from Reynolds Averaged Navier Stocks (RANS) equations instead of laminar boundary layer equations, although RANS equations need a much finer boundary-layer grid, and the linear stability analysis method is used to analyze the boundary layer and get the transition point. In this study, the influences of different actuation intensities and positions are investigated, and results show that if the actuation intensity is stronger and the actuation position is closer to the base transition point, more drag reduction can be obtained. However, the efficiency of plasma transition delay is really low. For example, when the actuation voltage is 16 kV, the actuation frequency is 1 kHz, and the main Mach number is 0.1, the saved power due to drag reduction is about 5.09 W, but the power consumed is about 32.61 W, and the efficiency is just 15.6%.     

基于动态非结构重叠网格法的直升机前飞非定常流场数值模拟研究

 发展了非结构重叠网格技术，提出适合于直升机旋翼运动的动态非结构重叠网格算法，该算法使得在采用非结构重叠网格模拟直升机流场时实现程序自动处理网格，不需任何人工干预。流动控制方程采用非定常可压缩Euler方程，空间离散采用Jameson有限体积法，时间推进采用双时间步长法，为加速收敛，在双时间步长内循环迭代中采用预处理技术。为验证本发展的预处理解算器以及动态非结构重叠网格的算法的可信性和有效性，开展一系列的数值实验，模拟了不同速度范围的NACA0012翼型绕流,Caradonna旋翼悬停和前飞流场以及Robin直升机机身-旋翼组合的前飞非定常流场，结果表明采用本文算法得到的结果是可信的，网格处理是自动高效的，可以应用于直升机全机流场模拟，为直升机设计提供可信指导数据。     

跨音速压气机级三维紊流流场数值模拟

在对跨音速压气机级的动／静双叶排中的三维定常紊流流场数值的模拟中，利用三阶高分辨率ＮＮＤ格式和ＬＵ－ＳＧＳ隐式推进迭代法，既保证了流场中激波的模拟质量和粘性流动特征的正确预估，也实现了求解过程的高效率。对某单级跨音速压气机的计算及其与实验数据的对比证实了本文方法的有效性。将计算推广至多级轴流叶轮机械的情况将会是十分简便的