Three-wave interactions of disturbances in a hypersonic boundary layer on a porous surface

Interactions of disturbances in a hypersonic boundary layer on a porous surface are considered within the framework of the weakly nonlinear stability theory. Acoustic and vortex waves in resonant three-wave systems are found to interact in the weak redistribution mode, which leads to weak decay of the acoustic component and weak amplification of the vortex component. Three-dimensional vortex waves are demonstrated to interact more intensively than two-dimensional waves. The feature responsible for attenuation of nonlinearity is the presence of a porous coating on the surface, which absorbs acoustic disturbances and amplifies vortex disturbances at high Mach numbers. Vanishing of the pumping wave, which corresponds to a plane acoustic wave on a solid surface, is found to assist in increasing the length of the regions of linear growth of disturbances and the laminar flow regime. In this case, the low-frequency spectrum of vortex modes can be filled owing to nonlinear processes that occur in vortex triplets.     

圆环旋转粘性液膜射流三维不稳定性研究

利用线性稳定性理论进行了射流液体粘性对圆环旋转液膜射流稳定性影响的研究,推导出了三维扰动下具有固体旋涡型速度分布的圆环旋转粘性液膜射流的色散方程;在此基础上进行了类反对称模式与类对称模式下的圆环旋转粘性液膜射流的三维不稳定性分析。研究结果表明,在类反对称模式下,液体粘性超过一定值后,射流最大扰动增长率随液体粘性的增加而迅速减小;轴对称模态的射流特征频率产生一个突降变化;随液体粘性增加,轴对称模态不稳定波数范围减小,非轴对称模态不稳定波数范围呈现出先减小后增大趋势。在类对称模式下,液体粘性对射流最大扰动增长率的影响主要体现在对非轴对称模态的影响上;液体粘性只在粘性较大时才会对非轴对称模态射流特征频率产生一定影响;液体粘性超过一定值后,轴对称模态与非轴对称模态的不稳定波数范围都会快速下降。     

Formation of radially expanding liquid sheet by impinging two round jets

A thin circular liquid sheet can be formed by impinging two identical round jets against each other. The liquid sheet expands to a certain critical radial distance and breaks. The unsteady process of the formation and breakup of the liquid sheet in the ambient gas is simulated numerically. Both liquid and gas are treated as incompressible Newtonian fluids. The flow considered is axisymmetric. The liquid-gas interface is modeled with a level set function. A finite difference scheme is used to solve the governing Navier-Stokes equations with physical boundary conditions. The numerical results show how a thin circular sheet can be formed and break at its circular edge in slow motion. The sheet continues to thin as it expands radially. Hence, the Weber number decreases radially. The Weber number is defined as ρu 2 h/σ, where ρ and σ are, respectively, the liquid density and the surface tension, and u and h are, respectively, the average velocity and the half sheet thickness at a local radial location in the liquid sheet. The numerical results show that the sheet indeed terminates at a radial location, where the Weber number reaches one as observed in experiments. The spatio-temporal linear theory predicts that the breakup is initiated by the sinuous mode at the critical Weber number We c =1, below which the absolute instability occurs. The other independent mode called the varicose mode grows more slowly than the sinuous mode according to the linear theory. However, our numerical results show that the varicose mode actually overtakes the sinuous mode during the nonlinear evolution, and is responsible for the final breakup. The linear theory predicts the nature of disturbance waves correctly only at the onset of the instability, but cannot predict the exact consequence of the instability.     

Linear stability of two-dimensional steady flow in wavy-walled channels

Linear stability of fully developed two-dimensional periodic steady flows in sinusoidal wavy-walled channels is investigated numerically. Two types of channels are considered: the geometry of wavy walls is identical and the location of the crest of the lower and upper walls coincides (symmetric channel) or the crest of the lower wall corresponds to the furrow of the upper wall (sinuous channel). It is found that the critical Reynolds number is substantially lower than that for plane channel flow and that when the non-dimensionalized wall variation amplitude is smaller than a critical value (about 0.26 for symmetric channel, 0.28 for sinuous channel), critical modes are three-dimensional stationary and for larger , two-dimensional oscillatory instabilities set in. Critical Reynolds numbers of sinuous channel flows are smaller for three-dimensional disturbances and larger for two-dimensional disturbances than those of symmetric channel flows. The disturbance velocity distribution obtained by the linear stability analysis suggests that the three-dimensional stationary instability is mainly caused by local concavity of basic flows near the reattachment point, while the critical two-dimensional mode resembles closely the Tollmien–Schlichting wave for plane Poiseuille flow.     

Absolute and convective instability of a liquid sheet with transverse temperature gradient

The spatial–temporal instability behavior of a viscous liquid sheet with temperature difference between the two surfaces was investigated theoretically. The practical situation motivating this investigation is liquid sheet heated by ambient gas, usually encountered in industrial heat transfer and liquid propellant rocket engines. The existing dispersion relation was used, to explore the spatial–temporal instability of viscous liquid sheets with a nonuniform temperature profile, by setting both the wave number and frequency complex. A parametric study was performed in both sinuous and varicose modes to test the influence of dimensionless numbers on the transition between absolute and convective instability of the flow. For a small value of liquid Weber number, or a great value of gas-to-liquid density ratio, the flow was found to be absolutely unstable. The absolute instability was enhanced by increasing the liquid viscosity. It was found that variation of the Marangoni number hardly influenced the absolute instability of the sinuous mode of oscillations; however it slightly affected the absolute instability in the varicose mode.     

Simple and efficient numerical methods for vortex sheet motion with surface tension

We present two simple and efficient explicit methods for the vortex sheet with surface tension. The first one is the standard point vortex method, which has been known to be unstable in the presence of surface tension, due to spurious growth of waves of high modes. We show, for the first time, that the standard point vortex method is able to calculate the vortex sheet motion with surface tension by employing a Fourier filtering. The second method is a modification of the Pullin method using central differences for numerical differentiations. This method is more convenient to implement than other spectral methods and is free from the aliasing instability. We give a linear stability analysis for the numerical methods and show results for the long‐time evolution of the vortex sheet. We also propose a new redistribution procedure to control point clustering, by setting limits of minimum and maximum distances between neighboring points. This procedure is found to be very efficient for long‐time computations of the explicit methods. Copyright © 2013 John Wiley & Sons, Ltd.     

Three-wave interactions between disturbances in the hypersonic boundary layer on impermeable and porous surfaces

The interaction between disturbances in the hypersonic boundary layer on impermeable and porous surfaces is considered within the framework of weakly-nonlinear stability theory. It is established that on the impermeable surface nonlinear interactions between different waves (acoustic and vortex) occur in the parametric resonance regime. The role of pumping wave is played by a plane acoustic wave. The nonlinear interactions take place over a wide frequency range and can lead to the packet growth of Tollmien-Schlichting waves. On the porous surface the analogous interactions are fairly weak and result in a slight decay of the acoustic mode and a slight amplification of the vortex mode. This leads to the dragging out of the laminar flow regime and the regions of linear disturbance growth. In this situation the low-frequency spectrum of the vortex modes may be filled on account of the nonlinear processes occurring in the three-wave systems between the vortex components.     

圆环旋转黏性液体射流空间不稳定性研究

利用线性稳定性理论, 进行了液体黏性对不同旋转强度下圆环旋转液体射流 空间不稳定性影响的研究. 在推导出的三维扰动下具有固体涡核型旋转速度分布的圆环旋转 黏性液体射流色散方程的基础上, 针对中低速射流, 进行了类反对称模式与类对称模式下圆 环旋转黏性液体射流的空间不稳定性分析. 研究结果表明, 对于旋转强度较大的圆环旋转液 体射流, 液体黏性的增加, 不利于射流的破碎; 随着液体黏性的增加, 射流的特征频率和最 不稳定波数减小. 然而, 对于旋转强度较小的圆环旋转液体射流, 液体黏性的增加, 有利于 射流的破碎; 随着液体黏性的增加, 类反对称模式下射流特征频率先减小后增大, 类对称模 式下射流特征频率增大; 随着液体黏性的增加, 类反对称模式下射流最不稳定波数先减小后 增大, 类对称模式下射流最不稳定波数增大.     

Spatial instability of a viscous liquid sheet

In the present study, the spatial instability for a two‐dimensional viscous liquid sheet, which is thinning with time, has been analysed. The study includes the derivation of a spatial dispersion equation, numerical solutions for the growth rate of sinuous disturbances, and parameter sensitivity studies. For a given wave number, the growth rate of the disturbance is essentially a function of Weber number, Reynolds number, and gas/liquid density ratio. The analysis indicates that the cut‐off wave number of the disturbance becomes larger with an increase in Weber number or gas/liquid density ratio. Thus, the liquid sheet should produce finer drops. When the Reynolds number decreases, the higher viscosity has a greater damping effect on shorter waves than longer waves. This could explain that only large drops and ligaments were observed in past measurements for the disintegration of a very viscous sheet. The spatial instability results of the present study were also compared with the temporal theory. The importance of spatial analysis was found and demonstrated for the cases of low Weber numbers. The temporal theory underestimates growth rates when the Weber number is less than 100. The discrepancy between the two theories increases as the Weber number further decreases. Copyright © 2005 John Wiley & Sons, Ltd.     

Instability of a viscoelastic liquid jet with axisymmetric and asymmetric disturbances

The temporal instability behavior of a viscoelastic liquid jet in the wind-induced regime with axisymmetric and asymmetric disturbances moving in an inviscid gaseous environment is investigated theoretically. The corresponding dispersion relation between the wave growth rate and the wavenumber is derived. The linear instability analysis shows that viscoelastic liquid jets are more unstable than their Newtonian counterparts, and less unstable than their inviscid counterparts, for both axisymmetric and asymmetric disturbances, respectively. The instability behavior of viscoelastic jets is influenced by the interaction of liquid viscosity and elasticity, in which the viscosity tends to dampen the instability, whereas the elasticity results in an enhancement of instability. Relatively, the effect of the ratio of deformation retardation to stress relaxation time on the instability of viscoelastic jets is weak. It is found that the liquid Weber number is a key measure that controls the viscoelastic jet instability behavior. At small Weber number, the axisymmetric disturbance dominates the instability of viscoelastic jets, i.e., the growth rate of an axisymmetric disturbance exceeds that of asymmetric disturbances. When the Weber number increases, both the growth rate and the instability range of disturbances increase drastically. The asymptotic analysis shows that at large Weber number, more asymmetric disturbance modes become unstable, and the growth rate of each asymmetric disturbance mode approaches that of the axisymmetric disturbance. Therefore, the asymmetric disturbances are more dangerous than that of axisymmetric disturbances for a viscoelastic jet at large Weber numbers. Similar to the liquid Weber number, the ratio of gas to liquid density is another key measure that affects the viscoelastic jet instability behavior substantially.     

Nonlinear instability of an annular liquid sheet exposed to gas flow

Nonlinear instability and breakup of an annular liquid sheet has been modeled in this paper. The liquid sheet is considered to move axially and is exposed to co-flowing inner and outer gas streams. Also, the effect of outer gas swirl on sheet breakup has been studied. In the developed model a perturbation expansion method has been used with the initial magnitude of the disturbance as the perturbation parameter. This is a comprehensive model in that other geometries of planar sheet and a coaxial jet can be obtained as limiting cases of very large inner radius and inner radius equal to zero, respectively. In this temporal analysis, the effect of liquid Weber number, initial disturbance amplitude, inner gas-to-liquid velocity ratio, outer gas-to-liquid velocity ratio and outer gas swirl strength on the breakup time is investigated. The model is validated by comparison with earlier analytical studies for the limiting case of a planar sheet as well as with experimental data of sheet breakup length available in literature. It is shown that the linear theory cannot predict breakup of an annular sheet and the developed nonlinear model is necessary to accurately determine the breakup length. In the limiting case of a coaxial jet, results show that gas swirl destabilizes the jet, makes helical modes dominant compared to the axisymmetric mode and decreases jet breakup length. These results contradict earlier linear analyses and agree with experimental observations. For an annular sheet, it is found that gas flow hastens the sheet breakup process and shorter breakup lengths are obtained by increasing the inner and the outer gas velocity. Axially moving inner gas stream is more effective in disintegrating the annular sheet compared to axially moving outer gas stream. When both gas streams are moving axially, the liquid sheet breakup is quicker compared to that with any one gas stream. In the absence of outer gas swirl, the axisymmetric mode is the dominant instability mode. However, when outer gas flow has a swirl component higher helical modes become dominant. With increasing outer gas swirl strength, the maximum disturbance growth rate increases and the most unstable circumferential wave number increases resulting in a highly asymmetric sheet breakup with shorter breakup lengths and thinner ligaments.     

Experimental and analytical investigation of liquid sheet breakup characteristics

The main objective of this research is to study analytically and experimentally the liquid sheet breakup of a flat fan jet nozzle resulting from pressure-swirling. In this study the effects of nozzle shape and spray pressure on the liquid sheet characteristics were investigated for four nozzles with different exit widths (1.0, 1.5, 2.0 and 2.5 mm). The length of liquid sheet breakup, liquid sheet velocity and the size of formed droplets were measured by a digital high speed camera. The breakup characteristics of plane liquid sheets in atmosphere are analytically investigated by means of linear and nonlinear hydrodynamic instability analyses. The liquid sheet breakup process was studied for initial sinuous and also varicose modes of disturbance. The results presented the effect of the nozzle width and the spray pressure on the breakup length and also on the size of the formed droplets. Comparing the experimental results with the theoretical ones for all the four types of nozzles, gives a good agreement with difference ranges from 4% to 12%. Also, the comparison between the obtained results and the results due to others shows a good agreement with difference ranged from 5% to 16%. Empirical correlations have been deduced describing the relation between the liquid sheet breakup characteristics and affecting parameters; liquid sheet Reynolds number, Weber number and the nozzle width.     

On the inviscid acoustic-mode instability of supersonic shear flows

The same methods used previously to study acoustic-mode instability in supersonic boundary layers are applied to free shear layers, and new calculations are made for boundary layers with cooling and suction. The objective is to obtain additional information about acoustic-mode instability, and to find what features of the instability are common to boundary layers and free shear flows. Acoustic modes exist whenever there is an embedded region of locally supersonic flow relative to the phase speed of the instability wave. Consequently, they can be found in boundary layers, wakes, and jets, but not in mixing layers unless the flow is confined. In this first part of a two-part paper, attention is directed principally to two-dimensional waves. The linear, inviscid stability theory is used to calculate spatial amplification rates at Mach number 3 for the sinuous and varicose modes of a single wake flow and a single jet flow, each made up of the same mixing-layer profile plus a central region of uniform flow. Along with sequences of sinuous and varicose unstable modes clearly identifiable as acoustic modes, both of these flows, unlike the boundary layer, have a lowest sinuous mode that is the most unstable. The unstable modes include both subsonic and radiating disturbances with large amplification rates. The latter phenomenon is also found for highly cooled boundary layers with suction. In these boundary layers, suction is generally stabilizing for nonradiating acoustic disturbances, but destabilizing for radiating disturbances.The work described in this paper was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration (NASA). Support from the Aerodynamics Division of the Office of Aeronautics and Exploration Technology is gratefully acknowledged. A preliminary version of this paper was presented at the Fourth Symposium on Numerical and Physical Aspects of Aerodynamic Flows, California State University, Long Beach, CA, 16–19 January 1989.     

Electrohydrodynamic instability of a dielectric compressible liquid sheet streaming into an ambient stationary compressible gas

The effect of compressibility of fluids on the linear electrohydrodynamic instability of a dielectric liquid sheet issued from a nozzle into an ambient dielectric stationary gas in the presence of a horizontal electric field is investigated. It is found that increasing the Mach number from subsonic to transonic causes the maximum growth rate and the dominant wavenumber of the disturbances to increase, and the increase is higher in the presence of the electric field. Liquid compressibility has been found to have a minimal effect on instability. At constant wavenumber and electric field values, the growth rate of disturbances increases as the gas Mach number tends to 1, and then begins to decrease with further increase in the gas Mach number. At small values of wavenumber, antisymmetrical disturbances grow faster than symmetrical ones, while the growth rate of both types of disturbances approach each other at large wavenumbers, which increases by increasing the electric field values. At small Weber numbers, antisymmetrical disturbances exhibit a higher maximum growth rate and a lower dominant wavenumber than symmetrical disturbances. However, the maximum growth rate and dominant wavenumber of the two types of disturbances are almost identical when both Weber number and electric field values become large. An increase in the gas to liquid density ratio enhances the instability, and this effect is enhanced for high electric field values. Surface tension and electric fields always oppose and increase the development of instability, respectively; and they have opposite effects for long wavelengths and high Weber numbers.     

Stability and capillary dynamics of circular vortex sheets

We investigate the motion of circular vortex sheets with surface tension. A linear stability analysis shows that high modes of the circular vortex sheet are stabilized by surface tension, and the sheet is stable if surface tension is larger than a critical value. The modes of perturbations, n = 1 and 2, are always stable, regardless of surface tension, and the mode n = 3 is also stable for large surface tension. The numerical results show that a stable vortex sheet rotates and oscillates weakly. The oscillations of each mode of the interface mainly consist of two travelling waves of different frequencies in time. The amplitude and the period of the oscillation depend on the mode of the perturbation and surface tension. We also perform long-time computations for the unstable evolution of circular sheets. For a high Weber number, ripples are produced on the sheets, as well as pinching and self-intersection. It is found that the appearance of ripples is associated with the growth of noise. For an intermediate Weber number, the sheet evolves to an exotic structure with small spikes on the fingers, while for a low Weber number, it is nonlinearly stable.     

Viscous and Inviscid Instabilities of Flow Along a Streamwise Corner

Here we consider the stability of flow along a streamwise corner formed by the intersection of two large flat plates held perpendicular to each other. Self-similar solutions for the steady laminar mean flow in the corner region have been obtained by solving the boundary layer equations for zero and nonzero streamwise pressure gradients. The stability of the mean flow is investigated using linear stability analysis. An eigensolver has been developed to solve the resulting linear eigenvalue problem either in a global mode to obtain an approximation to all the dominant eigenmodes or in a local mode to refine a particular eigenmode. The stability results indicate that the entire spectrum of two-dimensional and oblique viscous modes of a two-dimensional Blasius boundary layer is active in the case of a corner layer as well, but away from the cornerline. In a corner region of finite spanwise extent, the continuous spectrum of oblique modes degenerates to a discrete spectrum of modes of increasing spanwise wave number. The effect of the corner on the two-dimensional viscous instability is small and decreases the growth rate. The growth rate of outgoing oblique disturbances is observed to decrease, while the growth rate of incoming oblique disturbances is enhanced by the corner. This asymmetry between the outgoing and incoming viscous modes increases with increasing obliqueness of the disturbance. The instability of a zero pressure gradient corner layer is dominated by the viscous modes; however, an inviscid corner mode is also observed. The critical Reynolds number of the inviscid mode rapidly decreases with even a small adverse streamwise pressure gradient and the inviscid mode becomes the dominant one. Received 17 March 1998 and accepted 28 April 1999     

Three-dimensional instability analysis of boundary layers perturbed by streamwise vortices

A parametric study is presented for the incompressible, zero-pressure-gradient flat-plate boundary layer perturbed by streamwise vortices. The vortices are placed near the leading edge and model the vortices induced by miniature vortex generators (MVGs), which consist in a spanwise-periodic array of small winglet pairs. The introduction of MVGs has been experimentally proved to be a successful passive flow control strategy for delaying laminar-turbulent transition caused by Tollmien–Schlichting (TS) waves. The counter-rotating vortex pairs induce non-modal, transient growth that leads to a streaky boundary layer flow. The initial intensity of the vortices and their wall-normal distances to the plate wall are varied with the aim of finding the most effective location for streak generation and the effect on the instability characteristics of the perturbed flow. The study includes the solution of the three-dimensional, stationary, streaky boundary layer flows by using the boundary region equations, and the three-dimensional instability analysis of the resulting basic flows by using the plane-marching parabolized stability equations. Depending on the initial circulation and positioning of the vortices, planar TS waves are stabilized by the presence of the streaks, resulting in a reduction in the region of instability and shrink of the neutral stability curve. For a fixed maximum streak amplitude below the threshold for secondary instability (SI), the most effective wall-normal distance for the formation of the streaks is found to also offer the most stabilization of TS waves. By setting a maximum streak amplitude above the threshold for SI, sinuous shear layer modes become unstable, as well as another instability mode that is amplified in a narrow region near the vortex inlet position.     

Controlling supersonic boundary layer stability by means of distributed mass transfer through a porous wall

The interaction between disturbances in a compressible boundary layer in the presence of distributed mass transfer (injection or suction) through a permeable porous wall is considered in the linear and nonlinear approximations (weakly nonlinear stability theory). The regimes of moderate and high supersonic velocities (Mach numbers M = 2 and 5.35) are studied. The boundary conditions for the disturbances on a permeable wall are derived with account for the gas compressibility in pores and the presence of a suction chamber. Maximum pore dimensions, at which the surface properties have no effect on the disturbance characteristics, which are stabilized upon suction and destabilized upon injection, are determined. When the surface properties are taken into account, intense growth of the first-mode vortex disturbances occurs, which can completely undo the stabilizing effect of the suction. Injection leads to the vortex and acoustic mode destabilization on the linear range and the enhancement of the nonlinear processes on the transitional range.     

Experimental characterization of the instability of the vortex rings. Part II: Non-linear phase

The first stage of the instability of a vortex ring is linear and characterized by the growth of an azimuthal stationary wave which develops around the ring. Theoretical works predict its origin, shape, number of waves and growth rate. Apart for the growth rate, experimental and numerical results in viscous fluids fit well with the predictions based on an ideal fluid hypothesis. On the other hand, the next stages of the development of the instability (which are non-linear) are not well known. Only few phenomena are described, in an isolated way, in various partial contributions. The aim of this paper is to report on a complete experimental investigation of the non-linear phase of the instability of the vortex ring. The vortices were produced in water and their Reynolds number Re _p was varied from 2,650 to 6,100. Visualizations were performed using planar laser induced fluorescence and measurements with 2D2C and 2D3C particle image velocimetry. Based on a Fourier analysis of the results, it appears that the non-linear phase begins with the development of harmonics of the linear modes (first unstable modes). But the growth of those harmonics is rapidly stopped by the development of low order modes. Then appears an m=0 mode, which corresponds to a mean azimuthal velocity around the vortex. Simultaneously, secondary vortical structures develop all around the vortex in its peripheral zone. These vortical structures are linked with the ejection of vorticity in the wake of the ring and they appear just before the transition towards turbulence. A tentative is made here to place all these phenomena chronologically, in order to propose a scenario for the transition from the linear phase to turbulence.     

Thin‐tube vortex simulations for sinusoidal instability in a counter‐rotating vortex pair

A thin‐tube vortex method is developed to investigate the intrinsic instability within a counter‐rotating vortex pair system and the effects from the core size and the wavenumbers (or wavelengths). The numerical accuracy and the advantages of the scheme are theoretically estimated. A nearest‐neighbour‐image method is employed in this three‐dimensional vortex simulation. Agreement with Crow's instability analysis has been achieved numerically for the long‐wave cases. A short‐wave instability for the zeroth radial mode of bending instability has also been found using the thin‐tube vortex simulations. Then, the combinations of long‐ and short‐wave instability are investigated to elucidate the non‐linear effects due to the interactions of two different modes. It is shown that instability is enhanced if both long‐ and short‐wave instabilities occur simultaneously. Although the method used in the paper is not capable of including effects such as axial flow, vortex core deformation and other complicated viscous effects, it effectively predicts and clarifies the first‐order factor that dominates the sinusoidal instability behaviour in a vortex pair. Copyright © 2002 John Wiley & Sons, Ltd.