On the control of rapidly rotating convection by an axially varying magnetic field

The magnetic field in rapidly rotating dynamos is spatially inhomogeneous. The axial variation of the magnetic field is of particular importance because tall columnar vortices aligned with the rotation axis form at the onset of convection. The classical picture of magnetoconvection with constant or axially varying magnetic fields is that the Rayleigh number and wavenumber at onset decrease appreciably from their non-magnetic values. Nonlinear dynamo simulations show that the axial lengthscale of the self-generated azimuthal magnetic field becomes progressively smaller as we move towards a rapidly rotating regime. With a small-scale field, however, the magnetic control of convection is different from that in previous studies with a uniform or large-scale field. This study looks at the competing viscous and magnetic mode instabilities when the Ekman number E (ratio of viscous to Coriolis forces) is small. As the applied magnetic field strength (measured by the Elsasser number Lambda) increases, the critical Rayleigh number for onset of convection initially increases in a viscous branch, reaches an apex where both viscous and magnetic instabilities co-exist, and then falls in the magnetic branch. The magnetic mode of onset is notable for its dramatic suppression of convection in the bulk of the fluid layer where the field is weak. The viscous-magnetic mode transition occurs at Lambda similar to 1, which implies that small-scale convection can exist at field strengths higher than previously thought. In spherical shell dynamos with basal heating, convection near the tangent cylinder is likely to be in the magnetic mode. The wavenumber of convection is only slightly reduced by the self-generated magnetic field at Lambda similar to 1, in agreement with previous planetary dynamo models. The back reaction of the magnetic field on the flow is, however, visible in the difference in kinetic helicity between cyclonic and anticyclonic vortices.

Local structure of boundary layer transition in experiments with a single streamwise vortex

Transition induced by an isolated streamwise vortex embedded in a flat plate boundary layer was studied experimentally. The vortex was created by a gentle hill with a Gaussian profile that spanned on half of the width of a flat plate mounted in a low turbulence wind tunnel. PIV and hot-wire anemometry data were taken. Transition occurs as a non-inclined shear layer breaks up into a sequence of vortices, close to the boundary layer edge. The passing frequency of these vortices scales with square of the freestream velocity, similar to that in single-roughness induced transition. Quadrant analysis of streamwise and wall-normal velocity fluctuations show large ejection events in the outer layer. (C) 2015 Elsevier Inc. All rights reserved.

Intertwined vorticity and elastodynamics in flapping wing propulsion

We performed numerical experiments on a one-dimensional elastic solid oscillating in a two-dimensional viscous incompressible fluid with the intent of discerning the interplay of vorticity and elastodynamics in flapping wing propulsion. Perhaps for the first time, we have established the role of foil deflection topology and its influence on vorticity generation, through spatially and temporally evolving foil slope and curvature. Though the frequency of oscillation of the foil has a definite role, it is the phase relation between foil slope and pressure that determines thrust or drag. Similarly, the phase difference between flapping velocity, and pressure and inertial forces, determine the power input to the foil, and in turn drives propulsive efficiency. At low frequencies of oscillation, the sympathetic slope and curvature of deformation of the foil allow generation of leading-edge vortices that do not separate; they cause substantial rise in pressure between the leading edge and mid-chord. The circulatory component of pressure is determined primarily by the leading-edge vortex and therefore thrust too is predominantly circulatory in origin at low frequencies. In the intermediate and high-frequency range, thrust and drag on the foil spatially alternate and non-circulatory forces dominate over circulatory and viscous forces. For the mass ratios we simulated, thrust due to flapping varies quadratically as a function of Strouhal number or trailing-edge flapping velocity; further, the trailing edge flapping velocities peak at the same set of frequencies where the thrust is also a maximum. Propulsive efficiency, on the other hand, is roughly a mirror image of the thrust variation with respect to Strouhal number. Given that most instances of flapping propulsion in nature are primarily through distributed muscular actuation that enables precise control of deformation shape, leading to high thrust and efficiency, the results presented here are pointers towards understanding some of the mechanisms that drive thrust and propulsive efficiency.

Instabilities in unsteady boundary layers with reverse flow

Instabilities arising in unsteady boundary layers with reverse flow have been investigated experimentally. Experiments are conducted in a piston driven unsteady water tunnel with a shallow angle diffuser placed in the test section. The ratio of temporal (Pi(t)) to spatial (Pi(x)) component of the pressure gradient can be varied by a controlled motion of the piston. In all the experiments, the piston velocity variation with time is trapezoidal consisting of three phases: constant acceleration from rest, constant velocity and constant deceleration to rest. The adverse pressure gradient (and reverse flow) are due to a combination of spatial deceleration of the free stream in the diffuser and temporal deceleration of the free stream caused by the piston deceleration. The instability is usually initiated with the formation of one or more vortices. The onset of reverse flow in the boundary layer, location and time of formation of the first vortex and the subsequent flow evolution are studied for various values of the ratio Pi(x) (Pi(x) + Pi(t)) for the bottom and the top walls. Instability is due to the inflectional velocity profiles of the unsteady boundary layer. The instability is localized and spreads to the other regions at later times. At higher Reynolds numbers growth rate of instability is higher and localized transition to turbulence is observed. Scalings have been proposed for initial vortex formation time and wavelength of the instability vortices. Initial vortex formation time scales with convective time, delta/Delta U, where S is the boundary layer thickness and Delta U is the difference of maximum and minimum velocities in the boundary layer. Non-dimensional vortex formation time based on convective time scale for the bottom and the top walls are found to be 23 and 30 respectively. Wavelength of instability vortices scales with the time averaged boundary layer thickness. (C) 2015 Elsevier Masson SAS. All rights reserved.

Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.