Stability of the flow of a fluid through a flexible tube at high Reynolds number

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (ρVR / η), the ratio of the viscosities of the wall and fluid ηr = (ηs/η), the ratio of radii H and the dimensionless velocity Γ = (ρV2/G)1/2. Here ρ is the density of the fluid, G is the coefficient of elasticity of the wall and Vis the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter ε = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate do), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctruations due to the Reynolds strees. There is an O(ε1/2) correction to the growth rate, s(1), due to the presence of a wall layer of thickness ε1/2R where the viscous stresses are O(ε1/2) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Γ and wavenumber k where s(l) = 0. At these points, the wail layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(ε) correction to the growth rate s(2) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s(2) increases [is proportional to] (H − 1)−2 for (H − 1) [double less-than sign] 1 (thickness of wall much less than the tube radius), and decreases [is proportional to] (H−4 for H [dbl greater-than sign] 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube.

On cavity flow at high Reynolds numbers

The flow in a square cavity is studied by solving the full Navier–Stokes and energy equations numerically, employing finite-difference techniques. Solutions are obtained over a wide range of Reynolds numbers from 0 to 50000. The solutions show that only at very high Reynolds numbers (Re [gt-or-equal, slanted] 30000) does the flow in the cavity completely correspond to that assumed by Batchelor's model for separated flows. The flow and thermal fields at such high Reynolds numbers clearly exhibit a boundary-layer character. For the first time, it is demonstrated that the downstream secondary eddy grows and decays in a manner similar to the upstream one. The upstream and downstream secondary eddies remain completely viscous throughout the range of Reynolds numbers of their existence. It is suggested that the behaviour of the secondary eddies may be characteristic of internal separated flows.

Defining the ‘modified Griffin plot’ in vortex-induced vibration: revealing the effect of Reynolds number using controlled damping

In the present work, we study the transverse vortex-induced vibrations of an elastically mounted rigid cylinder in a fluid flow. We employ a technique to accurately control the structural damping, enabling the system to take on both negative and positive damping. This permits a systematic study of the effects of system mass and damping on the peak vibration response. Previous experiments over the last 30 years indicate a large scatter in peak-amplitude data ($A^*$) versus the product of mass–damping ($\alpha$), in the so-called ‘Griffin plot’. A principal result in the present work is the discovery that the data collapse very well if one takes into account the effect of Reynolds number ($\mbox{\textit{Re}}$), as an extra parameter in a modified Griffin plot. Peak amplitudes corresponding to zero damping ($A^*_{{\alpha}{=}0}$), for a compilation of experiments over a wide range of $\mbox{\textit{Re}}\,{=}\,500-33000$, are very well represented by the functional form $A^*_{\alpha{=}0} \,{=}\, f(\mbox{\textit{Re}}) \,{=}\, \log(0.41\,\mbox{\textit{Re}}^{0.36}$). For a given $\mbox{\textit{Re}}$, the amplitude $A^*$ appears to be proportional to a function of mass–damping, $A^*\propto g(\alpha)$, which is a similar function over all $\mbox{\textit{Re}}$. A good best-fit for a wide range of mass–damping and Reynolds number is thus given by the following simple expression, where $A^*\,{=}\, g(\alpha)\,f(\mbox{\textit{Re}})$: \[ A^* \,{=}\,(1 - 1.12\,\alpha + 0.30\,\alpha^2)\,\log (0.41\,\mbox{\textit{Re}}^{0.36}). \] In essence, by using a renormalized parameter, which we define as the ‘modified amplitude’, $A^*_M\,{=}\,A^*/A^*_{\alpha{=}0}$, the previously scattered data collapse very well onto a single curve, $g(\alpha)$, on what we refer to as the ‘modified Griffin plot’. There has also been much debate over the last three decades concerning the validity of using the product of mass and damping (such as $\alpha$) in these problems. Our results indicate that the combined mass–damping parameter ($\alpha$) does indeed collapse peak-amplitude data well, at a given $\mbox{\textit{Re}}$, independent of the precise mass and damping values, for mass ratios down to $m^*\,{=}\,1$.

Fluid Mechanics of Aquatic Locomotion at Large Reynolds Numbers

Abstract | There exist a huge range of fish species besides other aquatic organisms like squids and salps that locomote in water at large Reynolds numbers, a regime of flow where inertial forces dominate viscous forces. In the present review, we discuss the fluid mechanics governing the locomotion of such organisms. Most fishes propel themselves by periodic undulatory motions of the body and tail, and the typical classification of their swimming modes is based on the fraction of their body that undergoes such undulatory motions. In the angulliform mode, or the eel type, the entire body undergoes undulatory motions in the form of a travelling wave that goes from head to tail, while in the other extreme case, the thunniform mode, only the rear tail (caudal fin) undergoes lateral oscillations. The thunniform mode of swimming is essentially based on the lift force generated by the airfoil like crosssection of the fish tail as it moves laterally through the water, while the anguilliform mode may be understood using the “reactive theory” of Lighthill. In pulsed jet propulsion, adopted by squids and salps, there are two components to the thrust; the first due to the familiar ejection of momentum and the other due to an over-pressure at the exit plane caused by the unsteadiness of the jet. The flow immediately downstream of the body in all three modes consists of vortex rings; the differentiating point being the vastly different orientations of the vortex rings. However, since all the bodies are self-propelling, the thrust force must be equal to the drag force (at steady speed), implying no net force on the body, and hence the wake or flow downstream must be momentumless. For such bodies, where there is no net force, it is difficult to directly define a propulsion efficiency, although it is possible to use some other very different measures like “cost of transportation” to broadly judge performance.

Performance of cellular CDMA with voice/data traffic with an SIR based admission control

We analyze the performance of an SIR based admission control strategy in cellular CDMA systems with both voice and data traffic. Most studies In the current literature to estimate CDMA system capacity with both voice and data traf-Bc do not take signal-tlFlnterference ratio (SIR) based admission control into account In this paper, we present an analytical approach to evaluate the outage probability for voice trafllc, the average system throughput and the mean delay for data traffic for a volce/data CDMA system which employs an SIR based admission controL We show that for a dataaniy system, an improvement of about 25% In both the Erlang capacity as well as the mean delay performance is achieved with an SIR based admission control as compared to code availability based admission control. For a mixed voice/data srtem with 10 Erlangs of voice traffic, the Lmprovement in the mean delay performance for data Is about 40%.Ah, for a mean delay of 50 ms with 10 Erlangs voice traffic, the data Erlang capacity improves by about 9%.

Breakup and coalescence characteristics of a hollow cone swirling spray

This paper deals with an experimental study of the breakup characteristics of water emanating from hollow cone hydraulic injector nozzles induced by pressure-swirling. The experiments were conducted using two nozzles with different orifice diameters 0.3 mm and 0.5 mm and injection pressures (0.3-4 MPa) which correspond to Rep = 7000-26 000. Two types of laser diagnostic techniques were utilized: shadowgraph and phase Doppler particle anemometry for a complete study of the atomization process. Measurements that were made in the spray in both axial and radial directions indicate that both velocity and average droplet diameter profiles are highly dependent on the nozzle characteristics, Weber number and Reynolds number. The spatial variation of diameter and velocity arises principally due to primary breakup of liquid films and subsequent secondary breakup of large droplets due to aerodynamic shear. Downstream of the nozzle, coalescence of droplets due to collision was also found to be significant. Different types of liquid film breakup were considered and found to match well with the theory. Secondary breakup due to shear was also studied theoretically and compared to the experimental data. Coalescence probability at different axial and radial locations was computed to explain the experimental results. The spray is subdivided into three zones: near the nozzle, a zone consisting of film and ligament regime, where primary breakup and some secondary breakup take place; a second zone where the secondary breakup process continues, but weakens, and the centrifugal dispersion becomes dominant; and a third zone away from the spray where coalescence is dominant. Each regime has been analyzed in detail, characterized by timescale and Weber number and validated using experimental data. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4773065]

Experiments on thermal response of low aspect ratio packed beds at high Reynolds numbers with varying inflow temperature

An experimental study for transient temperature response of low aspect ratio packed beds at high Reynolds numbers for a free stream with varying inlet temperature is presented. The packed bed is used as a compact heat exchanger along with a solid propellant gas-generator, to generate room temperature gases for use in applications such as control actuation and air bottle pressurization. Packed beds of lengths similar to 200 mm and 300 mm were characterized for packing diameter based Reynolds numbers, Re-d ranging from 0.6 x 10(4) to 8.5 x 10(4). The solid packing used in the bed consisted of circle divide 9.5 mm and circle divide 5 mm steel spheres with suitable arrangements to eliminate flow entrance and exit effects. The ratios of packed bed diameter to packing diameter for 9.5 mm and 5 mm sphere packing were similar to 9.5 and 18 respectively, with the average packed bed porosities around 0.4. Gas temperatures were measured at the entry, exit and at three axial locations along centre-line in the packed beds. The solid packing temperature was measured at three axial locations in the packed bed. An average Nusselt number correlation of the form Nu(d) = 3.91Re(d)(05) for Re-d range of 10(4) is proposed. For engineering applications of packed beds such as pebble bed heaters, thermal storage systems, and compact heat exchangers a simple procedure is suggested for calculating unsteady gas temperature at packed bed exit for packing Biot number Bi-d < 0.1. (C) 2012 Elsevier Inc. All rights reserved.

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.

Stall flutter of NACA 0012 airfoil at low Reynolds numbers

In the present work, we experimentally study and demarcate the stall flutter boundaries of a NACA 0012 airfoil at low Reynolds numbers (Re similar to 10(4)) by measuring the forces and flow fields around the airfoil when it is forced to oscillate. The airfoil is placed at large mean angle of attack (alpha(m)), and is forced to undergo small amplitude pitch oscillations, the amplitude (Delta alpha) and frequency (f) of which are systematically varied. The unsteady loads on the oscillating airfoil are directly measured, and are used to calculate the energy transfer to the airfoil from the flow. These measurements indicate that for large mean angles of attack of the airfoil (alpha(m)), there is positive energy transfer to the airfoil over a range of reduced frequencies (k=pi fc/U), indicating that there is a possibility of airfoil excitation or stall flutter even at these low Re (c=chord length). Outside this range of reduced frequencies, the energy transfer is negative and under these conditions the oscillations would be damped. Particle Image Velocimetry (PIV) measurements of the flow around the oscillating airfoil show that the shear layer separates from the leading edge and forms a leading edge vortex, although it is not very clear and distinct due to the low oscillation amplitudes. On the other hand, the shear layer formed after separation is found to clearly move periodically away from the airfoil suction surface and towards it with a phase lag to the airfoil oscillations. The phase of the shear layer motion with respect to the airfoil motions shows a clear difference between the exciting and the damping case.

Effect of steady rotation on low Reynolds number vortex shedding behind a circular cylinder

In this paper control of oblique vortex shedding in the wake behind a straight circular cylinder is explored experimentally and computationally. Towards this, steady rotation of the cylinder about its axis is used as a control device. Some limited studies are also performed with a stepped circular cylinder, where at the step the flow is inevitably three-dimensional irrespective of the rotation rate. When there is no rotation, the vortex shedding pattern is three dimensional as described in many previous studies. With a non-zero rotation rate, it is demonstrated experimentally as well as numerically that the shedding pattern becomes more and more two-dimensional. At sufficiently high rotation rates, the vortex shedding is completely suppressed.

COMPARISON AND CROSS-VALIDATION OF OPTICAL TECHNIQUES IN DIFFERENT SWIRL SPRAY REGIMES

This paper deals with an experimental study of pressure-swirl hydraulic injector nozzles using non-intrusive optical techniques. Experiments were conducted to study atomization characteristics using two nozzles with different orifice diameters, 0.3 mm and 0.5 mm, and injection pressures, 0.3-3.5 Mpa, which correspond to Reynolds number (Re-p) = 7,000-45,000, depending on nozzle utilized. Three laser diagnostic techniques were utilized: Shadowgraph, PIV (Particle Image Velocimetry), and PDPA (Phase Doppler Particle Anemometry). Measurements made in the spray in both axial and radial directions indicate that velocity, average droplet diameter profiles, and spray dynamics are highly dependent on the nozzle characteristics and injection pressure. Limitations of these techniques in the different flow regimes, related to the primary and secondary breakups as well as coalescence, are provided. Results indicate that all three techniques provide similar results throughout the different regimes. Shadowgraph and PDPA were possible in the secondary atomization and coalescence regimes while PIV measurements could be made only at the end of secondary atomization and coalescence.

How to run a campaign: Optimal control of SIS and SIR information epidemics

Information spreading in a population can be modeled as an epidemic. Campaigners (e.g., election campaign managers, companies marketing products or movies) are interested in spreading a message by a given deadline, using limited resources. In this paper, we formulate the above situation as an optimal control problem and the solution (using Pontryagin's Maximum Principle) prescribes an optimal resource allocation over the time of the campaign. We consider two different scenarios-in the first, the campaigner can adjust a direct control (over time) which allows her to recruit individuals from the population (at some cost) to act as spreaders for the Susceptible-Infected-Susceptible (SIS) epidemic model. In the second case, we allow the campaigner to adjust the effective spreading rate by incentivizing the infected in the Susceptible-Infected-Recovered (SIR) model, in addition to the direct recruitment. We consider time varying information spreading rate in our formulation to model the changing interest level of individuals in the campaign, as the deadline is reached. In both the cases, we show the existence of a solution and its uniqueness for sufficiently small campaign deadlines. For the fixed spreading rate, we show the effectiveness of the optimal control strategy against the constant control strategy, a heuristic control strategy and no control. We show the sensitivity of the optimal control to the spreading rate profile when it is time varying. (C) 2014 Elsevier Inc. All rights reserved.

Approaches toward (eta(2)-HX) (X = H, SiR (R = Me-3 or Me2Ph), CH3) sigma-complexes of ruthenium

Two new Ru(II)-complexes RuH(Tpms)(PPh3)(2)] 1 (Tpms - (C3H3N2)(3)CSO3, tris-(pyrazolyl) methane sulfonate) and Ru(OTf)(Tpms)(PPh3)(2)] 2 (OTf = CF3SO3) have been synthesized and characterized wherein Ru-H and Ru-OTf are the key reactive centers. Reaction of 1 with HOTf results in the Ru(eta(2)-H-2)(Tpms)(PPh3)(2)]OTf] complex 3, whereas reaction of 1 with Me3SiOTf affords the dihydrogen complex 3 and complex 1 through an unobserved sigma-silane intermediate. In addition, an attempt to characterize the sigma methane complex via reaction of complex 1 with CH3OTf yields complex 2 and free methane. On the other hand, reaction of Ru(OTf)(Tpms)(PPh3)(2)] 2 with H-2 and PhMe2SiH at low temperature resulted in sigma-H-2, 3 and a probable sigma-silane complexes, respectively. However, no sigma-methane complex was observed for the reaction of complex 2 with methane even at low temperature. (C) 2014 Elsevier B. V. All rights reserved.