Computational fluid-structure interaction of DGB parachutes in compressible fluid flow

This paper describes large-scale simulations of compressible flows over a supersonic disk-gap-band parachute system. An adaptive mesh refinement method is used to resolve the coupled fluid-structure model. The fluid model employs large-eddy simulation to describe the turbulent wakes appearing upstream and downstream of the parachute canopy and the structural model employed a thin-shell finite element solver that allows large canopy deformations by using subdivision finite elements. The fluid-structure interaction is described by a variant of the Ghost-Fluid method. The simulation was carried out at Mach number 1.96 where strong nonlinear coupling between the system of bow shocks, turbulent wake and canopy is observed. It was found that the canopy oscillations were characterized by a breathing type motion due to the strong interaction of the turbulent wake and bow shock upstream of the flexible canopy. Copyright © 2010 by ASME.

Modelling of liquefaction-induced instability in pile

Centrifuge testing has been undertaken to investigate instability failure of pile groups during seismic liquefaction, with specific reference to the 'top-down' propagation of liquefaction during the earthquake and to account for initial imperfections in pile geometry. The results of these tests were used to validate numerical models within the finite element program ABAQUS, based on the popular p-y analysis method. Pseudostatic classical and post-buckling analyses were conducted to examine the collapse behaviour of the pile groups and were found to give reasonable predictions of collapse load and conservative predictions of the associated deflection conditions. This numerical model was compared to currently published methods which were found to over-predict collapse loads. The resulting insights into the collapse of axially loaded pile groups revealed that the failure load is strongly dependent on both the depth of liquefaction propagation and initial imperfections, which reduce the collapse load.

A unifying resolution-independent formulation for early vision

We present a model for early vision tasks such as denoising, super-resolution, deblurring, and demosaicing. The model provides a resolution-independent representation of discrete images which admits a truly rotationally invariant prior. The model generalizes several existing approaches: variational methods, finite element methods, and discrete random fields. The primary contribution is a novel energy functional which has not previously been written down, which combines the discrete measurements from pixels with a continuous-domain world viewed through continous-domain point-spread functions. The value of the functional is that simple priors (such as total variation and generalizations) on the continous-domain world become realistic priors on the sampled images. We show that despite its apparent complexity, optimization of this model depends on just a few computational primitives, which although tedious to derive, can now be reused in many domains. We define a set of optimization algorithms which greatly overcome the apparent complexity of this model, and make possible its practical application. New experimental results include infinite-resolution upsampling, and a method for obtaining subpixel superpixels. © 2012 IEEE.

Experimental and numerical study of grout injections in silty soils

Compensation grouting is increasingly employed as a mitigation technique of settlements induced by tunnelling and its effectiveness both in clayey and sandy soils is reported in a wide number of case histories. However, the results are highly dependent on grout properties, injection characteristics and soil properties. An experimental study was conducted to investigate the parameters that control grout injections in silty soils. The results from one injection test in a large sample of silty soil show that the compensation efficiency, defined as the ratio of the volume of heave obtained at ground surface and the injected grout volume, is much lower than one and tends to decrease with time, while the initial volume of grout lost due to pressure filtration is small. Finally, results from finite elements back analyses of the laboratory test show that a good agreement with the experimental data can be obtained if the development of large strains is taken into account. © 2012 Taylor & Francis Group.

Models for acoustic propagation through turbofan exhaust flows

Turbomachinery noise radiating into the rearward arc is an important problem. This noise is scattered by the trailing edges of the nacelle and the jet exhaust, and interacts with the shear layers between the external flow, bypass stream and jet, en route to the far field. In the past a range of relevant model problems involving semi-infinite cylinders have been solved. However, one limitation of previous solutions is that they do not allow for the jet nozzle to protrude a finite distance beyond the end of the nacelle (or in certain configurations being buried a finite distance upstream). In this paper we use the matrix Wiener-Hopf technique, which will allow precisely the finite nacelle-jet nozzle separation to be included. The crucial step in our work is to factorise a certain matrix as a product of terms analytic and invertible in the upper/lower halves of the complex plane. The way we do this matrix factorisation is quite different in the buried and protruding nozzle cases. In the buried case our solution method is the so-called pole-removal technique. In the technically more demanding protruding case, however, we must first use Pade approximants to generate a uniformly-valid, meromorphic representation of a certain function, before the same pole-removal method can be applied. Sample results are presented, investigating in particular the effects of exit plane stagger. © 2007 by B Veitch and N Peake.

Free energy surface reconstruction from umbrella samples using Gaussian process regression II: Multiple collective variables

We demonstrate how a prior assumption of smoothness can be used to enhance the reconstruction of free energy profiles from multiple umbrella sampling simulations using the Bayesian Gaussian process regression approach. The method we derive allows the concurrent use of histograms and free energy gradients and can easily be extended to include further data. In Part I we review the necessary theory and test the method for one collective variable. We demonstrate improved performance with respect to the weighted histogram analysis method and obtain meaningful error bars without any significant additional computation. In Part II we consider the case of multiple collective variables and compare to a reconstruction using least squares fitting of radial basis functions. We find substantial improvements in the regimes of spatially sparse data or short sampling trajectories. A software implementation is made available on www.libatoms.org.

Numerical modelling of hydrodynamic pressures on dams

This paper discusses several considerations related to appropriate numerical modelling of the reservoir hydrodynamic pressures on dams. The reservoir is modelled with 8-noded isoparametric displacement based solid finite elements. The study includes both stiff and flexible dams with vertical and sloped upstream faces under ramp, harmonic and random acceleration loads. The numerical results were compared and found to be in good agreement with available closed-form solutions. The same approach may be used in analyses of other waterfront structures such as quay walls. © 2013 Elsevier Ltd.

Anisotropy preserving DTI processing

Statistical analysis of diffusion tensor imaging (DTI) data requires a computational framework that is both numerically tractable (to account for the high dimensional nature of the data) and geometric (to account for the nonlinear nature of diffusion tensors). Building upon earlier studies exploiting a Riemannian framework to address these challenges, the present paper proposes a novel metric and an accompanying computational framework for DTI data processing. The proposed approach grounds the signal processing operations in interpolating curves. Well-chosen interpolating curves are shown to provide a computational framework that is at the same time tractable and information relevant for DTI processing. In addition, and in contrast to earlier methods, it provides an interpolation method which preserves anisotropy, a central information carried by diffusion tensor data. © 2013 Springer Science+Business Media New York.