Phase diagram of a two-species lattice model with a linear instability

We discuss the properties of a one-dimensional lattice model of a driven system with two species of particles in which the mobility of one species depends on the density of the other. This model was introduced by Lahiri and Ramaswamy (Phys. Rev. Lett., 79, 1150 (1997)) in the context of sedimenting colloidal crystals, and its continuum version was shown to exhibit an instability arising from linear gradient couplings. In this paper we review recent progress in understanding the full phase diagram of the model. There are three phases. In the first, the steady state can be determined exactly along a representative locus using the condition of detailed balance. The system shows phase separation of an exceptionally robust sort, termed strong phase separation, which survives at all temperatures. The second phase arises in the threshold case where the first species evolves independently of the second, but the fluctuations of the first influence the evolution of the second, as in the passive scalar problem. The second species then shows phase separation of a delicate sort, in which long-range order coexists with fluctuations which do not damp down in the large-size limit. This fluctuation-dominated phase ordering is associated with power law decays in cluster size distributions and a breakdown of the Porod law. The third phase is one with a uniform overall density, and along a representative locus the steady state is shown to have product measure form. Density fluctuations are transported by two kinematic waves, each involving both species and coupled at the nonlinear level. Their dissipation properties are governed by the symmetries of these couplings, which depend on the overall densities. In the most interesting case,, the dissipation of the two modes is characterized by different critical exponents, despite the nonlinear coupling.

Cold gas in cluster cores: global stability analysis and non-linear simulations of thermal instability

We perform global linear stability analysis and idealized numerical simulations in global thermal balance to understand the condensation of cold gas from hot/virial atmospheres (coronae), in particular the intracluster medium (ICM). We pay particular attention to geometry (e.g. spherical versus plane-parallel) and the nature of the gravitational potential. Global linear analysis gives a similar value for the fastest growing thermal instability modes in spherical and Cartesian geometries. Simulations and observations suggest that cooling in haloes critically depends on the ratio of the cooling time to the free-fall time (t(cool)/t(ff)). Extended cold gas condenses out of the ICM only if this ratio is smaller than a threshold value close to 10. Previous works highlighted the difference between the nature of cold gas condensation in spherical and plane-parallel atmospheres; namely, cold gas condensation appeared easier in spherical atmospheres. This apparent difference due to geometry arises because the previous plane-parallel simulations focused on in situ condensation of multiphase gas but spherical simulations studied condensation anywhere in the box. Unlike previous claims, our non-linear simulations show that there are only minor differences in cold gas condensation, either in situ or anywhere, for different geometries. The amount of cold gas depends on the shape of tcool/tff; gas has more time to condense if gravitational acceleration decreases towards the centre. In our idealized plane-parallel simulations with heating balancing cooling in each layer, there can be significant mass/energy/momentum transfer across layers that can trigger condensation and drive tcool/tff far beyond the critical value close to 10.

Non-linear resonance in strings under narrow band random excitation, part I: Planar response and stability

Non-linear planar response of a string to planar narrow band random excitation is investigated in this paper. A response equation for the mean square deflection σ2 is obtained under a single mode approximation by using the equivalent linearization technique. It is shown that the response is triple valued, as in the case of harmonic excitation, if the centre frequency of excitation Ω lies in a certain specified range. The triple valued response occurs only if the excitation bandwidth β is smaller than a critical value βcrit which is a monotonically increasing function of the intensity of excitation. An approximate method of investigating the almost sure asymptotic stability of the solution is presented and regions of instability in the Ω-σ2 plane have been charted. It is shown that planar response can become unstable either due to an unbounded growth of the in-plane component of motion or due to a spontaneous appearance of an out-of-plane component.

On the origin of the inflectional instability of a laminar separation bubble

This is an experimental and theoretical Study of a laminar separation bubble and the associated linear stability mechanisms. Experiments were performed over a flat plate kept in a wind tunnel, with an imposed pressure gradient typical of an aerofoil that would involve a laminar separation bubble. The separation bubble was characterized by measurement of surface-pressure distribution and streamwise velocity using hot-wire anemometry. Single component hot-wire anemometry was also used for a detailed study of the transition dynamics. It was foundthat the so-called dead-air region in the front portion of the bubble corresponded to a region of small disturbance amplitudes, with the amplitude reaching a maximum value close to the reattachment point. An exponential growth rate of the disturbance was seen in the region upstream of the mean maximum height of the bubble, and this was indicative of a linear instability mechanism at work. An infinitesimal disturbance was impulsively introduced into the boundary layer upstream of separation location, and the wave packet was tracked (in an ensemble-averaged sense) while it was getting advected downstream. The disturbance was found to be convective in nature. Linear stability analyses (both the Orr-Sommerfeld and Rayleigh calculations) were performed for mean velocity profiles, starting from an attached adverse-pressure-gradient boundary layer all the way up to the front portion of the separation-bubble region (i.e. up to the end of the dead-air region in which linear evolution of the disturbance could be expected). The conclusion from the present work is that the primary instability mechanism in a separation bubble is inflectional in nature, and its origin can be traced back to upstream of the separation location. In other words, the inviscid inflectional instability of the separated shear layer should be logically seen as an extension of the instability of the upstream attached adverse-pressure-gradient boundary layer. This modifies the traditional view that pegs the origin of the instability in a separation bubble to the detached shear layer Outside the bubble, with its associated Kelvin-Helmholtz mechanism. We contendthat only when the separated shear layer has moved considerably away from the wall (and this happens near the maximum-height location of the mean bubble), a description by the Kelvin-Helmholtz instability paradigm, with its associated scaling principles, Could become relevant. We also propose a new scaling for the most amplified frequency for a wall-bounded shear layer in terms of the inflection-point height and the vorticity thickness and show it to be universal.

Non-linear response of a gyrostabilized platform to transient input

In this paper the response of a gyrostabilized platform subjected to a transient torque has been analyzed by deliberately introducing non-linearity into the command of the servomotor. The resulting third-order non-linear differential equation has been solved by using a transformation technique involving the displacement variable. The condition under which platform oscillations may grow with time or die with time are important from the point of view of platform stabilization. The effect of deliberate addition of non-linearity with a view to achieving the ideal response—that is, to bring the platform back to its equilibrium position with as few oscillations as possible—has been investigated. The conditions under which instability may set in on account of the small transient input and small non-linearity has also been discussed. The analysis is illustrated by means of a numerical example. The results of analysis are compared with numerical solutions obtained on a digital computer.

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.

Non-Linear vibration of an elastic string

Equations proposed in previous work on the non-linear motion of a string show a basic disagreement, which is here traced to an assumption about the longitudinal displacement u. It is shown that it is neither necessary nor justifiable to assume that u is zero; and also that the velocity of propagation of u disturbances in a string is different from that in an infinite medium, although this difference is usually negligible. After formulating the exact equations of motion for the string, a systematic procedure is described for obtaining approximations to these equations to any order, making only the assumption that the strain in the material of the string is small. The lowest order equations in this scheme are non-linear, and are used to describe the response of a string near resonance. Finally, it is shown that in the absence of damping, planar motion of a string is always unstable at sufficiently high amplitudes, the critical amplitude falling to zero at the natural frequency and its subharmonics. The effect of slight damping on this instability is also discussed.

Instability and dewetting of ultrathin solid viscoelastic films on homogeneous and heterogeneous substrates

Instability and dewetting engendered by the van der Waals force in soft thin (<100 nm) linear viscoelastic solid (e. g., elastomeric gel) films on uniform and patterned surfaces are explored. Linear stability analysis shows that, although the elasticity of the film controls the onset of instability and the corresponding critical wavelength, the dominant length-scale remains invariant with the elastic modulus of the film. The unstable modes are found to be long-wave, for which a nonlinear long-wave analysis and simulations are performed to uncover the dynamics and morphology of dewetting. The stored elastic energy slows down the temporal growth of instability significantly. The simulations also show that a thermodynamically stable film with zero-frequency elasticity can be made unstable in the presence of physico-chemical defects on the substrate and can follow an entirely different pathway with far fewer holes as compared to the viscous films. Further, the elastic restoring force can retard the growth of a depression adjacent to the hole-rim and thus suppress the formation of satellite holes bordering the primary holes. These findings are in contrast to the dewetting of viscoelastic liquid films where nonzero frequency elasticity accelerates the film rupture and promotes the secondary instabilities. Thus, the zero-frequency elasticity can play a major role in imposing a better-defined long-range order to the dewetted structures by arresting the secondary instabilities. (C) 2011 American Institute of Physics. doi: 10.1063/1.3554748]

Prediction of nonlocal scaling parameter for armchair and zigzag single-walled carbon nanotubes based on molecular structural mechanics, nonlocal elasticity and wave propagation

Atomic vibration in the Carbon Nanotubes (CNTs) gives rise to non-local interactions. In this paper, an expression for the non-local scaling parameter is derived as a function of the geometric and electronic properties of the rolled graphene sheet in single-walled CNTs. A self-consistent method is developed for the linearization of the problem of ultrasonic wave propagation in CNTs. We show that (i) the general three-dimensional elastic problem leads to a single non-local scaling parameter (e(0)), (ii) e(0) is almost constant irrespective of chirality of CNT in the case of longitudinal wave propagation, (iii) e(0) is a linear function of diameter of CNT for the case of torsional mode of wave propagation, (iv) e(0) in the case of coupled longitudinal-torsional modes of wave propagation, is a function which exponentially converges to that of axial mode at large diameters and to torsional mode at smaller diameters. These results are valid in the long-wavelength limit. (C) 2011 Elsevier Ltd. All rights reserved.

Weakly non-linear stability of a hydrodynamic accretion disc

When the cold accretion disc coupling between neutral gas and a magnetic field is so weak that the magnetorotational instability is less effective or even stops working, it is of prime interest to investigate the pure hydrodynamic origin of turbulence and transport phenomena. As the Reynolds number increases, the relative importance of the non-linear term in the hydrodynamic equation increases. In an accretion disc where the molecular viscosity is too small, the Reynolds number is large enough for the non-linear term to have new effects. We investigate the scenario of the `weakly non-linear' evolution of the amplitude of the linear mode when the flow is bounded by two parallel walls. The unperturbed flow is similar to the plane Couette flow, but with the Coriolis force included in the hydrodynamic equation. Although there is no exponentially growing eigenmode, because of the self-interaction, the least stable eigenmode will grow in an intermediate phase. Later, this will lead to higher-order non-linearity and plausible turbulence. Although the non-linear term in the hydrodynamic equation is energy-conserving, within the weakly non-linear analysis it is possible to define a lower bound of the energy (alpha A(c)(2), where A(c) is the threshold amplitude) needed for the flow to transform to the turbulent phase. Such an unstable phase is possible only if the Reynolds number >= 10(3-4). The numerical difficulties in obtaining such a large Reynolds number might be the reason for the negative result of numerical simulations on a pure hydrodynamic Keplerian accretion disc.

Linear stability, transient growth and the role of viscosity stratification in compressible Couette flow

Linear stability and the nonmodal transient energy growth in compressible plane Couette flow are investigated for two prototype mean flows: (a) the uniform shear flow with constant viscosity, and (b) the nonuniform shear flow with stratified viscosity. Both mean flows are linearly unstable for a range of supersonic Mach numbers (M). For a given M, the critical Reynolds number (Re) is significantly smaller for the uniform shear flow than its nonuniform shear counterpart; for a given Re, the dominant instability (over all streamwise wave numbers, α) of each mean flow belongs to different modes for a range of supersonic M. An analysis of perturbation energy reveals that the instability is primarily caused by an excess transfer of energy from mean flow to perturbations. It is shown that the energy transfer from mean flow occurs close to the moving top wall for “mode I” instability, whereas it occurs in the bulk of the flow domain for “mode II.” For the nonmodal transient growth analysis, it is shown that the maximum temporal amplification of perturbation energy, Gmax, and the corresponding time scale are significantly larger for the uniform shear case compared to those for its nonuniform counterpart. For α=0, the linear stability operator can be partitioned into L∼L̅ +Re2 Lp, and the Re-dependent operator Lp is shown to have a negligibly small contribution to perturbation energy which is responsible for the validity of the well-known quadratic-scaling law in uniform shear flow: G(t∕Re)∼Re2. In contrast, the dominance of Lp is responsible for the invalidity of this scaling law in nonuniform shear flow. An inviscid reduced model, based on Ellingsen-Palm-type solution, has been shown to capture all salient features of transient energy growth of full viscous problem. For both modal and nonmodal instability, it is shown that the viscosity stratification of the underlying mean flow would lead to a delayed transition in compressible Couette flow.

Secondary Perturbation Effects in Keplerian Accretion Disks: Elliptical Instability

Origin of turbulence in cold accretion disks, particularly in 3D, which is expected to be hydrodynamic but not magnetohydrodynamic, is a big puzzle. While the flow must exhibit some turbulence in support of the transfer of mass inward and angular momentum outward, according to the linear perturbation theory it should always be stable. We demonstrate that the 3D secondary disturbance to the primarily perturbed disk which exhibits elliptical vortices into the system solves the problem. This result is essentially applicable to the outer region of accretion disks in active galactic nuclei where the gas is significantly cold and neutral in charge and the magnetic Reynolds number is smaller than 10^4.

Relevance of local parallel theory to the linear stability of laminar separation bubbles

Laminar separation bubbles are thought to be highly non-parallel, and hence global stability studies start from this premise. However, experimentalists have always realized that the flow is more parallel than is commonly believed, for pressure-gradient-induced bubbles, and this is why linear parallel stability theory has been successful in describing their early stages of transition. The present experimental/numerical study re-examines this important issue and finds that the base flow in such a separation bubble becomes nearly parallel due to a strong-interaction process between the separated boundary layer and the outer potential flow. The so-called dead-air region or the region of constant pressure is a simple consequence of this strong interaction. We use triple-deck theory to qualitatively explain these features. Next, the implications of global analysis for the linear stability of separation bubbles are considered. In particular we show that in the initial portion of the bubble, where the flow is nearly parallel, local stability analysis is sufficient to capture the essential physics. It appears that the real utility of the global analysis is perhaps in the rear portion of the bubble, where the flow is highly non-parallel, and where the secondary/nonlinear instability stages are likely to dominate the dynamics.

Stochastically driven instability in rotating shear flows

The origin of hydrodynamic turbulence in rotating shear flows is investigated, with particular emphasis on the flows whose angular velocity decreases but whose specific angular momentum increases with the increasing radial coordinate. Such flows are Rayleigh stable, but must be turbulent in order to explain the observed data. Such a mismatch between the linear theory and the observations/experiments is more severe when any hydromagnetic/magnetohydrodynamic instability and then the corresponding turbulence therein is ruled out. This work explores the effect of stochastic noise on such hydrodynamic flows. We essentially concentrate on a small section of such a flow, which is nothing but a plane shear flow supplemented by the Coriolis effect. This also mimics a small section of an astrophysical accretion disc. It is found that such stochastically driven flows exhibit large temporal and spatial correlations of perturbation velocities and hence large energy dissipations of perturbation, which presumably generate the instability. A range of angular velocity (Omega) profiles of the background flow, starting from that of a constant specific angular momentum (lambda = Omega r(2); r being the radial coordinate) to a constant circular velocity (v(phi) = Omega r), is explored. However, all the background angular velocities exhibit identical growth and roughness exponents of their perturbations, revealing a unique universality class for the stochastically forced hydrodynamics of rotating shear flows. This work, to the best of our knowledge, is the first attempt to understand the origin of instability and turbulence in three-dimensional Rayleigh stable rotating shear flows by introducing additive noise to the underlying linearized governing equations. This has important implications to resolve the turbulence problem in astrophysical hydrodynamic flows such as accretion discs.

A multifold reduction in the transition Reynolds number, and ultra-fast mixing, in a micro-channel due to a dynamical instability induced by a soft wall

A dynamical instability is observed in experimental studies on micro-channels of rectangular cross-section with smallest dimension 100 and 160 mu m in which one of the walls is made of soft gel. There is a spontaneous transition from an ordered, laminar flow to a chaotic and highly mixed flow state when the Reynolds number increases beyond a critical value. The critical Reynolds number, which decreases as the elasticity modulus of the soft wall is reduced, is as low as 200 for the softest wall used here (in contrast to 1200 for a rigid-walled channel) The instability onset is observed by the breakup of a dye-stream introduced in the centre of the micro-channel, as well as the onset of wall oscillations due to laser scattering from fluorescent beads embedded in the wall of the channel. The mixing time across a channel of width 1.5 mm, measured by dye-stream and outlet conductance experiments, is smaller by a factor of 10(5) than that for a laminar flow. The increased mixing rate comes at very little cost, because the pressure drop (energy requirement to drive the flow) increases continuously and modestly at transition. The deformed shape is reconstructed numerically, and computational fluid dynamics (CFD) simulations are carried out to obtain the pressure gradient and the velocity fields for different flow rates. The pressure difference across the channel predicted by simulations is in agreement with the experiments (within experimental errors) for flow rates where the dye stream is laminar, but the experimental pressure difference is higher than the simulation prediction after dye-stream breakup. A linear stability analysis is carried out using the parallel-flow approximation, in which the wall is modelled as a neo-Hookean elastic solid, and the simulation results for the mean velocity and pressure gradient from the CFD simulations are used as inputs. The stability analysis accurately predicts the Reynolds number (based on flow rate) at which an instability is observed in the dye stream, and it also predicts that the instability first takes place at the downstream converging section of the channel, and not at the upstream diverging section. The stability analysis also indicates that the destabilization is due to the modification of the flow and the local pressure gradient due to the wall deformation; if we assume a parabolic velocity profile with the pressure gradient given by the plane Poiseuille law, the flow is always found to be stable.