超声速涡轮叶型全局气动优化设计

针对涡轮叶型全局优化设计计算时间长、样本空间大等难点提出一种可行的优化设计方法,该方法将控制叶型的17个参数作为优化变量,采用第二代多目标遗传算法进行全局自动寻优。基于此方法,搭建了涡轮叶型全局优化设计平台。利用此平台,分别采用轴向稠度固定和自由优化两种方式对超声速涡轮叶型进行了优化设计。数值计算结果表明,两组优化设计叶型在设计工况下总压损失系数比参考叶型分别低19.5%和10.0%,流道中的激波强度更弱,且在变工况条件下都具有较好的气动性能。深入分析流场与激波结构后发现,外尾激波相比于内尾激波对总损失的影响更大,通过减小气流膨胀转折角或内尾激波气流转折角能够有效削弱外尾激波强度。     

一种高阶精度大涡模拟技术用于激波/火焰相互作用

基于Smagorinsky亚格子模型和增厚火焰技术,开发了一种高阶精度的反应流大涡模拟求解器,运用于数值研究边界层对激波与火焰相互作用的影响。该求解器的核心在于加入了基于超黏性的激波捕捉技术和时空三阶精度的两步Taylor-Galerkin紧致（TTGC）有限元格式,并通过对一维Shu-Osher问题和二维激波/气泡相互作用问题的计算,验证了求解器对激波、接触间断和湍流脉动等流动细节的捕捉精度,计算结果与实验数据吻合良好。通过对激波管内激波、火焰与边界层相互作用问题的数值模拟,发现由于激波与边界层的相互作用会产生不稳定的激波分叉现象,激波三分叉点传播速度的发展经历了水平匀速运动、小斜率线性增长和大斜率迅速增长3个阶段,由此揭示了激波分叉促进火焰加速的机理。当火焰面传入激波分叉区后,流场不均匀的回流区起到了稳定火焰的作用,一方面分叉结构内火焰面能够为激波的运动持续地供应热量,另一方面局部超声速区域为火焰的快速传播提供动力,使其能紧跟激波。通过对比相同条件下甲烷和乙烯燃烧的数值结果,发现两者爆震点触发的位置都出现在马赫杆后面,热量释放率的变化趋势也大致相同,但乙烯出现爆震的时间比甲烷早。     

新型桨尖抑制旋翼跨声速特性的影响分析

为研究不同形式的新型桨尖在抑制旋翼跨声速特性方面的作用,开展了多种桨尖对旋翼局部流动及气动特性影响的数值分析研究.发展了基于高效嵌套网格方法的旋翼流场高精度CFD求解方法.在此基础上,详细分析了桨尖外形对旋翼桨叶跨声速区域激波强度、激波诱导气流分离、桨尖涡尾迹及气动性能的影响.数值结果表明:桨尖的后掠和上反在缓解旋翼跨声速特性方面的作用相对较小;桨尖前掠和下反能更有效地减少桨尖外端跨声速区域,降低该位置激波强度并缓解激波-附面层干扰诱导的气流分离;后掠桨尖在减小旋翼反扭矩方面的整体效果良好,直线前掠桨尖在大桨盘拉力状态能够更有效降低旋翼扭矩(直线前掠30°时,扭矩降低达12.3%),桨尖下反可以有效抑制桨尖涡强度(抛物下反30°时,桨尖涡强度降低50%),并加快桨尖涡尾迹的耗散.      

基于耦合伴随方法的串/并行气动结构优化设计对比

采用RANS方程、线性有限元分析手段、映射点局部插值方法和逆距离动网格技术对气动和结构之间的耦合现象进行了分析。然后基于气动结构耦合伴随方法，对气动和结构设计变量梯度进行了高效求解。最后结合自由曲面变形技术、梯度优化方法构建气动结构优化设计框架。选取飞翼布局民用客机开展气动结构串行与并行优化设计研究。在相同阻力水平的前提下，并行优化结果比串行优化获得了8.4%航程收益和8.3%结构质量收益。同时，串行优化设计结果在中翼段上翼面仍然存在明显激波，而并行优化结果呈现光滑的等间距等压线分布特征。因此并行气动结构优化设计方法可以充分挖掘气动结构耦合下的设计潜力，更有利于对大柔性气动结构优化设计问题进行研究。      

热激励在超声速进气道内对激波诱导的边界层分离的控制机理

通过数值模拟的方法研究了马赫5的超声速进气道内,热激励对激波/边界层相互作用的控制机理。研究了热激励器放热功率(E)、热激励器展向放置数目(N)和热激励器到控制激波的距离(S)三个参数在超声速进气道内激波控制和边界层分离改善中的表现。分别针对以下四种条件进行了数值模拟:1)E=2k W、N=2、S=0.02m;2)E=3k W、N=2、S=0.02m;3)E=2k W、N=3、S=0.02m;4)E=2k W、N=2、S=0m。分析发现:在以上四种条件下,均可观察到热激励在控制激波和改善激波诱导边界层分离上有着显著的效果。热激励器的放热功率对激波的控制效果有着明显的影响,在本文所考虑的热激励能量范围内,放热功率越大,原激波角的改变越大,最终分离区的改变越明显;展向放置的热激励器数目N影响着输入到流场的能量密度,并且N越大,壁面附近的激波面越趋于平面。尽管N对上壁面沿展向的压力分布无明显影响,但对上壁面分离区大小有明显的影响;对比条件1和4下的计算结果,发现S=0.02m可以得到很好的控制效果,而S=0m时流场结构几乎没有变化,这就表明热激励器必须放置在控制激波上游的一定距离处才会有预期的效果。     

马赫数6柱-裙激波/边界层干扰直接模拟

高超声速激波/边界层干扰比超声速工况下具有更强的可压缩效应,再附之后会形成极高的局部力/热载荷,严重影响飞行器飞行安全。而激波/湍流边界层干扰区附近流动的三维特性使得流动更加复杂而难以预测。采用直接数值模拟对高超声速条件下的柱-裙激波/湍流边界层干扰进行了详细研究,特别是对Görtler涡结构对分离泡、物面压力和热流造成的展向差异开展了定性和定量分析。研究发现,干扰区附近的分离泡结构呈现出明显的三维效应,且Görtler涡展向分离位置截面的分离泡要明显小于再附位置,而这两个截面上分离泡的运动基本同步,没有明显的延迟或超前现象。物面压力和热流在展向出现显著的不均匀性,展向再附位置的平均压力和热流要比展向分离位置分别高13%和16.2%,脉动压力和热流比展向分离位置分别高28%和20%。截面流向速度特征正交分解结果显示两个位置上的能量都集中在剪切层附近,并且展向再附位置上低频模态占有更高的能量。在采用模态重构流场分析分离区面积发现,展向分离位置重构误差更小,而展向再附位置上的重构误差收敛更快。     

一种形状可控的激波增强管道型线设计新方法

激波管所产生的非定常运动激波,若强度和形状能够按照一定的设计要求进行可控条件下的调节,将可望为燃料点火燃烧试验等提供具有独到优势的研究手段。基于激波动力学理论,针对激波管中所产生的平面运动激波,通过设计特定的上下壁面收缩型线,使初始平面运动激波,经收缩段（包括光滑凹形曲线段、斜直线段和光滑凸形曲线段）的变形和强度增加,再以平面波面形状进入较小截面直管段的连续转变过渡,得到了强度增加的平面激波。进一步对所设计的典型型线分别采用数值计算和试验的方法,考核分析激波运动过程中的形状变化,验证了理论方法的可靠性。在此基础上,分析了型线设计的关键参数对激波增强幅度的影响,结果表明,相对于传统激波管方法,本文中所提出的收缩截面方法能更显著地增加平面激波强度;另外,还考察了初始入射激波马赫数对壁面型线和运动激波波面形状的影响,结果表明,对于较强的初始入射激波来说,壁面型线对入射激波强度依赖较小,也就是说,当实际入射激波马赫数即使稍偏离设计状态时,仍然能得到近乎完美的平面形状增强激波。     

跨声速风扇叶片的静态气动弹性问题

使用时域的流固耦合数值计算方法,研究了跨声速风扇叶片在气动力和离心力共同作用下的静态气动弹性问题,分析了叶片在不同工况下的变形规律及叶片变形对整体气动性能的影响.NASA rotor 67的静态气动弹性计算说明气动力对叶片最大变形的贡献达13.07%, 而且叶片变形明显地改变了通道激波的位置和强度.宽弦空心跨声速风扇叶片的静态气动弹性计算说明叶片变形对总体气动效率的影响为0.15%~ 0.5%,其中气动力对变形贡献在叶片尖部的前缘可达41%,考虑气动力引起的变形使得该风扇的流量增大,气动特性线整体向右偏移.计算结果说明:气动力的非线性对跨声速风扇叶片静态变形问题有显著的影响,工程实践中从设计叶型到制造叶型的反扭过程应该采用流固耦合方法以得到更准确的叶型.      

受限空间内激波与火焰作用的三维计算

激波与火焰的相互作用常发生在超声速燃烧和燃烧转爆轰过程中,为深入了解这一现象,采用带化学反应的三维Navier-Stokes方程,对平面入射激波及其反射激波与火焰的作用进行了计算研究,其中燃烧过程使用单步反应模型描述。研究结果显示:在受限空间内模拟激波与火焰作用,能更好地符合实验结果,从而体现出受限空间的三维效应;火焰在入射激波的作用下主要经Richtmyer-Meshkov不稳定而发生变形,此时火焰变形以物理作用为主,燃烧膨胀效果相对不明显;当反射激波与变形火焰再次作用后,火焰迅速膨胀变形,放热率维持在较高水平,此时化学反应过程起主要作用;在反射激波的作用下,变形火焰复杂三维涡结构的形成能强化热量与质量的输送,提高燃烧速率。     

冷气喷射法控制激波强度的数值研究

对冷气喷射时激波受到的影响进行了多方案数值研究.在高压级静叶吸力面反射激波生成点前后5个不同位置上采用相同总压及相同温度的冷气喷射,分析了不同位置的冷气喷射对激波强度和方向影响.结果表明:在吸力面激波折射点附近喷射冷气是一种行之有效的控制激波强度的方法.冷气喷射位置位于激波折射点附近时能够对激波的强度产生影响,在接近激波折射点前部位置注入冷气对流动有积极作用,减弱了激波.      

双斜切双压缩面进气道进口波系结构的数值研究

 双斜切双压缩面进气道具有良好的气动性能和隐身性能。生成了这种三维进气道的计算网格,采用数值计算的方法对双斜切双压缩面进气道进口波系结构进行了研究。通过数值计算,分析了双斜切双压缩面进气道进口的流场特性,总结了斜切角、压缩角、来流马赫数等多种影响参数对进气道进口波系结构的影响规律。     

具有分流叶片跨声速离心叶轮激波及能量损失分析

对具有分流叶片跨声速离心叶轮,采用激波基本关系式对数值信息进行分析和判定,理论证实了叶轮流道中激波的存在,并确定流道内高速气流易在主叶片前缘、分流叶片前缘以及主叶片的吸力面与分流叶片的压力面通道位置形成激波.进一步对比分析了叶顶间隙和运行工况对激波位置分布特征的影响及引起的能量损失.结果表明:固定叶顶间隙下,随着流量减小,激波位置向流道上游移动;叶顶间隙导致激波线位置向流道下游偏移;分流叶片前缘处有最强激波并引起较明显能量损失.      

Tracking characteristics of tracer particles for PIV measurements in supersonic flows

The tracking characteristics of tracer particles for particle image velocimetry (PIV) mea-surements in supersonic flows were investigated. The experimental tests were conducted at Mach number 4 in Multi-Mach Wind Tunnel (MMWT) of Shanghai Jiao Tong University. The motion of tracer particles carried by the supersonic flow across shockwaves was theoretically modelled, and then their aerodynamic characteristics with compressibility and rarefaction effects were evalu-ated. According to the proposed selection criterion of tracer particles, the PIV measured results clearly identified that the shockwave amplitude is in good agreement with theory and Schlieren visu-alizations. For the tracer particles in nanoscales, their effective aerodynamic sizes in the diagnostic zone can be faithfully estimated to characterize the tracking capability and dispersity performance based on their relaxation motion across oblique shockwaves. On the other hand, the seeding system enabled the tracer particles well-controlled and repeatable dispersity against the storage and humidity.     

超声叶栅前缘处的脱体激波预测

为了准确预测超声叶栅前缘处的脱体激波,以Moeckel法为基础,通过分析和公式推导,构造叶栅前缘处的脱体激波模型.首先对Moeckel法进行改进,提高均匀来流条件下的对称脱体激波逼近精度;然后再将Moeckel法推广到均匀来流条件下的非对称脱体激波逼近;最后结合超声叶栅流动特征,给出叶栅前缘处的脱体激波模型.将所得模型用于3个超声叶栅,预测叶栅前缘处的脱体激波形状和位置,并将预测结果与CFD软件求解结果进行比较.结果表明:在均匀来流条件下,改进后的Moeckel法能更准确地逼近对称脱体激波,并可用于逼近非对称脱体激波;由超声叶栅脱体激波模型确定的脱体激波形状和位置与CFD求解结果一致性很好.      

水下爆轰燃气泡形态与激波传播过程研究*(爆震专刊)

针对脉冲爆轰发动机在水下工作过程中形成的燃气射流问题,搭建了水下爆轰燃气实验系统,研究了第一个爆轰循环在水下的燃气泡发展变化过程。建立了基于气液两相双流体模型的脉冲爆轰发动机水下喷射模型,采用时-空守恒元和求解元方法,模拟了爆轰波退化为激波在水中的传播及衰减过程。研究结果表明：燃气泡前期受外界水环境阻滞作用呈现“豌豆状”形态,充分发展阶段气液交界面逐渐失稳,在达到最大尺寸后开始收缩并在中心轴线位置出现凸出的射流；水下爆轰燃气射流发展过程中同时存在脱离燃气泡的水中前导激波和管口燃气泡内的高压区两部分,水中前导激波在传播过程中压力迅速衰减至常压量级,而管口燃气泡内则一直保持较高压力；中心轴线区域气液交界面处反射激波的回传使管口附近出现回击现象,并导致前导激波波阵面上压力峰值逐渐出现在30°方向上。     

水中上升气泡体积变化率的图像分析技术

水中上升气泡的体积变化率是舰船自消隐特种气幕技术等诸多研究的重要基础。鉴于当前对这一体积变化率研究的紧迫需求,提出并较为深入地研究了水中上升气泡体积变化率的图像分析技术。首先,在理论研究的基础上,专门建立了分析计算的数学模型;进而给出了分析的实施方法,即利用摄像法获取水中上升气泡的图像序列,并从中得出所需图像的相关信息,再利用建立的模型即可求出其体积变化率;同时,设计了专门的实验,初步验证了这一分析技术的可行性。     

SPH方法对气液两相流自由界面运动的追踪模拟

气液两相流动是自然界中常见的流动现象,对其进行数值模拟时要求必须能够准确跟踪界面运动变形。本文利用光滑粒子流体动力学方法,结合微可压缩模型(SCM),引入界面控制方法XSPH速度修正以及Van Derwaals状态方程修正,对典型的气液两相流动如二维溃坝、气泡上浮等问题进行了数值模拟,分析了空气和水相互作用机理以及界面运动规律,并同实验结果和其他数值计算结果进行了比较。结果表明,该方法在模拟多相介质界面运动问题准确有效,可用于处理更为复杂的多相流动工程问题。     

Characteristics of pulsed plasma synthetic jet and its control effect on supersonic flow

The plasma synthetic jet is a novel flow control approach which is currently being studied. In this paper its characteristic and control effect on supersonic flow is investigated both experimentally and numerically. In the experiment, the formation of plasma synthetic jet and its propagation velocity in quiescent air are recorded and calculated with time resolved schlieren method. The jet velocity is up to 100 m/s and no remarkable difference has been found after changing discharge parameters. When applied in Mach 2 supersonic flow, an obvious shockwave can be observed. In the modeling of electrical heating, the arc domain is not defined as an initial condition with fixed temperature or pressure, but a source term with time-varying input power density, which is expected to better describe the influence of heating process. Velocity variation with different heating efficiencies is presented and discussed and a peak velocity of 850 m/s is achieved in still air with heating power density of 5.0 · 1012W/m3. For more details on the interaction between plasma synthetic jet and supersonic flow, the plasma synthetic jet induced shockwave and the disturbances in the boundary layer are numerically researched. All the results have demonstrated the control authority of plasma synthetic jet onto supersonic flow.     

Clutter metrics for target detection systems

The problem of calculating an effective contrast for diffused targets in visual scenery is addressed. The effect of local variations of intensity in the vicinity of the target is modeled by considering the distribution of the intensity of the edges. The effect of global clutter is discussed from the vantage point of selecting the best model for background scenery. Generalized contrast measures are described and applied to real scenes     

一种鲁棒的红外与可见光多级景象匹配算法

针对基于特征的景象匹配方法对图像噪声和局部形变适应性差以及对特征质量依赖性强等问题,从提高图像特征描述能力入手,将具有局部光照和对比度不变性的相位一致性变换引入到红外与可见光景象匹配的特征表征中,以获得对噪声和局部形变较强的适应性。在此基础上,首先采用抗旋转的圆投影实现图像粗略匹配,然后基于Zernike矩推导出图像的互相关匹配重构函数并用来剔除错误匹配点,最后基于曲面拟合方法获得亚像素定位,从而实现图像精确与鲁棒匹配。仿真实验验证了所提算法的有效性。