Stability and transition in the flow of polymer solutions

The linear stability analysis of a plane Couette flow of viscoelastic fluid have been studied with the emphasis on two dimensional disturbances with wave number k similar to Re-1/2, where Re is Reynolds number based on maximum velocity and channel width. We employ three models to represent the dilute polymer solution: the classical Oldroyd-B model, the Oldroyd-B model with artificial diffusivity and the non-homogeneous polymer model. The result of the linear stability analysis is found to be sensitive to the polymer model used. While the plane Couette flow is found to be stable to infinitesimal disturbances for the first two models, the last one exhibits a linear instability.

On the bursting criterion for laminar separation bubbles

Traditionally, laminar separation bubbles have been characterised as being 'long' or 'short' on the basis of a two parameter 'bursting' criterion involving a pressure gradient parameter and Reynolds Number at separation. In the present work we suggest a refined bursting criterion, which takes into account not just the length of the bubble but also the maximum height of the bubble, thereby shedding some light on the less understood phenomenon of 'bursting' in laminar separation bubbles.

Effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics of a sphere

Experiments were conducted in water and wind tunnels on spheres in the Reynolds number range 6 x 10(3) to 6.5 x 10(5) to study the effect of natural ventilation on the boundary layer separation and near-wake Vortex shedding characteristics. In the subcritical range of Re (<2 x 10(5)), ventilation caused a marginal downstream shift in the location of laminar boundary layer separation; there was only a small change in the vortex shedding frequency. In the supercritical range (Re > 4 x 10(5)), ventilation caused a downstream shift in the mean locations of boundary layer separation and reattachment; these lines showed significant axisymmetry in the presence of venting. No distinct vortex shedding frequency was found. Instead, a dramatic reduction occurred in the wake unsteadiness at all frequencies. The reduction of wake unsteadiness is consistent with the reduction in total drag already reported. Based on the present results and those reported earlier, the effects of natural ventilation on the flow past a sphere can be categorized in two broad regimes, viz., weak and strong interaction regimes. In the weak interaction regime (subcritical Re), the broad features of the basic sphere are largely unaltered despite the large addition of mass in the near wake. Strong interaction is promoted by the closer proximity of the inner and outer shear layers at supercritical Re. This results in a modified and steady near-wake flow, characterized by reduced unsteadiness and small drag.

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 study on boundary-layer transition induced by free-stream turbulence

Boundary-layer transition at different free-stream turbulence levels has been investigated using the particle-image velocimetry technique. The measurements show organized positive and negative fluctuations of the streamwise fluctuating velocity component, which resemble the forward and backward jet-like structures reported in the direct numerical simulation of bypass transition. These fluctuations are associated with unsteady streaky structures. Large inclined high shear-layer regions are also observed and the organized negative fluctuations are found to appear consistently with these inclined shear layers, along with highly inflectional instantaneous streamwise velocity profiles. These inflectional velocity profiles are similar to those in the ribbon-induced boundary-layer transition. An oscillating-inclined shear layer appears to be the turbulent spot-precursor. The measurements also enabled to compare the actual turbulent spot in bypass transition with the simulated one. A proper orthogonal decomposition analysis of the fluctuating velocity field is carried out. The dominant flow structures of the organized positive and negative fluctuations are captured by the first few eigenfunction modes carrying most of the fluctuating energy. The similarity in the dominant eigenfunctions at different Reynolds numbers suggests that the flow prevails its structural identity even in intermittent flows. This analysis also indicates the possibility of the existence of a spatio-temporal symmetry associated with a travelling wave in the flow.

MPD arcjet cathode surface current density

The radial current density distribution on the cathode longitudinal surface of magnetoplasmadynamic arcjets for axisymmetric geometries has been obtained by simultaneous solution of the electromagnetic equations for a given uniform gas dynamic field. The problem formulation permits a parametric study of the effects of the Hall parameter and the magnetic Reynolds number. The solution for the current density distribution displays current concentrations at two locations, that is, at the upstream and downstream ends of the cathode. This result is in conformity with known experimental data. The parameters responsible for these current concentrations are identified. It is shown that the effect of the magnetic Reynolds number on the current density distribution is different depending on whether or not the Hall effect is included. This result is also found to be consistent with experimental data.

Some observations of tip-vortex cavitation

Cavitation has been observed in the trailing vortex system of an elliptic planform hydrofoil. A complex dependence on Reynolds number and gas content is noted at inception. Some of the observations can be related to tension effects associated with the lack of sufficiently large-sized nuclei. Inception measurements are compared with estimates of pressure in the vortex obtained from LDV measurements of velocity within the vortex. It is concluded that a complete correlation is not possible without knowledge of the fluctuating levels of pressure in tip-vortex flows. When cavitation is fully developed, the observed tip-vortex trajectory shows a surprising lack of dependence on any of the physical parameters varied, such as angle of attack, Reynolds number, cavitation number, and dissolved gas content.

Steady incompressible flow of cohesionless granular materials through a wedge-shaped hopper: Frictional 14kinetic solution to the smooth wall, radial gravity problem

Hybrid frictional-kinetic equations are used to predict the velocity, grain temperature, and stress fields in hoppers. A suitable choice of dimensionless variables permits the pseudo-thermal energy balance to be decoupled from the momentum balance. These balances contain a small parameter, which is analogous to a reciprocal Reynolds number. Hence an approximate semi-analytical solution is constructed using perturbation methods. The energy balance is solved using the method of matched asymptotic expansions. The effect of heat conduction is confined to a very thin boundary layer near the exit, where it causes a marginal change in the temperature. Outside this layer, the temperature T increases rapidly as the radial coordinate r decreases. In particular, the conduction-free energy balance yields an asymptotic solution, valid for small values of r, of the form T proportional r-4. There is a corresponding increase in the kinetic stresses, which attain their maximum values at the hopper exit. The momentum balance is solved by a regular perturbation method. The contribution of the kinetic stresses is important only in a small region near the exit, where the frictional stresses tend to zero. Therefore, the discharge rate is only about 2.3% lower than the frictional value, for typical parameter values. As in the frictional case, the discharge rate for deep hoppers is found to be independent of the head of material.

Biorthogonal series method for Oseen type flows

A biorthogonal series method is developed to solve Oseen type flow problems. The theory leads to a new set of eigenfunctions for a specific class of linear non-selfadjoint operators containing the biharmonic one. These eigenfunctions differ from those given earlier in the literature for the biharmonic operator. The method is applied to the problem of thermocapillary flow in a cylindrical liquid bridge of finite length with axial through flow. Flow and temperature distributions are obtained at leading order of an expansion for small surface tension Reynolds number and Prandtl number. Another related problem considered is that of cylindrical cavity flow. Solutions for both cases are presented in terms of biorthogonal series. The effect of axial through flow on velocity and temperature fields is discussed by numerical evaluation of the truncated analytical series. The presence of axial through flow not only convectively shifts the vortices induced by surface forces in the direction of the through flow, but also moves their centers toward the outer cylindrical boundary. This process can lead to significantly asymmetric flow structures.

Desorption from concentrated polymer solutions in good and poor solvents

The objective of the study was to investigate the effects of the nature of solvent and polymer concentration on the mass-transfer coefficients in desorption of solvents and to develop a correlation to predict them. Desorption was experimentally studied in a Lewis cell with concentrated binary solutions of polymer in good and poor solvents. The range of parameters covered are polymer weight fraction between 0.25 and 0.6, Reynolds number between 3 and 100; Schmidt number between 1.4 X lo6 and 2.5 X lo8, and Sherwood number between 3.5 X lo2 and 1.2 X lo4. Desorption from moderately concentrated solutions (polymer weight fraction -0.25) is gas-phase controlled. Studies with more concentrated solutions showed that the effects of solvent and concentration were such that corrections due to concentration-dependent diffusivity and viscosity as well as high flux had to be applied to the mass-transfer coefficients before they could be correlated.

Exact N-Wave Solutions for the Non-Planar Burgers Equation

An exact representation of N-wave solutions for the non-planar Burgers equation u(t) + uu(x) + 1/2ju/t = 1/2deltau(xx), j = m/n, m < 2n, where m and n are positive integers with no common factors, is given. This solution is asymptotic to the inviscid solution for Absolute value of x < square-root (2Q0 t), where Q0 is a function of the initial lobe area, as lobe Reynolds number tends to infinity, and is also asymptotic to the old age linear solution, as t tends to infinity; the formulae for the lobe Reynolds numbers are shown to have the correct behaviour in these limits. The general results apply to all j = m/n, m < 2n, and are rather involved; explicit results are written out for j = 0, 1, 1/2, 1/3 and 1/4. The case of spherical symmetry j = 2 is found to be 'singular' and the general approach set forth here does not work; an alternative approach for this case gives the large time behaviour in two different time regimes. The results of this study are compared with those of Crighton & Scott (1979).

Nonlinear effects on rotating disk flow in a cylindrical casing

The flow due to a finite disk rotating in an incompressible viscous fluid has been studied. A modified Newton-gradient finite difference scheme is used to obtain the solution of full Navier-Stokes equations numerically for different disk and cylinder sizes for a wide range of Reynolds numbers. The introduction of the aspect ratio and the disk-shroud gap, significantly alters the flow characteristics in the region under consideration, The frictional torque calculated from the flow data reveals that the contribution due to nonlinear terms is not negligible even at a low Reynolds number. For large Reynolds numbers, the flow structure reveals a strong boundary layer character.

Wind loading over the top portions of tall stacks with and without external landing platforms

The variation of the drag force near the top portions of tall stacks with and without external landing platforms, and with the exit open and closed, has been examined by model studies in a wind tunnel at Reynolds numbers of about 10(5). Pressure measurements on three models of different height to diameter ratios have been supplemented by flow visualisation studies. Observations confirm that when there is no platform, significant load enhancement over the top three to four diameters occurs, due to the high suction caused by the sharp separation of the flow over the top from the rim, in the aft regions of the stack. The enhanced loading is found to be greater if the exit is closed. A platform at the top, of less than twice the exit diameter, further increases the drag force near the top, but a still larger platform at the top, of about three times the exit diameter, decreases the drag force to values less than those much further below, effectively nullifying the enhanced drag force. It was found that such a reduction of the enhanced drag force in the top regions can also be achieved by a smaller platform of 1.1 to 1.3 times the local diameter, located at about three to five diameters below the top.

Mixed convection heat transfer from the bottom tip of a cylinder spinning about a vertical axis in a saturated porous medium

The present work is a numerical study of heat transfer characteristics from the bottom tip of a cylinder spinning about a vertical axis in an infinitely saturated porous medium. The problem is axisymmetric. The non-dimensionalized governing equations are solved using the SIMPLER algorithm on a staggered grid. The influence of rotational Reynolds numbers and Darcy numbers on the heat transfer for a Grashof number of 104 and Prandtl number of 7.0 is studied. It is found that for very high Darcy numbers, over a wide range of rotational Reynolds numbers, the heat transfer takes place mainly due to conduction. The convective heat transfer takes place for lower Darcy numbers and for higher rotational Reynolds numbers. Moreover, there is a rapid increase in the overall Nusselt number below a certain Darcy number with increase in the rotational Reynolds numbers. The effect of the Darcy number and the rotational Reynolds number on the heat transfer and fluid flow in the porous medium is depicted in the form of streamline and isotherm plots. The variation of the overall Nusselt number with respect to the Darcy number for various rotational Reynolds numbers is plotted. The variation of the local Nusselt number with respect to the radial coordinate at the heated tip of the vertical cylinder is plotted for various Darcy and rotational Reynolds numbers.