Flexibility in flapping foil suppresses meandering of induced jet in absence of free stream

Thrust-generating flapping foils are known to produce jets inclined to the free stream at high Strouhal numbers St = fA/U-infinity, where f is the frequency and A is the amplitude of flapping and U-infinity is the free-stream velocity. Our experiments, in the limiting case of St —> infinity (zero free-stream speed), show that a purely oscillatory pitching motion of a chordwise flexible foil produces a coherent jet composed of a reverse Benard-Karman vortex street along the centreline, albeit over a specific range of effective flap stiffnesses. We obtain flexibility by attaching a thin flap to the trailing edge of a rigid NACA0015 foil; length of flap is 0.79 c where c is rigid foil chord length. It is the time-varying deflections of the flexible flap that suppress the meandering found in the jets produced by a pitching rigid foil for zero free-stream condition. Recent experiments (Marais et al., J. Fluid Mech., vol. 710, 2012, p. 659) have also shown that the flexibility increases the St at which non-deflected jets are obtained. Analysing the near-wake vortex dynamics from flow visualization and particle image velocimetry (PIV) measurements, we identify the mechanisms by which flexibility suppresses jet deflection and meandering. A convenient characterization of flap deformation, caused by fluid-flap interaction, is through a non-dimensional effective stiffness', EI* = 8 EI/(rho V-TEmax(2) s(f) c(f)(3)/2), representing the inverse of the flap deflection due to the fluid-dynamic loading; here, EI is the bending stiffness of flap, rho is fluid density, V-TEmax is the maximum velocity of rigid foil trailing edge, s(f) is span and c(f) is chord length of the flexible flap. By varying the amplitude and frequency of pitching, we obtain a variation in EI* over nearly two orders of magnitude and show that only moderate EI*. (0.1 less than or similar to EI * less than or similar to 1 generates a sustained, coherent, orderly jet. Relatively `stiff' flaps (EI* greater than or similar to 1), including the extreme case of no flap, produce meandering jets, whereas highly `flexible' flaps (EI* less than or similar to 0.1) produce spread-out jets. Obtained from the measured mean velocity fields, we present values of thrust coefficients for the cases for which orderly jets are observed.

The role of buoyancy in polarity reversals of the geodynamo

We investigate polarity reversals in the geodynamo using a rotating, convection-driven dynamo model. As the flow in rapidly rotating convection is dominated by columns aligned with the axis of rotation, the focus is on the dynamics of columnar vortices. By studying the growth of a seed magnetic field to a stable axial dipole field, we show that the magnetic field acts in ways that significantly enhance the relative helicity between cyclonic and anticyclonic vortices. This flow asymmetry is the hallmark of a dipolar dynamo. Strong buoyancy, on the other hand, offsets the effect of the magnetic field, establishing parity between positive and negative vortices. As the dipole field is deprived of the helicity required to support itself, the dynamo is pushed into a reversing state. This is a likely regime for polarity reversals in the Earth's core. The integral lengthscale at which buoyancy injects energy is not significantly different from the convective flow lengthscale, which implies that buoyancy does not feed vortices at the small scales where non-linear inertia is present. The lengthscale at which the Lorentz force acts in the reversing dynamo is small, which may allow the passive presence of non-linear inertia in the small scales.

Influence of melt pouring temperature and plate inclination on solidification and microstructure of A356 aluminum alloy produced using oblique plate

The preparation of semisolid slurry of A356 aluminum alloy using an oblique plate was investigated. A356 alloy melt undergoes partial solidification when it flows down on an oblique plate cooled from underneath by counter flowing water. It results in continuous formation of columnar dendrites on plate wall. Due to forced convection, these dendrites are sheared off into equiaxed/fragmented grains and then washed away continuously to produce semisolid slurry at plate exit. Melt pouring temperature provides required condition of solidification whereas plate inclination enables necessary shear for producing semisolid slurry of desired quality. Slurry obtained was solidified in metal mould to produce semisolid-cast billets of desired microstructure. Furthermore, semisolid-cast billets were heat treated to improve surface quality. Microstructures of both semisolid-cast and heat-treated billets were analyzed. Effects of melt pouring temperature and plate inclination on solidification and microstructure of billets produced using oblique plate were described. The investigations involved four different melt pouring temperatures (620, 625, 630 and 635 degrees C) associated with four different plate inclinations (30 degrees, 45 degrees, 60 degrees and 75 degrees). Melt pouring temperature of 625 degrees C with plate inclination of 60 degrees shows fine and globular microstructures and it is the optimum.

Flux Transport Dynamos: From Kinematics to Dynamics

Over the past several decades, Flux-Transport Dynamo (FTD) models have emerged as a popular paradigm for explaining the cyclic nature of solar magnetic activity. Their defining characteristic is the key role played by the mean meridional circulation in transporting magnetic flux and thereby regulating the cycle period. Most FTD models also incorporate the so-called Babcock-Leighton (BL) mechanism in which the mean poloidal field is produced by the emergence and subsequent dispersal of bipolar active regions. This feature is well grounded in solar observations and provides a means for assimilating observed surface flows and fields into the models in order to forecast future solar activity, to identify model biases, and to clarify the underlying physical processes. Furthermore, interpreting historical sunspot records within the context of FTD models can potentially provide insight into why cycle features such as amplitude and duration vary and what causes extreme events such as Grand Minima. Though they are generally robust in a modeling sense and make good contact with observed cycle features, FTD models rely on input physics that is only partially constrained by observation and that neglects the subtleties of convective transport, convective field generation, and nonlinear feedbacks. Here we review the formulation and application of FTD models and assess our current understanding of the input physics based largely on complementary 3D MHD simulations of solar convection, dynamo action, and flux emergence.

Agglomeration front dynamics: Drying in sessile nano-particle laden droplets

The drying of sessile, nano-silica laden water droplet is studied under ambient conditions, in the absence of any convection. The drying process can be divided into two distinct regimes. During regime 1, the outer edge of the droplet remains pinned and particles agglomerate at the droplet periphery similar to the traditional coffee ring. However in regime 2, with further evaporation, both the liquid contact line and the agglomeration front starts moving radially inwards from the initial contact edge. The contact between the liquid and the agglomerate is maintained throughout regime 2 and the vaporisation driven liquid edge recession essentially drives the inward growth of the particle deposition. Fast kinetics of particle aggregation results in rapid growth of this agglomeration front as seen from the experiments. A theoretical formulation involving a simplistic model of the agglomeration front growth based on particle mass balance has been proposed. (C) 2014 Elsevier Ltd. All rights reserved,

Cross-correlation-aided transport in stochastically driven accretion flows

The origin of linear instability resulting in rotating sheared accretion flows has remained a controversial subject for a long time. While some explanations of such non-normal transient growth of disturbances in the Rayleigh stable limit were available for magnetized accretion flows, similar instabilities in the absence of magnetic perturbations remained unexplained. This dichotomy was resolved in two recent publications by Chattopadhyay and co-workers Mukhopadhyay and Chattopadhyay, J. Phys. A 46, 035501 (2013); Nath et al., Phys. Rev. E 88, 013010 (2013)] where it was shown that such instabilities, especially for nonmagnetized accretion flows, were introduced through interaction of the inherent stochastic noise in the system (even a ``cold'' accretion flow at 3000Kis too ``hot'' in the statistical parlance and is capable of inducing strong thermal modes) with the underlying Taylor-Couette flow profiles. Both studies, however, excluded the additional energy influx (or efflux) that could result from nonzero cross correlation of a noise perturbing the velocity flow, say, with the noise that is driving the vorticity flow (or equivalently the magnetic field and magnetic vorticity flow dynamics). Through the introduction of such a time symmetry violating effect, in this article we show that nonzero noise cross correlations essentially renormalize the strength of temporal correlations. Apart from an overall boost in the energy rate (both for spatial and temporal correlations, and hence in the ensemble averaged energy spectra), this results in mutual competition in growth rates of affected variables often resulting in suppression of oscillating Alfven waves at small times while leading to faster saturations at relatively longer time scales. The effects are seen to be more pronounced with magnetic field fluxes where the noise cross correlation magnifies the strength of the field concerned. Another remarkable feature noted specifically for the autocorrelation functions is the removal of energy degeneracy in the temporal profiles of fast growing non-normal modes leading to faster saturation with minimum oscillations. These results, including those presented in the previous two publications, now convincingly explain subcritical transition to turbulence in the linear limit for all possible situations that could now serve as the benchmark for nonlinear stability studies in Keplerian accretion disks.

A robust layer adapted difference method for singularly perturbed two-parameter parabolic problems

A finite difference method for a time-dependent singularly perturbed convection-diffusion-reaction problem involving two small parameters in one space dimension is considered. We use the classical implicit Euler method for time discretization and upwind scheme on the Shishkin-Bakhvalov mesh for spatial discretization. The method is analysed for convergence and is shown to be uniform with respect to both the perturbation parameters. The use of the Shishkin-Bakhvalov mesh gives first-order convergence unlike the Shishkin mesh where convergence is deteriorated due to the presence of a logarithmic factor. Numerical results are presented to validate the theoretical estimates obtained.

Influence of plate length and plate cooling rate on solidification and microstructure of A356 alloy produced by oblique plate

A356 alloy melt solidifies partially when it flows down on an oblique plate cooled from bottom by counter flowing water. Columnar dendrites are continuously formed on the plate wall. Because of the forced convection, these dendrites are sheared off into equiaxed/fragmented grains and then washed away continuously by producing semisolid slurry at plate exit. Plate cooling rate provides required extent/amount of solidification whereas plate length enables necessary shear for producing semisolid slurry of desired quality. Slurry obtained is solidified in metal mould to produce semisolid-cast billets of desired microstructure. Furthermore, semisolid-cast billets are also heat-treated to improve surface quality. Microstructures of both semisolid-cast and heat-treated billets are compared. The effects of plate length and plate cooling rate on solidification and microstructure of billets produced by using oblique plate are illustrated. Three different plate lengths (200 mm, 250 mm, 300 mm) associated with three different heat transfer coefficients (1000, 2000 and 2500 W/(m(2).K)) are involved. Plate length of 250 mm with heat transfer coefficient of 2000 W/(m(2).K) gives fine and globular microstructures and is the optimum as there is absolutely no possibility of sticking of slurry to plate wall.

Improved design methodology for the development of electrically actuated MEMS structures

A comprehensive design flow is proposed for the design of Micro Electro Mechanical Systems that are fabricated using SOIMUMPs process. Many of the designers typically do not model the temperature dependency of electrical conductivity, thermal conductivity and convection coefficient, as it is very cumbersome to create/incorporate the same in the existing FEM simulators. Capturing these dependencies is very critical particularly for structures that are electrically actuated. Lookup tables that capture the temperature dependency of electrical conductivity, thermal conductivity and convection coefficient are created. These look up tables are taken as inputs for a commercially available FEM simulator to model the semiconductor behavior. It is demonstrated that when temperature dependency for all the above mentioned parameters is not captured, then the error in estimation of the maximum temperature (for a given structure) could be as high as 30%. Error in the estimated resistance value under the same conditions is as high as 40%. When temperature dependency of the above mentioned parameters is considered then error w.r.t the measured values is less than 5%. It is evident that error in temperature estimates leads to erroneous results from mechanical simulations as well.

First hyperpolarizability of Ru-half-sandwich complexes: The effect of halogen atom substitution on the ancillary ligand

The first hyperpolarizability (beta) of a series of half-sandwich Ru complexes with a mercaptobenzothiazole ligand bearing a halogen atom substitution in the para-position has been investigated by hyper-Rayleigh scattering and quantum chemical calculations. The heterocyclic ligand with a bromine atom in the para position makes it a very good donor and charge flows to the Ru center enhancing the beta value of the complex by a factor of 2 compared to the complex with the ligand without the halogen substitution. The resonance (+R) and the inductive (-I) effects exerted by the halogen atom in the para position push electrons in opposing directions in the complex. For the Br and Cl atoms the resonance effect dominates which enables the ligand to donate electrons to the metal center thereby increasing the hyperpolarizability whereas for the fluorine atom, the inductive effect is dominant which reduces the charge flow to the metal and the hyperpolarizability drops even below that of the unsubstituted ligand. This unprecedented halogen atom effect on beta of metal complexes is reported. (C) 2015 Elsevier By. All rights reserved.

Interaction of a vortex ring with a single bubble: bubble and vorticity dynamics

The interaction of a single bubble with a single vortex ring in water has been studied experimentally. Measurements of both the bubble dynamics and vorticity dynamics have been done to help understand the two-way coupled problem. The circulation strength of the vortex ring (Gamma) has been systematically varied, while keeping the bubble diameter (D-b) constant, with the bubble volume to vortex core volume ratio (V-R) also kept fixed at about 0.1. The other important parameter in the problem is a Weber number based on the vortex ring strength. (We = 0.87 rho(Gamma/2 pi a)(2)/(sigma/D-b); a = vortex core radius, sigma = surface tension), which is varied over a large range, We = 3-406. The interaction between the bubble and ring for each of the We cases broadly falls into four stages. Stage I is before capture of the bubble by the ring where the bubble is drawn into the low-pressure vortex core, while in stage II the bubble is stretched in the azimuthal direction within the ring and gradually broken up into a number of smaller bubbles. Following this, in stage III the bubble break-up is complete and the resulting smaller bubbles slowly move around the core, and finally in stage IV the bubbles escape. Apart from the effect of the ring on the bubble, the bubble is also shown to significantly affect the vortex ring, especially at low We (We similar to 3). In these low-We cases, the convection speed drops significantly compared to the base case without a bubble, while the core appears to fragment with a resultant large decrease in enstrophy by about 50 %. In the higher-We cases (We > 100), there are some differences in convection speed and enstrophy, but the effects are relatively small. The most dramatic effects of the bubble on the ring are found for thicker core rings at low We (We similar to 3) with the vortex ring almost stopping after interacting with the bubble, and the core fragmenting into two parts. The present idealized experiments exhibit many phenomena also seen in bubbly turbulent flows such as reduction in enstrophy, suppression of structures, enhancement of energy at small scales and reduction in energy at large scales. These similarities suggest that results from the present experiments can be helpful in better understanding interactions of bubbles with eddies in turbulent flows.

Surface motions and intraplate continental deformation in Alaska driven by mantle flow

The degree to which the lithosphere and mantle are coupled and contribute to surface deformation beneath continental regions remains a fundamental question in the field of geodynamics. Here we use a new approach with a surface deformation field constrained by GPS, geologic, and seismicity data, together with a lithospheric geodynamic model, to solve for tractions inferred to be generated by mantle convection that (1) drive extension within interior Alaska generating southward directed surface motions toward the southern convergent plate boundary, (2) result in accommodation of the relative motions between the Pacific and North America in a comparatively small zone near the plate boundary, and (3) generate the observed convergence within the North American plate interior in the Mackenzie mountains in northwestern Canada. The evidence for deeper mantle influence on surface deformation beneath a continental region suggests that this mechanism may be an important contributing driver to continental plate assemblage and breakup.

Effects of instabilities and coherent structures on the performance of a thermocline based thermal energy storage

We investigate the transient dynamics of disturbances inside a thermocline based molten salt thermal energy storage (TES). Numerical simulations were conducted with four inlet flow configurations. The disturbances introduced at the inlet grow via Rayleigh Taylor instability. The formed vortical motions inside the tank propagate downstream and destroy the thermocline. The vortex-thermocline interaction upsets the stratification inside the TES. The disturbance growth rate, penetration length and vortex Reynolds number are measured. The growth of penetration length prior to the vortex-thermocline interaction is quadratic. The vortex Reynolds number of the eddy which causes thermocline breakdown increases with increase in Atwood number. The impingement of vortex on thermocline is studied. (C) 2015 Elsevier Ltd. All rights reserved.

Simulations of impinging droplets with surfactant-dependent dynamic contact angle

An arbitrary Lagrangian-Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier-Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surface tension and surfactant-dependent dynamic contact angle. In particular, the dynamic contact angle (theta(d)) of the droplet is defined as a function of the surfactant concentration at the contact line and the equilibrium contact angle (theta(0)(e)) of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included into the model as surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effects of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when theta(0)(e) is less than 90 degrees, and increases it further when theta(0)(e) is greater than 90 degrees. Nevertheless, the presence of surfactants has no effect on the contact angle when theta(0)(e) = 90 degrees. The numerical study clearly demonstrates that the surfactant-dependent contact angle has to be considered, in addition to the Marangoni effect, in order to study the flow dynamics and the equilibrium states of surfactant droplet impingement accurately. The proposed numerical scheme guarantees the conservation of fluid mass and of the surfactant mass accurately. (C) 2015 Elsevier Inc. All rights reserved.

General framework for acoustic emission during plastic deformation

Despite the long history, so far there is no general theoretical framework for calculating the acoustic emission spectrum accompanying any plastic deformation. We set up a discrete wave equation with plastic strain rate as a source term and include the Rayleigh-dissipation function to represent dissipation accompanying acoustic emission. We devise a method of bridging the widely separated time scales of plastic deformation and elastic degrees of freedom. While this equation is applicable to any type of plastic deformation, it should be supplemented by evolution equations for the dislocation microstructure for calculating the plastic strain rate. The efficacy of the framework is illustrated by considering three distinct cases of plastic deformation. The first one is the acoustic emission during a typical continuous yield exhibiting a smooth stress-strain curve. We first construct an appropriate set of evolution equations for two types of dislocation densities and then show that the shape of the model stress-strain curve and accompanying acoustic emission spectrum match very well with experimental results. The second and the third are the more complex cases of the Portevin-Le Chatelier bands and the Luders band. These two cases are dealt with in the context of the Ananthakrishna model since the model predicts the three types of the Portevin-Le Chatelier bands and also Luders-like bands. Our results show that for the type-C bands where the serration amplitude is large, the acoustic emission spectrum consists of well-separated bursts of acoustic emission. At higher strain rates of hopping type-B bands, the burst-type acoustic emission spectrum tends to overlap, forming a nearly continuous background with some sharp acoustic emission bursts. The latter can be identified with the nucleation of new bands. The acoustic emission spectrum associated with the continuously propagating type-A band is continuous. These predictions are consistent with experimental results. More importantly, our study shows that the low-amplitude continuous acoustic emission spectrum seen in both the type-B and type-A band regimes is directly correlated to small-amplitude serrations induced by propagating bands. The acoustic emission spectrum of the Luders-like band matches with recent experiments as well. In all of these cases, acoustic emission signals are burstlike, reflecting the intermittent character of dislocation-mediated plastic flow.