A new method for the acquisition of ultrasonic strain image volumes.

This article presents a new method for acquiring three-dimensional (3-D) volumes of ultrasonic axial strain data. The method uses a mechanically-swept probe to sweep out a single volume while applying a continuously varying axial compression. Acquisition of a volume takes 15-20 s. A strain volume is then calculated by comparing frame pairs throughout the sequence. The method uses strain quality estimates to automatically pick out high quality frame pairs, and so does not require careful control of the axial compression. In a series of in vitro and in vivo experiments, we quantify the image quality of the new method and also assess its ease of use. Results are compared with those for the current best alternative, which calculates strain between two complete volumes. The volume pair approach can produce high quality data, but skillful scanning is required to acquire two volumes with appropriate relative strain. In the new method, the automatic quality-weighted selection of image pairs overcomes this difficulty and the method produces superior quality images with a relatively relaxed scanning technique.

Flow field measurements of a series of turbulent premixed and stratified methane/air flames

This paper presents flow field measurements for the turbulent stratified burner introduced in two previous publications in which high resolution scalar measurements were made by Sweeney et al. [1,2] for model validation. The flow fields of the series of premixed and stratified methane/air flames are investigated under turbulent, globally lean conditions (φg=0.75). Velocity data acquired with laser Doppler anemometry (LDA) and particle image velocimetry (PIV) are presented and discussed. Pairwise 2-component LDA measurements provide profiles of axial velocity, radial velocity, tangential velocity and corresponding fluctuating velocities. The LDA measurements of axial and tangential velocities enable the swirl number to be evaluated and the degree of swirl characterized. Power spectral density and autocorrelation functions derived from the LDA data acquired at 10kHz are optimized to calculate the integral time scales. Flow patterns are obtained using a 2-component PIV system operated at 7Hz. Velocity profiles and spatial correlations derived from the PIV and LDA measurements are shown to be in very good agreement, thus offering 3D mapping of the velocities. A strong correlation was observed between the shape of the recirculation zones above the central bluff body and the effects of heat release, stoichiometry and swirl. Detailed analyses of the LDA data further demonstrate that the flow behavior changes significantly with the levels of swirl and stratification, which combines the contributions of dilatation, recirculation and swirl. Key turbulence parameters are derived from the total velocity components, combining axial, radial and tangential velocities. © 2013 The Combustion Institute.

High resolution cortical bone thickness measurement from clinical CT data.

The distribution of cortical bone in the proximal femur is believed to be a critical component in determining fracture resistance. Current CT technology is limited in its ability to measure cortical thickness, especially in the sub-millimetre range which lies within the point spread function of today's clinical scanners. In this paper, we present a novel technique that is capable of producing unbiased thickness estimates down to 0.3mm. The technique relies on a mathematical model of the anatomy and the imaging system, which is fitted to the data at a large number of sites around the proximal femur, producing around 17,000 independent thickness estimates per specimen. In a series of experiments on 16 cadaveric femurs, estimation errors were measured as -0.01+/-0.58mm (mean+/-1std.dev.) for cortical thicknesses in the range 0.3-4mm. This compares with 0.25+/-0.69mm for simple thresholding and 0.90+/-0.92mm for a variant of the 50% relative threshold method. In the clinically relevant sub-millimetre range, thresholding increasingly fails to detect the cortex at all, whereas the new technique continues to perform well. The many cortical thickness estimates can be displayed as a colour map painted onto the femoral surface. Computation of the surfaces and colour maps is largely automatic, requiring around 15min on a modest laptop computer.

Image deconvolution using a Gaussian scale mixtures model to approximate the wavelet sparseness constraint

This paper proposes to use an extended Gaussian Scale Mixtures (GSM) model instead of the conventional ℓ1 norm to approximate the sparseness constraint in the wavelet domain. We combine this new constraint with subband-dependent minimization to formulate an iterative algorithm on two shift-invariant wavelet transforms, the Shannon wavelet transform and dual-tree complex wavelet transform (DTCWT). This extented GSM model introduces spatially varying information into the deconvolution process and thus enables the algorithm to achieve better results with fewer iterations in our experiments. ©2009 IEEE.

Reversible jump sampler for autoregressive time series

We use reversible jump Markov chain Monte Carlo (MCMC) methods to address the problem of model order uncertainty in autoregressive (AR) time series within a Bayesian framework. Efficient model jumping is achieved by proposing model space moves from the full conditional density for the AR parameters, which is obtained analytically. This is compared with an alternative method, for which the moves are cheaper to compute, in which proposals are made only for new parameters in each move. Results are presented for both synthetic and audio time series.

Bayesian model selection for time series using Markov chain Monte Carlo

We present a stochastic simulation technique for subset selection in time series models, based on the use of indicator variables with the Gibbs sampler within a hierarchical Bayesian framework. As an example, the method is applied to the selection of subset linear AR models, in which only significant lags are included. Joint sampling of the indicators and parameters is found to speed convergence. We discuss the possibility of model mixing where the model is not well determined by the data, and the extension of the approach to include non-linear model terms.

Microstructural modelling in friction stir welding of 2000 series aluminium alloys

Simple process models are applied to predict microstructural changes due to the thermal cycle imposed in friction stir welding. A softening model developed for heat-treatable aluminium alloys of the 6000 series is applied to the aerospace alloy 2014 in the peak-aged (T6) condition. It is found that the model is not readily applicable to alloy 2024 in the naturally aged (T3) temper, but the softening behaviour can still be described semi-empirically. Both analytical and numerical (finite element) thermal models are used to predict the thermal histories in trial welds. These are coupled to the microstructural model to investigate: (a) the hardness profile across the welded plate; (b) alloy softening ahead of the approaching welding tool. By incorporating the softening model applied to 6082-T6 alloy, the hardness profile of friction stir welds in dissimilar alloys is also predicted. © AFM, EDP Sciences 2005.

Uplift resistance of buried pipelines at low cover-diameter ratios

Reliable estimates for the maximum available uplift resistance from the backfill soil are essential to prevent upheaval buckling of buried pipelines. The current design code DNV RP F110 does not offer guidance on how to predict the uplift resistance when the cover:pipe diameter (H/D) ratio is less than 2. Hence the current industry practice is to discount the shear contribution from uplift resitance for design scenarios with H/D ratios less than 1. The necessity of this extra conservatism is assessed through a series of full-scale and centrifuge tests, 21 in total, at the Schofield Centre, University of Cambridge. Backfill types include saturated loose sand, saturated dense sand and dry gravel. Data revealed that the Vertical Slip Surface Model remains applicable for design scenarios in loose sand, dense sand and gravel with H/D ratios less than 1, and that there is no evidence that the contribution from shear should be ignored at these low H/D ratios. For uplift events in gravel, the shear component seems reliable if the cover is more than 1-2 times the average particle size (D50), and more research effort is currenty being carried out to verify this conclusion. Strain analysis from the Particle Image Velocimetry (PIV) technique proves that the Vertical Slip Surface Model is a good representation of the true uplift deformation mechanism in loose sand at H/D ratios between 0.5 and 3.5. At very low H/D ratios (H/D < 0.5), the deformation mechanism is more wedge-like, but the increased contribution from soil weight is likely to be compensated by the reduced shear contributions. Hence the design equation based on the Vertical Slip Surface Model still produces good estimates for the maximum available uplift resistance. The evolution of shear strain field from PIV analysis provides useful insight into how uplift resistance is mobilized as the uplift event progresses. Copyright 2010, Offshore Technology Conference.