ALBERT

Description: A new formulation of the stability of boundary-layer flows in pressure gradients is presented, taking into account the spatial development of the flow and utilizing a special coordinate transformation. The formulation assumes that disturbance wavelength and eigenfunction vary downstream no more rapidly than the boundary-layer thickness, and includes all terms nominally of order R-1in the boundary-layer Reynolds number R. In Blasius flow, the present approach is consistent with that of Bertolotti et al. (1992) to O(R-1) but simpler (i.e. has fewer terms), and may best be seen as providing a parametric differential equation which can be solved without having to march in space. The computed neutral boundaries depend strongly on distance from the surface, but the one corresponding to the inner maximum of the streamwise velocity perturbation happens to be close to the parallel flow (Orr-Sommerfeld) boundary. For this quantity, solutions for the Falkner-Skan flows show the effects of spatial growth to be striking only in the presence of strong adverse pressure gradients. As a rational analysis to O(R-1) demands inclusion of higher-order corrections on the mean flow, an illustrative calculation of one such correction, due to the displacement effect of the boundary layer, is made, and shown to have a significant destabilizing influence on the stability boundary in strong adverse pressure gradients. The effect of non-parallelism on the growth of relatively high frequencies can be significant at low Reynolds numbers, but is marginal in other cases. As an extension of the present approach, a method of dealing with non-similar flows is also presented and illustrated. However, inherent in the transformation underlying the present approach is a lower-order non-parallel theory, which is obtained by dropping all terms of nominal order R-1except those required for obtaining the lowest-order solution in the critical and wall layers. It is shown that a reduced Orr-Sommerfeld equation (in transformed coordinates) already contains the major effects of non-parallelism. © 1995, Cambridge University Press. All rights reserved.

Description: It has long been known from linear stability theory that heating a surface immersed in water flow tends to stabilize the boundary layer on the surface, suggesting that there may be a corresponding delay in transition. Experiments confirm the suggestion, but based on intermittency data on a heated body of revolution (Lauchle and Gurney 1984) it has been inferred that incremental changes in transition Reynolds number diminish as the overheat increases. The parameter chosen to locate transition in the analysis leading to this conclusion corresponds to the point where the intermittency is 0.5. However, intermittency distributions in the transition zone on an axisymmetric body may contain what have been called 'subtransitions' (Narasimha 1984). Taking this possibility into account, we formulate here a model for the variation of intermittency with flow Reynolds number at a fixed station on the body, as in the experiments. The rate at which turbulent spots merge with each other is shown to determine the location of subtransition. The transition onset Reynolds number (corresponding to the location where intermittency begins to depart from zero), inferred from the data on the basis of this model, shows a continuing increase with the temperature overheat, a trend in closer agreement with stability theory; but the axisymmetric body geometry results in a very short transition zone, countering in part the benefits of transition delay.

Description: We first demonstrate that, if the contributions of higher-order mean flow are ignored, the parabolized stability equations (Bertolotti et al. 1992) and the 'full' non-parallel equation of Govindarajan & Narasimha (1995, hereafter GN95) are both equivalent to order R-1 in the local Reynolds number R to Gaster's (1974) equation for the stability of spatially developing boundary layers. It is therefore of some concern that a detailed comparison between Gaster (1974) and GN95 reveals a small difference in the computed amplitude ratios. Although this difference is not significant in practical terms in Blasius flow, it is traced here to the approximation, in Gaster's method, of neglecting the change in eigenfunction shape due to flow non-parallelism. This approximation is not justified in the critical and the wall layers, where the neglected term is respectively O(R-2/3) and O(R-1) compared to the largest term. The excellent agreement of GN95 with exact numerical simulations, on the other hand, suggests that the effect of change in eigenfunction is accurately taken into account in that paper.

Description: Although the critical Reynolds number for linear instability of the laminar flow in a straight pipe is infinite, we show that it is finite for a divergent pipe, and approaches infinity as the inverse of the divergence angle. The velocity profile at the threshold of inviscid stability is obtained. A non-parallel analysis yields linear instability at surprisingly low Reynolds numbers, of about 150 for a divergence of 3°, which would suggest a role for such instabilities in the transition to turbulence. A multigrid Poisson equation solver is employed for the basic flow, and an extended eigenvalue method for the partial differential equations describing the stability. © 2005 Cambridge University Press.

Description: According to parallel-flow theory based on the Orr-Sommerfeld equation, a mixing layer is unstable at all Reynolds numbers. However this is untenable from energy considerations, which demand that there exist a non-zero Reynolds number below which disturbances cannot extract net energy from the mean flow. It is shown here that a linear stability analysis of similarity solutions of the plane mixing layer, including the effects of flow non-parallelism using the minimal composite theory and the properties of adjoints, following Govindarajan & Narasimha (Theor. Comput. Fluid Dyn. vol. 19, 2005, p. 229) resolves the issue by yielding a non-zero critical Reynolds number for co-flowing streams of any velocity ratio. The critical Reynolds number for the total disturbance kinetic energy is found to be close to 30 for all velocity ratios in the range from zero to unity. © 2006 Cambridge University Press.

Description: A comprehensive study of the effect of wall heating or cooling on the linear, transient and secondary growth of instability in channel flow is conducted. The effect of viscosity stratification, heat diffusivity and of buoyancy are estimated separately, with some unexpected results. From linear stability results, it has been accepted that heat diffusivity does not affect stability. However, we show that realistic Prandtl numbers cause a transient growth of disturbances that is an order of magnitude higher than at zero Prandtl number. Buoyancy, even at fairly low levels, gives rise to high levels of subcritical energy growth. Unusually for transient growth, both of these are spanwise-independent and not in the form of streamwise vortices. At moderate Grashof numbers, exponential growth dominates, with distinct Poiseuille-Rayleigh-Bénard and Tollmien-Schlichting modes for Grashof numbers up to ∼ 25 000, which merge thereafter. Wall heating has a converse effect on the secondary instability compared to the primary instability, destabilizing significantly when viscosity decreases towards the wall. It is hoped that the work will motivate experimental and numerical efforts to understand the role of wall heating in the control of channel and pipe flows. © 2007 Cambridge University Press.

Description: The linear stability of variable viscosity, miscible core-annular flows is investigated. Consistent with pipe flow of a single fluid, the flow is stable at any Reynolds number when the magnitude of the viscosity ratio is less than a critical value. This is in contrast to the immiscible case without interfacial tension, which is unstable at any viscosity ratio. Beyond the critical value of the viscosity ratio, the flow can be unstable even when the more viscous fluid is in the core. This is in contrast to plane channel flows with finite interface thickness, which are always stabilized relative to single fluid flow when the less viscous fluid is in contact with the wall. If the more viscous fluid occupies the core, the axisymmetric mode usually dominates over the corkscrew mode. It is demonstrated that, for a less viscous core, the corkscrew mode is inviscidly unstable, whereas the axisymmetric mode is unstable for small Reynolds numbers at high Schmidt numbers. For the parameters under consideration, the switchover occurs at an intermediate Schmidt number of about 500. The occurrence of inviscid instability for the corkscrew mode is shown to be consistent with the Rayleigh criterion for pipe flows. In some parameter ranges, the miscible flow is seen to be more unstable than its immiscible counterpart, and the physical reasons for this behaviour are discussed. A detailed parametric study shows that increasing the interface thickness has a uniformly stabilizing effect. The flow is least stable when the interface between the two fluids is located at approximately 0.6 times the tube radius. Unlike for channel flow, there is no sudden change in the stability with radial location of the interface. The instability originates mainly in the less viscous fluid, close to the interface. © 2007 Cambridge University Press.

Description: We study the stability of a vortex in an axisymmetric density distribution. It is shown that a light-cored vortex can be unstable in spite of the ‘stable stratification’ of density. Using a model flow consisting of step jumps in vorticity and density, we show that a wave interaction mediated by shear is the mechanism for the instability. The requirement is for the density gradient to be placed outside the vortex core but within the critical radius of the Kelvin mode. Conversely, a heavy-cored vortex, found in other studies to be unstable in the centrifugal Rayleigh–Taylor sense, is stabilized when the density jump is placed in this region. Asymptotic solutions at small Atwood numberAtshow growth rates scaling asAt1/3close to the critical radius, andAt1/2further away. By considering a family of vorticity and density profiles of progressively increasing smoothness, going from a step to a Gaussian, it is shown that sharp gradients are necessary for the instability of the light-cored vortex, consistent with recent work which found Gaussian profiles to be stable. For sharp gradients, it is argued that wave interaction can be supported due to the presence ofquasi-modes. Probably for the first time, a quasi-mode which decays exponentially is shown to interact with a neutral wave to give exponential growth in the combined case. We finally study the nonlinear stages using viscous direct numerical simulations. The initial exponential instability of light-cored vortices is arrested due to a restoring centrifugal buoyancy force, leading to stable non-axisymmetric structures, such as a tripolar state for an azimuthal wavenumber of 2. The study is restricted to two dimensions, and neglects gravity.