Effect of Martian atmosphere on aerodynamic performance of supersonic parachute two-body systems

Supersonic flows around parachute two-body systems are numerically investigated by solving the compressible Navier-Stokes equations. In the present study, both rigid and flexible parachute models are considered, which comprise a capsule and a canopy. The objective of the present study is to investigate the effects of the Martian atmosphere on the unsteady flows produced by these parachute two-body models and the structural behavior of the flexible canopy. It was found that in the Martian atmosph...     

压气机动静叶干扰的非定常流动特性

采用大涡模拟(large eddy simulation,LES)方法计算动叶尾迹对静叶干扰的流场信息。利用涡量分布揭示动叶尾迹在静叶通道内的演化过程,利用压力梯度识别激波结构及波振源,运用动力学模态分解(dynamic mode decomposition,DMD)方法对静叶通道内流场的时空结构进行模态分解。结果表明:流场中存在3处波振源,分别位于动叶尾缘、静叶前缘和静叶尾缘处;发现静叶通道内流场的频谱具有多峰值现象,模态分解的第1阶流动代表动叶尾迹在通道内随时间迁移,对应频率为动叶通过频率(blade passing frequency,BPF)是通道内旋涡非定常波动的主导频率;第2阶流动是动叶通过频率的2倍频流动,旋涡的空间尺度为1阶模态的1/2,为更小尺度的扰动。      

跨声速翼型上激波／边界层干扰的自适应控制计算

在激波区使用自适应壁对跨声速翼型上的激波/边界层干扰进行控制,可改变机翼的气动性能。这种被动控制可通过在翼型的激波区开一凹腔,其上覆盖一弹性橡胶膜柔壁来实现。本文给出用N-S方程数值模拟这一自适应控制翼型的跨声速粘性绕流,提出了一个适用于本特殊情况(物面边界局部地区在求解过程中有变化)的处理办法。并探讨了自适应柔壁对当代跨声速翼型绕流的影响。     

振动叶栅中尾迹/激波相干的数值研究

研究了尾迹/激波相干对振动叶栅非定常流动的影响.数值方法利用了先进的谐波方法.该方法基于傅立叶变换,把叶栅非定常流变量分解为一时间平均量和若干谐波分量,继而,时域非定常N-S方程被转化为频域定常傅立叶系数方程.定常傅立叶系数方程的求解采用了伪时间推进方法.为加速数值解收敛,采用了本地时间步长和多重网格技术.通过NASA67转子算例结果表明:上游尾迹对振动叶栅非定常流动具有显著影响.当存在尾迹相干时,叶片通道内的激波在轴向方向上表现出较大的漂移运动,叶片所受的非定常力幅值变大,同时失去了时间周期性.      

高压涡轮转子间隙泄漏流动的非定常特征研究

利用数值方法求解三维非定常雷诺平均Navier-Stokes方程模拟某跨声速高压涡轮流场,研究了某跨声速高压涡轮流场的非定常特征,通过详细分析动静干涉对间隙泄漏流动的影响,进一步明确了泄漏流周期性变化的规律和成因.研究发现:静子尾缘燕尾波的外侧分支外尾波是间隙内部流动结构变化的主要原因,间隙泄漏涡的周期性变化则受外尾波和尾迹的共同影响.外尾波深入转子通道内部周期性经过间隙,在间隙前缘附近产生很强的逆压梯度,使间隙前部流动方向明显改变而产生大范围分离.外尾波导致间隙泄漏流量明显增加并周期性震荡.在静子尾迹和外尾波共同作用下,泄漏涡强度出现波动且涡的位置前后移动,使泄漏涡呈现明显的非定常性.      

高压涡轮全环非定常流动数值模拟

为了进一步准确地掌握涡轮内部复杂流动特征,采用全环非定常数值模拟方法研究了某型高压涡轮内部非定常流动,重点分析了上游导叶尾迹和激波在动叶流动中的非定常作用机制,探讨了引起动叶通道内及叶片表面压力分布周期性变化的非定常扰动。数值结果表明:与尾迹作用相比,由导叶激波引起的非定常效应更为明显,尾迹引起的非定常扰动在动叶0.5倍轴向弦长位置下游较为突出。在动叶流场中及叶片表面上的静压变化显示,在动叶通道内还存在低频扰动,对局部流场的非定常扰动明显。最后指出,在非定常计算时采用区域缩放法对流场中非定常扰动的预测存在一定误差。     

低雷诺数NACA0012平面叶栅流场直接数值模拟

采用具有7阶精度的weighted essentially non-oscillatory(WENO)差分格式,直接求解可压缩二维非定常N-S方程组,研究了NACA0012翼型平面叶栅低雷诺数流动的特征.直接模拟及与文献对比的结果表明:叶栅尾缘涡脱落的形成过程与圆柱绕流涡脱落的形成过程非常相似.平面叶栅尾迹区的2阶统计量与孤立翼型尾迹区的2阶统计量具有相同的分布特征,但前者的强度显著大于后者.周期性的涡脱落不仅在上下翼面形成非定常分离,也使得尾迹区某点的总压发生准周期性的变化.随着栅距的减小,翼型上的平均分离位置向前缘移动;尾迹区某点的总压变化频率及其幅值均显著地增加;而且栅距越小,速度脉动2阶统计量反而越大.      

凸拐角附近激波与湍流边界层干扰的数值模拟研究

采用大涡模拟方法数值研究了来流马赫数Ma∞=3、雷诺数Reθ=2070的凸拐角附近激波与湍流边界层干扰问题。对于尖拐角和钝拐角两种情况,计算考察了激波入射到拐角下游时的流场结构、非定常分离流以及湍流统计特性。研究发现:尖拐角能够使流动产生局部集中的顺压梯度,导致壁面摩擦系数在拐角处出现峰值,因此分离点难以越过拐角向上游运动,而再附点发生前后运动。与尖拐角情况不同的是,钝拐角产生的局部顺压梯度不足以阻止分离点越过拐角,因此分离点和再附点都会发生运动。此外,尽管凸拐角产生的顺压梯度能抑制分离点附近的近壁小尺度涡结构,但在经过入射激波作用之后,近壁小尺度涡结构在再附点附近显著增强。     

高负荷跨声速涡轮激波损失机理及控制技术研究

为了指导高性能、高负荷跨声速涡轮的设计,对其叶栅内的流场结构、尾缘波系结构、减小激波损失的机理及其控制技术进行了分析研究。结果表明:跨声速涡轮尾缘流场结构复杂,存在分离膨胀波、分离激波、基底区、再附激波、尾迹、吸力面反射波甚至激波边界层相互干扰等流动现象。通过采用收缩-扩张通道,喉道后采用直线型吸力面,减小吸力面尾缘弯折角、尾缘厚度和尾缘附近局部修型等措施,从而减弱激波强度,减小激波损失。     

锯齿尾缘叶片气动特性和绕流流场的数值研究

以基于NACA 0018翼型的锯齿尾缘仿生叶片为研究对象，采用大涡模拟的方法研究锯齿相对齿宽与相对齿高对锯齿尾缘叶片的气动特性和非定常绕流流场的影响规律和机制.研究表明，尾缘锯齿参数对叶片气动性能的影响是复杂的非线性过程，在一定来流攻角范围内能提高升阻比，但失速提前.如在9.4°~14.8°来流攻角范围内，不同相对齿宽系列叶片的升阻比高于原始叶片，升阻比与锯齿相对齿宽之间没有线性关系.研究还表明，锯齿尾缘能延迟边界层分离，加速尾迹的流动掺混和能量扩散，改变非定常涡结构和涡脱落频率.相对齿高的变化对非定常流动特性的影响更为显著.尾缘锯齿诱导的二次湍流射流和吸力面侧反向涡对改变了原始叶片的绕翼环量，进而影响锯齿尾缘叶片的气动特性和绕流流场特性.      

基于DES方法的高超声速激波/边界层干扰的双微楔控制数值研究

高超声速飞行器在流场中通常会伴随激波/边界层干扰(SWBLI),其引发的流动分离将导致进气道性能下降。采用湍流离散涡模拟(DES)方法、结合有限体积离散方法与自适应网格加密(AMR)技术对来流马赫数为7.0的流场中SWBLI诱导的流动分离进行数值模拟,并分别采用单、双微楔对其进行控制。针对流场结构、近壁面流向速度、压力梯度及总压损失等参数,分析讨论了不同双微楔流向安装位置对SWBLI的控制效果。研究结果表明:双微楔产生的流向涡对与涡对之间的相互诱导促进了各自流向涡对之间的卷吸作用,使得双微楔对分离气泡的消除效果优于单只微楔;流动总压损失系数随着微楔后缘与分离气泡中心的距离的减小呈先减小后增加的趋势;综合讨论流向涡强度与形状阻力的影响,得到了双微楔最佳流向安装位置。     

机翼和叶栅非定常流的Kutta条件

本文在现有的实验研究和数值计算的基础上，对经典Ｋｕｔｔａ条件在非定常流动中的适用性作了进一步探讨，并对文献〔１，２〕中提出的一种新的Ｋｕｔｔａ条件－Ｇｉｅｓｉｎｇ－Ｍａｓｋｅｌｌ（Ｇ－Ｍ）模型进行了分析与比较，本文给出了二维振荡机翼和叶栅绕流的变域变分有限元解法中Ｋｕｔｔａ条件的处理方法。本文旨在为机翼和叶栅非定常绕流的数值计算提供更合理可靠的尾缘边界条件。     

Flow characteristics around airfoils near transonic buffet onset conditions

In transonic flow, buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation. The phenomenon is common in transonic flow, and it has serious impact on the structural strength and fatigue life of aircraft. In this paper, three typical airfoils: the supercritical OAT15A, the high-speed symmetrical NACA64A010, and the thin, transonic/supersonic NACA64A204 are selected as the research objects. The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method, and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD). Qualitative and quantitative analysis of the shock wave motion, shock wave intensity, shock foot bubble and trailing edge separation, and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions. The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack, while the growth rate of the primary mode increases inversely. It was found that at the onset of buffet, there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle, and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.     

Unsteady Flow Variability Driven by Rotor-stator Interaction at Rotor Exit

Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studied in the simulation. Particular attention is attached to the analysis of mechanisms involved in driving rotor wake oscillation, rotor wake skewing and flow angle fluctuation at rotor exit. The results show that smaller axial gap is favorable to enhance the interaction in the region between two adjacent rows, and the fluctuation of the static pressure difference between two sides of rotor wake is improved by potential field from down stator, which is the driving force for rotor wake oscillation. The interaction between rotor and stator is weakened by increasing axial distance, rotor wake shifts to suction side of rotor blade with 5%-10% of rotor pitch, the absolute value of flow angle at rotor exit is less than that in the case of close interspace for every time step, and the fluctuation amplitude is also decreased.     

BVF诊断在非定常转静匹配中的应用

对一台高负荷跨声速单级风扇进行了非定常数值模拟.并将BVF(boundary vorticity flux, 边界涡量流)诊断方法应用于非定常条件下的转静干涉研究中.结果表明BVF不仅能够清晰地捕捉到转静干涉对流场的影响, 指出转子叶排的尾迹干扰是引发下游叶排通道流场非定常特性的最主要因素, 还能够根据叶片表面出现的一些细节, 结合静压、熵等分析方法, 找出动叶的叶尖和静叶的叶根在转静匹配中存在的问题.表明BVF诊断方法能够反映出转静匹配的状态好坏, 并指出改进设计的方向, 是一种能够有效找出对流场产生负面影响的根源, 并指出削弱甚至消除这些负面影响办法的新型诊断方法.      

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

一种基于聚类分析的二维激波模式识别算法

在激波捕捉求解器计算的可压缩无黏流场基础上，提出了一种探测并识别二维激波干扰模式的新算法，从3个层面详细介绍了该算法的实施流程。首先，采用基于当地流场参数设计的传统激波探测方法，辨识出激波附近的一系列网格单元；其次，通过经典的K-means聚类算法将这些激波单元划分成许多簇，并根据簇的相邻信息定义每个簇的类别；最后，设定相关准则对某些紧邻的簇进行合并，进而确定各个激波干扰点的位置，记录各条激波分支所对应的簇，采用Bézier曲线拟合算法分别对其聚类中心进行拟合以获取更加光滑的激波线。数值试验表明，该算法不受网格类型的限制，不仅可以保证最终拟合的激波线具有较高的位置精度，还可以清晰地识别出流场中多激波干扰的模式，同时对分析非定常流场中激波的运动与演化过程也提供了一种有效的可视化手段。     

直升机主,尾桨自由尾迹相互干扰的数值模拟

本文给出了一种新的旋翼和尾桨干扰的数值计算了方法,通过对常值直升机肇翼和尾桨自由尾迹相互干扰的数值模拟,给出了干扰对旋翼和尾桨非定常气动力以及相应的非定常尾迹几何形状的影响,数值结果表明,悬停时尾桨叶面非定常压力脉动主要是由于在旋翼下洗流的诱导下,尾桨与其自身尾涡于干扰引起的,随着旋翼和尾桨间距的增加,旋翼对尾桨气运特性的影响逐渐减弱。     

喷气流与栅后流场掺混及干扰效应的探讨

应用有效的冷气喷射模型，结合涡轮叶栅三维流场计算程序，对叶片尾缘对开缝喷气进行了数值模拟，从栅后速度场、总压场以及气流角分布角度，研究了射流与主流的掺混及干扰过程，结果表明：采用对开缝喷气的涡轮叶栅，在尾缘附近射流与尾迹的掺混占主导地位，经过一段距离后，尾迹与主流的掺混起主要作用，而且后缘对开缝喷气对尾迹的流动结构及漂移过程有一定的影响。     

汽轮机末级三维非定常流动数值模拟

采用三维粘性非定常数值模拟方法,对某汽轮机设计状态末级非定常流动进行了模拟,并对设计状态下各因素引起末级动叶表面产生非定常力的大小和作用机理进行了分析。结果表明,本文中的汽轮机末级在设计条件下,上游静叶尾缘激波和位势场对下游转子的影响比尾迹作用的影响更大,静叶尾缘激波和位势场导致下游转子叶片表面非定常力大小呈周期性变化;此外上游静子尾迹以及动叶尾缘出口膨胀波等也是产生转子叶片表面非定常力的因素。