三维可压缩混合层中扰动演化的研究

为了利用被动控制增强可压缩混合层中混合效率,本文用数值模拟的方法研究了三维可压缩混合层中扰动向下游的演化.通过在入口处引入不同的不稳定T-S波,给出了扰动随空间的演化过程,及马赫数对扰动演化的影响.计算结果表明在马赫数较小时,扰动波增长很快,当马赫数增大时,扰动所形成的流向涡的强度会减弱.这一结论与流动稳定性分析结果是一致的.     

局部抽吸对边界层流动稳定性的影响和控制研究

在边界层壁面上,设计局部抽吸结构,采用直接数值模拟的方法,获得稳定的三维基本流.在此基础上,研究稳定及最不稳定的二维扰动T-S波的时、空演化机制;进一步探讨了局部抽吸的形式、强度大小及分布结构对二维T-S波的非线性演化影响及其对增长率的贡献大小.结果表明,局部抽吸结构诱导产生的三维基本流是扰动波得以快速增长的一个关键性因素,这是由于平均流剖面的改变及展向速度的出现,增强了流体运动中的不稳定性、扩大了中性曲线的不稳定区域范围.在最不稳定的二维扰动T-S波的非线性演化过程中,由于非线性作用的不断增强,逐渐激发产生出三维扰动波及高次谐波,其三维扰动波的流向波数和频率与二维扰动波的流向波数和频率相同;同时展向速度的大小对二维扰动波的增长、流动的失稳、流向涡的形成等方面都起着激励的作用.随着时、空的不断发展和非线性作用的迅速加强,正、负相间的流向涡逐渐形成,强度逐渐增大,流向涡的影响区域也在不断扩大,涡的形状逐渐拉伸变长,并出现强的剪切层,流动开始失稳等其它机制;这些结论与文献[4、5]的结果相吻合.     

对Klebanoff型转捩点实验测定方法的初步研究

近年来对边界层转捩的实验研究取得了许多进展,并已验证了最初的二维Tollmien Schlichting（T-S）波可以由不同的途径引起转捩。各种途径的区别在于叠加在二维T-S波上的三维扰动沿展向的不同分布。所谓Klebanoff型转捩是最常见到的一种型     

高精度数值方法在DNS中的应用

利用不同格式模拟扰动波在平板边界层中的传播,考察了格式对T-S波传播过程计算结果的影响。研究了不同计算格式对自由剪切层扰动波增长的影响,结果表明:对于二维扰动波模型,采用二阶精度格式,扰动波在平板边界层中只能出现衰减的结果;在自由剪切层中,虽然能模拟出扰动波在传播过程中幅值的增长,但却在相当长的计算范围内看不见三维波的出现。但对于同样的扰动模型,如果采用高精度格式,在边界层中能出现增长的T-S波,并与线性理论非常吻合;在自由剪切层中,扰动波不仅迅速增长,并能在不很长的距离内,明显的看到三维波的出现,并激发出了三维结构。      

激波边界层的相互作用对扰动波传播的影响

利用线性稳定性理论和直接数值模拟研究了带有入射斜激波的、来流马赫数Ma=4．5条件下的平板边界层的失稳特性。重点考察了在由于激波边界层相互干扰，平板边界层上形成分离区，又进而产生激波、膨胀波和旋涡等复杂流动现象的流场上游，引入小扰动的TS波后，扰动波传播通过带有这些复杂流动现象的流场时，扰动波的发展变化特点。通过对流场中扰动波（包括基本波和衍生波）演化特征的分析，研究分离和激波等复杂流动现象对平板边界层稳定性的影响特点。数值模拟发现，激波的出现不同于一般的压缩波，在亚声速区与超声速区对扰动波演化的影响是不同的，此外，分离对扰动有稳定作用。     

马赫数为6的高超声速钝锥湍流边界层空间演化的直接数值模拟

首先求解了来流马赫数为6的零攻角高超声速钝锥边界层的层流基本流场,在选定的计算域入口引入一组有限幅值的T-S波扰动,用高精度差分格式对流动进行了直接数值模拟.引入的扰动触发了转捩,从而得到了空间模式下的湍流边界层.研究了湍流平均场与脉动场的统计特性,给出了相干结构的流动显示图,并对强雷诺比拟的结论进行了检验.     

超声速边界层三维扰动引起小激波的数值研究

通过直接数值模拟的方法,研究了M=4.5的超声速边界层中三维扰动的演化。以某一剖面为入口,加入一个及一对三维T-S波,发现随展向波数的增加,扰动幅值的增长率逐渐减小,证实了M=4.5的超声速边界层中,当三维扰动达到一定幅值时会有小激波出现,为建立可压缩流动稳定性理论提供了依据。     

局部扰动对平板边界层流动稳定性影响的研究

在壁面上构建局部扰动,数值研究不同特征的局部扰动对平板边界层流动稳定性影响的作用机理,进一步探讨壁面扰动的分布、类型、强度大小对三维扰动波在平板边界层流动中的非线性演化规律以及失稳机制的问题.确定什么样的局部扰动能使三维扰动波快速增长或者是什么样的局部扰动能对流动起抑制的作用,同时寻求在三维扰动波向下游发展过程中的传播方向会发生怎样的变化等.对该问题的深入研究,将使人们能深入理解认识平板边界层流动过程中的转捩发生、流动失稳及湍流形成的理论机制.     

超声速小攻角下的圆锥转捩攻角效应

采用高阶差分格式求解三维可压缩滤波Navier-Stokes方程,对超声速零攻角和小攻角圆锥边界层转捩进行了空间大涡模拟研究.在三维Tollmien-Schlichting(T-S)波扰动作用下,计算得到了流场非定常转捩过程.计算结果清晰地展现了非定常大尺度结构演化特征和攻角效应引起的非对称转捩现象,流场发卡涡、横流涡等转捩主控大尺度结构与实验一致.     

钝体头部边界层“逾越”型转捩机理研究

利用数值模拟的方法,研究了再入条件下,钝头体边界层中一种新的失稳机制.这种失稳机制无法用线性稳定性理论解释.最近的研究工作表明,流动中存在着一种与局部有限幅值扰动相关的失稳机制.它与常规的T-S波或二次失稳的过程不同.数值模拟了受空间模式的局部扰动影响的钝体头部超声速边界层的失稳过程.在模拟中,明确了不同扰动模型非线性相互作用的物理机制.      

高精度数值方法与非线性扰动研究

通过数值模拟扰动波在平板边界层中的传播,考察了格式对T-S波传播过程计算的精确性。本文还研究了波-波与涡-涡两种典型的非线性扰动模型,结果表明:在自由剪切层中,波-波模型能加快扰动的增长速度,并激发出了λ结构,但仍然经过了线性失稳阶段;在钝体头部的超声速边界层中,涡-涡模型通过涡-涡相互作用,产生“发卡”结构,进而出现湍流斑。在此过程中,没有线性失稳的迹象。     

锥形液膜空间稳定性分析

对离心喷嘴喷出的锥形液膜做了空间不稳定性分析,导出了锥形液膜的色散方程,并对色散方程进行了数值求解和分析研究。研究结果表明,在液膜表面存在着对称扰动波和反对称扰动波。对于对称扰动波,液膜半锥角的增大使得液膜的不稳定范围和波增长率缩小,并且将扰动波的主波向长波的方向移动。而反对称扰动波既有类似于对称扰动波的变化规律,又有其自身的特点。仅考虑液膜曲率半径的变化时,半锥角的增大使得波增长率的减小趋缓。而同时考虑曲率半径和液膜形状的变化时,液膜的波增长率出现了两个局部峰值。这两个峰值分别对应长短不同的两个波长,并且短波的稳定范围和增长率均远小于长波。这说明锥形液膜的破碎受到两种长短不同的扰动波的共同影响,且在短波的影响下液膜更易破碎。      

Feedback control of boundary layer Tollmien-Schlichting waves using a simple model-based controller

The attenuation of spatially evolving instability Tollmien-Schlichting (T-S) waves in the boundary layer of a flat plate with zero pressure gradients using an active feedback control scheme is theoretically and numerically investigated. The boundary layer is excited artificially by various perturbations to create a three-dimensional field of instability waves. Arrays of actuators and sensors are distributed locally at the wall surface and connected together via a feedback controller. The key elements of this feedback control are the determination of the dynamic model of the flat plate boundary layer between the actuators and the sensors, and the design of the model-based feedback controller. The dynamic model is established based on the linear stability calculation which simulates the three-dimensional input-output behaviour of the boundary layer. To simplify the control problem, an uncoupled control mode of the dynamic model is made to capture only those dynamics that have greatest influences on the input-output behaviour. A Proportional-Integral-Derivative (PID) controller, i.e. a lead-lag compensator, combining with a standard Smith predictor is designed based on the system stability criterion and the specifications using frequency-response methods. Good performance of the feedback control with the uncoupled control mode is demonstrated by the large reduction of the three-dimensional disturbances in the boundary layer. This simple feedback control is realistic and competitive in a practical implementation of T-S wave cancellation using a limited number of localised sensors and actuators.     

恒热流斜板自然对流边界层的稳定性研究

对恒热流条件下的斜板自然对流边界层的稳定性进行了分析及数值模拟.研究发现在一定的倾角和Prandtl数范围内,三维波可存在两种不稳定模态,其临界频率分别接近于重力内波和二维T-S(Tol-mien-Schlichting)波的频率.更为特殊的是,在适当条件下该三维波模态可比二维T-S波更不稳定.这一看似违反Squire定理的结果是由于控制方程中的浮力项造成的.数值模拟初步验证了三维波的临界特性.由于稳定性总是同层流至湍流的转捩过程密切相关,计算结果可为高效换热翅片的设计提供理论依据.      

Breakdown过程长度及其前后边界层名义厚度的变化

采用数值模拟的方法,模拟了平板边界层的转捩过程及发展的湍流,在此基础上,研究了不同初始扰动对转捩持续时间和湍流边界层厚度影响的共性和差异.研究结果表明:在不同的波数情况下,计算得到不同的初始扰动引起的转捩的持续时间不同.在给定基本波数情况下,如果初始扰动在展向不是对称的,那么,它们引起的转捩后平均边界层厚度与转捩起始点处边界层厚度的比值δ_e/δ_b在相当长的时间内表现出相同的增长规律,这一增长规律与所加的初始扰动无关.然而,展向是非对称的表现出不同增长规律.      

富氢燃气旋转爆轰波传播特性实验研究

为研究富氢燃气旋转爆轰波传播特性,利用氢气与氧气预燃烧产生的富氢燃气作为燃料,空气为氧化剂,开展了旋转爆轰实验研究。对富氢燃气旋转爆轰压力变化、时频特性及传播速度等参数进行了分析,研究了不同传播模态下富氢燃气旋转爆轰波传播特性。研究表明：本文实验条件下,富氢燃气与空气旋转爆轰的传播模态主要受当量比影响,当量比高于1.06时呈现单波模态,随着当量比减小,旋转爆轰波呈现单波-双波过渡模态,即同一工况下,单波模态和双波模态交替出现,当量比减小到0.68左右时,基本呈现复杂的双波模态;在270 g/s的空气流量下,当量比增大,旋转爆轰波在环形燃烧室内的传播速度随之提高,但当量比到达临界点以后,传播速度提高不明显;在相同当量比下,当空气流量增大到370 g/s 时,旋转爆轰波的传播速度会进一步提高;空气流量越大,临界点对应的当量比越低,其中270 g/s空气流量对应临界当量比为1.32,370 g/s空气流量对应临界当量比1.16;达到临界当量比以后,传播速度受当量比和空气流量影响不大。