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.

Stability numbers for an unsupported vertical circular excavation in c-phi soil

By using an axisymmetric lower bound finite element limit analysis formulation, the stability numbers (gamma H/C) for an unsupported vertical circular excavation in a cohesive-frictional soil have been generated. The numerical results are obtained for values of normalized excavation height (H/b) and friction angle (phi) greater than those considered previously in the literature. The results compare well with those available in literature. The stability numbers presented in this note would be beneficial from a design point of view. (C) 2011 Elsevier Ltd. All rights reserved.

Numerical analysis of transpiration cooling in carbon dioxide atmosphere at hypersonic Mach numbers

The numerical solutions are obtained for skin friction, heat transfer to the wall and growth of boundary layer along the flat plate by employing two dimensional Navier-Stokes equations governing the hypersonic flow coupled with species continuity equations. Flow fields have been computed along the flat plate in CO2 atmosphere in the presence of transpiration cooling using air and carbon dioxide.

Wave propagation in single-walled carbon nanotube under longitudinal magnetic field using nonlocal Euler-Bernoulli beam theory

In the present work, the effect of longitudinal magnetic field on wave dispersion characteristics of equivalent continuum structure (ECS) of single-walled carbon nanotubes (SWCNT) embedded in elastic medium is studied. The ECS is modelled as an Euler-Bernoulli beam. The chemical bonds between a SWCNT and the elastic medium are assumed to be formed. The elastic matrix is described by Pasternak foundation model, which accounts for both normal pressure and the transverse shear deformation. The governing equations of motion for the ECS of SWCNT under a longitudinal magnetic field are derived by considering the Lorentz magnetic force obtained from Maxwell's relations within the frame work of nonlocal elasticity theory. The wave propagation analysis is performed using spectral analysis. The results obtained show that the velocity of flexural waves in SWCNTs increases with the increase of longitudinal magnetic field exerted on it in the frequency range: 0-20 THz. The present analysis also shows that the flexural wave dispersion in the ECS of SWCNT obtained by local and nonlocal elasticity theories differ. It is found that the nonlocality reduces the wave velocity irrespective of the presence of the magnetic field and does not influences it in the higher frequency region. Further it is found that the presence of elastic matrix introduces the frequency band gap in flexural wave mode. The band gap in the flexural wave is found to independent of strength of the longitudinal magnetic field. (C) 2011 Elsevier Inc. All rights reserved.

Badly approximable numbers and vectors in Cantor-like sets

We show that a large class of Cantor-like sets of R-d, d >= 1, contains uncountably many badly approximable numbers, respectively badly approximable vectors, when d >= 2. An analogous result is also proved for subsets of R-d arising in the study of geodesic flows corresponding to (d+1)-dimensional manifolds of constant negative curvature and finite volume, generalizing the set of badly approximable numbers in R. Furthermore, we describe a condition on sets, which is fulfilled by a large class, ensuring a large intersection with these Cantor-like sets.

Coupled polynomial field approach for elimination of flexure and torsion locking phenomena in the Timoshenko and Euler-Bernoulli curved beam elements

The curvature related locking phenomena in the out-of-plane deformation of Timoshenko and Euler-Bernoulli curved beam elements are demonstrated and a novel approach is proposed to circumvent them. Both flexure and Torsion locking phenomena are noticed in Timoshenko beam and torsion locking phenomenon alone in Euler-Bernoulli beam. Two locking-free curved beam finite element models are developed using coupled polynomial displacement field interpolations to eliminate these locking effects. The coupled polynomial interpolation fields are derived independently for Timoshenko and Euler-Bernoulli beam elements using the governing equations. The presented of penalty terms in the couple displacement fields incorporates the flexure-torsion coupling and flexure-shear coupling effects in an approximate manner and produce no spurious constraints in the extreme geometric limits of flexure, torsion and shear stiffness. the proposed couple polynomial finite element models, as special cases, reduce to the conventional Timoshenko beam element and Euler-Bernoulli beam element, respectively. These models are shown to perform consistently over a wide range of flexure-to-shear (EI/GA) and flexure-to-torsion (EI/GJ) stiffness ratios and are inherently devoid of flexure, torsion and shear locking phenomena. The efficacy, accuracy and reliability of the proposed models to straight and curved beam applications are demonstrated through numerical examples. (C) 2012 Elsevier B.V. All rights reserved.

Analogy Between Rotating Euler-Bernoulli and Timoshenko Beams and Stiff Strings

The governing differential equation of a rotating beam becomes the stiff-string equation if we assume uniform tension. We find the tension in the stiff string which yields the same frequency as a rotating cantilever beam with a prescribed rotating speed and identical uniform mass and stiffness. This tension varies for different modes and are found by solving a transcendental equation using bisection method. We also find the location along the rotating beam where equivalent constant tension for the stiff string acts for a given mode. Both Euler-Bernoulli and Timoshenko beams are considered for numerical results. The results provide physical insight into relation between rotating beams and stiff string which are useful for creating basis functions for approximate methods in vibration analysis of rotating beams.

Tribological response of soft materials sliding against hard surface textures at various numbers of cycles

In the present investigation, various kinds of textures were attained on the steel surfaces. Roughness of the textures was varied using different grits of emery papers or polishing powders. Pins made of pure Al, Al-4Mg alloy and pure Mg were then slid against prepared steel plate surfaces at various numbers of cycles using an inclined pin-on-plate sliding tester. Tests were conducted at a sliding velocity of 2mms(-1) in ambient conditions under both dry and lubricated conditions. Normal loads were increased up to 110N during the tests. The morphologies of the worn surfaces of the pins and the formation of transfer layer on the counter surfaces were observed using a scanning electron microscope. Surface roughness parameters of the plate were measured using an optical profilometer. In the experiments, it was observed that the coefficient of friction and formation of a transfer layer (under dry and lubricated conditions) only depended on surface texture during the first few sliding cycles. The steady-state variation in the coefficient of friction under both dry and lubrication conditions was attributed to the self-organisation of texture of the surfaces at the interface during sliding. Copyright (C) 2012 John Wiley & Sons, Ltd.

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.

Closed-form solutions for non-uniform Euler-Bernoulli free-free beams

In this paper, the free vibration of a non-uniform free-free Euler-Bernoulli beam is studied using an inverse problem approach. It is found that the fourth-order governing differential equation for such beams possess a fundamental closed-form solution for certain polynomial variations of the mass and stiffness. An infinite number of non-uniform free-free beams exist, with different mass and stiffness variations, but sharing the same fundamental frequency. A detailed study is conducted for linear, quadratic and cubic variations of mass, and on how to pre-select the internal nodes such that the closed-form solutions exist for the three cases. A special case is also considered where, at the internal nodes, external elastic constraints are present. The derived results are provided as benchmark solutions for the validation of non-uniform free-free beam numerical codes. (C) 2013 Elsevier Ltd. All rights reserved.

Modeling harvest rates and numbers from age and sex ratios: A demonstration for elephant populations

Illegal harvest rates of wildlife populations are often unknown or difficult to estimate from field data due to under-reporting or incomplete detection of carcasses. This is especially true for elephants that are killed for ivory or in conflicts with people. We describe a method to infer harvest rates from coarse field data of three population parameters, namely, adult female to male ratio, male old-adult to young-adult ratio, and proportion of adult males in the population using Jensen's (2000) 2-sex, density-dependent Leslie matrix model. The specific combination of male and female harvest rates and numbers can be determined from the history of harvest and estimate of population size. We applied this technique to two populations of elephants for which data on age structure and records of mortality were available-a forest-dwelling population of the Asian elephant (at Nagarahole, India) and an African savannah elephant population (at Samburu, Kenya) that had experienced male-biased harvest regimes over 2-3 decades. For the Nagarahole population, the recorded numbers of male and female elephants killed illegally during 1981-2000 were 64% and 88% of the values predicted by the model, respectively, implying some non-detection or incomplete reporting while for the Samburu population the recorded and modeled numbers of harvest during 1990-1999 closely matched. This technique, applicable to any animal population following logistic growth model, can be especially useful for inferring illegal harvest numbers of forest elephants in Africa and Asia.

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.

Globally coupled stochastic two-state oscillators: Fluctuations due to finite numbers

Infinite arrays of coupled two-state stochastic oscillators exhibit well-defined steady states. We study the fluctuations that occur when the number N of oscillators in the array is finite. We choose a particular form of global coupling that in the infinite array leads to a pitchfork bifurcation from a monostable to a bistable steady state, the latter with two equally probable stationary states. The control parameter for this bifurcation is the coupling strength. In finite arrays these states become metastable: The fluctuations lead to distributions around the most probable states, with one maximum in the monostable regime and two maxima in the bistable regime. In the latter regime, the fluctuations lead to transitions between the two peak regions of the distribution. Also, we find that the fluctuations break the symmetry in the bimodal regime, that is, one metastable state becomes more probable than the other, increasingly so with increasing array size. To arrive at these results, we start from microscopic dynamical evolution equations from which we derive a Langevin equation that exhibits an interesting multiplicative noise structure. We also present a master equation description of the dynamics. Both of these equations lead to the same Fokker-Planck equation, the master equation via a 1/N expansion and the Langevin equation via standard methods of Ito calculus for multiplicative noise. From the Fokker-Planck equation we obtain an effective potential that reflects the transition from the monomodal to the bimodal distribution as a function of a control parameter. We present a variety of numerical and analytic results that illustrate the strong effects of the fluctuations. We also show that the limits N -> infinity and t -> infinity(t is the time) do not commute. In fact, the two orders of implementation lead to drastically different results.

Effect of Rayleigh numbers on the evolution of double-diffusive salt fingers

This is a transient two-dimensional numerical study of double-diffusive salt fingers in a two-layer heat-salt system for a wide range of initial density stability ratio (R-rho 0) and thermal Rayleigh numbers (Ra-T similar to 10(3) - 10(11)). Salt fingers have been studied for several decades now, but several perplexing features of this rich and complex system remain unexplained. The work in question studies this problem and shows the morphological variation in fingers from low to high thermal Rayleigh numbers, which have been missed by the previous investigators. Considerable variations in convective structures and evolution pattern were observed in the range of Ra-T used in the simulation. Evolution of salt fingers was studied by monitoring the finger structures, kinetic energy, vertical profiles, velocity fields, and transient variation of R-rho(t). The results show that large scale convection that limits the finger length was observed only at high Rayleigh numbers. The transition from nonlinear to linear convection occurs at about Ra-T similar to 10(8). Contrary to the popular notion, R-rho(t) first decrease during diffusion before the onset time and then increase when convection begins at the interface. Decrease in R-rho(t) is substantial at low Ra-T and it decreases even below unity resulting in overturning of the system. Interestingly, all the finger system passes through the same state before the onset of convection irrespective of Rayleigh number and density stability ratio of the system. (C) 2014 AIP Publishing LLC.