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

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

Wave-flow theory for a thin viscous liquid layer

On the basis of a simplified system of equations we study the process of development and stability of wave flows in a thin layer of a viscous liquid. Any unstable disturbance of the laminar flow grows and leads to the establishment of the wave regime. The time to establish the flow changes little for large flow rates, but increases sharply with reduction of the flow rate. Given the same amplitudes of the initial disturbances, the optimum regimes which provide the greatest flow rate in a layer of given average thickness develop more rapidly than the other regimes. All the wave regimes are unstable to disturbances in the form of traveling waves. With moderate flow rates, the optimum regimes will be most stable to near-by disturbances.Strictly periodic wave flows in a thin layer of a viscous liquid under the influence of the gravity force were calculated in [1], Various flow wave regimes which differ in wavelength can theoretically be established for a given liquid flow rate. In particular, there is a wavelength for which the flowing layer exhibits minimum average thickness (and maximum flow rate for a given average thickness). These optimum regimes correspond closely to the experimental data [2]; however, with regard to calculation technique these regimes are no different from the others. In the following we make a comparison of the wave regimes on the basis of the nature of their development and stability.     

Transient stress and interfacial area generated during breakup of a Newtonian thread immersed in a Newtonian medium

The transient stress and the transient average orientation generated by the breakup process of a long Newtonian filament imbedded in a quiescent Newtonian viscous liquid have been calculated. Rayleigh disturbances were used to describe the relaxation of the filament and the variation of interfacial area in the absence of flow during the course of disintegration process. The effect of viscosity ratio and initial radius of the filament were discussed. It was demonstrated that the predictions of the model in terms of the time-evolution of interfacial area can be used to select the best conditions for carrying out the breaking thread experiments. The predictions of the proposed model were compared to some experimental data on polyamide/polyethylene system.     

Interaction Between a Time-Dependent Laminar Boundary Layer and an Inviscid Subsonic Flow Near a Region of Local Energy Supply

The results of numerical modeling of the time-dependent flows of a viscous heat-conducting gas occurring in the region of interaction between an external inviscid flow and a laminar boundary layer near a zone of local energy supply at high subcritical Reynolds numbers are presented. The solution of the Navier-Stokes equations is constructed on the basis of the method of matched asymptotic expansions. Numerical solutions of the nonlinear boundary-value problem describing the flow in the wall region of the boundary layer are given in similarity variables. It is shown that time- and space-localized energy supply results in the formation of a self-consistent flow disturbance, whose downstream propagation is accompanied by a disturbance amplitude growth during a short time interval, even after the energy supply has stopped. Calculations of the flows induced by two heat sources placed in tandem make it possible to conclude that the time lag for the second energy supply zone and the distance between the sources can be so chosen that superposition of the disturbances induced by the first and second sources leads, due to nonlinear effects, to a considerable increase in the amplitude of the total flow disturbance.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, 2005, pp. 64–75.Original Russian Text Copyright © 2005 by Kazakov.     

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.     

高Bond数下黏性液滴的Rayleigh-Taylor不稳定性

给出了高Bond数下黏性液滴表面Rayleigh-Taylor线性不稳定性的分析解，这种不稳定性对于超音速气流作用下液滴破碎的早期阶段起着至关重要的作用．基于稳定性分析的结果，导出了用于估算稳定液滴的最大直径及液滴无量纲初始破碎时间的计算式，这些计算式与相关文献给出的实验和分析结果比较显示了良好的一致．     

Development of three-dimensional disturbances in a boundary layer

In the region of transition from a two-dimensional laminar boundary layer to a turbulent one, three-dimensional flow occurs [1–3]. It has been proposed that this flow is formed as the result of nonlinear interaction of two-dimensional and three-dimensional disturbances predicted by linear hydrodynamic stability theory. Using many simplifications, [4, 5] performed a calculation of this interaction for a free boundary layer and a boundary layer on a wall with a very coarse approximation of the velocity profile. The results showed some argreement with experiment. On the other hand, it is known that disturbances of the Tollmin—Schlichting wave type can be observed at sufficiently high amplitude. This present study will use the method of successive linearization to calculate the primary two- and three-dimensional disturbances, and also the average secondary flow occurring because of nonlinear interaction of the primary disturbances. The method of calculation used is close to that of [4, 5], the disturbance parameters being calculated on the basis of a Blazius velocity profile. A detailed comparison of results with experimental data [1] is made. It developed that at large disturbance amplitude the amplitude growth rate differs from that of linear theory, while the spatial distribution of disturbances agree s well with the distribution given by the natural functions and their nonlinear interaction. In calculating the secondary flow an experimental correction was made to the amplitude growth rate.     

The breakup and atomization of a viscous liquid jet

Based on the linear analysis of stability, a dispersion equation is deduced which delineates the evolution of a general 3-dimensional disturbance on the free surface of an incompressible viscous liquid jet. With respect to the spatial growing disturbance mode, the numerical results obtained from the solution of the dispersion equation reveal that a dimensionless parameterJ _e exists. AsJ _e>1, the axisymmetric disturbance mode is most unstable; and whenJ _e<1, the asymmetric disturbances come into being, their growth rate increases with the decrease, ofJ _e, till one of them becomes the most unstable disturbance. The breakup of a low-speed liquid jet results from the developing of axisymmetric disturbances, whose instability is produced by the surface tension; while the atomization of a high-speed liquid jet is brought about by the evolution of nonaxisymmetric disturbance, whose instability is caused by the aerodynamic force on the interface between the jet and the ambient gas. The project supported by the National Natural Science Foundation of China     

Secondary atomization

When a drop is subjected to a surrounding dispersed phase that is moving at an initial relative velocity, aerodynamic forces will cause it to deform and fragment. This is referred to as secondary atomization. In this paper, the abundant literature on secondary atomization experimental methods, breakup morphology, breakup times, fragment size and velocity distributions, and modeling efforts is reviewed and discussed. Focus is placed on experimental and numerical results which clarify the physical processes that lead to breakup. From this, a consistent theory is presented which explains the observed behavior. It is concluded that viscous shear plays little role in the breakup of liquid drops in a gaseous environment. Correlations are given which will be useful to the designer, and a number of areas are highlighted where more work is needed. This material is based upon work supported under a National Science Foundation Graduate Research Fellowship.     

Technique for extensional rheology characterization of highly reactive viscoelastic liquids

This paper reports the elaboration and validation of a new technique for extensional rheological characterization of highly reactive viscous liquids contained in a transparent reaction bulb where considered liquids have been prepared. For measurements, a rod was pulled down to contact with a surface of the liquid and then instantly raised up resulting in formation of a liquid thread. During the test, the studied fluid is in an equilibrium state with an inherent vapor. The elasto-capillary self-thinning and ultimate breakup of threads was registered by a high-speed camera. From the fluid dynamics point of view, the main distinguishing feature of the proposed method from the standard capillary breakup extensional rheometer (CaBER) method is the presence of unlimited amount of a test liquid under the endplate of a moving rod. In the present work, the presumable effect of this distinguishing factor on the final test data was estimated. In the frame of this research, the simulations of Newtonian and non-Newtonian fluids were conducted for the tests performed with the newly elaborated and CaBER-like configurations of devices. The test and computed data show some difference in the results obtained only at low values of Hencky strain (at the initial period of thread evolution). Increase in Hencky strain leads to a good agreement between results obtained with the help of CaBER-like devices as well as the proposed one.     

Flow and breakup characteristics of elliptical liquid jets

This paper presents the results of an experimental study on liquid jets discharging from elliptical orifices into still ambient air. The experiments were conducted with a set of elliptical orifices of approximately same area of cross section but varying orifice aspect ratio using water and water–glycerol mixture as experimental fluids. The flow behavior of liquid jets was analyzed using their photographs captured by an imaging system. The measurements obtained for the elliptical liquid jets were compared with the circular liquid jets discharging from a circular orifice of the same area of cross section. Elliptical geometry of the orifice results in a flow process by which the emanating liquid jet periodically switches its major and minor axes as it flows downstream of the orifice. In this paper, we attempt to characterize the axis-switching process through its wavelength and amplitude. For a given elliptical orifice, the axis-switching process is dominantly seen in a particular range of flow conditions. The effects of the orifice aspect ratio and liquid viscosity on the axis-switching process are revealed through this study. The experimental results on jet breakup show that axis-switching process has a destabilizing effect on elliptical liquid jets within a particular range of flow conditions and it results in shorter breakup lengths compared to the circular jet. The extent to which axis-switching destabilizes the jet is dictated by the viscosity of liquid. An increase in orifice aspect ratio destabilizes elliptical liquid jets with low viscosity like water; however, this behavior seems to get obscured in water–glycerol mixture elliptical jets due to high viscosity.     

Numerical simulation of the nonuniform breakup of capillary jets

The problem of the evolution of the surface of a jet up to the stage at which it breaks up into droplets is solved numerically for two initial wave disturbances. The wave number of one of these coincides with the wave number of the disturbance that grows most strongly according to the linear theory, while the wave number of the other is varied. The effect of the wave numbers and the amplitude ratio of the initial disturbances on the breakup time and the appearance of nonuniformity is investigated.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.2, pp. 12–17, March–April, 1993.     

Influence of periodic disturbances on inception of cavitation in sharp-edged orifices

The influence of periodic flow disturbances on the onset of cavitation in sharp-edged orifices was determined through a series of experiments. Different amplitudes and frequencies of disturbance were introduced in the flow. The threshold value of cavitation parameter for which transition from attached to cavitating flow took place was shown to depend on the amplitude of the disturbance. Cavitating flow was initiated more spontaneously as the amplitude of disturbance increased. Hysteresis in the flow was, however, observed to decrease. Variations in the frequency of the disturbance did not bring about any significant change in the cavitation parameter at transition. The results are interpreted in terms of the available dwell time for cavitation and the characteristic residence time of the fluid in the sharp-edged orifice and in the pipeline upstream of the orifice.     

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.     

Simulation of the primary breakup of a high-viscosity liquid jet by a coaxial annular gas flow

The conversion of low-grade fossil and biogenic energy resources (petcoke, biomass) to a synthesis gas in a high pressure entrained flow gasification process opens a wide spectrum for high efficient energy conversion processes. The synthesis gas can be used for production of methane (SNG), liquid fuels (BtL, CtL) or as fuel for operation of a gas turbine in a combined cycle power plant (IGCC). The production of a tar free high quality syngas is a challenging objective especially due to the fact that typical liquid or suspension fuels for entrained flow gasifiers feature viscosities up to 1000 mPas. Fuel droplet conversion at typical entrained flow gasification conditions is characterized by heat up, evaporation and subsequent degradation of the vapour phase. To guarantee a high fuel conversion rate in the gasifier an efficient atomization of the fuel is required. Mainly twin-fluid burner nozzles are used for atomization of those typically high viscous fuels. The present study is focused on the assessment of the accuracy of CFD computations for the primary breakup of high-viscosity liquids using an external mixing twin fluid nozzle. In a first step experiments were performed with a Newtonian glycerol-water-mixture featuring a liquid viscosity of 400 mPas. Jet breakup was investigated using a high speed camera as well as PIV and LDA-System for a detailed investigation of the flow field. In a second step the experimental results serve as reference data to assess the accuracy of CFD computations. Compressible large eddy simulations (LES) were performed to capture the morphology of the primary breakup as well as the important flow field characteristics. A Volume of Fluid (VOF) approach was used to track the unsteady evolution and breakup of the liquid jet. Comparison of experimental and numerical results showed good agreement with respect to breakup frequency, velocity fields and morphology.     

The effect of surfactants on the instability of a rotating liquid jet

It has long been known that the presence of surfactants on the free surface of a liquid jet can create surface tension gradients along the interface. The resulting formation of tangential stresses along the surface lead to Marangoni type flows and greatly affect the resulting dynamics of rupture. In this way surfactants can be used to manipulate the breakup of a liquid jet and control the size of droplets produced. In this paper we investigate the effects of insoluble surfactants on the breakup of rotating liquid jets with applications to industrial prilling. Using a long wavelength approximation we reduce the governing equations into a set of one-dimensional equations. We use an asymptotic theory to find steady solutions and then carry out a linear instability analysis on these solutions. We show that steady state centreline solutions are independent of viscosity to leading order and that the most unstable wavenumber and growth rate of disturbances decrease as the effectiveness of surfactants is increased. We also numerically solve these equations using a finite difference scheme to investigate the effects of changing the initial surfactant concentration and other fluid parameters. Our results show that differences in breakup lengths between rotating surfactant-laden jets and surfactant-free jets increase with the rate of rotation. Moreover, we find that satellite droplet sizes decrease as the rate of rotation is decreased with the effect of surfactants amplifying the reduction in sizes. Furthermore, the presence of surfactants at fixed rotation rates is shown to produce larger main droplets at low disturbance wavenumbers whilst satellite droplets are smaller for moderate disturbance wavenumbers κ≈0.7.     

Stability of a viscous liquid film flowing down a periodic surface

The paper is devoted to a theoretical analysis of linear stability of the viscous liquid film flowing down a wavy surface. The study is based on the Navier–Stokes equations in their full statement. The developed numerical algorithm allows us to obtain pioneer results in the stability of the film flow down a corrugated surface without asymptotic approximations in a wide range over Reynolds and Kapitsa’s numbers. It is shown that in the case of moderate Reynolds numbers there is a region of the corrugation parameters (amplitude and period) where all disturbances decay in time and the wall corrugation demonstrates a stabilizing effect. At the same time, there exist corrugation parameters at which the steady-state solution is unstable with respect to perturbations of the same period as the period of corrugation. In this case the waveless solution cannot be observed in reality and the wall corrugation demonstrates a destabilizing effect.     

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.     

Stability of a laminar rivulet liquid flow in a cylindrical duct in the approximation of one-dimensional waves

The problem of the stability of a viscous laminar liquid flow with a liquid free surface in an inclined duct is theoretically considered. Since the dependence of the flow rate on the free-surface height is not monotonic (the highest flow rate in a cylindrical duct is observed at H_*=1.7R), primary attention is given to the region H>H_*. It is proved that there is aw region of instability: for an arbitrarity low Reynolds number, there is a free-surface level above which the flow becomes unstable against one-dimensional disturbances. When the height of the liquid layer is close to the vertical dimension of the duct, the one-dimensional disturbances propagate mainly upstream (for moderate Reynolds numbers). Hence it follows that there is not steady regime of liquid flow from a fully filled duct with an open end. Kutateladze Institute of Thermal Physics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 3, pp. 90–96, May–June, 1999.     

Stability Analysis of a Falling Film Flow Down a Plane with Sinusoidal Corrugations

The liquid viscous film falling down a vertical wall with sinusoidal relief is considered. The linear stability of steady-state flow with respect to time-periodic disturbances is studied using the Floquet theory. It is shown that in the case of applying corrugations the variation in the disturbance growth rate is proportional to the second power of their undulations. Depending on the relief parameters there exist two possibilities: the instability domain can expand or certain disturbances can be stabilized. The growth rates are obtained numerically and analytically in the approximation of low-amplitude corrugations. The development of waves from small disturbances is simulated within the framework of nonlinear equations and the formation of structures whose wavelength is significantly greater than the space relief period is found out.