Comparative investigations of surface instabilities (“sharkskin”) of a linear and a long-chain branched polyethylene

An experimental study of the physical origin and the mechanisms of the sharkskin instability is presented. Extrusion flows through a slit die are studied for two materials: a linear low density polyethylene (LLDPE) which exhibits sharkskin instability for flow rates larger than an onset value and a low density polyethylene (LDPE) which does not show any instability over a broad range of flow rates. By combining laser-Doppler velocimetry (LDV) with rheological measurements in both uniaxial extension and shear, the distributions of tensile and shear stresses in extrusion flows are measured for both materials. The experimentally measured flow fields appear to be qualitatively similar for both the unstable (LLDPE) and stable case (LDPE): around the die exit the flow accelerates near the boundaries and decelerates around the flow axis. The fields of the axial gradients of the axial velocity component are, however, quite different in the two cases. In the unstable case there exists a strongly non-uniform transversal distribution of velocity gradients near the die exit. This non-uniformity of the distribution of gradients is significantly smaller in the stable case. Significant differences in the extensional rheological properties of the two materials are found as well. Due to its branched structure, the LDPE is able to sustain higher tensile stresses prior to failure. Measurements of the distributions of tensile stresses around the die exit reveal a stress boundary layer and a stress imbalance between the boundaries and the bulk. The magnitude of the stress imbalance exceeds the melt strength in the experiments with the LLDPE which causes the failure of the material in the superficial layers and results in the emergence of the sharkskin instability. In the experiments with the LDPE, the magnitude of the stress imbalance remains smaller than the melt strength which explains the lack of an instability. The measured shear stresses around the die exit are significantly smaller than the tensile stresses, suggesting that the shear component of the flow plays no significant role in the emergence of the sharkskin instability.     

Vortical Structures in Rotating Uniformly Sheared Turbulence

Rotating homogeneous turbulence with and without mean uniform shear is investigated numerically. It is found that in the shearless case the two-dimensionalization process is most effective when the initial small-scale Rossby number is around unity and the resonant triad interactions play a central role in the process. The vortical structures are studied systematically by changing the relative strength of the mean shear and the system rotation as well as the sense of rotation. (The system is called cyclonic (or anti-cyclonic) when the direction of the vorticity associated with the rotation is the same as (or opposite to) that of the mean shear.) A distinct coherent structure appears in the anti-cyclonic system when the vorticities associated with the rotation and the mean shear cancel out, i.e. the absolute vorticity of the mean shear vanishes. For the linearly most unstable case in the anti-cyclonic system, the vortex tubes develop in the sheared direction, which is caused by instability of vortex layers. For linearly stable cases in both the cyclonic and the anti-cyclonic systems, there appear three typical structures, that is, the oblique vortex tubes, the pancake-like structures and the ribbon-like structures. It is interesting that the flow behaves quite differently between the cyclonic and anti-cyclonic systems even at the same Bradshaw number.     

Linear instability characteristics of incompressible coaxial jets

The present study deals with the local linear instability of axisymmetric coaxial jets with a duct wall separating the two streams. The flow is assumed to be locally parallel, inviscid and incompressible. The objective of the work is to understand how the various parameters describing this flow geometry (i.e. boundary layers thicknesses at the exit, velocity ratio, wall thickness) may influence the instability of the flow and, in particular, the convective/absolute instability transition. A specific family of profiles is chosen for the modelling of the mean undisturbed flow and a spatial stability analysis is performed in order to identify the unstable modes and to assess how they are affected by the wake region behind the wall. An absolutely unstable mode is found, and its characteristics, depending on the velocity ratio and shear layers thicknesses, are determined. Results show that the absolute unstable mode is present only for a limited range of velocity ratios and that the corresponding frequency is almost constant if normalized with the mean velocity and wake thickness. This frequency value and the extension of the range of velocity ratios is similar to those found in the experiments on a similar geometry. Finally, a specific velocity ratio is found that maximizes the region at the jet exit for which an absolute instability behind the wall is present. This may increase the possibility for the onset of a global mode that may sustain the instability of the whole jet, enhancing considerably the mixing and entrainment characteristics between the two streams.     

Shallow water modeling of rolling pad instability in liquid metal batteries

Magnetohydrodynamically induced interface instability in liquid metal batteries is analyzed. The batteries are represented by a simplified system in the form of a rectangular cell, in which strong vertical electric current flows through three horizontal layers: the layer of a heavy metal at the bottom, the layer of a light metal at the top, and the layer of electrolyte in the middle. A new two-dimensional nonlinear model based on the conservative shallow water approximation is derived and utilized in a numerical study. It is found that in the case of small density difference between the electrolyte and one of the metals, the instability closely resembles the rolling pad instability observed earlier in the aluminum reduction cells. When the two electrolyte-metal density differences are comparable, the dynamics of unstable systems is more complex and characterized by interaction between two nearly synchronized or nearly anti-synchronized interfacial waves.     

Dissipative instability of fluid flows with piecewise linear velocity profiles

The viscous dissipative instability of two flows with continuous spectrum of neutrally-stable perturbations in the absence of dissipation is investigated. Ranges of wave numbers in which viscosity leads to flow destabilization are determined for a shear discontinuity in a smoothly-stratified fluid. A shear flow with a velocity in the transition layer that depends linearly on the coordinate has a continuum of neutral modes even in the case of an unstratified fluid. When viscosity is present in one of the layers with constant velocity, one of the branches of the spectrum becomes unstable. When the viscosity is the same above and below the shear layer, dissipation only leads to the damping of the perturbations. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 14–19, November–December, 1986.     

Suppression or enhancement of interfacial instability in two-layer plane Couette flow of FENE-P fluids past a deformable solid layer

The linear stability of two-layer plane Couette flow of FENE-P fluids past a deformable solid layer is analyzed in order to examine the effect of solid deformability on the interfacial instability due to elasticity and viscosity stratification at the two-fluid interface. The solid layer is modeled using both linear viscoelastic and neo-Hookean constitutive equations. The limiting case of two-layer flow of upper-convected Maxwell (UCM) fluids is used as a starting point, and results for the FENE-P case are obtained by numerically continuing the UCM results for the interfacial mode to finite values of the chain extensibility parameter. For the case of two-layer plane Couette flow past a rigid solid surface, our results show that the finite extensibility of the polymer chain significantly alters the neutral stability boundaries of the interfacial instability. In particular, the two-layer Couette flow of FENE-P fluids is found to be unstable in a larger range of nondimensional parameters when compared to two-layer flow of UCM fluids. The presence of the deformable solid layer is shown to completely suppress the interfacial instability in most of the parameter regimes where the interfacial mode is unstable, while it could have a completely destabilizing effect in other parameter regimes even when the interfacial mode is stable in rigid channels. When compared with two-layer UCM flow, the two-layer FENE-P case is found in general to require solid layers with relatively lower shear modulii in order to suppress the interfacial instability. The results from the linear elastic solid model are compared with those obtained using the (more rigorous) neo-Hookean model for the solid, and good agreement is found between the two models for neutral stability curves pertaining to the two-fluid interfacial mode. The present study thus provides an important extension of the earlier analysis of two-layer UCM flow [V. Shankar, Stability of two-layer viscoelastic plane Couette flow past a deformable solid layer: implications of fluid viscosity stratification, J. Non-Newtonian Fluid Mech. 125 (2005) 143–158] to more accurate constitutive models for the fluid and solid layers, and reaffirms the central conclusion of instability suppression in two-layer flows of viscoelastic fluids by soft elastomeric coatings in more realistic settings.     

Crack growth instability studied by the strain energy density theory

The strain energy density theory has successfully been used to address the problem of material damage and structural failure in problems of engineering interest. The theory makes use of the strain energy density function, dW/dV, and focuses attention in its stationary values. The directions of crack growth and yielding are determined from the minimum and maximum values of dW/dV, respectively, along the circumference of a circle centered at the point of failure initiation. Failure by crack growth or yielding takes place when these values of dW/dV become equal to their critical values which are material constants. In the present work the basic principles of the strain energy density theory were reviewed. Furthermore, this theory was used to study three problems of structural failure, namely the problem of slow stable growth of an inclined crack in a plate subjected to uniaxial tension, the problem of fracture instability of a plate with a central crack and two notches, and the problem of unstable crack growth in a circular disc subjected to two equal and opposite forces. The results of stress analysis were combined with the strain energy density theory to obtain the whole history of crack growth from initiation to instability. A length parameter was introduced to define the fracture instability of a mechanical system. Fracture trajectories were obtained for fast unstable crack propagation.     

Thermal instability of a viscoelastic fluid in a fluid-porous system with a plane Poiseuille flow

The thermal convection of a Jeffreys fluid subjected to a plane Poiseuille flow in a fluid-porous system composed of a fluid layer and a porous layer is studied in the paper. A linear stability analysis and a Chebyshev τ-QZ algorithm are employed to solve the thermal mixed convection. Unlike the case in a single layer, the neutral curves of the two-layer system may be bi-modal in the proper depth ratio of the two layers. We find that the longitudinal rolls(LRs) only depend on the depth ratio. Wi...     

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.     

A numerical study of an annular liquid jet in a compressible gas medium

An annular liquid jet in a compressible gas medium has been examined using an Eulerian approach with mixed-fluid treatment. The governing equations have been solved by using highly accurate numerical methods. An adapted volume of fluid method combined with a continuum surface force model was used to capture the gas–liquid interface dynamics. The numerical simulations showed the existence of a recirculation zone adjacent to the nozzle exit and unsteady large vortical structures at downstream locations, which lead to significant velocity reversals in the flow field. It was found that the annular jet flow is highly unstable because of the existence of two adjacent shear layers in the annular configuration. The large vortical structures developed naturally in the flow field without external perturbations. Surface tension tends to promote the Kelvin–Helmholtz instability and the development of vortical structures that leads to an increased liquid dispersion. A decrease in the liquid sheet thickness resulted in a reduced liquid dispersion. It was identified that the liquid-to-gas density and viscosity ratios have opposite effects on the flow field with the reduced liquid-to-gas density ratio demoting the instability and the reduced liquid-to-gas viscosity ratio promoting the instability characteristics.     

On Stabilization of Multi-Layer Hele-Shaw and Porous Media Flows in the Presence of Gravity

Stabilization of multi-layer Hele-Shaw flows is studied here by including the influence of Rayleigh?CTaylor instability in our earlier work (Daripa, J. Stat. Mech. 12:28, 2008a) on stabilization of multi-layer Saffman?CTaylor instability. Furthermore, this article goes beyond our previous work with few extensions, improvements, new interpretations, and clarifications on the use of some terminologies. Results of two complete studies have been presented: the first investigates the effect of individually unstable interfaces on the overall stability of the flow, and the second studies the cumulative effect of unstable interfaces as well as unstable internal viscous layers. In each case, modal and absolute upper bounds on the growth rate are reported. Next, these bounds are used to investigate (i) stabilization of long waves on various interfaces; (ii) stabilization of all waves on all interfaces in comparison to pure Taylor instability; (iii) stabilization of disturbances on interior interfaces instead of exterior interfaces. In the first study, notions of partial and total stabilization with respect to the pure Taylor growth rate are introduced. Then necessary and sufficient conditions for partial and total stabilizations are found. Proof of stabilization of long waves on one of the two external interfaces in multi-layer flows is also proved. In the second study, an absolute upper bound is obtained in the presence of stabilizing density stratification across each internal interface even though all interfaces and layers have unstable viscous profiles. Exact results on the upper bounds, and necessary and sufficient conditions for control of instabilities driven by stable/unstable density stratification, unstable viscous layers and unstable interfaces are new and may be relevant to explain observed phenomena in many complex flows generating these kinds of viscous profiles and density stratification as they evolve. The present work builds upon and goes much further in details and new results than our previous work. The gravity effect included here brings with it restrictions which have not been addressed before in this multi-layer context.     

Shear layer instabilities in low-density transverse jets

The near-field shear layer instabilities forming in round transverse jets of variable (reduced) densities relative to the crossflow are investigated through gas-phase experiments. Jets composed of helium and nitrogen mixtures are injected from a converging nozzle mounted flush with an injection wall into air crossflow, allowing the jet-to-crossflow density ratio S to be varied between 1.00, the equidensity case, and 0.14, at constant jet Reynolds number Re _j ?=?1,800. Jet-to-crossflow momentum flux ratios J are examined in the range $\infty>J\geq5$ at incremental values of the density ratio S. The results of single-component hotwire measurements in the jet shear layer indicate that a transition to global instability likely occurs as J is brought below approximately 10, and/or as S is brought below approximately 0.45?C0.40. This observation appears to link many previous independent studies of both equidensity transverse jets and low-density free jets, which may become globally unstable under alteration of J and S, respectively. However, the dynamical character of the transition to global instability in the low-density transverse jet displays differences under independent variation of J and S, which may indicate the predominance of different modes.     

Instability of a vertical columnar vortex in a stratified fluid

In this article, we study experimentally the evolution of a vertical columnar vortex in a stratified fluid. Three different measurement techniques are used. Particle image velocimetry allows us to monitor the time evolution of the characteristics of the vortex (Froude and Reynolds numbers). Dye visualizations reveal the existence of an instability for Froude numbers smaller than one, which creates an undulation of the vortex centerline. Synthetic schlieren visualization shows that the density structure of the unstable mode is very similar to the structure found recently numerically for the radiative instability of a Lamb–Oseen vortex (Riedinger et al. in J Fluid Mech, 2010). The experimental stability diagram and unstable wavelengths are compared with these numerical results. A secondary instability associated with the presence of critical layers is also observed for Froude numbers larger than one.     

On the stability of plane parallel viscoelastic shear flows in the limit of infinite Weissenberg and Reynolds numbers

Elastic effects on the hydrodynamic instability of inviscid parallel shear flows are investigated through a linear stability analysis. We focus on the upper convected Maxwell model in the limit of infinite Weissenberg and Reynolds numbers. We study the effects of elasticity on the instability of a few classes of simple parallel flows, specifically plane Poiseuille and Couette flows, the hyperbolic-tangent shear layer and the Bickley jet.The equation for stability is derived and solved numerically using the spectral Chebyshev collocation method. This algorithm is computationally efficient and accurate in reproducing the eigenvalues. We consider flows bounded by walls as well as flows bounded by free surfaces. In the inviscid, nonelastic case all the flows we study are unstable for free surfaces. In the case of wall bounded flow, there are instabilities in the shear layer and Bickley jet flows. In all cases, the effect of elasticity is to reduce and ultimately suppress the inviscid instability.     

Optimal disturbances of a dusty-gas plane-channel flow with a nonuniform distribution of particles

The classical stability theory for multiphase flows, based on an analysis of one (most unstable) mode, is generalized. A method for studying an algebraic (non-modal) instability of a disperse medium, which consists in examining the energy of linear combinations of three-dimensional modes with given wave vectors, is proposed. An algebraic instability of a dusty-gas flow in a plane channel with a nonuniform particle distribution in the form of two layers arranged symmetrically with respect to the flow axis is investigated. For all possible values of governing parameters, the optimal disturbances of the disperse flow have zero wavenumber in the flow direction, which indicates their banded structure (“streaks”). The presence of dispersed particles in the flow increases the algebraic instability, since the energy of optimal disturbances in the disperse medium exceeds that for the pure-fluid flow. It is found that for a homogeneous particle distribution the increase in the energy of optimal perturbations is proportional to the square of the sum of unity and the particle mass concentration and is almost independent of particle inertia. For a non-uniform distribution of the dispersed phase, the largest increase in the initial energy of disturbances is achieved in the case when the dust layers are located in the middle between the center line of the flow and the walls.     

Investigation the effects of shear rate on stationary droplets coalescence by lattice Boltzmann

Coalescence is the process by which two or more droplets merge to form a single droplet. It can take place in many processes, ranging from meteorology to astrophysics. When two stationary droplets are suspended in a bulk vapor, if the gap between the droplets has been smaller than a critical value, the two droplets will coalesce. In this paper, single component, two phase flow is modeled under shear flow using a free energy lattice Boltzmann approach and the coalescence of stationary droplets are investigated for different cases of radius and shear rate. The results show that there is a critical gap between droplets and for the values larger than that they will not coalesce. Also in the case of constant thermophysical properties, this critical gap is a function of droplet radius and shear rate.     

Experiments on entrainment in an annulus with and without velocity gradient across the density interface

The entrainment process in a two layer density stratified fluid column was studied experimentally by imposing external shear stress on one or both layers. The experiments have been conducted in an annular tank containing two water layers of different salt concentration and the shear stress was applied by means of rotating screens. The following quantities were measured: the screen velocity (which was kept constant during each experiment), the stress at the upper screen, and vertical profiles of circumferential velocity and density at different radial locations. When equal stress was imposed at the surface of the upper layer and at the bottom of the lower layer, entrainment took place from the two sides of the density interface at equal rate so that the interface was stationary in the central position between the two screens and there was no velocity gradient across the interface. The dependence of the entrainment coefficient on Richardson number obtained in these experiments was similar in form to that obtained in the shear-free experiments with an oscillating grid (e.g. Nokes 1988). When a shear stress was applied at the upper surface only, the upper layer depth increased with time and a velocity gradient existed at the interface. The influence of the interfacial velocity gradient on the entrainment rate was studied by comparing the rates obtained with and without this velocity gradient. The entrainment rates were approximately the same for high values of the Richardson number while at low Richardson number the entrainment rate was much larger when a velocity gradient existed across the interface. The main results of this work are as follows:

1. Despite the curved geometry of the annular system, the dependence of the entrainment coefficient on Richardson number for shear-free interface experiments was found to be similar in form to that obtained for oscillating grid experiments.
2. The entrainment across the interface is due to turbulent energy generated at some distance from the interface by an external source (i.e. shear stress induced by a screen) and due to turbulence produced locally at the interface by a velocity gradient. The relative contribution of each turbulence source to the total entrainment was found to depend on the stability of the interface.

     

Linear Dynamics of Perturbations in Flows with Constant Horizontal Shear

The dynamics of perturbations in shallow water and incompressible stratified fluid flows with constant horizontal shear are described using the nonmodal analysis. It is shown that the shear flow perturbations can be divided into two classes on the basis of the potential vorticity: rapidly oscillating wave perturbations with zero potential vorticity and slow vortex perturbations with nonzero potential vorticity. In the cases of weak and strong shear the main features of the dynamics of wave and vortex perturbations are studied analytically (using the WKBG method) and numerically. It is shown that for large times the wave perturbation energy increases linearly, i.e., the shear flow is algebraically unstable due to the growth of rapid wave perturbations. This instability can be of importance in processes of turbulence development and surface and internal wave generation.     

Instability theory of shock wave in a channel

The instability theory of shock wave was extended from the case with an infinitefront to the case of a channel with a rectangular cross section.First,themathematical formulation of the problem was given which included a system ofdisturbed equations and three kinds of boundary conditions.Then,the general solutionsof the equations upstream and downstream were given and each contained fiveconstants to be determined.Thirdly,under one boundary condition and oneassumption,it was proved that all of the disturbances in front of the shock front andone of the two acoustic disturbances behind the shock front should be zero.Theboundary condition was that all of the disturbed physical quantities should approach tozero at infinity.The assumption was that only the unstable shock wave was concernedhere.So it was reasonable to assumeω=iγ,γwas the instability growth rate andwas a positive real number.Another kind of boundary conditions was that the normaldisturbed velocities should be zero at the solid wall of the cha