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.

Some characteristics of thin shear layers in homogeneous turbulent flow.

Tomographic particle image velocimetry measurements of homogeneous isotropic turbulence that have been made in a large mixing tank facility at Cambridge are analysed in order to characterize thin highly sheared regions that have been observed. The results indicate that such regions coincide with regions of high enstrophy, dissipation and stretching. Large velocity jumps are observed across the width of these regions. The thickness of the shear layers seems to scale with the Taylor microscale, as has been suggested previously. The results discussed here concentrate on examining individual realizations rather than statistics of these regions.

Modelling of waterfront cantilever sheet pile walls subjected to earthquake loading

The seismic performance of waterfront cantilever sheet pile retaining walls is of continuing interest to geotechnical engineers as these structures suffer severe damage and even complete failure during earthquakes. This is often precipitated by liquefaction of the surrounding soil, either in the backfill or in front of the wall. This paper presents results from a series of small-scale plane strain models that were tested on a 1-g shaking table and recorded using a high-speed, high-resolution digital camera. The technique of Particle Image Velocimetry (PIV) was applied in order to allow the failure mechanisms to be visualised. It is shown that using PIV analyses it is possible to obtain failure mechanisms for a cantilever wall in liquefiable soil. These failure mechanisms are compared with those obtained for a cantilever wall in dry soil, previously carried out at a similar scale. It was observed that seismic liquefaction causes significant displacement in much larger zones of soil near the retaining wall compared to an equivalent dry case. The failure mechanism for a cantilever wall with liquefiable backfill, but with a remediated zone designed not to liquefy, is also presented and compared to the unremediated case.

Three-dimensional effects on sliding and waving wings

Like large insects, micro air vehicles operate at low Reynolds numbers O(1; 000 - 10; 000) in a regime characterized by separated flow and strong vortices. The leading-edge vortex has been identified as a significant source of high lift on insect wings, but the conditions required for the formation of a stably attached leading-edge vortex are not yet known. The waving wing is designed to model the translational phase of an insect wing stroke by preserving the unsteady starting and stopping motion as well as three-dimensionality in both wing geometry (via a finite-span wing) and kinematics (via wing rotation). The current study examines the effect of the spanwise velocity gradient on the development of the leading-edge vortex along the wing as well as the effects of increasing threedimensionalityby decreasing wing aspect ratio from four to two. Dye flow visualization and particle image velocimetry reveal that the leading-edge vortices that form on a sliding or waving wing have a very high aspect ratio. The structure of the flow is largely two-dimensional on both sliding and waving wings and there is minimal interaction between the leading-edge vortices and the tip vortex. Significant spanwise flow was observed on the waving wing but not on the sliding wing. Despite the increased three-dimensionality on the aspect ratio 2 waving wing, there is no evidence of an attached leading-edge vortex and the structure of the flow is very similar to that on the higher-aspect-ratio wing and sliding wing. © Copyright 2010.

Pipeline uplift response in thawed sandy backfills: Preliminary findings

Offshore and onshore buried pipelines under high operating temperature and pressures may lead to upheaval buckling (UHB) if sufficient soil cover is not present to prevent the upward movement of the pipeline. In regions where seasonal changes involve ground soil undergoing freezing-thawing cycles, the uplift resistance from soil cover may be minimum when the soil is undergoing thawing. This paper presents the results from 2 directly-comparable minidrum centrifuge tests conducted at the Schofield Centre, University of Cambridge, to investigate the difference in uplift resistance responses between fully-saturated and thawed sandy backfill conditions. Both tests were conducted drained at 30g using an 8.6 mm diameter aluminium model pipe, corresponding to a prototype pipe diameter of 258 mm. The soil cover/pipe diameter ratio, H/D, was kept at 1. Fraction E fine silica sand was used as the backfill. Preliminary experimental results indicated that the ultimate uplift resistance of a thawing sand backfill to be lower than that of a fully saturated sand backfill. This suggests that in regions where backfill soil undergoes freeze-thaw cycles, the thawing backfill may be more critical than fully saturated backfill for uplift resistance. The 2-dimensional displacement field during the experiment was accurately measured and analysed using the Particle Image Velocimetry technique. Copyright © 2011 by the International Society of Offshore and Polar Engineers (ISOPE).

Instabilities in the oscillatory flow of a complex fluid.

The dynamics of a fluid in a vertical tube, subjected to an oscillatory pressure gradient, is studied experimentally for both a Newtonian and a viscoelastic shear-thinning fluid. Particle image velocimetry is used to determine the two-dimensional velocity fields in the vertical plane of the tube axis, in a range of driving amplitudes from 0.8 to 2.5 mm and of driving frequencies from 2.0 to 11.5 Hz. The Newtonian fluid exhibits a laminar flow regime, independent of the axial position, in the whole range of drivings. For the complex fluid, instead, the parallel shear flow regime exhibited at low amplitudes [Torralba, Phys. Rev. E 72, 016308 (2005)] becomes unstable at higher drivings against the formation of symmetric vortices, equally spaced along the tube. At even higher drivings the vortex structure itself becomes unstable, and complex nonsymmetric structures develop. Given that inertial effects remain negligible even at the hardest drivings (Re < 10(-1)), it is the complex rheology of the fluid that is responsible for the instabilities observed. The system studied represents an interesting example of the development of shear-induced instabilities in nonlinear complex fluids in purely parallel shear flow.

Experimental demonstration of the Rayleigh acoustic viscous boundary layer theory.

Amplitude and phase velocity measurements on the laminar oscillatory viscous boundary layer produced by acoustic waves are presented. The measurements were carried out in acoustic standing waves in air with frequencies of 68.5 and 114.5 Hz using laser Doppler anemometry and particle image velocimetry. The results obtained by these two techniques are in good agreement with the predictions made by the Rayleigh viscous boundary layer theory and confirm the existence of a local maximum of the velocity amplitude and its expected location.

2D PIV measurements in the near field of grid turbulence using stitched fields from multiple cameras

We present measurements of grid turbulence using 2D particle image velocimetry taken immediately downstream from the grid at a Reynolds number of Re M = 16500 where M is the rod spacing. A long field of view of 14M x 4M in the down- and cross-stream directions was achieved by stitching multiple cameras together. Two uniform biplanar grids were selected to have the same M and pressure drop but different rod diameter D and crosssection. A large data set (10 4 vector fields) was obtained to ensure good convergence of second-order statistics. Estimations of the dissipation rate ε of turbulent kinetic energy (TKE) were found to be sensitive to the number of meansquared velocity gradient terms included and not whether the turbulence was assumed to adhere to isotropy or axisymmetry. The resolution dependency of different turbulence statistics was assessed with a procedure that does not rely on the dissipation scale η. The streamwise evolution of the TKE components and ε was found to collapse across grids when the rod diameter was included in the normalisation. We argue that this should be the case between all regular grids when the other relevant dimensionless quantities are matched and the flow has become homogeneous across the stream. Two-point space correlation functions at x/M = 1 show evidence of complex wake interactions which exhibit a strong Reynolds number dependence. However, these changes in initial conditions disappear indicating rapid cross-stream homogenisation. On the other hand, isotropy was, as expected, not found to be established by x/M = 12 for any case studied. © Springer-Verlag 2012.

Statistics and dynamics of turbulence-flame alignment in premixed combustion

The geometric alignment of turbulent strain-rate structures with premixed flames greatly influences the results of the turbulence-flame interaction. Here, the statistics and dynamics of this alignment are experimentally investigated in turbulent premixed Bunsen flames using high-repetition-rate stereoscopic particle image velocimetry. In all cases, the statistics showed that the most extensive principal strain-rate associated with the turbulence preferentially aligned such that it was more perpendicular than parallel to the flame surface normal direction. The mean turbulence-flame alignment differed between the flames, with the stronger flames (higher laminar flame speed) exhibiting stronger preferential alignment. Furthermore, the preferential alignment was greatest on the reactant side of the mean flame brush. To understand these differences, individual structures of fluid-dynamic strain-rate were tracked through time in a Lagrangian manner (i.e., by following the fluid elements). It was found that the flame surface affected the orientation of the turbulence structures, with the majority of structures rotating as they approached the flame such that their most extensive principal strain-rate was perpendicular to the flame normal. The maximum change in turbulent structure orientation was found to decrease with the strength of the structure, increase with the strength of the flame, and exhibit similar trends when the structure strength and flame strength were represented by a Karlovitz number. The mean change in orientation decreased from the unburnt to burnt side of the flame brush and appears to be influenced by the overall flame shape. © 2011 The Combustion Institute.

A new apparatus for modelling excavations

Underground space is commonly exploited both to maximise the utility of costly land in urban development and to reduce the vertical load acting on the ground. Deep excavations are carried out to construct various types of underground infrastructure such as deep basements, subways and service tunnels. Although the soil response to excavation is known in principle, designers lack practical calculation methods for predicting both short- and long-term ground movements. As the understanding of how soil behaves around an excavation in both the short and long term is insufficient and usually empirical, the judgements used in design are also empirical and serious accidents are common. To gain a better understanding of the mechanisms involved in soil excavation, a new apparatus for the centrifuge model testing of deep excavations in soft clay has been developed. This apparatus simulates the field construction sequence of a multi-propped retaining wall during centrifuge flight. A comparison is given between the new technique and the previously used method of draining heavy fluid to simulate excavation in a centrifuge model. The new system has the benefit of giving the correct initial ground conditions before excavation and the proper earth pressure distribution on the retaining structures during excavation, whereas heavy fluid only gives an earth pressure coefficient of unity and is unable to capture any changes in the earth pressure coefficient of soil inside the zone of excavation, for example owing to wall movements. Settlements of the ground surface, changes in pore water pressure, variations in earth pressure, prop forces and bending moments in the retaining wall are all monitored during excavation. Furthermore, digital images taken of a cross-section during the test are analysed using particle image velocimetry to illustrate ground deformation and soil-structure interaction mechanisms. The significance of these observations is discussed.

The structure of turbulent stratified and premixed methane/air flames I: Non-swirling flows

A series of flames in a turbulent methane/air stratified swirl burner is presented. The degree of stratification and swirl are systematically varied to generate a matrix of experimental conditions, allowing their separate and combined effects to be investigated. Non-swirling flows are considered in the present paper, and the effects of swirl are considered in a companion paper (Part II). A mean equivalence ratio of φ=0.75 is used, with φ for the highest level of stratification spanning 0.375-1.125. The burner features a central bluff-body to aid flame stabilization, and the influence of the induced recirculation zone is also considered. The current work focuses on non-swirling flows where two-component particle image velocimetry (PIV) measurements are sufficient to characterize the main features of the flow field. Scalar data obtained from Rayleigh/Raman/CO laser induced fluorescence (CO-LIF) line measurements at 103μm resolution allow the behavior of key combustion species-CH 4, CO 2, CO, H 2, H 2O and O 2-to be probed within the instantaneous flame front. Simultaneous cross-planar OH-PLIF is used to determine the orientation of the instantaneous flame normal in the scalar measurement window, allowing gradients in temperature and progress variable to be angle corrected to their three dimensional values. The relationship between curvature and flame thickness is investigated using the OH-PLIF images, as well as the effect of stratification on curvature.The main findings are that the behavior of the key combustion species in temperature space is well captured on the mean by laminar flame calculations regardless of the level of stratification. H 2 and CO are significant exceptions, both appearing at elevated levels in the stratified flames. Values for surface density function and by extension thermal scalar dissipation rate are found to be substantially lower than laminar values, as the thickening of the flame due to turbulence dominates the effect of increased strain. These findings hold for both premixed and stratified flames. The current series of flames is proposed as an interesting if challenging set of test cases for existing and emerging turbulent flame models, and data are available on request. © 2012 The Combustion Institute.

Development of a new apparatus for modeling deep excavation related problems in geotechnical centrifuge

Underground space is commonly exploited both to maximise the utility of costly land in urban development and to reduce the vertical load acting on the ground. Deep excavations are carried out to construct various types of underground infrastructure such as deep basements, subways and service tunnels. Although the soil response to excavation is known in principle, designers lack practical calculation methods for predicting both short- and long-term ground movements. As the understanding of how soil behaves around an excavation in both the short and long term is insufficient and usually empirical, the judgements used in design are also empirical and serious accidents are common. To gain a better understanding of the mechanisms involved in soil excavation, a new apparatus for the centrifuge model testing of deep excavations in soft clay has been developed. This apparatus simulates the field construction sequence of a multi-propped retaining wall during centrifuge flight. A comparison is given between the new technique and the previously used method of draining heavy fluid to simulate excavation in a centrifuge model. The new system has the benefit of giving the correct initial ground conditions before excavation and the proper earth pressure distribution on the retaining structures during excavation, whereas heavy fluid only gives an earth pressure coefficient of unity and is unable to capture any changes in the earth pressure coefficient of soil inside the zone of excavation, for example owing to wall movements. Settlements of the ground surface, changes in pore water pressure, variations in earth pressure, prop forces and bending moments in the retaining wall are all monitored during excavation. Furthermore, digital images taken of a cross-section during the test are analysed using particle image velocimetry to illustrate ground deformation and soil–structure interaction mechanisms. The significance of these observations is discussed.

Large-eddy simulations and measurements of a small-scale high-speed coflowing jet

Measurements and predictions are made of a short-cowl coflowing jet with a bypass ratio of 8:1. The Reynolds number is 300,000, and the inlet Mach numbers are representative of aeroengine conditions. The low Reynolds number of the measurements makes the case well suited to the assessment of large-eddy-simulation-related strategies. The nozzle concentricity is carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both laser Doppler anemometry and particle image velocimetry. The simulations are completed on 6× 106, 12× 106, and 50 × 106 cell meshes. To overcome near-wall modeling problems, a hybrid large-eddy-simulation-Reynolds-averaged-Navier-Stokesrelated method is used. The near-wall Reynolds-averaged-Navier-Stokes layer is helpful in preventing nonphysical separation from the nozzle wall.Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

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.

Flow field measurements of pulverized coal combustion using optical diagnostic techniques

This paper demonstrates the application of laser Doppler velocimetry (LDV) and particle image velocimetry (PIV) techniques to a particle-laden reacting flow of pulverized coal. A laboratory-scale open-type annular burner is utilized to generate velocity profiles of coal particles and micrometric alumina particles. Pair-wise two-component LDV measurements and high-speed stereo PIV measurements provide three-dimensional velocity components of the flow field. A detailed comparison of velocities for alumina and coal particle seeding revealed differences attributed to the wide size distribution of coal particles. In addition, the non-spherical shape and high flame luminosity associated with coal particle combustion introduces noise to the Mie scatter images. The comparison of mean and RMS velocities measured by LDV and PIV techniques showed that PIV measurements are affected by the wide size distribution of coal particles, whereas LDV measurements become biased toward the velocity of small particles, as signals from large particles are rejected. This small-particle bias is also reflected in the spectral characteristics for both techniques, which are in good agreement within the range of frequencies accessible. PIV measurements showed an expected lack of response of large coal particles to the turbulence fluctuations. The overall good agreement between LDV and PIV measurements demonstrates the applicability of the high-speed PIV technique to a particle-laden, high luminosity coal flame while highlighting some of its limitations. © 2013 Springer-Verlag Berlin Heidelberg.