Bistable Vibration Energy Harvesters: a review

Powering electronics without depending on batteries is an open research field. Mechanical vibrations prove to be a reliable energy source, but low-frequency broadband vibrations cannot be harvested effectively using linear oscillators. This article discusses an alternative for harvesting such vibrations, with energy harvesters with two stable configurations. The challenges related to nonlinear dynamics are briefly discussed. Different existing designs of bistable energy harvesters are presented and classified, according to their feasibility for miniaturization. A general dynamic model for those designs is described. Finally, an extensive discussion on quantitative measures of evaluating the effectiveness of energy harvesters is accomplished, resulting in the proposition of a new dimensionless metric suited for a broadband analysis.

Density stratified environments: The double-tank method

We present an alternative method of producing density stratifications in the laboratory based on the 'double-tank' method proposed by Oster (Sci Am 213:70-76, 1965). We refer to Oster's method as the 'forced-drain' approach, as the volume flow rates between connecting tanks are controlled by mechanical pumps. We first determine the range of density profiles that may be established with the forced-drain approach other than the linear stratification predicted by Oster. The dimensionless density stratification is expressed analytically as a function of three ratios: the volume flow rate ratio n, the ratio of the initial liquid volumes λ and the ratio of the initial densities ψ. We then propose a method which does not require pumps to control the volume flow rates but instead allows the connecting tanks to drain freely under gravity. This is referred to as the 'free-drain' approach. We derive an expression for the density stratification produced and compare our predictions with saline stratifications established in the laboratory using the 'free-drain' extension of Oster's method. To assist in the practical application of our results we plot the region of parameter space that yield concave/convex or linear density profiles for both forced-drain and free-drain approaches. The free-drain approach allows the experimentalist to produce a broad range of density profiles by varying the initial liquid depths, cross-sectional and drain opening areas of the tanks. One advantage over the original Oster approach is that density profiles with an inflexion point can now be established. © 2008 Springer-Verlag.

Pollutant flushing with natural displacement ventilation

We examine the time taken to flush pollutants from a naturally ventilated room. A simple theoretical model is developed to predict the time taken for neutrally-buoyant pollutants to be removed from a room by a flow driven by localised heat inputs; both line and point heat sources are considered. We show that the rate of flushing is a function of the room volume, vent areas ( A) and the distribution, number (n) and strength (B) of the heat sources. We also show that the entire problem can be reduced to a single parameter ( μ) that is a measure of the vent areas, and a dimensionless time ( τ) that is a function of B, V and μ. Small-scale salt-bath experiments were conducted to measure the flushing rates in order to validate our modelling assumptions and predictions. The predicted flushing times show good agreement with the experiments over a wide range of μ. We apply our model to a typical open plan office and lecture theatre and discuss some of the implications of our results. © 2005 Elsevier Ltd. All rights reserved.

The fluid mechanics of a two-dimensional Aaberg exhaust hood

In this paper a study of the air flow pattern created by a two-dimensional Aaberg exhaust hood local ventilation system is presented. A mathematical model of the flow, in terms of the stream function ψ, is derived analytically for both laminar and turbulent injections of fluid. Streamlines and lines of constant speed deduced from the model are examined for various values of the governing dimensionless operating parameter and predictions are given as to the area in front of the hood from which the air can be sampled. The effect of the injection of fluid on the centre-line velocity of the flow is examined and a comparison of the results with the available experimental data is given. © 1992.

Design guidelines for granular particles in a conical centrifugal filter

Essential design criteria for successful drying of granular particles in a conical continuous centrifugal filter are developed in a dimensionless fashion. Four criteria are considered: minimum flow thickness (to ensure sliding bulk flow rather than particulate flow), desaturation position, output dryness and basket failure. The criteria are based on idealised physical models of the machine operation and are written explicitly as functions of the basket size lout, spin velocity Ω and input flow rate of powder Qp. The separation of sucrose crystals from liquid molasses is taken as a case study and the successful regime of potential operating points (lout, Ω) is plotted for a wide range of selected values of flow rate Qp. Analytical expressions are given for minimum and maximum values of the three independent parameters (lout, Ω, Qp) as a function of the slurry and basket properties. The viable operating regime for a conical centrifugal filter is thereby obtained as a function of the slurry and basket properties. © 2012 The Institution of Chemical Engineers.

Blade Loading and Its Application in the Mean-Line Design of Low Pressure Turbines

In order to minimize the number of iterations to a turbine design, reasonable choices of the key parameters must be made at the preliminary design stage. The choice of blade loading is of particular concern in the low pressure (LP) turbine of civil aero engines, where the use of high-lift blades is widespread. This paper considers how blade loading should be measured, compares the performance of various loss correlations, and explores the impact of blade lift on performance and lapse rates. To these ends, an analytical design study is presented for a repeating-stage, axial-flow LP turbine. It is demonstrated that the long-established Zweifel lift coefficient (Zweifel, 1945, "The Spacing of Turbomachine Blading, Especially with Large Angular Deflection" Brown Boveri Rev., 32(1), pp. 436-444) is flawed because it does not account for the blade camber. As a result the Zweifel coefficient is only meaningful for a fixed set of flow angles and cannot be used as an absolute measure of blade loading. A lift coefficient based on circulation is instead proposed that accounts for the blade curvature and is independent of the flow angles. Various existing profile and secondary loss correlations are examined for their suitability to preliminary design. A largely qualitative comparison demonstrates that the loss correlations based on Ainley and Mathieson (Ainley and Mathieson, 1957, "A Method of Performance Estimation for Axial-Flow Turbines," ARC Reports and Memoranda No. 2974; Dunham and Came, 1970, "Improvements to the Ainley-Mathieson Method of Turbine Performance Prediction," Trans. ASME: J. Eng. Gas Turbines Power, July, pp. 252-256; Kacker and Okapuu, 1982, "A Mean Line Performance Method for Axial Flow Turbine Efficiency," J. Eng. Power, 104, pp. 111-119). are not realistic, while the profile loss model of Coull and Hodson (Coull and Hodson, 2011, "Predicting the Profile Loss of High-Lift Low Pressure Turbines," J. Turbomach., 134(2), pp. 021002) and the secondary loss model of (Traupel, W, 1977, Thermische Turbomaschinen, Springer-Verlag, Berlin) are arguably the most reasonable. A quantitative comparison with multistage rig data indicates that, together, these methods over-predict lapse rates by around 30%, highlighting the need for improved loss models and a better understanding of the multistage environment. By examining the influence of blade lift across the Smith efficiency chart, the analysis demonstrates that designs with higher flow turning will tend to be less sensitive to increases in blade loading. © 2013 American Society of Mechanical Engineers.

Probabilistic soil identification based on cone penetration tests

In geotechnical engineering, soil classification is an essential component in the design process. Field methods such as the cone penetration test (CPT) can be used as less expensive and faster alternatives to sample retrieval and testing. Unfortunately, current soil classification charts based on CPT data and laboratory measurements are too generic, and may not provide an accurate prediction of the soil type. A probabilistic approach is proposed here to update and modify soil identification charts based on site-specific CPT data. The probability that a soil is correctly classified is also estimated. The updated identification chart can be used for a more accurate prediction of the classification of the soil, and can account for prior information available before conducting the tests, site-specific data, and measurement errors. As an illustration, the proposed approach is implemented using CPT data from the Treporti Test Site (TTS) near Venice (Italy) and the National Geotechnical Experimentation Sites (NGES) at Texas A&M University. The applicability of the site-specific chart for other sites in Venice Lagoon is assessed using data from the Malamocco test site, approximately 20 km from TTS.

Simulations of thermal upgrading methods for oil shale reservoirs

This paper analyzes reaction and thermal front development in porous reservoirs with reacting flows, such as those encountered in shale oil extraction. A set of dimensionless parameters and a 3D code are developed in order to investigate the important physical and chemical variables of such reservoirs when heated by in situ methods. This contribution builds on a 1D model developed for the precursor study to this work. Theory necessary for this study is presented, namely shale decomposition chemical mechanisms, governing equations for multiphase flow in porous media and necessary closure models. Plotting the ratio of the thermal wave speed to the fluid speed allows one to infer that the reaction wave front ends where this ratio is at a minimum. The reaction front follows the thermal front closely, thus allowing assumptions to be made about the extent of decomposition solely by looking at thermal wave progression. Furthermore, this sensitivity analysis showed that a certain minimum permeability is required in order to ensure the formation of a traveling thermal wave. It was found that by studying the non-dimensional governing parameters of the system one can ascribe characteristic values for these parameters for given initial and boundary conditions. This allows one to roughly predict the performance of a particular method on a particular reservoir given approximate values for initial and boundary conditions. Channelling and flow blockage due to carbon residue buildup impeded each method's performance. Blockage was found to be a result of imbalanced heating. Copyright 2012, Society of Petroleum Engineers.

Nonlinear drillstring dynamics analysis

This paper studies the dynamical response of a rotary drilling system with a drag bit, using a lumped parameter model that takes into consideration the axial and torsional vibration modes of the bit. These vibrations are coupled through a bit-rock interaction law. At the bit-rock interface, the cutting process introduces a state-dependent delay, while the frictional process is responsible for discontinuous right-hand sides in the equations governing the motion of the bit. This complex system is characterized by a fast axial dynamics compared to the slow torsional dynamics. A dimensionless formulation exhibits a large parameter in the axial equation, enabling a two-time-scales analysis that uses a combination of averaging methods and a singular perturbation approach. An approximate model of the decoupled axial dynamics permits us to derive a pseudoanalytical expression of the solution of the axial equation. Its averaged behavior influences the slow torsional dynamics by generating an apparent velocity weakening friction law that has been proposed empirically in earlier work. The analytical expression of the solution of the axial dynamics is used to derive an approximate analytical expression of the velocity weakening friction law related to the physical parameters of the system. This expression can be used to provide recommendations on the operating parameters and the drillstring or the bit design in order to reduce the amplitude of the torsional vibrations. Moreover, it is an appropriate candidate model to replace empirical friction laws encountered in torsional models used for control. © 2009 Society for Industrial and Applied Mathematics.

The rhythm of fountains: The length and time scales of rise height fluctuations at low and high Froude numbers

The magnitude and frequency of vertical fluctuations of the top of an axisymmetric miscible Boussinesq fountain forms the focus of this work. We present measurements of these quantities for saline-aqueous fountains in uniform quiescent surroundings. Our results span source Froude numbers 0.3 ≤ Fr 0 ≤ 40 and, thereby, encompass very weak, weak, intermediate and forced classes of fountain. We identify distinct scalings, based on known quantities at the fountain source, for the frequency of fountain height fluctuations which collapse our data within bands of Fr0. Notably, our scalings reveal that the (dimensionless) frequency takes a constant value within each band. These results highlight characteristic time scales for the fluctuations which we decompose into a single, physically apparent, length scale and velocity scale within each band. Moreover, within one particular band, spanning source Froude numbers towards the lower end of the full range considered, we identify unexpectedly long-period fluctuations indicating a near balance of inertia and (opposing) buoyancy at the source. Our analysis identifies four distinct classes of fluctuation behaviour (four bands of Fr 0) and this classification matches well with existing classifications of fountains based on rise heights. As such, we show that an analysis of the behaviour of the fountain top alone, rather than the entire fountain, provides an alternative approach to classifying fountains. The similarity of classifications based on the two different methods confirms that the boundaries between classes mark tangible changes in the physics of fountains. For high Fr0 we show that the dominant fluctuations occur at the scale of the largest eddies which can be contained within the fountain near its top. Extending this, we develop a Strouhal number, Strtop, based on experimental measures of the fountain top, defined such that Strtop = 1 would suggest the dominant fluctuations are caused by a continual cycle of eddies forming and collapsing at this largest physical scale. For high- Fr 0 fountains we find Strtop ≈ 0. 9. © 2013 Cambridge University Press.

Direct numerical simulation of turbulent flows over superhydrophobic surfaces with gas pockets using linearized boundary conditions

Superhydrophobic surfaces are shown to be effective for surface drag reduction under laminar regime by both experiments and simulations (see for example, Ou and Rothstein, Phys. Fluids 17:103606, 2005). However, such drag reduction for fully developed turbulent flow maintaining the Cassie-Baxter state remains an open problem due to high shear rates and flow unsteadiness of turbulent boundary layer. Our work aims to develop an understanding of mechanisms leading to interface breaking and loss of gas pockets due to interactions with turbulent boundary layers. We take advantage of direct numerical simulation of turbulence with slip and no-slip patterned boundary conditions mimicking the superhydrophobic surface. In addition, we capture the dynamics of gas-water interface, by deriving a proper linearized boundary condition taking into account the surface tension of the interface and kinematic matching of interface deformation and normal velocity conditions on the wall. We will show results from our simulations predicting the dynamical behavior of gas pocket interfaces over a wide range of dimensionless surface tensions.

Swimming dynamics of bidirectional artificial flagella

We study magnetic artificial flagella whose swimming speed and direction can be controlled using light and magnetic field as external triggers. The dependence of the swimming velocity on the system parameters (e.g., length, stiffness, fluid viscosity, and magnetic field) is explored using a computational framework in which the magnetostatic, fluid dynamic, and solid mechanics equations are solved simultaneously. A dimensionless analysis is carried out to obtain an optimal combination of system parameters for which the swimming velocity is maximal. The swimming direction reversal is addressed by incorporating photoresponsive materials, which in the photoactuated state can mimic natural mastigonemes. © 2013 American Physical Society.

Swimming dynamics of bidirectional artificial flagella.

We study magnetic artificial flagella whose swimming speed and direction can be controlled using light and magnetic field as external triggers. The dependence of the swimming velocity on the system parameters (e.g., length, stiffness, fluid viscosity, and magnetic field) is explored using a computational framework in which the magnetostatic, fluid dynamic, and solid mechanics equations are solved simultaneously. A dimensionless analysis is carried out to obtain an optimal combination of system parameters for which the swimming velocity is maximal. The swimming direction reversal is addressed by incorporating photoresponsive materials, which in the photoactuated state can mimic natural mastigonemes.

Investigation of a radial-inflow bleed as a potential for compressor clearance control

The mismatch in thermal response between a High Pressure Compressor (HPC) drum and casing is a limiting factor in the reduction of compressor clearance. An experimental test rig has been used to demonstrate the concept of radial inflow to reduce the thermal time constant of HPC discs. The testing uses a simulated idle - Maximum Take Off (MTO) - idle transient in order to measure the thermal response directly. The testing is fully scaled in the dimensionless sense to engine conditions. A simple closure model based on lumped capacitance is used to illustrate the scope of potential benefits. The proof-of-concept testing shows that HPC disc time constant reductions of the order 2 are feasible with a radial-inflow bleed of only 4% of bore flow at scaled MTO conditions. Using the experimental results, the simple closure modelling suggests that for a stage with a significant mismatch in thermal response, reductions in 2D axis-symmetric clearance of as much as 50% at MTO conditions may be possible along with significant scope for improvements at cruise conditions. Copyright © 2013 by ASME.

Transient ventilation dynamics induced by heat sources of unequal strength

We examine theoretically the transient displacement flow and density stratification that develops within a ventilated box after two localized floor-level heat sources of unequal strengths are activated. The heat input is represented by two non-interacting turbulent axisymmetric plumes of constant buoyancy fluxes B1 and B2 > B1. The box connects to an unbounded quiescent external environment of uniform density via openings at the top and base. A theoretical model is developed to predict the time evolution of the dimensionless depths λj and mean buoyancies δj of the 'intermediate' (j = 1) and 'top' (j = 2) layers leading to steady state. The flow behaviour is classified in terms of a stratification parameter S, a dimensionless measure of the relative forcing strengths of the two buoyant layers that drive the flow. We find that dδ1/dτ α 1/λ1 and dδ2/dτ α 1/λ2, where τ is a dimensionless time. When S 1, the intermediate layer is shallow (small λ1), whereas the top layer is relatively deep (large λ2) and, in this limit, δ1 and δ2 evolve on two characteristically different time scales. This produces a time lag and gives rise to a 'thermal overshoot', during which δ1 exceeds its steady value and attains a maximum during the transients; a flow feature we refer to, in the context of a ventilated room, as 'localized overheating'. For a given source strength ratio ψ = B1/B2, we show that thermal overshoots are realized for dimensionless opening areas A < Aoh and are strongly dependent on the time history of the flow. We establish the region of {A, ψ} space where rapid development of δ1 results in δ1 > δ2, giving rise to a bulk overturning of the buoyant layers. Finally, some implications of these results, specifically to the ventilation of a room, are discussed. © Cambridge University Press 2013.