Stability enhancement by casing grooves: The importance of stall inception mechanism and solidity

This paper concerns the optimisation of casing grooves and the important influence of stall inception mechanism on groove performance. Installing casing grooves is a well known technique for improving the stable operating range of a compressor, but the wide-spread use of grooves is restricted by the loss of efficiency and flow capacity. In this paper, laboratory tests are used to examine the conditions under which casing treatment can be used to greatest effect. The use of a single casing groove was investigated in a recently published companion paper. The current work extends this to multiple-groove treatments and considers their performance in relation to stall inception mechanisms. Here it is shown that the stall margin gain from multiple grooves is less than the sum of the gains if the grooves were used individually. By contrast, the loss of efficiency is additive as the number of grooves increases. It is then shown that casing grooves give the greatest stall margin improvement when used in a compressor which exhibits spike-type stall inception, while modal activity before stall can dramatically reduce the effectiveness of the grooves. This finding highlights the importance of being able to predict the stall inception mechanism which might occur in a given compressor before and after grooves are added. Some published prediction techniques are therefore examined, but found wanting. Lastly, it is shown that casing grooves can, in some cases, be used to remove rotor blades and produce a more efficient, stable and light-weight rotor. © 2010 by ASME.

Additive Gaussian processes

We introduce a Gaussian process model of functions which are additive. An additive function is one which decomposes into a sum of low-dimensional functions, each depending on only a subset of the input variables. Additive GPs generalize both Generalized Additive Models, and the standard GP models which use squared-exponential kernels. Hyperparameter learning in this model can be seen as Bayesian Hierarchical Kernel Learning (HKL). We introduce an expressive but tractable parameterization of the kernel function, which allows efficient evaluation of all input interaction terms, whose number is exponential in the input dimension. The additional structure discoverable by this model results in increased interpretability, as well as state-of-the-art predictive power in regression tasks.

Linear gaussian computations for near-exact Bayesian Monte Carlo inference in skewed alpha-stable time series models

In this paper we study parameter estimation for time series with asymmetric α-stable innovations. The proposed methods use a Poisson sum series representation (PSSR) for the asymmetric α-stable noise to express the process in a conditionally Gaussian framework. That allows us to implement Bayesian parameter estimation using Markov chain Monte Carlo (MCMC) methods. We further enhance the series representation by introducing a novel approximation of the series residual terms in which we are able to characterise the mean and variance of the approximation. Simulations illustrate the proposed framework applied to linear time series, estimating the model parameter values and model order P for an autoregressive (AR(P)) model driven by asymmetric α-stable innovations. © 2012 IEEE.

Verification of periodically controlled hybrid systems: Application to an autonomous vehicle

This article introduces Periodically Controlled Hybrid Automata (PCHA) for modular specification of embedded control systems. In a PCHA, control actions that change the control input to the plant occur roughly periodically, while other actions that update the state of the controller may occur in the interim. Such actions could model, for example, sensor updates and information received from higher-level planning modules that change the set point of the controller. Based on periodicity and subtangential conditions, a new sufficient condition for verifying invariant properties of PCHAs is presented. For PCHAs with polynomial continuous vector fields, it is possible to check these conditions automatically using, for example, quantifier elimination or sum of squares decomposition. We examine the feasibility of this automatic approach on a small example. The proposed technique is also used to manually verify safety and progress properties of a fairly complex planner-controller subsystem of an autonomous ground vehicle. Geometric properties of planner-generated paths are derived which guarantee that such paths can be safely followed by the controller. © 2012 ACM.

ℓasso MPC: Smart regulation of over-actuated systems

In this paper, a novel MPC strategy is proposed, and referred to as asso MPC. The new paradigm features an 1-regularised least squares loss function, in which the control error variance competes with the sum of input channels magnitude (or slew rate) over the whole horizon length. This cost choice is motivated by the successful development of LASSO theory in signal processing and machine learning. In the latter fields, sum-of-norms regularisation have shown a strong capability to provide robust and sparse solutions for system identification and feature selection. In this paper, a discrete-time dual-mode asso MPC is formulated, and its stability is proven by application of standard MPC arguments. The controller is then tested for the problem of ship course keeping and roll reduction with rudder and fins, in a directional stochastic sea. Simulations show the asso MPC to inherit positive features from its corresponding regressor: extreme reduction of decision variables' magnitude, namely, actuators' magnitude (or variations), with a finite energy error, being particularly promising for over-actuated systems. © 2012 AACC American Automatic Control Council).

Geometry and interaction of structures in homogeneous isotropic turbulence

A strategy to extract turbulence structures from direct numerical simulation (DNS) data is described along with a systematic analysis of geometry and spatial distribution of the educed structures. A DNS dataset of decaying homogeneous isotropic turbulence at Reynolds number Reλ = 141 is considered. A bandpass filtering procedure is shown to be effective in extracting enstrophy and dissipation structures with their smallest scales matching the filter width, L. The geometry of these educed structures is characterized and classified through the use of two non-dimensional quantities, planarity' and filamentarity', obtained using the Minkowski functionals. The planarity increases gradually by a small amount as L is decreased, and its narrow variation suggests a nearly circular cross-section for the educed structures. The filamentarity increases significantly as L decreases demonstrating that the educed structures become progressively more tubular. An analysis of the preferential alignment between the filtered strain and vorticity fields reveals that vortical structures of a given scale L are most likely to align with the largest extensional strain at a scale 3-5 times larger than L. This is consistent with the classical energy cascade picture, in which vortices of a given scale are stretched by and absorb energy from structures of a somewhat larger scale. The spatial distribution of the educed structures shows that the enstrophy structures at the 5η scale (where η is the Kolmogorov scale) are more concentrated near the ones that are 3-5 times larger, which gives further support to the classical picture. Finally, it is shown by analysing the volume fraction of the educed enstrophy structures that there is a tendency for them to cluster around a larger structure or clusters of larger structures. Copyright © 2012 Cambridge University Press.

Flame interactions in turbulent premixed twin V-flames

Multiple flame-flame interactions in premixed combustion are investigated using direct numerical simulations of twin turbulent V-flames for a range of turbulence intensities and length scales. Interactions are identified using a novel automatic feature extraction (AFE) technique, based on data registration using the dual-tree complex wavelet transform. Information on the time, position, and type of interactions, and their influence on the flame area is extracted using AFE. Characteristic length and time scales for the interactions are identified. The effect of interactions on the flame brush is quantified through a global stretch rate, defined as the sum of flamelet stretch and interaction stretch contributions. The effects of each interaction type are discussed. It is found that the magnitude of the fluctuations in flamelet and interaction stretch are comparable, and a qualitative sensitivity to turbulence length scale is found for one interaction type. Implications for modeling are discussed. © 2013 Copyright Taylor and Francis Group, LLC.