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.     

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.     

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.     

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.     

Instability of a planar liquid sheet with surrounding fluids between two parallel walls

We investigate the linear and nonlinear instability of a planar liquid sheet with surrounding fluids between two parallel plane solid walls. Linear analysis shows that the maximum temporal growth rate and unstable wave number region of disturbances increase for the dilational and sinuous modes when the gap between the sheet and the wall decreases. The walls have more influence on the instability when the density ratio of the surrounding fluid to the sheet and/or the Weber number decrease. On the other hand, nonlinear analysis is performed by means of the discrete vortex method, where double vortex rows and their mirror images are placed so as to satisfy the boundary condition on the walls. Numerical results show that the walls enhance nonlinearity, which causes deformation and distortion of the sheet, whereas the nonlinearity diminishes linear growth rates except for long dilational disturbances. In particular, as the walls are placed more closely to the sheet, local sheet thinning becomes more pronounced in the long dilational mode, while the dilational mode is more strongly induced from the sinuous mode through monotonic or periodic energy exchanges between the two modes.     

Liquid jet in a high Mach number air stream

The instability of circular liquid jet immersed in a coflowing high velocity air stream is studied assuming that the flow of the viscous gas and liquid is irrotational. The basic velocity profiles are uniform and different. The instabilities are driven by Kelvin–Helmholtz instability due to a velocity difference and neckdown due to capillary instability. Capillary instabilities dominate for large Weber numbers. Kelvin–Helmholtz instability dominates for small Weber numbers. The wavelength for the most unstable wave decreases strongly with the Mach number and attains a very small minimum when the Mach number is somewhat larger than one. The peak growth rates are attained for axisymmetric disturbances (n = 0) when the viscosity of the liquid is not too large. The peak growth rates for the first asymmetric mode (n = 1) and the associated wavelength are very close to the n = 0 mode; the peak growth rate for n = 1 modes exceeds n = 0 when the viscosity of the liquid jet is large. The effects of viscosity on the irrotational instabilities are very strong. The analysis predicts that breakup fragments of liquids in high speed air streams may be exceedingly small, especially in the transonic range of Mach numbers.     

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

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

Temporal analysis of breakup for a power law liquid jet in a swirling gas

The breakup mechanism and instability of a power law liquid jet are investigated in this study. The power law model is used to account for the non-Newtonian behavior of the liquid fluid. A new theoretical model is established to explain the breakup of a power law liquid jet with axisymmetric and asymmetric disturbances, which moves in a swirling gas. The corresponding dispersion relation is derived by a linear approximation, and it is applicable for both shear-thinning and shear-thickening liquid jets. Analysis results are calculated based on the temporal mode. The analysis includes the effects of the generalized Reynolds number, the Weber number, the power law exponent, and the air swirl strength on the breakup of the jet. Results show that the shear-thickening liquid jet is more unstable than its Newtonian and shear-thinning counterparts when the effect of the air swirl is taken into account. The axisymmetric mode can be the dominant mode on the power law jet breakup when the air swirl strength is strong enough, while the non-axisymmetric mode is the domination on the instability of the power liquid jet with a high We and a low Re _n . It is also found that the air swirl is a stabilizing factor on the breakup of the power law liquid jet. Furthermore, the instability characteristics are different for different power law exponents. The amplitude of the power law liquid jet surface on the temporal mode is also discussed under different air swirl strengths.     

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.     

Acoustically excited air-assisted liquid sheets

The effect of acoustic excitation on the disintegration characteristics of air-assisted liquid sheets, which utilize water at ambient temperature, and for velocities up to 1.8 m/s, is investigated. The study using high-speed imaging techniques revealed that optimum frequency modulation of the perturbation generator has a pronounced influence on the associated surface waves and the subsequent breakup of the liquid sheet. The analysis includes characterization of critical wave amplitude, breakup length, and breakup frequency, for Weber numbers in the range 0.30<We_abs<0.44, which are compared with flow features in the absence of acoustic excitation. The results show that acoustic perturbation can effectively suppress the dominance of gravitational and surface tension effects. As a consequence, for low Weber number flows, the interfacial waves exhibit regularity, and thus a better control of primary breakup processes of liquid sheet may be accomplished.     

Experimental investigation on splashing and nonlinear fingerlike instability of large water drops

The fluid physics of the splashing and spreading of a large-scale water drop is experimentally observed and investigated. New phenomena of drop impact that differ from the conventional Rayleigh–Taylor instability theory are reported. Our experimental data shows good agreement with previous work at low Weber number but the number of fingers or instabilities begins to deviate from the R–T equation of Allen at high Weber numbers. Also observed were multiple waves (or rings) on the spreading liquid surface induced from pressure bouncing (or pulsation) within the impacting liquid. The first ring is transformed into a radially ejecting spray whose initial speed is accelerated to a velocity of 4–5 times that of the impacting drop. This first ring is said to be “splashing,” and its structure is somewhat chaotic and turbulent, similar to a columnar liquid jet surrounded by neighboring gas jets at relatively high impact speed. At lower impact speeds, splashing occurs as a crown-shaped cylindrical sheet. A second spreading ring is observed that transforms into fingers in the circumferential direction during spreading. At higher Weber number, the spreading of a third ring follows that of the second. This third ring, induced by the pressure pulsation, overruns and has fewer fingers than the second, which is still in a transitional spreading stage. Several important relationships between the drop impact speed, the spray ejection speed of the first ring, and the number of fingers of the second and third rings are presented, based on data acquired during a set of drop impact experiments. Issues related to the traditional use of the R–T instability are also addressed.     

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.     

Numerical study of flow and heat transfer during a high-speed micro-drop impact on thin liquid films

Numerical simulation of high-speed micro-droplet impingement on thin liquid film covering a heated solid surface has been carried out. Effect of droplet Weber number and liquid film thickness on the characteristics of flow and heat transfer has been investigated using the coupled level set and volume of fluid method. The code is validated against both the experimental and numerical results from the literature. Results show that the crown dynamics is mostly affected by variations in the initial film thickness but is weakly influenced by changes in the Weber number. The liquid within the film can be categorized as three regions based on the heat transfer distribution: the static film region, the transition region, and the impact region. The transient local wall temperature shows three stages: first stage when the temperature decreases rapidly, followed by a second stage in which the temperature starts to rise and then becomes almost constant in the third stage. After drop impact, the local Nusselt number continuously increases until reaching a maximum value, and then decreases approaching the initial impact stage. Our analysis of the change in Weber number shows that larger Weber number contributes to intense temperature variation at the crater core relative to other radial locations. Lastly, the results reveal that the thinner liquid film leads to lower wall temperature and hence, higher average Nusselt number.     

Numerical studies of the instability and breakdown of a boundary-layer low-speed streak

The experimental configuration in [M. Asai, M. Minagawa, M. Nishioka, The instability and breakdown of a near-wall low-speed streak, J. Fluid Mech. 455 (2002) 289–314] is numerically reproduced in order to examine the instability of a single low-speed streak in a laminar boundary layer and to investigate the resulting generation of coherent structures. Such a configuration is chosen since the experimental data show that the two instability modes, varicose and sinuous, are of comparable strength. The instability characteristics are retrieved from the simulation of the flow impulse response. The varicose instability is associated to higher frequencies and lower group velocities than those of the sinuous modes. The latter are less affected by the diffusion of the streak mean shear and are amplified for a longer streamwise distance. Analysis of the perturbation kinetic energy production reveals that both the varicose and the sinuous instability are driven by the work of the Reynolds stress against the wall-normal shear of the streak. The base flow considered here therefore presents an exception to the common knowledge, supported by several previous studies, that the sinuous instability is associated to the streak spanwise shear. The vortical structures at the late stage of the varicose breakdown are identified from the numerical data. By comparing them with those pertaining to other transition scenarios, it is confirmed that streaks and streamwise vortices are universal features of boundary layer transition.     

Linear Stability Analysis of a Viscous Liquid Sheet in a High-Speed Viscous Gas

Air-assisted atomizers in which a thin liquid sheet is deformed under the action of a high-speed air flow are extensively used in industrial applications, e.g., in aircraft turbojet injectors. Primary atomization in these devices is a consequence of the onset and growth of instabilities on the air/liquid interfaces. To better understand this process, a temporal linear instability analysis is applied to a thin planar liquid sheet flowing between two semi-infinite streams of a high-speed viscous gas. This study includes the full viscous effects both in the liquid and gas basic states and perturbations. The relevant dimensionless groups entering the non-dimensional Orr–Sommerfeld equations and boundary conditions are the liquid and gas stream Reynolds numbers, the gas to liquid momentum flux ratio, the gas/liquid velocity ratio, the Weber number and the equivalent gas boundary layer to liquid sheet thickness ratio. Growth rates and temporal frequencies as a function of the wave number, varying the different dimensionless parameters are presented, together with neutral stability curves. From the results of this parametric study it is concluded that when the physical properties of gas and liquid are fixed, the momentum flux ratio is especially relevant to determine the instability conditions. It is also observed that the gas boundary layer thickness strongly influences the wave propagation, and acts by damping sheet oscillation frequency and growth. This is especially important because viscosity in the basic gas velocity profile has always been ignored in instability analysis applied to the geometry under study. This revised version was published online in July 2006 with corrections to the Cover Date.     

Full-field spreading velocity measurement inside droplets impinging on a dry solid surface

Liquid droplet impacts onto solid surfaces have attracted enormous amount of attention from wide range of research fields including experimental and numerical investigations. Unlike experimental efforts, numerical and analytical studies generated various sets of data. In this study, we investigated the spreading velocities inside the water droplets impinging onto a dry glass substrate using time-resolved PIV. The method, together with the high spatiotemporal resolution and the additional treatments improving the robustness, allowed us to resolve the radial velocity profiles efficiently in the spreading phase. Several impact velocity cases ranging from 0.40 to 0.96 m/s were studied. They correspond to low and moderate level Weber numbers (4.9–27.6). We observed that instantaneous radial velocity distributions exhibit linear and nonlinear modes. The nonlinearity is caused by the vortical flows formed at outer regions of the spreading liquid lamella. We demonstrated that even at low impact velocities the linear parts of the profiles obey a quasi-one-dimensional theory proposed in the literature. The comparison of obtained results with a literature-based numerical study, performed for high range of Weber numbers, confirmed the simultaneous existence of linear and nonlinear parts in the radial velocity profiles. In spite of the scale differences in terms of Weber number, the agreements in the tendencies of the profiles imply the validity of the mechanism considered in the numerical study even at low and moderate level range of Weber numbers.     

POD analysis of the instability mode of a low-speed streak in a laminar boundary layer

The instability of one single low-speed streak in a zero-pressure-gradient laminar boundary layer is investigated experimentally via both hydrogen bubble visualization and planar particle image velocimetry (PIV) measurement. A single low-speed streak is generated and destabilized by the wake of an interference wire positioned normal to the wall and in the upstream. The downstream development of the streak includes secondary instability and self-reproduction process, which leads to the generation of two additional streaks appearing on either side of the primary one. A proper orthogonal decomposition (POD) analysis of PIV measured velocity field is used to identify the components of the streak instability in the POD mode space: for a sinuous/varicose type of POD mode, its basis functions present anti-symmetric/symmetric distributions about the streak centerline in the streamwise component, and the symmetry condition reverses in the spanwise component. It is further shown that sinuous mode dominates the turbulent kinematic energy (TKE) through the whole streak evolution process, the TKE content first increases along the streamwise direction to a saturation value and then decays slowly. In contrast, varicose mode exhibits a sustained growth of the TKE content, suggesting an increasing competition of varicose instability against sinuous instability.     

Experiments on the instability modes of buoyant diffusion flames and effects of ambient atmosphere on the instabilities

Large scale dynamic behavior of buoyant diffusion flames were studied experimentally. It was found that buoyant diffusion flames originating from circular nozzles exhibit two different modes of flame instabilities. The first mode results in a sinuous meandering of the diffusion flame, characteristic of flames originating from small diameter nozzles. This instability originates at some distance downstream of the nozzle exit in the contraction region of the buoyant flame envelope and develops into a sinuous motion of the flame. The second mode is the varicose mode which develops very close to the nozzle exit as axisymmetric perturbations of a contracting flame surface. In this mode, flame oscillations result in the formation of toroidal vortical structures that convect through the flame and cause periodic burn out at the flame top resulting in the observed flame height fluctuations. The average flame heights are found to be typically shorter for these flames. The oscillation frequencies and their scaling for the two modes are also different with the sinuous mode having higher frequencies than the varicose mode. It was also observed that the instability can switch from one mode to the other and the probability of observing the varicose mode appears to increase with increasing Richardson number. Additionally, the feasibility of altering the behavior of buoyant diffusion flames was explored through variation of the oxidizer medium density. It was found that the flame oscillations can be completely suppressed for flames burning in helium rich helium–oxygen mixtures. At lower helium concentrations, the oscillation frequency can be significantly reduced. In order to enhance the buoyancy effect, CO₂–O₂ mixtures were also studied. However, the density increase and its effects on flame oscillation frequency were found to be small compared to those flames burning in air. These experiments point towards the feasibility of altering buoyant flame behavior under earth gravity and studying the large scale dynamical aspects of buoyant flames without the need of variable gravity environment. Received: 2 March 1999/Accepted: 6 August 1999     

A study on the aerodynamic instability of attenuating liquid sheets

We describe the characteristics of a radially spreading unstable liquid sheet in quiescent air via optical measurement techniques and linear instability theory. A high speed CCD camera system and a complimentary laser refraction method were employed to measure the intact sheet diameter, unstable wave lengths, wave speed, wave frequency spectrum and spatial wave growth rates. Linear instability models for thinning, viscous and inviscid liquid sheets, which are available from the literature, allow for a comparison of experimental data and predicted sheet behaviour. The last section evaluates the differences and similarities between the current liquid sheet experiment and industrial spray applications such as fuel atomisation via pressure-swirl nozzles.     

电场作用下导电射流稳定性的理论与实验研究

建立导电射流在径向电场作用下的线性稳定性粘性模型,通过正则模方法,推导了轴对称和非轴对称模态下的色散关系,通过计算求得增长率随波数及电欧拉数的变化,并在理论上预测了最有可能波长.选用酒精和酒精甘油混合物作实验液体,观察了径向电场对射流不稳定性行为的影响规律,并测量射流表面波的波长.实验结果和理论结果在定性方面取得了较好的一致.但通过与实验比较,理论预测的最有可能波长在非轴对称模态出现较大偏差,普遍比实验结果小.而且,实验表明,最大增长率并不是判断主导模态的好标准,因为在非轴对称的最大增长率小于轴对称的最大增长率情况下,实验显示非轴对称模态要比轴对称模态明显了.因此,对于非轴对称的不稳定机理,需要进一步研究.对轴对称模态,理论给出了较好的预测.