Instability of Combustion

THE following equations governing the phenomenon of intrinsic instability of combustion, leading to low frequency oscillations in a rocket motor using a single liquid propellant, were derived and investigated by L. Crocco.

Nonlinear saturation of ion-acoustic instability

It is proposed that the wave mediated indirect wave-particle interaction may be responsible for nonlinear saturation of current driven low frequency ion-acoustic turbulence. This process decreases the growth rate and increases the damping rate of the wave. Comparison has been made with some experiments.

Global lopsided instability in a purely stellar galactic disc

It is shown that pure exponential discs in spiral galaxies are capable of supporting slowly varying discrete global lopsided modes, which can explain the observed features of lopsidedness in the stellar discs. Using linearized fluid dynamical equations with the softened self-gravity and pressure of the perturbation as the collective effect, we derive self-consistently a quadratic eigenvalue equation for the lopsided perturbation in the galactic disc. On solving this, we find that the ground-state mode shows the observed characteristics of the lopsidedness in a galactic disc, namely the fractional Fourier amplitude A(1), increases smoothly with the radius. These lopsided patterns precess in the disc with a very slow pattern speed with no preferred sense of precession. We show that the lopsided modes in the stellar disc are long-lived because of a substantial reduction (approximately a factor of 10 compared to the local free precession rate) in the differential precession. The numerical solution of the equations shows that the groundstate lopsided modes are either very slowly precessing stationary normal mode oscillations of the disc or growing modes with a slow growth rate depending on the relative importance of the collective effect of the self-gravity. N-body simulations are performed to test the spontaneous growth of lopsidedness in a pure stellar disc. Both approaches are then compared and interpreted in terms of long-lived global m = 1 instabilities, with almost zero pattern speed.

Surface instability of confined elastic bilayers: Theory and simulations

The surface of a soft elastic film becomes unstable and forms a self-organized undulating pattern because of adhesive interactions when it comes in contact proximity with a rigid surface. For a single film, the pattern length scale lambda, which is governed by the minimization of the elastic stored energy, gives lambda similar to 3h, where h is the film thickness. Based on a linear stability analysis and simulations of adhesion and debonding, we consider the contact instability of an elastic bilayer, which provides greater flexibility in the morphological control of interfacial instability. Unlike the case of a single film, the morphology of the contact instability patterns, debonding distance, and debonding force in a bilayer can be controlled in a nonlinear way by varying the thicknesses and shear moduli of the films. Interestingly, the pattern wavelength in a bilayer can be greatly increased or decreased compared to a single film when the adhesive contact is formed by the stiffer or the softer of the two films, respectively. In particular, lambda as small as 0.5h can be obtained. This indicates a new strategy for pattern miniaturization in elastic contact lithography.

Instability of Coupled Longitudinal and Transverse Electromagnetic Modes of Wave‐Propagation in Bounded Streaming Plasmas

The instability of coupled longitudinal and transverse electromagnetic modes associated with long wavelengths is studied in bounded streaming plasmas. The main conclusions are as follows: (i) For long waves for which O (k 2)=0, in the absence of relative streaming motion of electrons and ions and aωp/c<0.66, the whole spectrum of harmonic waves is excited due to finite temperature and boundary effects consisting of two subseries. One of these subseries can be identified with Tonks-Dattner resonance oscillations for the electrons, and arises primarily due to the electrons with frequencies greater than the electrostatic plasma frequency corresponding to the electron density in the midplane in the undisturbed state. The other series arises primarily due to ion motion. When aωp/c>0.66, in addition to the above spectrum of harmonic waves, the system admits an infinite number of growing and decaying waves. The instability associated with these modes is found to arise due to the interaction of the waves inside the plasma with the external electromagnetic field. (ii) For modes with comparatively shorter wavelengths for which O (k3)=0, the coupling due to finite temperature sets in, and it is found that the two series of harmonic waves obtained in (i) deriving energy from the transverse modes also become unstable. Thus, for these wavelengths the system admits three sets of growing and decaying modes, first two for all values of aωp/c and the third for (aωp/c) > 0.66. (iii) The presence of streaming velocities introduces various other coupling mechanisms, and we find that even for the wavelengths for which O (k2)=0, we get three sets of growing and decaying waves. The numerical values for the growth rates show that the streaming velocities enhance the growth rates of instability significantly.

Two-stream instability in plasma

The problem of two-stream instability in plasma is studied by specifying the importance of initial magnetic field associated with the motion of the charged particles and the boundary effects. In Part I the accurate initial steady state is studied when the streams of electrons and ions move with different uniform speeds in plasmas with plane and cylindrical geometry. In Part II, in order to show the effects of finiteness and inhomogeneity of the system, small transverse plasma oscillations are studied in the case of plane plasmas. The role of plasma-sheath oscillations at the boundaries is found to be very important in driving the instabilities associated with the electromagnetic modes. The numerical estimates of the growth rates of the instability are given for the specific case of the physical data in discharge tubes.

Instability of the self-similar flow into a cavity

In this paper we have investigated the instability of the self-similar flow behind the boundary of a collapsing cavity. The similarity solutions for the flow into a cavity in a fluid obeying a gas law p = Kργ, K = constant and 7 ≥ γ > 1 has been solved by Hunter, who finds that for the same value of γ there are two self-similar flows, one with accelerating cavity boundary and other with constant velocity cavity boundary. We find here that the first of these two flows is unstable. We arrive at this result only by studying the propagation of disturbances in the neighbourhood of the singular point.

The effect of viscosity and surface tension on the Taylor instability of the surface of a cavitation bubble*

This paper is devoted to a consideration of the following problem: A spherical mass of fluid of density varrho1, viscosity μ1 and external radius R is surrounded by a fluid of density varrho2 and viscosity μ2.The fluids are immiscible and incompressible. The interface is accelerated radially by g1: to study the effect of viscosity and surface tension on the stability of the interface. By analyzing the problem in spherical harmonics the mathematical problem is reduced to one of solution of the characteristic determinant equation. The particular case of a cavity bubble, where the viscosity μ1 of the fluid inside the bubble is negligible in comparison with the viscosity μ2 of the fluid outside the bubble, is considered in some detail. It is shown that viscosity has a stabilizing role on the interface; and when g1 > T(n − 1) (n + 2)/R2(varrho2 − varrho1) the stabilizing role of both viscosity and surface tension is more pronounced than would result when either of them is taken individually.

Localised instability in unsteady separation bubbles

We propose in experimental method to study the instability of thin unsteady separation bubbles, i.e. of unsteady boundary layers with reverse flow. The unsteady boundary layer is created by controlled temporal and spatial variations of the velocity external to the boundary layer. We present results of the evolution of instability in different temporally varying flows in a shallow angle diffuser. Depending on the extent of reverse flow in the boundary we observe that instability can be spatially localised.

Nonlinear intrinsic instability of solid propellant combustion including gas-phase thermal inertia

The problem of homogeneous solid propellant combustion instability is studied with a one-dimensional flame model, including the effects of gas-phase thermal inertia and nonlinearity. Computational results presented in this paper show nonlinear instabilities inherent in the equations, due to which periodic burning is found even under steady ambient conditions such as pressure. The stability boundary is obtained in terms of Denison-Baum parameters. It is found that inclusion of gas-phase thermal inertia stabilizes the combustion. Also, the effect of a distributed heat release in the gas phase, compared to the flame sheet model, is to destabilize the burning. Direct calculations for finite amplitude pressure disturbances show that two distinct resonant modes exist, the first one near the natural frequency as obtained from intrinsic instability analysis and a second mode occurring at a much higher driving frequency. It is found that er rn in the low frequency region, the response of the propellant is significantly affected by the specific type of gas-phase chemical heat-release model employed. Examination of frequency response function reveals that the role of gas-phase thermal inertia is to stabilize the burning near the first resonant mode. Calculations made for different amplitudes of driving pressure show that the mean burning rate decreases with increasing amplitude. Also, with an increase in the driving amplitude, higher harmonics are generated in the burning rate.

Adhesion-Induced Instability in Asperities

Adhesive forces between two approaching asperities will deform the asperities, and under certain conditions this will result in a sudden runaway deformations leading to a jump-to-contact instability. We present finite element-based numerical studies on adhesion-induced deformation and instability in asperities. We consider the adhesive force acting on an asperity, when it is brought near a rigid half-space, due to van der Waals interaction between the asperity and the half-space. The adhesive force is considered to be distributed over the volume of the asperity (body force), thus resulting in more realistic simulations for the length scales considered. Iteration scheme based on a ``residual stress update'' algorithm is used to capture the effect of deformation on the adhesion force, and thereby the equilibrium configuration and the corresponding force. The numerical results are compared with the previous approximate analytical solutions for adhesion force, deformation of the asperity and adhesion-induced mechanical instability (jump-to-contact). It is observed that the instability can occur at separations much higher,and could possibly explain the higher value of instability separation observed in experiments. The stresses in asperities, particularly in case of small ones, are found to be high enough to cause yielding before jump -to-contact. The effect of roughness is considered by modeling a spherical protrusion on the hemispherical asperity.This small-scale roughness at the tip of the asperities is found to control the deformation behavior at small separations, and hence are important in determining the friction and wear due to the jump-to-contact instability.

Karyotype instability in the ponerine ant genus Diacamma

The queenless ponerine ant Diacamma ceylonense and a population of Diacamma from the Nilgiri hills which we refer to as `nilgiri', exhibit interesting similarities as well as dissimilarities. Molecular phylogenetic study of these morphologically almost similar taxa has shown that D ceylonense is closely related to `nilgiri' and indicates that `nilgiri' is a recent diversion in the Diacamma phylogenetic tree. However, there is a striking behavioural difference in the way reproductive monopoly is maintained by the respective gamergates (mated egg laying workers), and there is evidence that they are genetically differentiated, suggesting a lack of gene flow To develop a better understanding of the mechanism involved in speciation of Diacamma, we have analysed karyotypes of D. ceylonense and `nilgiri' In both, we found surprising inter-individual and intra-individual karyotypic mosaicism. The observed numerical variability, both at intra-individual and inter-individual levels, does not appear to have hampered the sustainability of the chromosomal diversity in each population under study Since the related D. indicum, displays no such intra-individual or inter-Individual variability whatsoever under identical experimental conditions, these results are unlikely to he artifacts. Although no known mechanisms can account for the observed karyotypic variability of this nature, we believe that the present findings on the ants under study would provide opportunities for exciting new discoveries concerning the origin, maintenance and significance of intra-individual and inter-individual karyotypic mosaicism.

A sufficient criterion for Rayleigh-Taylor instability of incompressible viscous three-layer flow

Using normal mode analysis Rayleigh-Taylor instability is investigated for three-layer viscous stratified incompressible steady flow, when the top 3rd and bottom 1st layers extend up to infinity, the middle layer has a small thickness δ. The wave Reynolds number in the middle layer is assumed to be sufficiently small. A dispersion relation (a seventh degree polynomial in wave frequency ω) valid up to the order of the maximal value of all possible Kj (j less-than-or-equals, slant 0, K is the wave number) in each coefficient of the polynomial is obtained. A sufficient condition for instability is found out for the first time, pursuing a medium wavelength analysis. It depends on ratios (α and β) of the coefficients of viscosity, the thickness of the middle layer δ, surface tension ratio T and wave number K. This is a new analytical criterion for Rayleigh-Taylor instability of three-layer fluids. It recovers the results of the corresponding problem for two-layer fluids. Among the results obtained, it is observed that taking the coefficients of viscosity of 2nd and 3rd layers same can inhibit the effect of surface tension completely. For large wave number K, the thickness of the middle layer should be correspondingly small to keep the domain of dependence of the threshold wave number Kc constant for fixed α, β and T.

A sufficient criterion for Kelvin–Helmholtz instability in the magnetopause boundary‐layer region

The Kelvin–Helmholtz instability has been investigated for the magnetopause boundary‐layer region by the linearized method. The plasma in magnetosheath and magnetopause is assumed to be semi‐infinitely extended homogeneous, nondissipative, and incompressible. It is observed that, if one relation of two plasma speeds on the two sides of the magnetopause, wave number, and boundary‐layer thickness exceeds a certain threshold, the instability sets in. This new analytically sufficient criterion for excitation of instability in the three‐layer plasma flow generalizes the corresponding Chandrasekhar’s instability criterion for two‐layer plasma flow. The known results have been recovered and modified, the new results have been discovered. It is proved that the velocity threshold for the onset of instability is low when the magnitude of the magnetosheath and boundary‐layer region magnetic field and the angle between them are small. Also the threshold depends on the direction of plasma flow. The following results are observed numerically. The growth of the instability is sensitive to the magnetic field direction in the magnetosheath. A slight variation in the magnetic field direction in the second region can substantially change the relative velocity threshold for instability. When the ratio of the density of the second and third layer (magnetosphere) increases or that of the first and third layer decreases, the threshold decreases. Apart from this a necessary criterion for instability is obtained for a particular case.