Imaging the femoral cortex: Thickness, density and mass from clinical CT

There is growing evidence that focal thinning of cortical bone in the proximal femur may predispose a hip to fracture. Detecting such defects in clinical CT is challenging, since cortices may be significantly thinner than the imaging system's point spread function. We recently proposed a model-fitting technique to measure sub-millimetre cortices, an ill-posed problem which was regularized by assuming a specific, fixed value for the cortical density. In this paper, we develop the work further by proposing and evaluating a more rigorous method for estimating the constant cortical density, and extend the paradigm to encompass the mapping of cortical mass (mineral mg/cm 2) in addition to thickness. Density, thickness and mass estimates are evaluated on sixteen cadaveric femurs, with high resolution measurements from a micro-CT scanner providing the gold standard. The results demonstrate robust, accurate measurement of peak cortical density and cortical mass. Cortical thickness errors are confined to regions of thin cortex and are bounded by the extent to which the local density deviates from the peak, averaging 20% for 0.5mm cortex. © 2012 Elsevier B.V.

Interpreting the input admittance of violins and guitars

The most widespread vibration measurement on musical instrument bodies is of the point mobility at the bridge. Analysis of such measurements is presented, with a view to assessing what range of information could feasibly be extracted from the corpus of data. Analysis approaches include (1) pole-residue extraction; (2) damping trend analysis based on time decay information; (3) statistical estimates based on SEA power-balance and variance theory. Comparative results are shown for some key quantities. Damping trends with frequency are shown to have unexpectedly different forms for violins and for guitars. Linear averaging to estimate the "direct field" component gives a simple and clear visualisation of any local resonance behaviour near the bridge, such as the "bridge hill", and reveals some violins that show a double hill, while viols show only weak hills, and guitars none at all. © S. Hirzel Verlag · EAA.

Effects of multi-scattering on the performance of a single-beam acoustic manipulation device.

The effects of multiple scattering on acoustic manipulation of spherical particles using helicoidal Bessel-beams are discussed. A closed-form analytical solution is developed to calculate the acoustic radiation force resulting from a Bessel-beam on an acoustically reflective sphere, in the presence of an adjacent spherical particle, immersed in an unbounded fluid medium. The solution is based on the standard Fourier decomposition method and the effect of multi-scattering is taken into account using the addition theorem for spherical coordinates. Of particular interest here is the investigation of the effects of multiple scattering on the emergence of negative axial forces. To investigate the effects, the radiation force applied on the target particle resulting from a helicoidal Bessel-beam of different azimuthal indexes (m = 1 to 4), at different conical angles, is computed. Results are presented for soft and rigid spheres of various sizes, separated by a finite distance. Results have shown that the emergence of negative force regions is very sensitive to the level of cross-scattering between the particles. It has also been shown that in multiple scattering media, the negative axial force may occur at much smaller conical angles than previously reported for single particles, and that acoustic manipulation of soft spheres in such media may also become possible.

Decentralised minimal-time dynamic consensus

This paper considers a group of agents that aim to reach an agreement on individually received time-varying signals by local communication. In contrast to static network averaging problem, the consensus considered in this paper is reached in a dynamic sense. A discrete-time dynamic average consensus protocol can be designed to allow all the agents tracking the average of their reference inputs asymptotically. We propose a minimal-time dynamic consensus algorithm, which only utilises a minimal number of local observations of a randomly picked node in a network to compute the final consensus signal. Our results illustrate that with memory and computational ability, the running time of distributed averaging algorithms can be indeed improved dramatically as suggested by Olshevsky and Tsitsiklis. © 2012 AACC American Automatic Control Council).

RFID Tag Location Systems

We describe a methods of locating an RFID tag. One method comprises: transmitting tag location signals from a plurality of different transmit antennas, wherein said antennas are spaced apart by more than a near field limit distance at a frequency of a said signal; receiving a corresponding plurality of receiving return signals from said tag; and processing said tag return signals to determine a range to said tag; wherein said transmitting comprises transmitting at a plurality of different frequencies; wherein said processing comprises determining a phase difference at said plurality of different frequencies to determine said range, and wherein said determining of said phase difference determines a phase difference between either i) two or more of said transmit signals resulting in a maxima in the returned signal RSSI or ii) a first transmit signal and its corresponding return signal; and wherein said determining of said range to said tag uses said return signals weighted responsive to a respective received signal strength of the return signal. Further data which may be used for averaging may be generated by using the above techniques along with changes in the polarisation state of the transmit and receive antennas and/or physical reconfiguration of the antennas (e.g. switch the transmit and receive elements).

A neurophysiologically plausible population code model for feature integration explains visual crowding.

An object in the peripheral visual field is more difficult to recognize when surrounded by other objects. This phenomenon is called "crowding". Crowding places a fundamental constraint on human vision that limits performance on numerous tasks. It has been suggested that crowding results from spatial feature integration necessary for object recognition. However, in the absence of convincing models, this theory has remained controversial. Here, we present a quantitative and physiologically plausible model for spatial integration of orientation signals, based on the principles of population coding. Using simulations, we demonstrate that this model coherently accounts for fundamental properties of crowding, including critical spacing, "compulsory averaging", and a foveal-peripheral anisotropy. Moreover, we show that the model predicts increased responses to correlated visual stimuli. Altogether, these results suggest that crowding has little immediate bearing on object recognition but is a by-product of a general, elementary integration mechanism in early vision aimed at improving signal quality.

Lift and the leading edge vortex

Leading edge vortices are considered to be important in generating the high lift coefficients observed in insect flight and may therefore be relevant to micro-air vehicles. A potential flow model of an impulsively started flat plate, featuring a leading edge vortex (LEV) and a trailing edge vortex (TEV) is fitted to experimental data in order to provide insight into the mechanisms that influence the convection of the LEV and to study how the LEV contributes to lift. The potential flow model fits the experimental data best with no bound circulation, which is in accordance with Kelvin's circulation theorem. The lift-to-drag ratio is well approximated by the function 'cot α' for α > 15°, which supports the tentative conclusion that shortly after an impulsive start, at post-stall angles of attack, lift is caused non-circulatory forces and by the action of the LEV as opposed to bound circulation. Copyright © 2012 by C. W. Pitt Ford.

Closely correlating lower and upper bound plastic analysis of real slabs

An engineer assessing the load-carrying capacity of an existing reinforced concrete slab is likely to use elastic analysis to check the load at which the structure might be expected to fail in flexure or in shear. In practice, many reinforced concrete slabs are highly ductile in flexure, so an elastic analysis greatly underestimates the loads at which they fail in this mode. The use of conservative elastic analysis has led engineers to incorrectly condemn many slabs and therefore to specify unnecessary and wasteful flexural strengthening or replacement. The lower bound theorem is based on the same principles as the upper bound theorem used in yield line analysis, but any solution that rigorously satisfies the lower bound theorem is guaranteed to be a safe underestimate of the collapse load. Jackson presented a rigorous lower bound method that obtains very accurate results for complex real slabs.

Lift and the leading-edge vortex

Flapping wings often feature a leading-edge vortex (LEV) that is thought to enhance the lift generated by the wing. Here the lift on a wing featuring a leading-edge vortex is considered by performing experiments on a translating flat-plate aerofoil that is accelerated from rest in a water towing tank at a fixed angle of attack of 15°. The unsteady flow is investigated with dye flow visualization, particle image velocimetry (PIV) and force measurements. Leading-and trailing-edge vortex circulation and position are calculated directly from the velocity vectors obtained using PIV. In order to determine the most appropriate value of bound circulation, a two-dimensional potential flow model is employed and flow fields are calculated for a range of values of bound circulation. In this way, the value of bound circulation is selected to give the best fit between the experimental velocity field and the potential flow field. Early in the trajectory, the value of bound circulation calculated using this potential flow method is in accordance with Kelvin's circulation theorem, but differs from the values predicted by Wagner's growth of bound circulation and the Kutta condition. Later the Kutta condition is established but the bound circulation remains small; most of the circulation is contained instead in the LEVs. The growth of wake circulation can be approximated by Wagner's circulation curve. Superimposing the non-circulatory lift, approximated from the potential flow model, and Wagner's lift curve gives a first-order approximation of the measured lift. Lift is generated by inertial effects and the slow buildup of circulation, which is contained in shed vortices rather than bound circulation. © 2013 Cambridge University Press.

A six sensor method for measuring acoustic properties in ducts.

An accurate description of sound propagation in a duct is important to obtain the sound power radiating from a source in both near and far fields. A technique has been developed and applied to decompose higher-order modes of sound emitted into a duct. Traditional experiments and theory based on two-sensor methods are limited to the plane-wave contribution to the sound field at low frequency. Due to the increase in independent measurements required, a computational method has been developed to simulate sensitivities of real measurements (e.g., noise) and optimize the set-up. An experimental rig has been constructed to decompose the first two modes using six independent measurements from surface, flush-mounted microphones. Experiments were initially performed using a loudspeaker as the source for validation. Subsequently, the sound emitted by a mixed-flow fan has been investigated and compared to measurements made in accordance with the internationally standardized in-duct fan measurement method. This method utilizes large anechoic terminations and a procedure involving averaging over measurements in space and time to account for the contribution from higher-order modes. The new method does not require either of these added complications and gives detail about the underlying modal content of the emitted sound.

An experimental study on engine dynamics model based in-cylinder pressure estimation

The information provided by the in-cylinder pressure signal is of great importance for modern engine management systems. The obtained information is implemented to improve the control and diagnostics of the combustion process in order to meet the stringent emission regulations and to improve vehicle reliability and drivability. The work presented in this paper covers the experimental study and proposes a comprehensive and practical solution for the estimation of the in-cylinder pressure from the crankshaft speed fluctuation. Also, the paper emphasizes the feasibility and practicality aspects of the estimation techniques, for the real-time online application. In this study an engine dynamics model based estimation method is proposed. A discrete-time transformed form of a rigid-body crankshaft dynamics model is constructed based on the kinetic energy theorem, as the basis expression for total torque estimation. The major difficulties, including load torque estimation and separation of pressure profile from adjacent-firing cylinders, are addressed in this work and solutions to each problem are given respectively. The experimental results conducted on a multi-cylinder diesel engine have shown that the proposed method successfully estimate a more accurate cylinder pressure over a wider range of crankshaft angles. Copyright © 2012 SAE International.

Contraction and observer design on cones

We consider the problem of positive observer design for positive systems defined on solid cones in Banach spaces. The design is based on the Hilbert metric and convergence properties are analyzed in the light of the Birkhoff theorem. Two main applications are discussed: positive observers for systems defined in the positive orthant, and positive observers on the cone of positive semi-definite matrices with a view on quantum systems. © 2011 IEEE.

Consensus in non-commutative spaces

Convergence analysis of consensus algorithms is revisited in the light of the Hilbert distance. The Lyapunov function used in the early analysis by Tsitsiklis is shown to be the Hilbert distance to consensus in log coordinates. Birkhoff theorem, which proves contraction of the Hilbert metric for any positive homogeneous monotone map, provides an early yet general convergence result for consensus algorithms. Because Birkhoff theorem holds in arbitrary cones, we extend consensus algorithms to the cone of positive definite matrices. The proposed generalization finds applications in the convergence analysis of quantum stochastic maps, which are a generalization of stochastic maps to non-commutative probability spaces. ©2010 IEEE.

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.

Large eddy simulations in low-pressure turbines: Effect of wakes at elevated free-stream turbulence

The transition of a separated shear layer over a flat plate, in the presence of periodic wakes and elevated free-stream turbulence (FST), is numerically investigated using Large Eddy Simulation (LES). The upper wall of the test section is inviscid and specifically contoured to impose a streamwise pressure distribution over the flat plate to simulate the suction surface of a low-pressure turbine (LPT) blade. Two different distributions representative of a 'high-lift' and an 'ultra high-lift' turbine blade are examined. Results obtained from the current LES compare favourably with the extensive experimental data previously obtained for these configurations. The LES results are then used to further investigate the flow physics involved in the transition process.In line with experimental experience, the benefit of wakes and FST obtained by suppressing the separation bubble, is more pronounced in 'ultra high-lift' design when compared to the 'high-lift' design. Stronger 'Klebanoff streaks' are formed in the presence of wakes when compared to the streaks due to FST alone. These streaks promoted much early transition. The weak Klebanoff streaks due to FST continued to trigger transition in between the wake passing cycles.The experimental inference regarding the origin of Klebanoff streaks at the leading edge has been confirmed by the current simulations. While the wake convects at local free-stream velocity, its impression in the boundary layer in the form of streaks convects much slowly. The 'part-span' Kelvin-Helmholtz structures, which were observed in the experiments when the wake passes over the separation bubble, are also captured. The non-phase averaged space-time plots manifest that reattachment is a localized process across the span unlike the impression of global reattachment portrayed by phase averaging. © 2013 Elsevier Inc.