Spectral response of a droplet in pulsating external flow field

A droplet introduced in an external convective flow field exhibits significant multimodal shape oscillations depending upon the intensity of the aerodynamic forcing. In this paper, a theoretical model describing the temporal evolution of normal modes of the droplet shape is developed. The fluid is assumed to be weakly viscous and Newtonian. The convective flow velocity, which is assumed to be incompressible and inviscid, is incorporated in the model through the normal stress condition at the droplet surface and the equation of motion governing the dynamics of each mode is derived. The coupling between the external flow and the droplet is approximated to be a one-way process, i.e., the external flow perturbations effect the droplet shape oscillations and the droplet oscillation itself does not influence the external flow characteristics. The shape oscillations of the droplet with different fluid properties under different unsteady flow fields were simulated. For a pulsatile external flow, the frequency spectra of the normal modes of the droplet revealed a dominant response at the resonant frequency, in addition to the driving frequency and the corresponding harmonics. At driving frequencies sufficiently different from the resonant frequency of the prolate-oblate oscillation mode of the droplet, the oscillations are stable. But at resonance the oscillation amplitude grows in time leading to breakup depending upon the fluid viscosity. A line vortex advecting past the droplet, simulated as an isotropic jump in the far field velocity, leads to the resonant excitation of the droplet shape modes if and only if the time taken by the vortex to cross the droplet is less than the resonant period of the P-2 mode of the droplet. A train of two vortices interacting with the droplet is also analysed. It shows clearly that the time instant of introduction of the second vortex with respect to the droplet shape oscillation cycle is crucial in determining the amplitude of oscillation. (C) 2014 AIP Publishing LLC.

Steady mixed convection flow of Maxwell fluid over an exponentially stretching vertical surface with magnetic field and viscous dissipation

The steady mixed convection flow and heat transfer from an exponentially stretching vertical surface in a quiescent Maxwell fluid in the presence of magnetic field, viscous dissipation and Joule heating have been studied. The stretching velocity, surface temperature and magnetic field are assumed to have specific exponential function forms for the existence of the local similarity solution. The coupled nonlinear ordinary differential equations governing the local similarity flow and heat transfer have been solved numerically by Chebyshev finite difference method. The influence of the buoyancy parameter, viscous dissipation, relaxation parameter of Maxwell fluid, magnetic field and Prandtl number on the flow and heat transfer has been considered in detail. The Nusselt number increases significantly with the Prandtl number, but the skin friction coefficient decreases. The Nusselt number slightly decreases with increasing viscous dissipation parameter, but the skin friction coefficient slightly increases. Maxwell fluid reduces both skin friction coefficient and Nusselt number, whereas buoyancy force enhances them.

Deformation field heterogeneity in punch indentation

Plastic heterogeneity in indentation is fundamental for understanding mechanics of hardness testing and impression-based deformation processing methods. The heterogeneous deformation underlying plane-strain indentation was investigated in plastic loading of copper by a flat punch. Deformation parameters were measured, in situ, by tracking the motion of asperities in high-speed optical imaging. These measurements were coupled with multi-scale analyses of strength, microstructure and crystallographic texture in the vicinity of the indentation. Self-consistency is demonstrated in description of the deformation field using the in situ mechanics-based measurements and post-mortem materials characterization. Salient features of the punch indentation process elucidated include, among others, the presence of a dead-metal zone underneath the indenter, regions of intense strain rate (e. g. slip lines) and extent of the plastic flow field. Perhaps more intriguing are the transitions between shear-type and compression-type deformation modes over the indentation region that were quantified by the high-resolution crystallographic texture measurements. The evolution of the field concomitant to the progress of indentation is discussed and primary differences between the mechanics of indentation for a rigid perfectly plastic material and a strain-hardening material are described.

Influence of O-2 flow rate on HfO2 gate dielectrics for back-gated graphene transistors

HfO2 thin films deposited on Si substrate using electron beam evaporation, are evaluated for back-gated graphene transistors. The amount of O-2 flow rate, during vaporation is optimized for 35 nm thick HfO2 films, to achieve the best optical, chemical and electrical properties. It has been observed that with increasing oxygen flow rate, thickness of the films increased and refractive index decreased due to increase in porosity resulting from the scattering of the evaporant. The films deposited at low O-2 flow rates (1 and 3 SCCM) show better optical and compositional properties. The effects of post-deposition annealing and post-metallization annealing in forming gas ambience (FGA) on the optical and electrical properties of the films have been analyzed. The film deposited at 3 SCCM O-2 flow rate shows the best properties as measured on MOS capacitors. To evaluate the performance of device properties, back-gated bilayer graphene transistors on HfO2 films deposited at two O-2 flow rates of 3 and 20 SCCM have been fabricated and characterized. The transistor with HfO2 film deposited at 3 SCCM O-2 flow rate shows better electrical properties consistent with the observations on MOS capacitor structures. This suggests that an optimum oxygen pressure is necessary to get good quality films for high performance devices.

Effect of angle of incidence on mixed convective wake dynamics and heat transfer past a square cylinder in cross flow at Re=100

In this paper, a numerical investigation is performed to study the mixed convective flow and heat transfer characteristics past a square cylinder in cross flow at incidence. Utilizing air (Pr = 0.71) as an operating fluid, computations are carried out at a representative Reynolds number (Re) of 100. Angles of incidences are varied as, 0 degrees <= alpha <= 45 degrees. Effect of superimposed positive and negative cross-flow buoyancy is brought about by varying the Richardson number (RI) in the range -1.0 <= Ri <= 1.0. The detail features of flow topology and heat transport are analyzed critically for different angles of incidences. The thermo fluidic forces acting on the cylinder during mixed convection are captured in terms of the drag (C-D), lift (C-L), and moment (C-M) coefficients. The results show that the lateral width of the cylinder wake reduces with increasing alpha and the isotherms spread out far wide. In the range 0 degrees < alpha < 45 degrees, C-D reduces with increasing Ri. The functional dependence of C-M with Ri reveals a linear relationship. Thermal boundary layer thickness reduces with increasing angle of incidences. The global rate of heat transfer from the cylinder increases with increasing alpha. (C) 2014 Elsevier Ltd. All rights reserved.

Film thickness mediated transition in the kinetics of electric current induced flow of thin liquid metal films

Application of high electric-field between two points in a thin metallic film results in liquefaction and subsequent flow of the liquid-film from one electrode to another in a radially symmetric fashion. Here, we report the transition of the flow kinetics driven by the liquid film thickness varying from 3 to 100 nm. The mechanism of the flow behavior is observed to be independent of the film thickness; however, the kinetics of the flow depends on the film thickness and the applied voltage. An analytical model, incorporating viscosity and varying electrical resistivity with film thickness, is developed to explain the experimental observations. (C) 2014 AIP Publishing LLC.

Energy based approach for the evaluation of damage under partial slip and gross sliding condition

Three possible contact conditions may prevail at a contact interface depending on the magnitude of normal and tangential loads, that is, stick condition, partial slip condition or gross sliding condition. Numerical techniques have been used to evaluate the stress field under partial slip and gross sliding condition. Cattaneo and Mindlin approach has been adapted to model partial slip condition. Shear strain energy density and normalized strain energy release rate have been evaluated at the surface and in the subsurface region. It is apparent from the present study that the shear strain energy density gives a fair prediction for the nucleation of damage, whereas the propagation of the crack is controlled by normalized strain energy release rate. Further, it has been observed that the intensity of damage strongly depends on coefficient of friction and contact conditions prevailing at the contact interface. (C) 2014 Elsevier B.V. All rights reserved.

Slow Slip Acceleration beneath Andaman Islands Triggered by the 11 April 2012 Indian Ocean Earthquakes

The M-w 8.6 and 8.2 strike-slip earthquakes that struck the northeast Indian Ocean on 11 April 2012 resulted in coseismic deformation both at near and distant sites. The slip distribution, deduced using seismic-wave analysis for the orthogonal faults that ruptured during these earthquakes, is sufficient to predict the coseismic displacements at the Global Positioning System (GPS) sites, such as NTUS, PALK, and CUSV, but fall short at four continuous sites in the Andaman Islands region. Slip modeling, for times prior to the events, suggests that the lower portion of the thrust fault beneath the Andaman Islands has been slipping at least at the rate of 40 cm/yr, in response to the 2004 Sumatra-Andaman coseismic stress change. Modeling of GPS displacements suggests that the en echelon and orthogonal fault ruptures of the 2012 intraplate oceanic earthquakes could have possibly accelerated the ongoing slow slip, along the lower portion of the thrust fault beneath the islands with a month-long slip of 4-10 cm. The misfit to the coseismic GPS displacements along the Andaman Islands could be improved with a better source model, assuming that no local process contributed to this anomaly.

Parallel Flow-Sensitive Pointer Analysis by Graph-Rewriting

Precise pointer analysis is a problem of interest to both the compiler and the program verification community. Flow-sensitivity is an important dimension of pointer analysis that affects the precision of the final result computed. Scaling flow-sensitive pointer analysis to millions of lines of code is a major challenge. Recently, staged flow-sensitive pointer analysis has been proposed, which exploits a sparse representation of program code created by staged analysis. In this paper we formulate the staged flow-sensitive pointer analysis as a graph-rewriting problem. Graph-rewriting has already been used for flow-insensitive analysis. However, formulating flow-sensitive pointer analysis as a graph-rewriting problem adds additional challenges due to the nature of flow-sensitivity. We implement our parallel algorithm using Intel Threading Building Blocks and demonstrate considerable scaling (upto 2.6x) for 8 threads on a set of 10 benchmarks. Compared to the sequential implementation of staged flow-sensitive analysis, a single threaded execution of our implementation performs better in 8 of the benchmarks.

TRANSIENT NATURAL CONVECTION FLOW IN A RECTANGULAR CAVITY FILLED WITH A POROUS MATERIAL WITH LOCALIZED HEATING FROM BELOW AND THERMAL STRATIFICATION

The transient natural convection flow with thermal stratification in a rectangular cavity filled with fluid saturated porous medium obeying Darcy's law has been studied. Prior to the time t* = 0, the flow in the cavity is assumed to be motionless and all four walls of the cavity are at the same constant temperature. At time t* = 0, the temperatures of the vertical walls are suddenly increased which vary linearly with the distance y and at the same time on the bottom wall an isothermal heat source is placed centrally. This sudden change in the wall temperatures gives rise to unsteadiness in the problem. The horizontal temperature difference induces and sustains a buoyancy driven flow in the cavity which is then controlled by the vertical temperature difference. The partial differential equations governing the transient natural convection flow have been solved numerically. The local and average Nusselt numbers decrease rapidly in a small time interval after the start of the impulsive change in the wall temperatures and the steady state is reached quickly. The time required to reach the steady state depends on the Rayleigh number and the thermal stratification parameter.

Dense shallow granular flows

Simplified equations are derived for a granular flow in the `dense' limit where the volume fraction is close to that for dynamical arrest, and the `shallow' limit where the stream-wise length for flow development (L) is large compared with the cross-stream height (h). The mass and diameter of the particles are set equal to 1 in the analysis without loss of generality. In the dense limit, the equations are simplified by taking advantage of the power-law divergence of the pair distribution function chi proportional to (phi(ad) - phi)(-alpha), and a faster divergence of the derivativ rho(d chi/d rho) similar to (d chi/d phi), where rho and phi are the density and volume fraction, and phi(ad) is the volume fraction for arrested dynamics. When the height h is much larger than the conduction length, the energy equation reduces to an algebraic balance between the rates of production and dissipation of energy, and the stress is proportional to the square of the strain rate (Bagnold law). In the shallow limit, the stress reduces to a simplified Bagnold stress, where all components of the stress are proportional to (partial derivative u(x)/partial derivative y)(2), which is the cross-stream (y) derivative of the stream-wise (x) velocity. In the simplified equations for dense shallow flows, the inertial terms are neglected in the y momentum equation in the shallow limit because the are O(h/L) smaller than the divergence of the stress. The resulting model contains two equations, a mass conservation equations which reduces to a solenoidal condition on the velocity in the incompressible limit, and a stream-wise momentum equation which contains just one parameter B which is a combination of the Bagnold coefficients and their derivatives with respect to volume fraction. The leading-order dense shallow flow equations, as well as the first correction due to density variations, are analysed for two representative flows. The first is the development from a plug flow to a fully developed Bagnold profile for the flow down an inclined plane. The analysis shows that the flow development length is ((rho) over barh(3)/B) , where (rho) over bar is the mean density, and this length is numerically estimated from previous simulation results. The second example is the development of the boundary layer at the base of the flow when a plug flow (with a slip condition at the base) encounters a rough base, in the limit where the momentum boundary layer thickness is small compared with the flow height. Analytical solutions can be found only when the stream-wise velocity far from the surface varies as x(F), where x is the stream-wise distance from the start of the rough base and F is an exponent. The boundary layer thickness increases as (l(2)x)(1/3) for all values of F, where the length scale l = root 2B/(rho) over bar. The analysis reveals important differences between granular flows and the flows of Newtonian fluids. The Reynolds number (ratio of inertial and viscous terms) turns out to depend only on the layer height and Bagnold coefficients, and is independent of the flow velocity, because both the inertial terms in the conservation equations and the divergence of the stress depend on the square of the velocity/velocity gradients. The compressibility number (ratio of the variation in volume fraction and mean volume fraction) is independent of the flow velocity and layer height, and depends only on the volume fraction and Bagnold coefficients.

Modelling of Effect of Non-Uniform Current Density on the Performance of Soluble Lead Redox Flow Batteries

Soluble lead acid redox flow battery (SLRFB) offers a number of advantages. These advantages can be harnessed after problems associated with buildup of active material on. electrodes (residue) are resolved. A mathematical model is developed to understand residue formation in SLRFB. The model incorporates fluid flow, ion transport, electrode reactions, and non-uniform current distribution on electrode surfaces. A number of limiting cases are studied to conclude that ion transport and electrode reaction on anode simultaneously control battery performance. The model fits the reported cell voltage vs. time profiles very well. During the discharge cycle, the model predicts complete dissolution of deposited material from trailing edge side of the electrodes. With time, the active surface area of electrodes decreases rapidly. The corresponding increase in current density leads to precipitous decrease in cell potential before all the deposited material is dissolved. The successive charge-discharge cycles add to the residue. The model correctly captures the marginal effect of flow rate on cell voltage profiles, and identifies flow rate and flow direction as new variables for controlling residue buildup. Simulations carried out with alternating flow direction and a SLRFB with cylindrical electrodes show improved performance with respect to energy efficiency and residue buildup. (C) 2014 The Electrochemical Society. All rights reserved.

Breakup processes of pressure swirl spray in gaseous cross-flow

Injection of liquid fuel in cross flowing air has been a strategy for future aircraft engines in order to control the emissions. In this context, breakup of a pressure swirl spray in gaseous cross-flow is investigated experimentally. The atomizer discharges a conical swirling sheet of liquid that interacts with cross-flowing air. This complex interaction and the resulting spray structures at various flow conditions are studied through flow visualization using still as well as high speed photography. Experiments are performed over a wide range of aerodynamic Weber number (2-300) and liquid-to-air momentum flux ratio (5-150). Various breakup regimes exhibiting different breakup processes are mapped on a parameter space based on flow conditions. This map shows significant variations from breakup regime map for a plain liquid jet in cross-flow. It is observed that the breakup of leeward side of the sheet is dominated by bag breakup and the windward side of the sheet undergoes breakup through surface waves. Similarities and differences between bag breakup present in plain liquid jet in cross-flow and swirl spray in cross-flow are explained. Multimodal drop size distribution from bag breakup, frequency of bag breakup, wavelength of surface waves and trajectory of spray in cross-flow are measured by analyzing the spray images and parametric study of their variations is also presented. (C) 2014 Elsevier Ltd. All rights reserved.

Effect of surfactants on the instability of a two-layer film flow down an inclined plane

The effect of insoluble surfactants on the instability of a two-layer film flow down an inclined plane is investigated based on the Orr-Sommerfeld boundary value problem. The study, focusing on Stokes flow P. Gao and X.-Y. Lu, ``Effect of surfactants on the inertialess instability of a two-layer film flow,'' J. Fluid Mech. 591, 495-507 (2007)], is further extended by including the inertial effect. The surface mode is recognized along with the interface mode. The initial growth rate corresponding to the interface mode accelerates at sufficiently long-wave regime in the presence of surface surfactant. However, the maximum growth rate corresponding to both interface and surface modes decelerates in the presence of surface surfactant when the upper layer is more viscous than the lower layer. On the other hand, when the upper layer is less viscous than the lower layer, a new interfacial instability develops due to the inertial effect and becomes weaker in the presence of interfacial surfactant. In the limit of negligible surface and interfacial tensions, respectively, two successive peaks of temporal growth rate appear in the long-wave and short-wave regimes when the interface mode is analyzed. However, in the case of the surface mode, only the long-wave peak appears. (C) 2014 AIP Publishing LLC.

Effect of shrinkage induced flow on binary alloy dendrite growth: An equivalent undercooling model

This paper presents a theoretical model for studying the effects of shrinkage induced flow on the growth rate of binary alloy dendrites. An equivalent undercooling of the melt is defined in terms of ratio of the phase densities to represent the change in dendrite growth rate due to variation in solutal and thermal transport resulting from shrinkage induced flow. Subsequently, results for dendrite growth rate predicted by the equivalent undercooling model is compared with the corresponding predictions obtained using an enthalpy based numerical method for dendrite growth with shrinkage. The agreement is found to be good. Published by Elsevier Ltd.