逆流矢量喷管流动自身导致的非定常现象

基于数值方法,研究了逆流推力矢量喷管中由于流动自身导致的非定常现象。结果表明,当抽吸背压下降至某个范围时(50.66~60.8 kPa),主流在接近外套壁面的某个位置处呈现出周期性振荡,流动参数(抽吸二次流量、推力矢量角、合成推力系数)也以小振幅振荡,但是同向二次流量却没有振荡。深入的流动机理分析表明,此时在逆流剪切层中,分离涡周期性地形成和发展,没有呈现出流动失稳,同时剪切层中的速度比又远大于定常流动下的值,因此流动表现为自激振荡特征。     

超声速对流推力矢量喷管的数值模拟

利用二次流对流剪切激励主射流方式来产生射流对流推力矢量方法是一种效率很高的矢量化喷流新方法,又称为对流推力矢量法。本文运用三阶精度MUSCL TVD有限体积格式求解二维非定常、可压缩、雷诺平均N-S方程和Spalart-Allmaras湍流模型,计算了对流推力矢量喷管的超声速射流问题,获得几种不同二次流吸气压强下对狭缝和喷管肩部壁面压强分布、剪切层诱导旋涡的结构特征和主射流矢量角等结果,计算结果与实验数据对比取得良好的一致性,并对计算结果进行了详细分析和讨论。     

矩形射流矢量喷管数值模拟研究

根据流场特性设计了1个矩形射流矢量喷管。提出了在喷管内加入射流,促使喷管主流速度增高、和形成出口矢量角的设想;采用FLUENT流体计算软件,对射流矢量喷管内的二维内流场进行了模拟,分析了喷管出口位置对射流矢量喷管内流场、推力、出口气流方向和总压损失等气动参数的影响,并分析了其在不对称射流情况下,对喷管射流位置与气流方向、推力和总压损失等气动参数的影响,确定了第2道射流的位置。     

轴对称矢量喷管外流对内流干扰研究

应用有限体积法在同位网格中对控制方程进行离散化,对某一轴对称矢量喷管的内外流场进行了数值模拟,数值模拟和试验结果进行了比较,分析了该喷管在矢量偏角20°下,不同的落压比、不同的外流速度对轴对称矢量喷管内流的影响。研究结果表明,喷管内的分离激波是影响轴对称矢量喷管性能的重要因素。      

逆流推力矢量喷管基本流动特征的数值研究

 利用数值模拟的方法,通过对逆流方案中喷管气动性能的研究验证了逆流推力矢量方案的可行性。在非矢量状态下主喷管出口截面上的流量系数和推力系数分别达到99.2%和98.8%;矢量化状态下最大推力矢量角超过了20°,而推力系数与非矢量状态下的比较下降不超过3.7%,且最大抽吸二次流量比仅为2.1%。此外,对该方案中一些基本的流动特征进行了分析,得到了抽吸二次流量比与推力矢量角的变化关系所揭示的流场结构,并对此进行了详细解释,同时揭示了逆流剪切层强烈的湍动特性和大涡结构的特点。     

增强型双喉道射流推力矢量喷管的流动特性试验

对一种增强型双喉道射流推力矢量喷管开展了内部流动特性的试验研究,获得了其在不同次流压比状态下的内流结构和沿程静压分布.试验结果显示:在基准双喉道矢量喷管尾部附加扩张段后,能够以2.8%的次流消耗率获得超过20°的平均气流偏角,这表明通过附加扩张段来增加喷管矢量角的设计概念是可行的.在凹腔内,增强型双喉道射流推力矢量喷管的静压分布规律与基准双喉道矢量喷管一致,但在附加的扩张段内,下壁面的压强要明显高于上壁面,这正是其推力矢量角得到显著增大的原因.随着次流压比的增加,喷管获得的推力矢量角单调增加,但是喷管附加扩张段的矢量增强效果基本维持不变.      

双喉道推力矢量喷管的内流特性研究

为了研究双喉道推力矢量喷管(DTN)在非推力矢量和推力矢量情况下的内流特性,基于数值模拟的方法,计算分析了不同几何参数和气动参数对DTN的影响。结果表明,DTN在非推力矢量时,仅在落压比(ZNPR)为3～4之间才具有较高的内流性能(推力系数达0.97,流量系数为0.94),当落压比增加时,推力系数迅速下降。在推力矢量时,DTN可以获得很大的推力矢量效率(当落压比为4,引射量为3%时达到4),且推力系数也较高(0.94以上),其综合性能优于单喉道偏移和激波操纵式矢量喷管。二次流量、落压比、凹腔扩张角和收敛角、引射角度都对推力矢量状态下的DTN内流性能有着不同的影响。为了实现DTN在推力矢量和非推力矢量下都有较好的内流综合性能,所建议的设计参数为:落压比为3～4,引射量为3%,凹腔扩张角为10°左右,收敛角在20°～30°,引射角度为30°逆流引射角(β=30°)。     

轴对称激波诱导矢量喷管三维两相流数值模拟

运用颗粒轨道模型对轴对称激波诱导矢量喷管内流场进行了数值模拟,并与纯气相条件下的计算结果进行了对比分析,研究了1～30μm不同直径颗粒的运动轨迹和内流场参数的分布特征,并对矢量喷管的推力性能进行了研究。结果表明,颗粒相的加入对射流推力矢量喷管内流场和推力性能都带来很大影响,影响规律与颗粒直径密切相关。为两相流条件下射流矢量喷管的试验研究提供了一定依据。     

外流马赫数对逆流推力矢量喷管性能的影响

为了研究外流马赫数变化对逆流推力矢量喷管内部流场结构和性能特点的影响,参考实验,数值模拟了不同马赫数下的流场结构、喷管矢量角和逆向二次流质量流量系数。结果表明:外流马赫数从0增加到0.8的过程中,喷管内流场结构会随着外流马赫数的增大而不断变化从而改变主流两侧压力差,进而导致矢量角逐渐减小,与此同时二次流质量流量系数快速增大;外流马赫数从0.8增大到1.5的过程中,流场结构和矢量角随马赫数变化没有明显变化,二次流质量流量系数缓慢增大。     

轴对称矢量喷管三维传热计算研究

采用封闭腔理论算法，利用矢量喷管的三维内流场数值计算结果，对轴对称矢量喷管的壁温分布进行了数值计算。计算结果表明，在非矢量状态，喷管同一轴向位置处的周向壁温基本相同，壁温分布可按一维处理。在矢量状态，喷管同一轴向位置处的周向壁温相对差值达15％，壁温分布呈现复杂的形式，必须采用三维计算方法计算壁温分布；喷管扩张段壁温明显高于非矢量状态，所开发的计算方法和程序既可用于轴对称矢量喷管的工程设计，也可为试验验证提供理论依据。     

倾转涵道倾转过渡阶段的非定常气动力

为研究倾转涵道动力装置在倾转过渡阶段的非定常气动力,使用基于滑移网格技术的非定常计算方法,利用基于内场信息和叶素理论的压力盘模型模拟风扇螺旋桨,通过求解Navier-Stokes(N-S)方程,对涵道风扇俯仰拉起过程进行数值模拟.结果表明:在倾转过程中涵道风扇非定常气动力的迟滞特性明显,俯仰角速度延迟了流场分离,增加了涵道风扇的升力和阻力;低速飞行时,涵道风扇在整个倾转过渡阶段气动性能优良;在高速大迎角飞行时,涵道风扇气动性能恶化,俯仰力矩曲线紊乱,不利于进行倾转过渡飞行.      

变展长战术导弹非定常气动特性研究

以某轴对称基准弹为模型,通过数值模拟对其进行了准定常和非定常气动特性研究,探索可变翼飞行器在变展长过程中的非定常气动特性。结果表明:在弹翼连续变展长的过程中,非定常气动特性曲线在对应的准定常气动特性曲线周围形成滞回环,若来流马赫数增大,则滞回环变小。当来流为超声速时,轴对称基准弹飞行迎角和变形周期对其非定常气动特性影响较小,而非定常气动特性并不显著,接近于准定常气动特性,在可变翼飞行器设计过程中可以忽略不计;当来流为亚声速时,展长变化对非定常气动特性有一定的影响,不能忽略。     

超声速进气道反压引起的不起动/起动现象研究

为了探索节流反压变化对进气道起动性能的影响,在来流马赫数为2.0的情况下,采用定常/非定常数值仿真方法对出口可移动壁面运动速率为12,24,36mm/s的进气道进行了研究,分析了不同可移动壁面运动速率引起的堵塞度变化对进气道流场特性及气动性能的影响。结果表明:采用非定常数值仿真获得的进气道性能参数和流场特性变化规律与定常数值仿真时一致,定常/非定常数值仿真时的进气道性能参数均存在"迟滞回路"现象;进气道出口堵塞度变化速率越大,"迟滞回路"现象越明显;堵塞度变化速率越大,进气道进入不起动状态时的堵塞度越高,进气道再次进入起动状态时的堵塞度越低。     

Mechanism analysis of hysteretic aerodynamic characteristics on variable-sweep wings

Variable-sweep wings have large shape-changing capabilities and wide flight envelops, which are considered as one of the most promising directions for intelligent morphing UAVs. Aerodynamic investigations always focus on several static states in the varying sweep process, which ignore the unsteady aerodynamic characteristics. However, deviations to static aerodynamic forces are inevitably caused by dynamic sweep motion. In this work, first, unsteady aerodynamic characteristics on a typical variable-sweep UAV with large aspect ratio were analyzed. Then, deep mechanism of unsteady aerodynamic characteristics in the varying sweep process was studied. Finally, numerical simulation method integrated with structured moving overset grids was applied to solve the unsteady fluid of varying sweep process. The simulation results of a sweep forward-backward circle show a distinct dynamic hysteresis loop surrounding the static data for the aerodynamic forces. Compared with the static lift coefficients , at the same sweep angles, dynamic  lift coefficient in sweep forward process are all smaller, while dynamic sweep backward  lift coefficient are all larger. In addition, dynamic deviations to static  lift coefficient are positively related with the varying sweep speeds. Mechanism study on the unsteady aerodynamic characteristics indicates that three key factors lead to the dynamic hysteresis loop in varying sweep process. They are the effects of additional velocity caused by varying sweep motion, the effects of flow hysteresis and viscosity. The additional velocity induced by sweep motion affects the transversal flow direction along the wing and the effective angle of attack at the airfoil profile. The physical properties of flow, the hysteresis and viscosity affect the unsteady aerodynamic characteristics by flow separation and induced vortexes.     

双旋翼微型飞行器的非定常气动特性

发展了一种基于计算流体力学(CFD)的非定常气动特性预测方法,计算方法包括了动量源模型、预处理方法、非结构嵌套网格和Spalart-Allmaras(S-A)湍流模型等技术.通过计算悬停Caradonna-Tung算例和俯仰振荡NACA0012算例,验证了计算方法模拟双旋翼微型飞行器动态流场的有效性.数值模拟了双旋翼微型飞行器动态流场,给出了非定常气动系数的迟滞曲线,分析了缩减频率、前飞速度和螺旋桨转速对非定常气动特性的影响.计算结果表明:力矩系数迟滞效应随缩减频率、前飞速度和螺旋桨转速增大而增大,升力系数迟滞效应随缩减频率和前飞速度增大而增大,但基本不随螺旋桨转速变化而变化.      

三角翼俯仰振荡中若干参数影响的实验研究

本文研究一个 6 0°三角翼以不同振荡频率 ,不同振幅从不同起始迎角进行俯仰振荡的纵向气动力特性 ,研究发现跨越具有不同时间尺度的流态的振荡会产生明显的气动力迟滞现象 ,这一迟滞现象随减缩频率的增大有一个极限值 ;在同一种流态下振荡 ,气动力迟滞现象不明显 ,这一结果对充分利用非定常气动力有指导意义。     

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

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

离心叶轮和扩压器相互作用

采用定常和非定常数值模拟技术,研究了某跨声离心压气机级内的复杂流动,分析了叶轮与径向扩压器之间的非定常相互作用.结果表明,扩压器的势干扰使叶轮出口静压发生周期性地变化,引起叶轮后部负荷、叶尖泄漏和损失等也产生周期性地变化,影响范围达30%流向弦长.叶轮出口流场的周向不均匀性及非定常性使扩压器进口参数发生周期性的变化,导致扩压器负荷、总压恢复系数分布呈现出较强的非定常性,影响范围贯穿整个扩压器通道.      

一种定几何混压式二元进气道的再起动特性研究

针对一种设计飞行马赫数范围为2.25～4.0的定几何混压式二元超声速进气道由于反压引起的不起动开展了再起动特性风洞实验研究.根据实验结果和二维非定常数值仿真结果分析表明:在Ma=2.51和3.01下进气道性能在进锥/退锥实验中规律一致;按照Ma=2.25的起动面积收缩比确定的喉道面积,使进气道在来流Ma≥2.51时具有自起动能力;而在稳定亚临界状态下具有高的总压恢复系数及形成类似于外压式的气动通道是进气道无回路迟滞现象的主要原因.      

耦合升沉运动的七叶侧斜螺旋桨空泡性能研究

为建立有效的预报耦合升沉运动的螺旋桨空泡性能的方法,并基于该方法分析升沉运动对螺旋桨空泡性能的影响,基于RANS方法对耦合升沉运动的七叶侧斜螺旋桨的空泡性能进行了数值模拟,螺旋桨的升沉和旋转运动耦合采用自定义运动方程实现,其中升沉运动简化为简单正弦函数。非定常流场中网格和物理信息的传递采用重叠网格技术实现。将非定常流动中螺旋桨(DTMB 4381)的空泡计算结果与试验结果及现象进行了对比,验证了数值计算模型的准确性。计算中对不同升沉运动周期下的螺旋桨非定常推力及扭矩系数进行了分析,同时对螺旋桨的空泡性能进行了实时监测。计算结果显示,升沉运动使螺旋桨的推力和扭矩系数的非定常特性更加明显,同时导致螺旋桨各叶片上的空泡分布不均匀;升沉运动周期越小,上述现象越明显。同时较小升沉运动周期时,推进系数和空泡面积均为周期性变化,其变化周期均为升沉运动的1/2,螺旋桨空泡面积最小值较无升沉运动时增大约17%,空泡面积最大值则增大约57%。