自激振荡脉冲射流喷嘴装置系统频率特性理论研究

根据相似系统原理和流体网络理论建立了自激振荡脉冲射流喷嘴装置的等效网络模型,利用系统传递函数推导了系统频率特性方程并进行了数值计算。结果表明:喷嘴装置的固有频率主要由喷嘴形状、结构参数、入口流速、射流中压力扰动波波速决定;自激振荡腔腔径、自激振荡腔腔长、上喷嘴直径、下喷嘴直径都对系统频率特性影响很大。提出了相应的自激振荡脉冲射流喷嘴设计准则,即喷嘴装置在最佳阻尼比下产生谐波共振。     

带喷流激波针流动特性实验研究

采用动态测力、动态测压和纹影等风洞实验技术,对加装了带喷流激波针的钝头体的绕流特性、稳定和非稳模态的形成条件和机理进行了研究.结果表明:带喷流激波针流场存在稳态和非稳态两种模态,超声速喷流的压比大于临界压比时流动处于稳定模态,反之则为非稳模态;增大激波针长度可减小钝头体阻力,但达到一定长度后,进一步减阻的效果不再显著;增大喷流压比能够有效减弱再附激波强度,有利于缓解单独激波针的肩部热斑问题;非稳模态下波系自激振荡对再附激波在钝头体表面所围的区域影响剧烈,振荡是周期性的,且存在确定的主导频率,主导频率随喷流压力比增大而减小;自激振荡的产生是由于喷流出口周围的反压在喷流压比小于临界压比时无法获得持续的平衡而导致.      

自振脉冲气液射流振动分析

从提高射流利用率出发, 在自振脉冲水射流喷嘴振荡腔上开孔自吸入空气, 以形成一种自振脉冲气液射流, 研 究其振荡频率以及打击力的实际影响因素.基于水声学与流体动力学原理建立自振脉冲气液射流频率模型, 并以悬臂梁为靶物, 通过考察射流冲击悬臂梁时的振动特性获得其时域及频域特征.结果表明, 自振脉冲气液射流的振荡频率随腔长增加单调递减, 随含气率变化存在极小值, 且在初始段发生骤变, 实测值与理论值相吻合;对应自振脉冲气液射流产生最大谐振峰值的最优腔长略小于自振脉冲水射流, 且前者最大谐振峰值较后者明显提高.     

Flow induced oscillation of two rigid rectangular plates in a side-by-side configuration

The present work is an experimental study of two oscillating rigid plates placed in side-by-side configuration, hinged at their leading edges, subjected to low subsonic flow. This problem is investigated using smoke-wire flow visualization, hot-wire anemometry, and time resolved particle image velocimetry. It is found that beyond a critical Reynolds number, the plates set into oscillatory motion. This critical Reynolds number depends on the gap between the plates. It is also seen that this value of Reynolds number, at lower values of gap to thickness ratio (<7) is significantly higher than that of the single plate configuration value. The frequency and amplitude of the oscillating plates at various gaps and Reynolds numbers have been studied and compared with the characteristics of an oscillating single plate. It is also found that depending on the gap and acceleration of the free-stream, there exist two modes of oscillation - (i) in-phase and (ii) out-of-phase. For gap to thickness ratio less than 10, only in-phase oscillations take place for all values of free-stream velocity considered in the present work, whereas, when this ratio is greater than 10, the mode of oscillation depends on the initial conditions up to a certain free-stream velocity, beyond which the plates switch to in-phase mode. Smoke wire flow visualization technique along with time resolved particle image velocimetry reveal that the vorticity distributions around the plates are responsible for the initiation of the two modes of oscillations.     

Coherent structures in countercurrent axisymmetric shear flows

The dynamical behaviors of coherent structures in countercurrent axisymmetric shear flows are experimentally studied.The forward velocity U1 and the velocity ratio R=(U1-U2)/(U1+U2),where U2 denotes the suction velocity,are consldered as the control parameters.Two kinds of vortex structures,i.e.,axisymmetric and helical structures,were discovered with respect to different reginmes in the R versus U1 diagram .In the case of U1 rangjing from 3 to 20m/s and R from 1 to 3,the axisymmetric structures plan an important role.Based on the dynamical behaviors of axisymmetric structures,a critical forward velocity U1cr=6.8m/s was defined,subsequently,the subcritical velocity regime:U1&gt;U1cr and the supercritical velocity regime:U1&lt;U1er,In the subcritical velocity regine,the flow system contains shear layer self-excited oscillations in a certain range of the velocity ratio with respect to any forward velocity.In the supercritical velocity regime,the effect of the velocity ratio could be explained by the relative movement and the spatial evolution of the axisymmetric structure undergoes the following stages:(1) Kelvin-Helmholtz instability leading to vortex rolling up,(2) first time vortex agglomeration.(3) jet colunn self-excited oscillation,(4) shear layer self-excited oscillation,(5)“ordered tearing“,(6) turbulence in the case of U1&lt;4m/s (the “ordered tearing“ does not exist when U1&gt;4m/s),correspondingly,the spatial evolution of the temporal asymptotic behavior of a dynamical system can be described as follows:(1) Hopt bifurcation,(5) chaos(“weak turbulence“)in the case of U1&lt;4m/s(superharmonic bifurcation does not exist when U1&gt;4m/s).The proposed new terms,superharmonic and reversed superbarmonic bifurcations,are characterized of the frequency doubling rather than the period doubling.A kind of unfamiliar vortices referred to as the helical structure was discovered experimentally when the forward velocity around 2m/s and the velocity range from 1.1 to 2.3,There are two base frequencies contained in the flow system and they could coexist as indicated by the Wigner-Ville-Distribution and the temporal asymptotic behavior of the dynamical system corresponding to the helical vortex could be described as 2-torus as indicted by the 3D reconstructed phase trajectory and correlation dimension.The scenario of the spatial evolution of helical structures could be described as follows:the jet column is separated into two parts at a certain spatial location and they entangle each other to form the helical vortex until the occurrence of those separated jet columns to reconnect further downstream with the result that the flow system evolves into turbulence in a catastrophic form.Correspondingly,the dynamical system evolves directly into 2-tiorus through the supercritical Hopf bifurcation followed by a transition from a quasi-periodic attractor to a strange attractor.In the case of U1=2m/s,the parametric evolution of the temporal asymptotic behavior of the dynamical system as the velocity ratio increases from 1 to 3 could be described as follows:(1)2-torus(Hopf bifurcation),(2) limit cycle(reversed Hopf bifurcation),(3) strange attractor (subbarmonic bifurcation).     

开放剪切流空间全局动力学与空间相位斑图间的关系

开放流动空间动力学可基于两类全局能量关系式进行研究;而空间相位斑图则可通过互谱空间演化加以测定。全局能量关系式以时间Fourier系数的形式建立流场任意两点问速度脉动能量间的关系,籍此可定义全局意义上的线性、非线性和线性一非线性机制。基于轴对称剪切流、变密度轴对称圆射流以及平面对称剪切流的实验发现：轴对称旋涡结构的配对由线性、线性一非线性机制表征,对应有序空间相位斑图;并且能量可通过线性一非线性机制在具有相同相速度的扰动间传递。螺旋结构由线性机制表征,对应有序相位斑图。全局自激励振荡由非线性的能量共振表征,对应无序相位斑图。籍此,有序空间相位斑图对应线性和线性一非线性机制;而混沌相位斑图则对应非线性机制。     

THE OSCILLATION MODE OF INCOMPRESSIBLE CAVITY FLOW

采用二维大涡模拟方法进行了空腔水流场的数值计算, 考察空腔前缘动量损失厚度及来流速度等因素如何影响空腔流的振荡, 同时考察了空腔长深比与空腔流振荡模式的关系. 用空腔水流场的粒子图像测速测量结果验证了数值计算的可信性.结果表明, 空腔水流是否发生振荡取决于壁面摩擦速度.瞬时涡结构和空腔阻力系数2个方面的特征显示空腔水流场有2种典型的振荡模式, 剪切层模式与尾流模式, 确定振荡模式的关键因素是空腔长深比.     

Three-dimensional numerical simulation of the wake flow of an afterbody at subsonic speeds

We numerically investigate the wake flow of an afterbody at low Reynolds number in the incompressible and compressible regimes. We found that, with increasing Reynolds number, the initially stable and axisymmetric base flow undergoes a first stationary bifurcation which breaks the axisymmetry and develops two parallel steady counter-rotating vortices. The critical Reynolds number (Re _cs) for the loss of the flow axisymmetry reported here is in excellent agreement with previous axisymmetric BiGlobal linear stability (BiGLS) results. As the Reynolds number increases above a second threshold, Re _co, we report a second instability defined as a three-dimensional peristaltic oscillation which modulates the vortices, similar to the sphere wake, sharing many points in common with long-wavelength symmetric Crow instability. Both the critical Reynolds number for the onset of oscillation, Re _co, and the Strouhal number of the time-periodic limit cycle, St_sat, are substantially shifted with respect to previous axisymmetric BiGLS predictions neglecting the first bifurcation. For slightly larger Reynolds numbers, the wake oscillations are stronger and vortices are shed close to the afterbody base. In the compressible regime, no fundamental changes are observed in the bifurcation process. It is shown that the steady state planar-symmetric solution is almost equal to the incompressible case and that the break of planar symmetry in the vortex shedding regime is retarded due to compressibility effects. Finally, we report the developments of a low frequency which depends on the afterbody aspect ratio, as well as on the Reynolds and on the Mach number, prior to the loss of the planar symmetry of the wake.     

圆形月池流激振荡实验研究

采用基于经验模态分解(EMD)的本征模态函数(IMF)分析方法和对比研究的实验方 法,实验研究了圆形月池在均匀流条件下的流激振荡特性. 通过实验研究发现,尽管腔内的 流体在不同来流的条件下运动趋势有所不同,但归纳起来对于月池结构,其内部存在两类振 荡源,一类是腔内流体的固有振荡,该频率与腔体结构形式、尺寸、腔内流体深度有关,与 来流无关;另一类则是流体流经腔口时产生的流体动力振荡,该频率与来流特性、速度有关. 当二者频率相近时将会激起腔内流体的共振,即会有``活塞(piston)'现象产生. 实验中,还在腔体侧壁上布置了压力传感器,用来分析腔体侧壁受到的流体压力特性及 其变化趋势.     

Mass transfer characteristics in an axisymmetric wavy-walled tube for pulsatile flow with backward flow

Under the pulsatile flow with backward flow (PFBF) conditions, flow mixing and mass transfer characteristics are experimentally investigated in an axisymmetric wavy-walled tube at a net flow Reynolds number from 50 to 1,000. An electrochemical technique is employed to measure the mass transfer rate. An optimal Strouhal number corresponding to the peak value of the mass transfer enhancement factor is observed, which is independent of the oscillatory fraction of the flow rate, but decreases with the increasing net flow Reynolds number. It was found that the mass transfer enhancement under PFBF has the similar characters of resonant enhancement in two-dimensional (2-D) channels, but there also exists an essential difference since no self-sustained oscillation occurs in the wavy-walled tube.     

Flow instability of nanofuilds in jet

The flow instability of nanofluids in a jet is studied numerically under various shape factors of the velocity profile, Reynolds numbers, nanoparticle mass loadings,Knudsen numbers, and Stokes numbers. The numerical results are compared with the available theoretical results for validation. The results show that the presence of nanoparticles enhances the flow stability, and there exists a critical particle mass loading beyond which the flow is stable. As the shape factor of the velocity profile and the Reynolds number increase, the flow becomes more unstable. However, the flow becomes more stable with the increase of the particle mass loading. The wavenumber corresponding to the maximum of wave amplification becomes large with the increase of the shape factor of the velocity profile, and with the decrease of the particle mass loading and the Reynolds number. The variations of wave amplification with the Stokes number and the Knudsen number are not monotonic increasing or decreasing, and there exists a critical Stokes number and a Knudsen number with which the flow is relatively stable and most unstable,respectively, when other parameters remain unchanged. The perturbation with the first azimuthal mode makes the flow unstable more easily than that with the axisymmetric azimuthal mode. The wavenumbers corresponding to the maximum of wave amplification are more concentrated for the perturbation with the axisymmetric azimuthal mode.     

Flow due to parallel disks rotating about non-coincident axis with one of them oscillating in its plane

An exact solution of the Navier–Stokes equations is obtained for the flow between two eccentric disks rotating with the same angular velocity and one of them executing non-torsional oscillations. An analytical solution describing the flow at large and small times after the start is given. The solutions depend on the ratio of the frequency of oscillation to the angular velocity of the disks and the ratio of the amplitude of oscillation to the angular velocity of the disks and to the distance between the axes of rotation, and the Reynolds number based on the distance between the disks and the angular velocity of the disks. The solutions for three cases when the angular velocity is greater than the frequency of oscillation or it is smaller than the frequency or it is equal to the frequency are discussed.     

Oscillation of shallow flow past a cavity: Resonant coupling with a gravity wave

Self-excited oscillations of flow past a cavity are generated in a shallow free-surface system. The shear layer past the cavity opening has two basic forms: a separated free-shear flow; and a shear flow along a slotted plate. Instabilities of these classes of shear flows can couple with the fundamental gravity-wave mode of the adjacent cavity. The dimensionless frequencies of both types of oscillations scale on the length of the cavity opening, rather than the gap distance between the slats, i.e., a large-scale instability is always prevalent. A technique of high-image-density particle image velocimetry allows acquisition and interpretation of global, instantaneous images of the flow pattern, including patterns of vorticity and Reynolds stress correlation. Use of a cinema approach provides representations of the timewise evolution of the global, instantaneous flow structure, and thereby definition of the amplitude peaks and phase angles of the coupled fluctuations via auto- and cross-spectral techniques. These methods, along with global, averaged representations of the fluctuating flow field, provide insight into the onset of fully coupled (phase-locked) oscillations of the shear flow past the resonator cavity. The common, as well as the distinctive, features of the resonant-coupled instability of the shear flow past the slotted plate are characterized, relative to the corresponding coupled instability of the free-shear layer. Varying degrees of resonant coupling between the unstable shear layer and the adjacent resonator are attained by variations of the inflow velocity, which yield changes of the predominant oscillation frequency, relative to the resonant frequency of the adjacent cavity. Well-defined, coherent oscillations are indeed attainable for the case of the shear flow along the slotted plate, though their amplitude is significantly mitigated relative to the case of a free-shear layer. The degree of organization of the self-excited, resonant-coupled oscillation and the manner in which it varies with open area ratio and geometry of the plate are interpreted in terms of the flow structure on either side of, and within, the slotted plate; these features are compared with the corresponding structure of the free-shear layer oscillations.     

来流脉动对自激振荡脉冲流的影响

基于瞬变流和流体网络理论建立了由往复柱塞泵,管路和自激振荡呀嘴组成的装置系统的数学模型,分析了泵源脉动对自激振荡脉冲射流振荡特性的影响,表明当来流脉动频率小于或等于嘴嘴装置的固有频率时,能提高射流的振荡,冲蚀效果,理论分析与实验结果吻合,其结论对合理设计自激振荡射流喷嘴有一定的指导意义。     

Asymptotic solution of the two-dimensional equations for a thin viscous rarefied shock layer on a cold surface

Hypersonic rarefied flow past blunt bodies is studied in the continuum-free-molecular transition regime. On the basis of an asymptotic analysis three rarefied gas flow patterns are established depending on the relation between the relevant parameters of the problem. In the first regime corresponding to a cold surface asymptotic solutions of the equations of a thin viscous shock layer are derived at low Reynolds numbers in the axisymmetric and plane cases. Simple analytical expressions for the pressure and the heat transfer and friction coefficients are obtained as functions of the freestream parameters and the body geometry. With decrease in the Reynolds number the coefficients approach the values corresponding to free-molecular flow. In this regime a similarity parameter for the hypersonic rarefied flow past bodies is determined. The asymptotic solutions are compared with numerical solutions and the results of direct statistical simulation by the Monte Carlo method.     

Three-dimensional vortex wake structure of a flapping-wing micro aerial vehicle in forward flight configuration

The unsteady flow field past a backward-facing step in a rectangular duct is investigated by adopting time-resolved particle image velocimetry (PIV) in the Reynolds number range of 2,640–9,880 based on step height and the inlet average velocity. The PIV realizations are subjected to post-processing techniques, namely, proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD). At low Reynolds numbers, the second spatial POD modes indicate the presence of the shear layer mode, whereas this feature shifts to higher modes at higher Reynolds numbers. The corresponding temporal modes are Fourier-transformed to obtain the dominant frequency, whose Strouhal number corroborates the above observation. Short-time windows in the transverse velocity component along the shear layer are selected to investigate the temporal stability of the flow field by DMD to quantify the growth rate of the shear layer mode. The higher harmonics of this mode are also observed to grow, albeit at lesser rate. By relating to POD analysis, the most energetic structures were found to correspond to the unstable modes. The correlation between these unstable DMD modes and the Fourier-filtered flow fields for the same frequencies indicate better match for the lower operating Reynolds number case as compared to higher ones. The spatial stability analysis demonstrates the growth of the shear layer vortices, which is combined with the temporal stability analysis to evaluate the phase velocity of the identified shear layer structures. The calculated phase velocity magnitude of the shear layer is found to be reasonably below the local velocity as expected.     

同轴圆柱间滑移旋转流动问题的同伦解析解

从理论上研究了具有延伸柱面的同轴圆柱间滑移流动问题.通过引入适当的相似变换将控制方程组转化为一类非线性边值问题, 利用同伦分析方法首次获得问题的近似解析解, 分析讨论滑移参数、雷诺数和内外筒半径比对流动的影响.结果表明: 滑移边界参数对剪切应力有较大的影响, 滑移边界小的流体对壁面和边界层内流场施加了更大的剪切力; 增大雷诺数Re, 能增大内筒的壁面剪切力和扭矩; 内外筒半径比对流动结构也有较大的影响, 内筒半径固定, 增大外筒半径能够减少内筒的纵向剪切力.      

Effect of nozzle thickness on the self-excited impinging planar jet

The self-excited oscillation of a large aspect ratio planar jet impinging on a flat plate is investigated experimentally at a single transonic jet velocity to clarify the effect of varying the jet thickness on pattern of jet oscillation and frequency of resulting acoustic tone. The study has been performed for a series of jet thicknesses, 1 mm to 4 mm, each of which is tested for the complete range of plate position, i.e. impingement distance, over which acoustic tones are generated. The results reveal that the jet oscillation is controlled by a fluid-dynamic mechanism for small impingement distances, where the hydrodynamic flow instability controls the jet oscillation without any coupling with local acoustic resonances. At larger impingement distances, a fluid-resonant mechanism becomes dominant, in which one of the various hydrodynamic modes of the jet couples with one of the resonant acoustic modes occurring between the jet nozzle and the impingement plate. Within the fluid-resonant regime, the acoustic tones are found to be controlled by the impingement distance, which is the length scale of the acoustic mode, with the jet thickness having only minor effects on the tone frequency. Flow visualization images of the jet oscillation pattern at a constant impingement distance show that the oscillation occurs at the same hydrodynamic mode of the jet despite a four-fold increase in its thickness. Finally, a feedback model has been developed to predict the frequency of acoustic tones, and has been found to yield reasonable predictions over the tested range of impingement distance and nozzle thickness.     

Dynamic fluid–structure interaction of an elastic capsule in a viscous shear flow at moderate Reynolds number

The unsteady behavior of a 2-D circular elastic capsule was investigated in three viscous shear flows. An immersed boundary method (IBM) has been used to solve the dynamic fluid-structure interaction of the capsule. Computations were carried out in finite parameter ranges where the Reynolds number is Re=1-40 and the capillary number is Ca=0.0005-0.05, which is the ratio of the external viscous shear stress to the resistant elastic tensions of the membrane. For the simple shear flow, the effect of inertia on the transient behavior of the capsule was studied. For the pulsatile shear flow, two values of the peak fluid strain, T_f=1 and 5, were considered for the quasi-steady capsule mechanics. The capsule shows a cyclic structural response that includes subharmonics as the Reynolds number is elevated to 10 and 40. The capsule dynamic response includes a phase lag, which is a function of the capillary number, the Reynolds number, and the peak fluid strain. Finally, the capsule flowing in the Couette flow shows lateral migration due to the transient lift force, which is higher for lower Ca and higher Re. When capsules with diverse elasticity are dispersed along the velocity gradient, the capsule with a hard membrane experienced greater lift than the one with a soft membrane.     

The interaction between a plane shear layer and a slender body

The interaction of a plane shear layer with a thin flat plate located in the nonlinear region of the shear layer has been investigated experimentally. The shear layer's velocity ratio is 0.375, and its Reynolds number is ΔUθ₀/μ = 625. It is found that for different angles of attack of the plate, the mixing layer is deflected toward the slower stream. In addition, the plate attenuates turbulent fluctuations of the shear layer structures.