非预混条件下的旋转爆轰燃烧室双波头演化过程数值模拟

针对旋转爆轰燃烧室双波头演化过程中流场结构变化的问题,对非预混条件下的旋转爆轰燃烧室从起爆到形成稳定的双波头过程进行了数值模拟研究。研究结果表明,从起爆到形成稳定爆轰过程,燃烧室主要经历了起爆、爆轰波对撞和稳定爆轰三个阶段；在爆轰波对撞阶段,首次对撞是两个爆轰波间的对撞,由于对撞点处缺少新鲜混合气,从而在对撞结束后衰减为两个压力波。第二次对撞是两个压力波间的对撞,因为在第二次对撞点附近存在新鲜混合气来支撑爆轰波的持续传播,故对撞结束后产生了一个爆轰波和一个较弱的压力波；第二次对撞发生后,燃烧室内的压力波反射叠加并形成局部高压区,此高压区压缩气体使气体温度升高,高温气体引燃混合气后,最终发展成为第二个爆轰波；稳定阶段,两个爆轰波均能稳定自持传播,爆轰波峰面压力可达1.45MPa,波后温度为2500K,爆轰波速度稳定在1738m/s,产生的推力与比冲分别为79.76N和2312.15s；斜激波的存在使燃烧室出口平面流场产生了较大波动。      

Stall flutter prediction based on multi-layer GRU neural network

The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model (ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit (GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation (LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.     

考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。     

基于旋转高频信号注入法的内置式永磁同步电机低速段转子位置检测及其误差补偿

旋转高频信号注入法注入信号较为稳定,且位置估计过程不依赖电机参数,因而十分适用于内置式永磁同步电机(IPMSM)的零、低速转子位置检测。针对传统高频信号注入法无法辨别磁极的问题,用电压方波注入法检测磁极,结合有限元软件仿真,来合理选取方波电压幅值和时长,有效缩短了磁极判断耗时。分析了滤波器和信号离散化对位置估计精度的影响,提出在低速段可用线段拟合带通滤波器中心频率处的相频特性曲线,推导所需补偿角度与电机转速的关系。在理论分析的基础上,采用基于DSP28335的样机平台进行试验,结果表明磁极判断过程稳定,耗时较短,补偿后的位置估计值相比补偿前有明显改善,调速过程中系统动态性能良好。     

跨声速轴流压气机的失速发展机理

为了研究跨声速轴流压气机由近失速工况发展为失速工况的动态演化机理,对跨声速轴流压气机转子NASA Rotor 37进行了多通道全三维数值模拟,着重分析了激波及前缘溢流对失速发展的影响。结果显示:在峰值效率工况下,叶片通道内存在一道斜激波;在失速工况下,斜激波演化为脱体激波。泄漏涡通过脱体激波后发生破碎,破碎的泄漏涡在向通道下游发展的过程中,受逆压梯度的影响,在通道中部形成一个明显的涡结构。在失速工况初期,由于泄漏涡的自维持现象,前缘溢流现象随着叶顶阻塞区的周期性发展而间歇性出现;随着流量的降低,通道阻塞程度逐渐增加,会出现前缘溢流一直存在的现象,这一特性可以作为流场开始急剧恶化的标志。     

当量比对连续旋转爆轰燃烧室特性的影响

为了探究当量比对甲烷-空气连续旋转爆轰燃烧室（CRDC）特性的影响,利用二维可压缩欧拉方程对CRDC进行了数值研究,分析了爆轰波的发展过程和贫燃熄火过程,对比了不同工况下CRDC特性参数的变化情况。结果表明:CRDC起爆后燃烧场在由不稳定状态到相对稳定状态的过程中发生了2次碰撞,当进气当量比较低时,CRDC未能完全发生2次碰撞过程就已经熄爆。随着进气当量比的降低,爆轰波传播速度、轴向平均速度、出口平均温度、出口平均总压均呈下降趋势;增压比随当量比降低而减小的根本原因在于旋转爆轰燃烧过程和等压燃烧过程的熵增差减小,使吉布斯自由能增量差减小。CRDC的燃料驻留时间处于亚毫秒量级,燃烧热效率保持在99%以上。      

旋转整流罩积冰过程中水膜流动和脱离的数值研究

为了更准确地模拟旋转整流罩积冰过程,在现有三维积冰与冰层表面薄水膜流动耦合模型基础上,基于功平衡分析的方法引入了旋转部件表面水膜脱离模型,并发展了相应的计算方法,给出了水膜脱离的判定依据:当气流曳力做的功和由于离心力使水膜增加的潜能之和大于黏附功时整流罩表面的水膜会发生脱离。对旋转整流罩积冰进行数值模拟,计算结果与实验结果吻合得较好,验证了该模型的合理性和计算方法的可行性。之后分析了转速和来流速度对整流罩表面水膜脱离和积冰的影响,结果表明:转速和来流速度越大,水膜发生脱离的比例越大。在研究范围内,转速为3000r/min和6000r/min时,因水膜脱离导致积冰总量分别减少13.4%和15.8%;来流速度为40、50m/s和60m/s时,因水膜脱离导致积冰总量分别减少为12.2%、13.4%、14.2%。      

离心压气机失速模式及自循环机匣处理的作用机制

通过对离心压气机不同转速下实壁机匣和处理机匣两种状态流场的数值模拟,研究了失速模式随转速的变化及不同失速模式下机匣处理作用机制的转变.结果表明:随转速的降低,失速模式由100%设计转速时的扩压器失速逐步向导风轮失速过渡,导风轮失速也经历了由80%设计转速时的端壁失速和叶片前缘失速向50%设计转速的通道大范围完全失速的发展过程.自循环机匣处理对离心压气机失速裕度和效率的影响都与失速模式紧密相关,其扩稳作用仅对80%设计转速时的叶片前缘失速有效.在80%左右设计转速,机匣处理使主流效率有朝大流量工况点变化的趋势,从而使效率降低.当转速低至50%设计转速时,主流通道形成大范围回流区,机匣处理的二次回流可以减少叶轮对通道回流区气流的做功量,从而提高效率.      

低速轴流压气机转子尾迹三维动态测量

对一台低速轴流压气机的转子尾迹进行了三维动态测量,采用锁相平均、方均根和功率谱密度等分析方法对低速轴流压气机转子在节流过程中尾迹流场的稳态和动态特征进行了研究.结果显示:在绝对坐标系下,转子尾迹有高总压和高方均根特征,且在不同叶高的强度不同;转子尾迹与主流的掺混作用导致平均总压升系数的快速衰减;功率谱密度分析表明,在接近失速点时,转子尾迹90%叶高以上出现0.5倍叶片通过频率的特征频带且其功率随着流量减小而增大.      

翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。