Opengazer: open-source gaze tracker for ordinary webcams

Opengazer is an open source application that uses an ordinary webcam to estimate head pose, facial gestures, or the direction of your gaze. This information can then be passed to other applications. For example, used in conjunction with Dasher, opengazer allows you to write with your eyes. Opengazer aims to be a low-cost software alternative to commercial hardware-based eye trackers. The first version of Opengazer was developed by Piotr Zieliński, supported by Samsung and the Gatsby Charitable Foundation. Research and development for Opengazer has been continued by Emli-Mari Nel, and was supported until 2012 by the European Commission in the context of the AEGIS project, and also by the Gatsby Charitable Foundation.

Surface-waves, stability, and scattering for a lined duct with flow

We consider a straight cylindrical duct with a steady subsonic axial flow and a reacting boundary (e.g. an acoustic lining). The wave modes are separated into ordinary acoustic duct modes, and surface modes confined to a small neighbourhood of the boundary. Many researchers have used a mass-spring-damper boundary model, for which one surface mode has previously been identified as a convective instability; however, we show the stability analysis used in such cases to be questionable. We investigate instead the stability of the surface modes using the Briggs-Bers criterion for a Flügge thin-shell boundary model. For modest frequencies and wavenumbers the thin-shell has an impedance which is effectively that of a mass-spring-damper, although for the large wavenumbers needed for the stability analysis the thin-shell and mass-spring-damper impedances diverge, owing to the thin shell's bending stiffness. The thin shell model may therefore be viewed as a regularization of the mass-spring-damper model which accounts for nonlocally-reacting effects. We find all modes to be stable for realistic thin-shell parameters, while absolute instabilities are demonstrated for extremely thin boundary thicknesses. The limit of vanishing bending stiffness is found to be a singular limit, yielding absolute instabilities of arbitrarily large temporal growth rate. We propose that the problems with previous stability analyses are due to the neglect of something akin to bending stiffness in the boundary model. Our conclusion is that the surface mode previously identified as a convective instability may well be stable in reality. Finally, inspired by Rienstra's recent analysis, we investigate the scattering of an acoustic mode as it encounters a sudden change from a hard-wall to a thin-shell boundary, using a Wiener-Hopf technique. The thin-shell is considered to be clamped to the hard-wall. The acoustic mode is found to scatter into transmitted and reflected acoustic modes, and surface modes strongly linked to the solid waves in the boundary, although no longitudinal or transverse waves within the boundary are excited. Examples are provided that demonstrate total transmission, total reflection, and a combination of the two. This thin-shell scattering problem is preferable to the mass-spring-damper scattering problem presented by Rienstra, since the thin-shell problem is fully determined and does not need to appeal to a Kutta-like condition or the inclusion of an instability in order to avoid a surface-streamline cusp at the boundary change.

Long-term changes and variability of monthly extreme temperatures in India

Extreme temperatures are changing worldwide together with changes in the mean temperatures. This study investigates the long-term trends and variations of the monthly maximum and minimum temperatures and their effects on seasonal fluctuations in various climatological regions in India. The magnitude of the trends and their statistical significance were determined by parametric ordinary least square regression techniques and the variations were determined by the respective coefficient of variations. The results showed that the monthly maximum temperature increased, though unevenly, over the last century. Minimum temperature changes were more variable than maximum temperature changes, both temporally and spatially, with results of lesser significance. The results of this study are good indicators of Indian climate variability and its changes over the last century. © Springer-Verlag 2009.

Assessing seasonal precipitation trends in India using parametric and non-parametric statistical techniques

This paper provides an insight into the long-term trends of the four seasonal and annual precipitations in various climatological regions and sub-regions in India. The trends were useful to investigate whether Indian seasonal rainfall is changing in terms of magnitude or location-wise. Trends were assessed over the period of 1954-2003 using parametric ordinary least square fits and non-parametric Mann-Kendall technique. The trend significance was tested at the 95% confidence level. Apart from the trends for individual climatological regions in India and the average for the whole of India, trends were also specifically determined for the possible smaller geographical areas in order to understand how different the trends would be from the bigger spatial scales. The smaller geographical regions consist of the whole southwestern continental state of Kerala. It was shown that there are decreasing trends in the spring and monsoon rainfall and increasing trends in the autumn and winter rainfalls. These changes are not always homogeneous over various regions, even in the very short scales implying a careful regional analysis would be necessary for drawing conclusions regarding agro-ecological or other local projects requiring change in rainfall information. Furthermore, the differences between the trend magnitudes and directions from the two different methods are significantly small and fall well within the significance limit for all the cases investigated in Indian regions (except where noted). © 2010 Springer-Verlag.

The post-processed Galerkin method applied to non-linear shell vibrations

We present the results of a computational study of the post-processed Galerkin methods put forward by Garcia-Archilla et al. applied to the non-linear von Karman equations governing the dynamic response of a thin cylindrical panel periodically forced by a transverse point load. We spatially discretize the shell using finite differences to produce a large system of ordinary differential equations (ODEs). By analogy with spectral non-linear Galerkin methods we split this large system into a 'slowly' contracting subsystem and a 'quickly' contracting subsystem. We then compare the accuracy and efficiency of (i) ignoring the dynamics of the 'quick' system (analogous to a traditional spectral Galerkin truncation and sometimes referred to as 'subspace dynamics' in the finite element community when applied to numerical eigenvectors), (ii) slaving the dynamics of the quick system to the slow system during numerical integration (analogous to a non-linear Galerkin method), and (iii) ignoring the influence of the dynamics of the quick system on the evolution of the slow system until we require some output, when we 'lift' the variables from the slow system to the quick using the same slaving rule as in (ii). This corresponds to the post-processing of Garcia-Archilla et al. We find that method (iii) produces essentially the same accuracy as method (ii) but requires only the computational power of method (i) and is thus more efficient than either. In contrast with spectral methods, this type of finite-difference technique can be applied to irregularly shaped domains. We feel that post-processing of this form is a valuable method that can be implemented in computational schemes for a wide variety of partial differential equations (PDEs) of practical importance.

Adaptive lenses based on polarization modulation

We present and demonstrate a technique for producing a high-speed variable focus lens using a fixed birefringent lens and a ferroelectric liquid crystal cell as a polarization switch. A calcite lenses with ordinary and extraordinary focal lengths of 109mm and 88mm respectively, was used to demonstrate focus switching at frequencies of up to 3kHz. Two identical lenses and a single liquid crystal were also used to demonstrate zoom.

Cortical thickness mapping to identify focal osteoporosis in patients with hip fracture.

BACKGROUND: Individuals with osteoporosis are predisposed to hip fracture during trips, stumbles or falls, but half of all hip fractures occur in those without generalised osteoporosis. By analysing ordinary clinical CT scans using a novel cortical thickness mapping technique, we discovered patches of markedly thinner bone at fracture-prone regions in the femurs of women with acute hip fracture compared with controls. METHODS: We analysed CT scans from 75 female volunteers with acute fracture and 75 age- and sex-matched controls. We classified the fracture location as femoral neck or trochanteric before creating bone thickness maps of the outer 'cortical' shell of the intact contra-lateral hip. After registration of each bone to an average femur shape and statistical parametric mapping, we were able to visualise and quantify statistically significant foci of thinner cortical bone associated with each fracture type, assuming good symmetry of bone structure between the intact and fractured hip. The technique allowed us to pinpoint systematic differences and display the results on a 3D average femur shape model. FINDINGS: The cortex was generally thinner in femoral neck fracture cases than controls. More striking were several discrete patches of statistically significant thinner bone of up to 30%, which coincided with common sites of fracture initiation (femoral neck or trochanteric). INTERPRETATION: Femoral neck fracture patients had a thumbnail-sized patch of focal osteoporosis at the upper head-neck junction. This region coincided with a weak part of the femur, prone to both spontaneous 'tensile' fractures of the femoral neck, and as a site of crack initiation when falling sideways. Current hip fracture prevention strategies are based on case finding: they involve clinical risk factor estimation to determine the need for single-plane bone density measurement within a standard region of interest (ROI) of the femoral neck. The precise sites of focal osteoporosis that we have identified are overlooked by current 2D bone densitometry methods.

PH-dependent leaching behaviour and other performance properties of cement-treated mixed contaminated soil

Portland cement has been widely used for stabilisation/solidification (S/S) treatment of contaminated soils. However, there is a dearth of literature on pH-dependent leaching of contaminants from cement-treated soils. This study investigates the leachability of Cu, Pb, Ni, Zn and total petroleum hydrocarbons (TPH) from a mixed contaminated soil. A sandy soil was spiked with 3000 mg/kg each of Cd, Cu, Pb, Ni and Zn, and 10,000 mg/kg of diesel, and treated with ordinary Portland cement (CEM I). Four different binder dosages, 5%, 10%, 15% and 20% (m/m) and different water contents ranging from 13%-19% dry weight were used in order to find a safe operating envelope for the treatment process. The pH-dependent leaching behaviour of the treated soil was monitored over an 84-day period using a 3-point acid neutralisation capacity (ANC) test. The monolithic leaching test was also conducted. Geotechnical properties such as unconfined compressive strength (UCS), hydraulic conductivity and porosity were assessed over time. The treated soils recorded lower leachate concentrations of Ni and Zn compared to the untreated soil at the same pH depending on binder dosage. The binder had problems with Pb stabilisation and TPH leachability was independent of pH and binder dosage. The hydraulic conductivity of the mixes was generally of the order, 10-8 m/sec, while the porosity ranged from 26%-44%. The results of selected performance properties are compared with regulatory limits and the range of operating variables that lead to acceptable performance described. © 2012 The Research Centre for Eco-Environmental Sciences, Chinese Academy of Sciences.

Classical plume theory: 1937-2010 and beyond

Developing a theoretical description of turbulent plumes, the likes of which may be seen rising above industrial chimneys, is a daunting thought. Plumes are ubiquitous on a wide range of scales in both the natural and the man-made environments. Examples that immediately come to mind are the vapour plumes above industrial smoke stacks or the ash plumes forming particle-laden clouds above an erupting volcano. However, plumes also occur where they are less visually apparent, such as the rising stream of warmair above a domestic radiator, of oil from a subsea blowout or, at a larger scale, of air above the so-called urban heat island. In many instances, not only the plume itself is of interest but also its influence on the environment as a whole through the process of entrainment. Zeldovich (1937, The asymptotic laws of freely-ascending convective flows. Zh. Eksp. Teor. Fiz., 7, 1463-1465 (in Russian)), Batchelor (1954, Heat convection and buoyancy effects in fluids. Q. J. R. Meteor. Soc., 80, 339-358) and Morton et al. (1956, Turbulent gravitational convection from maintained and instantaneous sources. Proc. R. Soc. Lond. A, 234, 1-23) laid the foundations for classical plume theory, a theoretical description that is elegant in its simplicity and yet encapsulates the complex turbulent engulfment of ambient fluid into the plume. Testament to the insight and approach developed in these early models of plumes is that the essential theory remains unchanged and is widely applied today. We describe the foundations of plume theory and link the theoretical developments with the measurements made in experiments necessary to close these models before discussing some recent developments in plume theory, including an approach which generalizes results obtained separately for the Boussinesq and the non-Boussinesq plume cases. The theory presented - despite its simplicity - has been very successful at describing and explaining the behaviour of plumes across the wide range of scales they are observed. We present solutions to the coupled set of ordinary differential equations (the plume conservation equations) that Morton et al. (1956) derived from the Navier-Stokes equations which govern fluid motion. In order to describe and contrast the bulk behaviour of rising plumes from general area sources, we present closed-form solutions to the plume conservation equations that were achieved by solving for the variation with height of Morton's non-dimensional flux parameter Γ - this single flux parameter gives a unique representation of the behaviour of steady plumes and enables a characterization of the different types of plume. We discuss advantages of solutions in this form before describing extensions to plume theory and suggesting directions for new research. © 2010 The Author. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.

Reconstruction of arbitrary biochemical reaction networks: A compressive sensing approach

Reconstruction of biochemical reaction networks (BRN) and genetic regulatory networks (GRN) in particular is a central topic in systems biology which raises crucial theoretical challenges in system identification. Nonlinear Ordinary Differential Equations (ODEs) that involve polynomial and rational functions are typically used to model biochemical reaction networks. Such nonlinear models make the problem of determining the connectivity of biochemical networks from time-series experimental data quite difficult. In this paper, we present a network reconstruction algorithm that can deal with ODE model descriptions containing polynomial and rational functions. Rather than identifying the parameters of linear or nonlinear ODEs characterised by pre-defined equation structures, our methodology allows us to determine the nonlinear ODEs structure together with their associated parameters. To solve the network reconstruction problem, we cast it as a compressive sensing (CS) problem and use sparse Bayesian learning (SBL) algorithms as a computationally efficient and robust way to obtain its solution. © 2012 IEEE.

Seafloor interaction with steel catenary risers

This paper describes the key features of a seafloor-riser interaction model. The soil is represented in terms of non-linear load-deflection (P- y) relationships, which are also able to account for soil stiffness degradation due to cyclic loading. The analytical framework considers the riser-seafloor interaction problem in terms of a pipe resting on a bed of springs, and requires the iterative solution of a fourth-order ordinary differential equation. A series of simulations is used to illustrate the capabilities of the model. Thanks to the non-linear soil springs with stiffness degradation it is possible to simulate the trench formation process and estimate moments in a riser. Copyright © 2008 by The International Society of Offshore and Polar Engineers (ISOPE).

Optimization algorithms and ODE's in MDO

There is a need for a stronger theoretical understanding of Multidisciplinary Design Optimization (MDO) within the field. Having developed a differential geometry framework in response to this need, we consider how standard optimization algorithms can be modeled using systems of ordinary differential equations (ODEs) while also reviewing optimization algorithms which have been derived from ODE solution methods. We then use some of the framework's tools to show how our resultant systems of ODEs can be analyzed and their behaviour quantitatively evaluated. In doing so, we demonstrate the power and scope of our differential geometry framework, we provide new tools for analyzing MDO systems and their behaviour, and we suggest hitherto neglected optimization methods which may prove particularly useful within the MDO context. Copyright © 2013 by ASME.